PHP36 - プログラム倉庫

<?php
require("calldata.php");

$sql = "SELECT * FROM eli16";
$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;
}
$maxn=$n;
for ($n=1; $n < $maxn; $n++) {
$byear =$matrix[$n][1];
$age =$matrix[$n][2];
$m2by[$byear][$age] = $matrix[$n][3];
$f2by[$byear][$age] = $matrix[$n][4];
}

$sql = "SELECT * FROM data14";
$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 data151";
$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 data152";
$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;
}
$maxn=$n;
for ($n=1; $n < $maxn; $n++) {
$age =$matrix[$n][1];
$theta[$age] = $matrix[$n][2];
}



For ($byear=1990;$byear<2035;$byear++){

For ($age = 16;$age<65;$age++){
For ($car=1;$car<50;$car++){  
$mis[$age][$car] = 0;
$mos[$age][$car] = 0;
}
}
$age = 15;
$mis[15][1] = $m2by[$byear][15];
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[$age][1] = $z3;
$z5 = 0;
For ($car = 1;$car<50;$car++){  
$z5 = $z5 + $mos[$age-1][$car];
}
If ($z5 == 0)$z5 = 1;  
$v[$age] = $z4 / $z5;
If ($v[$age]> 1)$v[$age] = 1;  
For ($car=2;$car<50;$car++){  
$mis[$age][$car] = (1-$gamma[$age-1]) * $mis[$age-1][$car-1] + $v[$age] * $mos[$age-1][$car-1];
}
For ($car=1;$car<50;$car++){  
$mos[$age][$car]  = ($gamma[$age-1] - $alpha[$age-1] - $beta[$age]) * $mis[$age-1][$car]  + (1 - $v[$age]) * $mos[$age-1][$car];
}
}

For ($age = 16;$age<65;$age++){
For ($car=1;$car<50;$car++){  
$wmis[$age][$car] = 0;
$wmos[$age][$car] = 0;
}
}
$wmis[15][1]=$mw[15];
For ($age = 16;$age<65;$age++){
$wmis[$age][1]=$mw[$age];
For ($car=2;$car<50;$car++){  
$w1=$wmis[$age-1][$car-1]*(1-$gamma[$age-1]) * $mis[$age-1][$car-1] + $wmos[$age-1][$car-1]*$v[$age] * $mos[$age-1][$car-1];
$w2=$mis[$age][$car]+0.00001;
$w3=$w1/$w2;
$wmis[$age][$car]=($mw[$age]+($car-1)*$w3)/$car;
}
For ($car=1;$car<50;$car++){  
$w1=$wmis[$age-1][$car]*($gamma[$age-1] - $alpha[$age-1] - $beta[$age]) * $mis[$age-1][$car]  + $wmos[$age-1][$car]*(1 - $v[$age]) * $mos[$age-1][$car];
$w2=$mos[$age][$car]+0.00001;
$wmos[$age][$car]=$w1/$w2;
}
}
$z1=0.00001;
$z2=0;
For ($car=25;$car<50;$car++){  
$z1 = $z1 + $mis[64][$car]+$mos[64][$car];
$z2 = $z2 +$wmis[64][$car]*$mis[64][$car]+$wmos[64][$car]*$mos[64][$car];
}
$wage[$byear] = $z2/$z1;
}

$sql = "truncate eli36";
$exe= mysql_query($sql,$link);
for ($byear=1990;$byear<2035;$byear++){
$x=$wage[$byear];
$sql = "insert into eli36 values($byear,$x)";
$exe= mysql_query($sql,$link);
}
$close_flag = mysql_close($link);
print("36 男子賃金１ 計算終了");
?>