ｔｇｆ４５ - プログラム倉庫

package dsge;
 
 
public class pro {
 
    double[] cn=new double[101];
    double[] ln=new double[101];
    double[] k=new double[101];
    double a,beta,ep;
 
    double ks,ls,h;
    int n,mode,time;
 
 
    public static void main(String[] args) {
        pro test=new pro();
    }
pro(){
 
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);  
} 
 
 
 
 
ep=100;
 
mode=0;
 
while(mode<100){
 
simple();
 
 
if(ep<0.001)mode=1000;
 
mode=mode+1;
 
System.out.println(ep);
 
}
 
 
 
 
 
}
void simple(){ 
 
     double[] cns=new double[101];
    double[] lns=new double[101];
    double k1,c1,l1,n1,r1,uc,w1;
 
 
 for (n=10;n<91;n++){ 
k1=k[n]+Math.pow(k[n],a)*Math.pow(ln[n],1-a)-cn[n]; 
n1=k1/h; 
 
c1=con(n1); 
l1=lab(n1); 
r1=a*Math.pow(k1,a-1)*Math.pow(l1,1-a); 
uc=(beta*(1+r1))/c1; 
cns[n]=1/uc; 
w1=(1-a)*Math.pow(k[n],a)*Math.pow(ln[n],-a); 
lns[n]=1-cn[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]; 
} 
 
 
 
} 
double con(double n1){
int n2,n3;
double c1;
n2=(int)n1; 
n3=n2+1;      
c1=cn[n2]+(n1-n2)*(cn[n3]-cn[n2]);     
return c1;    
}
 
double lab(double n1){
int n2,n3;
double c1;
n2=(int)n1; 
n3=n2+1;      
c1=ln[n2]+(n1-n2)*(ln[n3]-ln[n2]);     
return c1;    
}
 
}