dsge01 x - プログラム倉庫

var a,beta,h;
k = new Array(101);
cx = new Array(101);
c = new Array(101);
 
function sample(){
 
a=0.33;
beta=0.95;
 
faststep();
 
dsge();
 
path();
 
line();
 
}
 
function faststep(){
 
var s;
var ks;
 
ks=Math.pow((1/ beta - 1)/a ,1/(a - 1)); 
h=2*ks/100;
 
for(s=1;s<101;s++){
k[s]=s*h;
}
 
for(s=1;s<101;s++){
cx[s]=f(k[s]);
}
 
}
 
function dsge(){
 
cxs = new Array(101);
var mode;
var n1;
var uc,c1,k1,r1;
var ep; 
var s;
 
mode=0;
 
while(mode<100){
 
for(s=10;s<91;s++){
k1=k[s]+f(k[s])-cx[s];
n1=k1/h;
c1=con(n1);
r1=a*Math.pow(k1,a-1);
uc=(beta*(1+r1))/c1;
cxs[s]=1/uc;
}
 
ep=0;
for(s=10;s<91;s++){
ep=ep+(cx[s]-cxs[s])*(cx[s]-cxs[s]);
}
 
for(s=10;s<91;s++){
cx[s]=cxs[s];
}
 
if(ep<0.0001)mode=1000;
 
mode=mode+1;
}
 
}
 
function path(){
 
var t;
var k1,c1; 
var n1;
 
k1=k[60];
 
for(t=1;t<101;t++){
n1=k1/h;
c1=con(n1);
k1=k1+Math.pow(k1,a)-c1;
c[t]=c1;
}
 
}
 
function con(n1){
var c1;
var n2,n3;
n2=Math.floor(n1);
n3=n2+1;
c1=cx[n2]+(n1-n2)*(cx[n3]-cx[n2]);
return c1;
}
 
 
function f(k1){
var f1;
f1=Math.pow(k1,a);
return f1;
}