ｋｊｒ４５ - プログラム倉庫

package dsge2;
 
 
class pro{ 
double maxc,minc,maxk,mink;   
double uc,w1,a,beta,h,ks,k1,ls,c1,r1,l1; 
int cc1,cc2;   
double th[]=new double[11]; 
double k[]=new double[101]; 
double cx[]=new double[101]; 
double cxs[]=new double[101]; 
double lx[]=new double[101]; 
double lxs[]=new double[101]; 
double cp[]=new double[11]; 
double lp[]=new double[11];
double cps[]=new double[11]; 
double lps[]=new double[11];
 
double px[]=new double[101]; 
double pxs[]=new double[101]; 
double pp[]=new double[11]; 
double pps[]=new double[11]; 
 
 
double n1; 
int t,m,n,mx,n2,n3;   
double ep;         
 
double ms;
 
public static void main(String [] args) { 
        pro f=new pro(); 
 
    } 
 
 
pro(){ 
 
for (m=1;m<11;m++){ 
th[m]=0.95+0.01*m; 
} 
 
 
for (m=1;m<11;m++){ 
lp[m]=0;
cp[m]=0;
}
 
for (n=10;n<91;n++){ 
px[n]=1; 
} 
 
for (m=1;m<11;m++){ 
pp[m]=0; 
} 
 
 
 
 
 
ms=20;
 
m=5;
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; 
lx[n]=ls-lp[m]; 
cx[n]=th[m]*Math.pow(k[n],a)*Math.pow(lx[n]+lp[m],1-a)-cp[m];  
} 
 
double ep1,ep2;
int t,t1,t2;
int mode;
 
mode=0;
 
while (mode<50){
 
    fast();
 
if(ep<0.00001)mode=100;
 
 
    mode=mode+1;
 
}
 
mode=0;
 
while (mode<50){
 
    sec();
 
if(ep<0.00001)mode=100;
 
 
    mode=mode+1;
 
}
 
mode=0;
 
while (mode<100){
 
 
  fast();
 
  ep1=ep;
 
   sec();
 
  ep2=ep;
 
  h=0;
  if(ep1<0.0001)h=h+1;
  if(ep2<0.0001)h=h+1;
if(h==2)mode=100;
 
    mode=mode+1;
 
    System.out.println(ep1+","+ep2);
 
 
}
 
mode=0;
 
while (mode<50){
 
    money1();
 
if(ep<0.00001)mode=100;
 
 
    mode=mode+1;
 
}
 
 
mode=0;
 
while (mode<50){
 
    money2();
 
if(ep<0.00001)mode=100;
 
 
    mode=mode+1;
 
}
mode=0;
 
while (mode<100){
 
 
  money1();
 
  ep1=ep;
 
   money2();
 
  ep2=ep;
 
  h=0;
  if(ep1<0.0001)h=h+1;
  if(ep2<0.0001)h=h+1;
if(h==2)mode=100;
 
    mode=mode+1;
 
    System.out.println(ep1+","+ep2);
 
 
}
 
 
 
 
 
} 
 
double seekdp(double pp){ 
double dp,l1,r1,c1,pc,pi,i1; 
dp=0; 
for (mx=1;mx<11;mx++){ 
l1=lab(mx,n2); 
r1=th[mx]*a*Math.pow(k1,a-1)*Math.pow(l1,1-a); 
c1=con(mx,n2); 
pc=price(mx,n2); 
pi=pc/pp-1; 
i1=(1+r1)*(1+pi)-1; 
dp=dp+i1*ms*beta/(c1*(1+pi)); 
} 
dp=0.1*dp; 
return dp; 
} 
 
