QB18 - プログラム倉庫

package eli;

class make18{  
int byear,car;
double mis[][]=new double[100][51];
double mos[][]=new double[100][51];
double fnew1[]=new double[200];
double f2by[][]=new double[200][100];
double falpha[]=new double[100];
double fbeta[]=new double[100];
double fgamma[]=new double[100];
double ftheta[]=new double[100];
int age;
int cho;
double z1,z2,z3,z4,z5,z6,z7;
double v;

void makedata(){

make16 mk16=new make16();
mk16.makedata();
f2by=mk16.f2by;

datainput18 di18=new datainput18();
di18.makedata();
falpha=di18.falpha;
fbeta=di18.fbeta;
fgamma=di18.fgamma;

datainput152 di152=new datainput152();
di152.makedata();
ftheta=di152.ftheta;

cho=1900;

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] = f2by[byear-cho][15];

for(age=16;age<65;age++){  
z1= (1 - fgamma[age - 1])* f2by[byear-cho][age-1];
z2=f2by[byear-cho][age] - z1;
if(z2 < 0)z2 = 0;
z3 =(1 - ftheta[age]) * z2;
z4 = ftheta[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=(double)z4/z5;
if(v > 1)v = 1;
for(car=2;car<50;car++){  
z6 =(1 - fgamma[age - 1]) *mis[age-1][car-1];
mis[age][car]=z6+v*mos[age-1][car-1];
}
for(car=1;car<50;car++){
z7=(fgamma[age - 1] - falpha[age - 1] - fbeta[age-1]) * mis[age-1][car];
mos[age][car] = z7+(1 - falpha[age - 1] - v) * mos[age-1][car];
}
}  
z1 = 0;
for(car=25;car<50;car++){
z1 = z1 + mis[64][car] + mos[64][car];
}
fnew1[byear-cho] = z1;
}

}

}