make89 - プログラム倉庫

class ex89{

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 ftuwage2[]=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];
double fw[]=new double[100];
double v[]=new double[100];
double fwinx[][]=new double[100][100];
double fwout[][]=new double[100][100];
double finx[][]=new double[100][100];
double fout[][]=new double[100][100];

int cho;

void makedata(String file){

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

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

data18 di18=new data18();
di18.makedata();
falpha=di18.alpha;
fbeta=di18.beta;
fgamma=di18.gamma;

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

ex20 mk20=new ex20();
mk20.makedata();
finx=mk20.mint;

ex22 mk22=new ex22();
mk22.makedata();
fout=mk22.mint;

ex83 mk83=new ex83();
mk83.makedata();
fwinx=mk83.fwage;

ex84 mk84=new ex84();
mk84.makedata();
fwout=mk84.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]=finx[stage][car];
mos[stage][car]=fout[stage][car];
}

for(age=stage+1;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[age]=(double)z4/z5;
if(v[age] > 1)v[age] = 1;
for(car=2;car<50;car++){  
z6 =(1 - fgamma[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=(fgamma[age - 1] - falpha[age - 1] - fbeta[age-1]) * mis[age-1][car];
mos[age][car] = z6+(1 - falpha[age - 1] - v[age]) * mos[age-1][car];
}

}

for (car=1;car<50;car++){  
wmis[stage][car] =fwinx[stage][car];
wmos[stage][car] =fwout[stage][car];
}
  
for (age =stage+1;age<65;age++){
wmis[age][1]=fw[age];
for (car=2;car<50;car++){  
w1=wmis[age-1][car-1]*(1-fgamma[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]=(fw[age]+(car-1)*w3)/car;
}
for (car=1;car<50;car++){  
w1=wmis[age-1][car]*(fgamma[age-1] - falpha[age-1] - fbeta[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=1;car<25;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]);
}
ftuwage2[byear-cho] = z2/z1;
}

}
}