void money1(){ 
 
int t1,t2; 
double p1,p2,p3; 
double dp; 
double z1,z2; 
 
m=5;    
 
for (n=10;n<91;n++){ 
p1=1.05*px[n]+pp[m]; 
p2=0.95*px[n]+pp[m]; 
k1=k[n]+th[m]*Math.pow(k[n],a)*Math.pow(lab(m,n),1-a)-con(m,n); 
n1=k1/h; 
dp=seekdp(p1); 
z1=dp-p1; 
t1=0; 
while (t1<100){ 
dp=seekdp(p2); 
z2=dp-p2; 
p3=p2-z2*(p2-p1)/(z2-z1); 
z1=z2; 
p1=p2; 
p2=p3; 
if (z2*z2<0.001)t1=1000;  
t1=t1+1; 
} 
pxs[n]=p2-pp[m]; 
} 
 
 
ep=0; 
for (n=10;n<90;n++){ 
ep=ep+Math.pow(px[n]-pxs[n],2); 
} 
 
for (n=10;n<90;n++){ 
px[n]=pxs[n]; 
} 
 
 
} 
 
void money2(){ 
 
int t1,t2; 
double p1,p2,p3; 
double dp; 
double z1,z2; 
 
n=50;    
 
for (m=1;m<11;m++){ 
p1=1.05*px[n]+pp[m]; 
p2=0.95*px[n]+pp[m]; 
k1=k[n]+th[m]*Math.pow(k[n],a)*Math.pow(lab(m,n),1-a)-con(m,n); 
n1=k1/h; 
dp=seekdp(p1); 
z1=dp-p1; 
t1=0; 
while (t1<100){ 
dp=seekdp(p2); 
z2=dp-p2; 
p3=p2-z2*(p2-p1)/(z2-z1); 
z1=z2; 
p1=p2; 
p2=p3; 
if (z2*z2<0.001)t1=1000;  
t1=t1+1; 
} 
pps[m]=p2-pp[m]; 
} 
 
ep=0; 
for (m=1;m<11;m++){ 
ep=ep+Math.pow(pp[m]-pps[m],2); 
} 
 
 
 
for (m=1;m<11;m++){ 
pp[m]=pps[m]; 
}   
 
} 
 
void fast(){
 
  m=5;
for (n=10;n<90;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<11;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.1*uc; 
cxs[n]=1/uc; 
w1=th[m]*(1-a)*Math.pow(k[n],a)*Math.pow(lab(m,n),-a); 
lxs[n]=1-con(m,n)/w1; 
} 
 
ep=0; 
for (n=10;n<90;n++){ 
ep=ep+Math.pow(cx[n]-cxs[n],2)+Math.pow(lx[n]-lxs[n],2); 
} 
for (n=10;n<90;n++){ 
cx[n]=cxs[n]; 
lx[n]=lxs[n]; 
} 
 
 
}
 
void sec(){
 
  n=50;
for (m=1;m<11;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<11;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.1*uc; 
cps[m]=1/uc-con(m,n); 
w1=th[m]*(1-a)*Math.pow(k[n],a)*Math.pow(lab(m,n),-a); 
lps[m]=1-con(m,n)/w1-lx[n]; 
} 
 
ep=0; 
for (m=1;m<11;m++){ 
ep=ep+Math.pow(cp[m]-cps[m],2)+Math.pow(lp[m]-lps[m],2); 
} 
for (m=1;m<11;m++){ 
cp[m]=cps[m]; 
lp[m]=lps[m]; 
}   
 
 
}
double price(int m,double nx){
double c1;
int n1,n2;
n1=(int)nx;
n2=n1+1;
if(n2>90)n2=90;
if(n1<10)n2=10;
if(n1<10)n1=10;
c1=px[n1]+(nx-n1)*(px[n2]-px[n1])+pp[m];
return c1;
}
 
double lab(int m,double nx){
double c1;
int n1,n2;
n1=(int)nx;
n2=n1+1;
if(n2>90)n2=90;
if(n1<10)n2=10;
if(n1<10)n1=10;
c1=lx[n1]+(nx-n1)*(lx[n2]-lx[n1])+lp[m];
return c1;
}
double con(int m,double nx){
double c1;
int n1,n2;
n1=(int)nx;
n2=n1+1;
if(n1>90)n2=90;
if(n1>90)n1=90;
if(n1<10)n2=10;
if(n1<10)n1=10;
c1=cx[n1]+(nx-n1)*(cx[n2]-cx[n1])+cp[m];
return c1;
}
 
 
}