ｄｓげｊｊ - 経済統計練習帳

var k=new Array();
var c=new Array();
var cs=new Array();
var con=new Array();
 
function sample() {
 
a=0.33;
beta=0.95;
 
k1=(1/beta)-1
k1=k1/a;
ks=Math.pow(k1,1/(a-1))
h=ks/50;
for(s=1;s<101;s++){
k[s]=s*h;    
c[s]=Math.pow(k[s],a);
}
 
mode=0;
while(mode<50){
 
for(s=10;s<91;s++){
k1=k[s]+Math.pow(k[s],a)-c[s];
n1=k1/h;
n2=Math.floor(n1);
n3=n2+1;
c1=c[n2]+(n1-n2)*(c[n3]-c[n2]);
cs[s]=c1/((1+a*Math.pow(k1,a-1))*beta);
 
}
 
e=0;
for(s=10;s<91;s++){
e=e+(c[s]-cs[s])*(c[s]-cs[s]);
}  
 
if(e<0.01)mode=100;
 
for(s=10;s<91;s++){
c[s]=cs[s];
}  
 
}
 
k1=k[40];
for(t=1;t<101;t++){
n1=k1/h;
n2=Math.floor(n1);
n3=n2+1;
con[t]=c[n2]+(n1-n2)*(c[n3]-c[n2]);
k1=k1+Math.pow(k1,a)-con[t];
}
 
 
paint();
 
}
 
function paint(){
 
var sub=new svgpaint();
 
str="<svg width=\"500\" height=\"500\">";
str=str+sub.line(0,0,0,500,5,"black");    
str=str+sub.line(0,500,500,500,5,"black");    
 
maxc=con[100];
minc=con[1];
 
for(s=1;s<90;s++){
x1=5*s-5;
y1=500*(con[s]-minc)/(maxc-minc);
y2=500*(con[s+1]-minc)/(maxc-minc);  
str=str+sub.line(x1,500-y1,x1+5,500-y2,3,"red");
}
 
str=str+"</svg>";
 
$("#memo").html(str);
 
 
}