いれ９８７６ - プログラム倉庫

package dsge;
 
 
public class pro {
 
    double[] chcn=new double[101];
    double[] chln=new double[101];
    double[] chcm=new double[101];
    double[] chlm=new double[101];
 
    double[] cn=new double[101];
    double[] ln=new double[101];
    double[] cm=new double[101];
    double[] lm=new double[101];
    double[] k=new double[101];
    double a,beta,ep;
double[] th=new double[101];
 
    double ks,ls,h;
    int n,m,mode,time;
 
 
    public static void main(String[] args) {
        pro test=new pro();
    }
pro(){
 
 
    for (m=1;m<101;m++){ 
chlm[m]=0;
chcm[m]=0;
}
 
for (n=1;n<101;n++){ 
chln[n]=0;
chcn[n]=0;
}  
 
 
  for (m=1;m<101;m++){   
   th[m]=0.95+0.001*m; 
  }
 
 
for (m=1;m<101;m++){ 
lm[m]=0;
cm[m]=0;
}
 
a=0.33;
beta=0.95;
ls=(1-a)/(2-a);
 
ks=ls*Math.pow((1 / beta - 1) / a , 1 / (a - 1)); 
h=2*ks/100; 
for (n=1;n<101;n++){ 
k[n]=n*h; 
ln[n]=ls; 
cn[n]=Math.pow(k[n],a)*Math.pow(ln[n],1-a);  
} 
 
int card;
 
card=0;
 
while(card<100){
 
mode=0;
 
while(mode<100){
 
fast();
 
if(ep<0.0001)mode=1000;
 
mode=mode+1;
 
}
 
mode=0;
 
while(mode<100){
 
sec();
 
 
if(ep<0.0001)mode=1000;
 
mode=mode+1;
 
}
 
 
ep=0; 
for (n=10;n<91;n++){ 
ep=ep+Math.pow(cn[n]-chcn[n],2)+Math.pow(ln[n]-chln[n],2); 
} 
 
for (n=10;n<91;n++){ 
chcn[n]=cn[n]; 
chln[n]=ln[n]; 
} 
 
 
for (m=1;m<101;m++){ 
ep=ep+Math.pow(cm[m]-chcm[m],2)+Math.pow(chlm[n]-chlm[n],2); 
} 
 
for (m=10;m<91;m++){ 
chcm[n]=cm[n]; 
chlm[n]=lm[n]; 
} 
 
card=card+1;
 
if(ep<0.0001)card=100;
 
System.out.println(ep);
 
 
 
 
}
 
}
 
 
 
 
void fast(){
 
int m,n,mx;    
double k1,uc,c1,l1,r1,w1;
double lns[]=new double[101]; 
double cns[]=new double[101]; 
double n1;  
 
m=50;
for (n=10;n<91;n++){ 
k1=k[n]+th[m]*Math.pow(k[n],a)*Math.pow(lab(m,n),1-a)-con(m,n); 
n1=k1/h; 
uc=0; 
for (mx=1;mx<101;mx++){ 
c1=con(mx,n1); 
l1=con(mx,n1); 
r1=th[mx]*a*Math.pow(k1,a-1)*Math.pow(l1,1-a); 
uc=uc+(beta*(1+r1))/c1; 
} 
uc=0.01*uc; 
cns[n]=1/uc; 
w1=th[m]*(1-a)*Math.pow(k[n],a)*Math.pow(lab(m,n),-a); 
lns[n]=1-con(m,n)/w1; 
} 
 
ep=0; 
for (n=10;n<91;n++){ 
ep=ep+Math.pow(cn[n]-cns[n],2)+Math.pow(ln[n]-lns[n],2); 
} 
 
for (n=10;n<91;n++){ 
cn[n]=cns[n]; 
ln[n]=lns[n]; 
} 
 
 
}
 
void sec(){
 
int m,n,mx;    
double k1,uc,c1,l1,r1,w1;
double lms[]=new double[101]; 
double cms[]=new double[101]; 
double n1;  
 
n=50;
for (m=1;m<101;m++){ 
k1=k[n]+th[m]*Math.pow(k[n],a)*Math.pow(lab(m,n),1-a)-con(m,n); 
n1=k1/h; 
uc=0; 
for (mx=1;mx<101;mx++){ 
c1=con(mx,n1); 
l1=con(mx,n1); 
r1=th[mx]*a*Math.pow(k1,a-1)*Math.pow(l1,1-a); 
uc=uc+(beta*(1+r1))/c1; 
} 
uc=0.01*uc; 
cms[n]=1/uc-cn[n]; 
w1=th[m]*(1-a)*Math.pow(k[n],a)*Math.pow(lab(m,n),-a); 
lms[m]=1-con(m,n)/w1-ln[n]; 
} 
 
ep=0; 
for (m=1;m<101;m++){ 
ep=ep+Math.pow(cm[m]-cms[m],2)+Math.pow(lm[m]-lms[m],2); 
} 
 
 
for (m=1;m<101;m++){ 
cm[m]=cms[m]; 
lm[m]=lms[m]; 
} 
 
 
}
 
double lab(int m,double nx){
double c1;
int n1,n2;
n1=(int)nx;
n2=n1+1;
c1=ln[n1]+(nx-n1)*(ln[n2]-ln[n1])+lm[m];
return c1;
 
}   
 
double con(int m,double nx){
double c1;
int n1,n2;
n1=(int)nx;
n2=n1+1;
c1=cn[n1]+(nx-n1)*(cn[n2]-cn[n1])+cm[m];
return c1;
 
}   
 
}