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

package memo2;
 
 
public class pro {
 
   double[][] x=new double[3][3];
   double[] d=new double[3];
   double[] w=new double[3];
  double[] y=new double[3];
   double[] ex=new double[3];
   double[] im=new double[3];
   double[] c=new double[3];
  int s,sx,number; 
  double[][] io=new double[3][3];
  double[] alpha=new double[3];
   double[] beta=new double[3];
   double[] theta=new double[3];
 
   double[] cs=new double[3];
   double[] ys=new double[3];
   double[] ims=new double[3];
   double[][] xs=new double[3][3];
 
    public static void main(String[] args) {
        pro test=new pro();
    }
 pro(){   
 
     number=2;
  y[1]=100;
  x[1][1]=20;
  x[1][2]=30;
  d[1]=20;
  w[1]=30;
  y[2]=200;
  x[2][1]=60;
  x[2][2]=30;
  d[2]=70;
  w[2]=40;
  ex[1]=10;
  im[1]=30;
  ex[2]=6;
  im[2]=40;
 
  for(s=1;s<number+1;s++){   
  c[s]=y[s]+im[s]-ex[s];
  }
 
  for(s=1;s<number+1;s++){   
  for(sx=1;sx<number+1;sx++){ 
  io[s][sx]=x[s][sx]/y[s];
  }
  }
 
   for(s=1;s<number+1;s++){   
   alpha[s]=w[s]/y[s];
   beta[s]=d[s]/y[s];
   }
 
   for(s=1;s<number+1;s++){   
   theta[s]=im[s]/(y[s]+im[s]-ex[s]);
   }
 
 
  cs[1]=c[1]+10;
  cs[2]=c[2];
  ys=y;
  xs=x;
  ims=im;
 
  int t;
  t=0;
  while(t<30){
  step();
  t=t+1;
  }
 
 
 
 
 
 }
void step(){
 
double[] yp=new double[3];
double[] imp=new double[3];
double[][] xp=new double[3][3];
 
for(s=1;s<number+1;s++){
for(sx=1;sx<number+1;sx++){
xp[s][sx]=io[s][sx]*ys[s];
}
}  
 
for(s=1;s<number+1;s++){
imp[s]=theta[s]*(xs[1][s]+xs[2][s]+cs[s]);
}
 
double y1;
 
for(s=1;s<number+1;s++){
y1=ex[s]+cs[s]-ims[s];
for(sx=1;sx<number+1;sx++){
y1=y1+x[sx][s];
}
yp[s]=y1;
}
 
ys=yp;
ims=imp;
xs=xp;
 
}
}