ｐｈｐ９３ - プログラム倉庫

<?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 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;
}
$maxn=$n;
for ($n=1; $n < $maxn; $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];
}
$sql = "SELECT * FROM data385";
$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 < 116; $n++) {
$age =$matrix[$n][1];  
$m2rate[$age]=$matrix[$n][2];
}

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

for ($age = 15;$age<65;$age++){
for ($car=1;$car<50;$car++){  
$mis[$age][$car]=0;
$mos[$age][$car]=0;
}
}

$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-1]) * $mis[$age-1][$car]  + (1 - $v[$age]) * $mos[$age-1][$car];
}
}
for ($age=15;$age<65;$age++){
$z6 = 0;  
for ($car=1;$car<25;$car++){  
$z6 =$z6+$m2rate[$age]*$alpha[$age]*$mis[$age][$car];
}
for ($car=25;$car<50;$car++){  
$z6 =$z6+$m2rate[$age]*$alpha[$age]*($mis[$age][$car]+$mos[$age][$car]);
}
$mgnew[$byear][$age] = $z6;
}
}

$sql = "truncate eli93";
$exe= mysql_query($sql,$link);
for ($byear=1990;$byear<2035;$byear++){
for ($age=20;$age<65;$age++){
$x=$mgnew[$byear][$age];
$sql = "insert into eli93 values($byear,$age,$x)";
$exe= mysql_query($sql,$link);
}
}
$close_flag = mysql_close($link);
print("93 男子現役遺族新規裁定者1 終了");  
?>