ｈｄｓふぇ - プログラム倉庫

package dsge01;
 
class pro{ 
 
double a,beta; 
double ks; 
double h; 
double k[]=new double[101]; 
double cx[]=new double[101]; 
double ct[]=new double[101]; 
double k1; 
double n1; 
int n2,n3; 
int n;   
int t; 
 
public static void main(String [] args) { 
        pro test=new pro(); 
 
    } 
 
pro(){ 
 
a=0.33; 
beta=0.95; 
ks=Math.pow((1 / beta - 1) / a , 1 / (a - 1)); 
h=2*(double)ks/100; 
 
for (n=1;n<101;n++){ 
k[n]=n*h; 
cx[n]=Math.pow(k[n],a); 
} 
 
simple(); 
 
 
} 
 
void simple(){ 
int t; 
double k1,n1,c1,r1,ep; 
int n2,n3; 
double cxs[]=new double[101]; 
 
t=0; 
while(t<100){ 
for (n=10;n<91;n++){ 
k1=k[n]+Math.pow(k[n],a)-cx[n]; 
n1=(double)k1/h; 
n2=(int)n1; 
n3=n2+1; 
c1=cx[n2]+(n1-n2)*(cx[n3]-cx[n2]); 
r1=a*Math.pow(k1,a-1); 
cxs[n]=c1/(beta*(1+r1)); 
} 
 
ep=0; 
for (n=10;n<91;n++){ 
ep=ep+Math.pow(cx[n]-cxs[n],2); 
} 
 
for (n=10;n<91;n++){ 
cx[n]=cxs[n]; 
} 
if(ep<0.0001)t=1000;   
t=t+1; 
} 
 
 
} 
 
 
 
 
}