まけ３３ - プログラム倉庫

class make33{
int byear,cho,car,stage,age;
double z1,z2,z3,z4,z5,z6,z7;
double v;
double mtime2[]=new double[200];
double mis[][]=new double[100][51];
double mos[][]=new double[100][51];
double m2by[][]=new double[200][100];
double malpha[]=new double[100];
double mbeta[]=new double[100];
double mgamma[]=new double[100];
double mtheta[]=new double[100];
double minx[][]=new double[100][100];
double mout[][]=new double[100][100];

void makedata(){

cho=1900;

make19 mk19=new make19();
mk19.makedata();
minx=mk19.minx;

make21 mk21=new make21();
mk21.makedata();
mout=mk21.mout;

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

datainput151 di151=new datainput151();
di151.makedata();
malpha=di151.malpha;
mbeta=di151.mbeta;
mgamma=di151.mgamma;

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

for(byear=1950; byear<1990; byear++){

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

stage =2008-byear;

for(car=1;car<50;car++){
mis[stage][car] = minx[stage][car];
mos[stage][car] = mout[stage][car];
}

for(age=stage+1;age<65;age++){  
z1= (1 - mgamma[age - 1])* m2by[byear-cho][age-1];
z2=m2by[byear-cho][age] - z1;
if(z2 < 0)z2 = 0;
z3 =(1 - mtheta[age]) * z2;
z4 = mtheta[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 - mgamma[age - 1]) *mis[age-1][car-1];
mis[age][car]=z6+v * mos[age-1][car-1];
}
for(car=1;car<50;car++){
z7=(mgamma[age - 1] - malpha[age - 1] - mbeta[age-1]) * mis[age-1][car];
mos[age][car] = z7+(1 - malpha[age - 1] - v) * mos[age-1][car];
}
}
z1=0;
z2=0;
for(car=25;car<50;car++){
z1 = z1 + mis[64][car] + mos[64][car];
z2 = z2 +car*( mis[64][car] + mos[64][car]);

}
mtime2[byear-cho] = (double)z2/z1;
}

}

}