産業連関３ - 経済統計練習帳

package iox;
 
import Jama.*; 
 
public class pro {
 
    double[] y=new double[100];
 double[][] a=new double[3][3];
 double[][] io=new double[3][3];
  double[] v=new double[100];
 double[] d=new double[100];
 
    public static void main(String[] args) {
       pro t=new pro();
    }
 
 
    pro(){
 
 
        a[1][1]=20;
        a[2][1]=30;
        a[1][2]=25;
        a[2][2]=25;
 
        y[1]=120;
        y[2]=180;
 
 
 
        int s,sx;
 
        for(s=1;s<3;s++){
        for(sx=1;sx<3;sx++){    
        io[s][sx]=a[s][sx]/y[sx];
        }}
 
        double y1;
 
        for(s=1;s<3;s++){
        y1=0;
        for(sx=1;sx<3;sx++){  
        y1=y1+io[s][sx]*y[sx];
        }
        d[s]=y[s]-y1;
        }
 
        int number=2;
 
       double[][] beta=new double[number][1];
 
 double[][] alpha=new double[number][number];
 double[][] ra=new double[500][500];
 
  double[][] i=new double[500][500];
 
 
 
for(s=0;s<number;s++){
for(sx=0;sx<number;sx++){
i[s][sx]=0;    
}}
 
for(s=0;s<number;s++){
i[s][s]=1;
}
 
for(s=0;s<number;s++){
for(sx=0;sx<number;sx++){
alpha[s][sx]=i[s][sx]-io[s+1][sx+1];
}    
}  
 
  for(s=0;s<number;s++){
  beta[s][0]=d[s+1];
  }
 
Matrix Ma = new Matrix(alpha);
 
Matrix Mb=new Matrix(beta);
 
Matrix x = Ma.solve(Mb);
 
double[] ys=new double[500];
 
for(s=0;s<number;s++){
ys[s+1]=x.get(s, 0);
System.out.println(ys[s+1]) ;
}
 
 
 
 
 
 
 
    }
 
 
 
}