netu - 経済統計練習帳

var a=new Array();
var b=new Array();
var sol=new Array();
var number;
var h;
var f;
 
function sample(){
 
var s,sx;
 
f=5;
h=1/100;
 
number=102;
 
for(s=1;s<number+1;s++){
a[s]=new Array();
}
 
for(s=1;s<number+1;s++){
for(sx=1;sx<number+1;sx++){
a[s][sx]=0.0;
}}
 
 
for(s=1;s<101;s++){
a[s][s]=1.0;
a[s][s+1]=-2.0;
a[s][s+2]=1.0;
b[s]=-f*h*h;
}
 
a[101][1]=1.0;
b[101]=0.0;
a[102][102]=1.0;
b[102]=0.0;
 
 
mat();    
 
paint();
 
 
}
 
function paint(){
 
var sub=new svgpaint();
 
str="<svg width=\"700\" height=\"700\">"; 
for(s=1;s<number+1;s++){
str=str+sub.rect(5*s,500-500*sol[s],5,500*sol[s],"black");
}
str=str+"</svg>";
 
$("#memo").html(str);
 
 
}
 
 
function mat(){
 
var s,sx,tx;    
var z;
 
for(s=1;s<number+1;s++){    
for(sx=1;sx<number+1;sx++){
z=0;
z=a[sx][s]/a[s][s];
if(s==sx)z=0;
for(tx=1;tx<number+1;tx++){
a[sx][tx]=a[sx][tx]-z*a[s][tx];
}
b[sx]=b[sx]-z*b[s];
}
}
 
for(s=1;s<number+1;s++){
sol[s]=b[s]/a[s][s];
}
 
}