９６ ｊ８６ - プログラム倉庫

class ex86{

double z1,z2,z3,z4,z5,z6;
double w1,w2,w3,w4,w5,w6;
int byear,age,car,stage;
double mis[][]=new double[100][51];
double mos[][]=new double[100][51];
double wmis[][]=new double[100][51];
double wmos[][]=new double[100][51];
double mwage2[]=new double[200];
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 mw[]=new double[100];
double v[]=new double[100];
double mwinx[][]=new double[100][100];
double mwout[][]=new double[100][100];
double minx[][]=new double[100][100];
double mout[][]=new double[100][100];

int cho;

void makedata(String file){

data14 di14=new data14();
di14.makedata();
mw=di14.mw;

ex16 mk16=new ex16();
mk16.makedata(file);
m2by=mk16.m2by;

data151 di151=new data151();
di151.makedata();
malpha=di151.alpha;
mbeta=di151.beta;
mgamma=di151.gamma;

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

ex19 mk19=new ex19();
mk19.makedata();
minx=mk19.mint;

ex21 mk21=new ex21();
mk21.makedata();
mout=mk21.mint;

ex82 mk82=new ex82();
mk82.makedata();
mwinx=mk82.mwage;

ex85 mk85=new ex85();
mk85.makedata();
mwout=mk85.mwage;

cho=1900;

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;
wmis[age][car] =0;
wmos[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]=(double)z4/z5;
if(v[age] > 1)v[age] = 1;
for(car=2;car<50;car++){  
z6 =(1 - mgamma[age - 1]) *mis[age-1][car-1];
mis[age][car]=z6+v[age] * mos[age-1][car-1];
}
for(car=1;car<49;car++){
z6=(mgamma[age - 1] - malpha[age - 1] - mbeta[age-1]) * mis[age-1][car];
mos[age][car] = z6+(1 - malpha[age - 1] - v[age]) * mos[age-1][car];
}

}

for (car=1;car<50;car++){  
wmis[stage][car] =mwinx[stage][car];
wmos[stage][car] =mwout[stage][car];
}
  
for (age =stage+1;age<65;age++){
wmis[age][1]=mw[age];
for (car=2;car<50;car++){  
w1=wmis[age-1][car-1]*(1-mgamma[age-1]) * mis[age-1][car-1];
w1=w1+wmos[age-1][car-1]*v[age] * mos[age-1][car-1];
w2=mis[age][car]+0.00001;
w3=(double)w1/w2;
wmis[age][car]=(mw[age]+(car-1)*w3)/car;
}
for (car=1;car<50;car++){  
w1=wmis[age-1][car]*(mgamma[age-1] - malpha[age-1] - mbeta[age-1]) * mis[age-1][car];  
w1=w1+wmos[age-1][car]*(1 - v[age]) * mos[age-1][car];
w2=mos[age][car]+0.00001;
wmos[age][car]=(double)w1/w2;
}

}

z1=0.000001;
z2=0;
for(car=25;car<50;car++){
z1 = z1 + car*(mis[64][car] + mos[64][car]);
z2 = z2 + car*(wmis[64][car]*mis[64][car] +wmos[64][car]* mos[64][car]);
}
mwage2[byear-cho] = z2/z1;
}

}
}