行列２ - 経済統計練習帳

var number;
var a=new Array();
var b=new Array();
var sol=new Array();
var y=new Array();
var c=new Array();
 
 
function sample(){
 
for(s=1;s<101;s++){
y[s]=100+s;
c[s]=50+0.7*y[s]+Math.random();
}
 
number=2;    
 
for(s=1;s<number+1;s++){
a[s]=new Array();
}
 
d1=0;    
for(s=1;s<101;s++){
d1=d1+c[s];
}
b[1]=d1;
 
 
d1=0;    
for(s=1;s<101;s++){
d1=d1+c[s]*y[s];
}
b[2]=d1;   
 
a[1][1]=100;
 
d1=0;    
for(s=1;s<101;s++){
d1=d1+y[s]*y[s];
}
a[2][2]=d1;     
 
d1=0;    
for(s=1;s<101;s++){
d1=d1+y[s];
}
a[1][2]=d1;    
a[2][1]=d1;
 
for(s=1;s<number+1;s++){
for(sx=1;sx<number+1;sx++){
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];
}
 
$("#memo").html(sol[1]+","+sol[2]);
 
 
}