まけ９５ - プログラム倉庫

class make95{

double z1,z2,z3,z4,z5,z6;
int byear,car,cho,age,stage;
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 mgnew2[][]=new double[200][100];
double v[]=new double[100];
double m2rate[]=new double[200];
double minx[][]=new double[100][100];
double mout[][]=new double[100][100];

void makedata(){

cho=1900;

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;

datainput401 di401=new datainput401();
di401.makedata();
m2rate=di401.m2rate;

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

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

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

for (age = 15;age<65;age++){
for (car=1;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[age]=z4 /z5;
if (v[age] > 1)v[age]=1;  
for (car=2;car<50;car++){  
mis[age][car] = (1-mgamma[age-1]) * mis[age-1][car-1] + v[age] * mos[age-1][car-1];
}
for (car=1;car<50;car++){  
mos[age][car]  = (mgamma[age-1] - malpha[age-1] - mbeta[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]*malpha[age]*mis[age][car];
}
for (car=25;car<50;car++){  
z6 =z6+m2rate[age]*malpha[age]*(mis[age][car]+mos[age][car]);
}
mgnew2[byear-cho][age] = z6;
}
}


}
}