産業連関表１０ - プログラム倉庫

var x=new Array();
var y=new Array();
var ex=new Array();
var im=new Array();
var s,number;
 
var c=new Array();
var rate=new Array();
var alpha=new Array();
 
var ys=new Array();
var ims=new Array();
var cs=new Array();
var delta=new Array();
 
var gamma=new Array();
 
function sample(){
 
 
number=2;
 
for(s=1;s<number+1;s++){    
x[s]=new Array();
rate[s]=new Array();
gamma[s]=new Array();
}
 
y[1]=100;
y[2]=200;
 
 
x[1][1]=50;
x[1][2]=10;
x[2][1]=30;
x[2][2]=20;
 
ex[1]=10;
ex[2]=20;
im[1]=1;
im[2]=2;
 
faststep();    
 
for(n=1;n<number+1;n++){
 
sol(n);
 
for(s=1;s<number+1;s++){
gamma[n][s]=delta[s];
}
 
}    
 
 
 
}
 
function sol(page){
 
cs=c;
ys=y;
ims=im;
 
cs[page]=cs[page]+1;
 
t=0;
while(t<50){    
t=t+1;
cal();    
}
 
for(s=1;s<number+1;s++){    
delta[s]=ys[s]-y[s];
}
 
}
 
function cal(){
 
var yp=new Array();
var imp=new Array();  
 
for(s=1;s<number+1;s++){    
y1=cs[s]+ex[s]-ims[s];
x1=0;
for(ou=1;ou<number+1;ou++){       
x1=x1+rate[ou][s]*ys[ou];
}
yp[s]=y1+x1;
}    
 
for(s=1;s<number+1;s++){    
x1=cs[s];
for(ou=1;ou<number+1;ou++){       
x1=x1+rate[ou][s]*ys[ou];
}
imp[s]=alpha[s]*x1;
}      
 
ims=imp;
ys=yp;    
 
}
 
function faststep(){
 
 
var c1,x1,ou,ip;    
 
for(s=1;s<number+1;s++){    
c1=y[s]-ex[s]+im[s];
x1=0;
for(ou=1;ou<number+1;ou++){       
x1=x1+x[ou][s];
}
c[s]=c1-x1;
}
 
 
 
for(ou=1;ou<number+1;ou++){   
  for(ip=1;ip<number+1;ip++){  
  rate[ou][ip]=x[ou][ip]/y[ou];
  }
 }
 
for(s=1;s<number+1;s++){   
x1=c[s];  
for(ou=1;ou<number+1;ou++){       
x1=x1+x[ou][s];
}    
alpha[s]=im[s]/x1;    
}
 
 
}