panda28 女子通算平均賃金 - プログラム倉庫

<?php
require("calldata.php");
$sql = "SELECT * FROM pgdata05";
$result = mysql_query($sql, $link);
$n=1;
while( $row = mysql_fetch_row( $result ) ){
$matrix[$n][1] = $row[0];
$matrix[$n][2]  =$row[1];
$matrix[$n][3]  =$row[2];
$matrix[$n][4]  =$row[3];
$n=$n+1;
}
for ($n=1; $n < 5501; $n++) {
$year =$matrix[$n][1];
$age =$matrix[$n][2];
$m2[$year][$age] = $matrix[$n][3];
$f2[$year][$age] = $matrix[$n][4];
}

$sql = "SELECT * FROM pdata09";
$result = mysql_query($sql, $link);
$n=1;
while( $row = mysql_fetch_row( $result ) ){
$matrix[$n][1] = $row[0];
$matrix[$n][2]  =$row[1];
$matrix[$n][3]  =$row[2];
$n=$n+1;
}
for ($n=1; $n < 56; $n++) {
$age =$matrix[$n][1];
$mw[$age] = $matrix[$n][2];
$fw[$age] = $matrix[$n][3];
}

$sql = "SELECT * FROM pdata12";
$result = mysql_query($sql, $link);
$n=1;
while( $row = mysql_fetch_row( $result ) ){
$matrix[$n][1] = $row[0];
$matrix[$n][2]  =$row[1];
$matrix[$n][3]  =$row[2];
$matrix[$n][4]  =$row[3];
$n=$n+1;
}
for ($n=1; $n < 56; $n++) {
$age =$matrix[$n][1];
$gamma[$age] = $matrix[$n][2];
$alpha[$age] = $matrix[$n][3];
$beta[$age] = $matrix[$n][4];
}

$sql = "SELECT * FROM pdata11";
$result = mysql_query($sql, $link);
$n=1;
while( $row = mysql_fetch_row( $result ) ){
$matrix[$n][1] = $row[0];
$matrix[$n][2]  =$row[1];
$matrix[$n][3]  =$row[2];
$n=$n+1;
}
for ($n=1; $n < 46; $n++) {
$age =$matrix[$n][1];
$theta[$age] = $matrix[$n][3];
}

For ($year = 5;$year<100;$year++){
For ($age = 15;$age<70;$age++){
$byear=$year-$age+2000;
$m2by[$byear][$age]=$m2[$year][$age];
}
}
For ($byear=1990;$byear<2035;$byear++){
$age = 15;
$mis[$byear][$age][1] = $f2by[$byear][$age];
For ($age = 16;$age<65;$age++){
$z1 =(1-$gamma[$age-1])*$m2by[$byear][$age-1];
$z2 = $m2by[$byear] [$age]-$z1;
If ($z2 < 0) {
$z2 = 0;
}
$z3 = (1 - $theta[$age]) * $z2;
$z4 = $theta[$age] * $z2;
$mis[$byear][$age][1] = $z3;
$z5 = 0;
For ($car = 1;$car<50;$car++){  
$z5 = $z5 + $mos[$byear][$age-1][$car];
}
If ($z5 == 0){
$z5 = 1;
}
$v = $z4 / $z5;
If ($v > 1){
$v = 1;
}
For ($car=2;$car<50;$car++){  
$mis[$byear][$age][$car] = (1-$gamma[$age-1]) * $mis[$byear][$age-1][$car-1] + $v * $mos[$byear][$age-1][$car-1];
}
For ($car=1;$car<50;$car++){  
$mos[$byear][$age][$car]  = ($gamma[$age-1] - $alpha[$age-1] - $beta[$age]) * $mis[$byear][$age-1][$car]  + (1 - $v) * $mos[$byear][$age-1][$car];
}
$wmis[$byear][$age][1]=$fw[$age];
For ($car=2;$car<50;$car++){  
$w1=$wmis[$byear][$age-1][$car-1]*(1-$gamma[$age-1]) * $mis[$byear][$age-1][$car-1] + $wmos[$byear][$age-1][$car-1]*$v * $mos[$byear][$age-1][$car-1];
$w2=$mis[$byear][$age][$car]+0.00001;
$w3=$w1/$w2;
$wmis[$byear][$age][$car]=($fw[$age]+($car-1)*$w3)/$car;
}
For ($car=1;$car<50;$car++){  
$w1=$wmis[$byear][$age-1][$car]*($gamma[$age-1] - $alpha[$age-1] - $beta[$age]) * $mis[$byear][$age-1][$car]  + $wmos[$byear][$age-1][$car]*(1 - $v) * $mos[$byear][$age-1][$car];
$w2=$mos[$byear][$age][$car]+0.00001;
$wmos[$byear][$age][$car]=$w1/$w2;
}
}
$z1=0.00001;
$z2=0;
For ($car=1;$car<25;$car++){  
$z1 = $z1 + $mis[$byear][64][$car]+$mos[$byear][64][$car];
$z2 = $z2 +$wmis[$byear][64][$car]*$mis[$byear][64][$car]+$wmos[$byear][64][$car]*$mos[$byear][64][$car];
}
$wage[$byear] = $z2/$z1;
}
$sql = "truncate pgdata28";
$exe= mysql_query($sql,$link);

for ($byear=1990;$byear<2035;$byear++){
$x=$wage[$byear];
$sql = "insert into pgdata28 values($byear,$x)";
$exe= mysql_query($sql,$link);
}
$close_flag = mysql_close($link);
for ($byear=1990; $byear < 2034; $byear++) {
print($wage[$byear]);
print(",");
}
$year=2034;
print($wage[$byear]);
?>