産業５ - 経済統計練習帳

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];
   double[] ys=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;
 
 
        mat sub=new mat();
        sub.number=number;
        sub.a=a;
        sub.d=d;
        sub.y=y;
        sub.makedata();
        ys=sub.ys;
 
        System.out.println(ys[2]);
 
 
 
 
 
 
 
 
 
    }
 
 
 
}
 
 
 
package iox;
 
import Jama.*; 
 
public class mat{
 
    double[] y=new double[100];
 double[][] a=new double[100][100];
 double[][] io=new double[100][100];
  double[] v=new double[100];
 double[] d=new double[100];
  int number;
     double[] ys=new double[100];
 
 
    void makedata(){
 
 
 
 
 
        int s,sx;
 
        for(s=1;s<number+1;s++){
        for(sx=1;sx<number+1;sx++){    
        io[s][sx]=a[s][sx]/y[sx];
        }}
 
 
 
 
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);
 
for(s=0;s<number;s++){
ys[s+1]=x.get(s, 0);
}
 
 
 
 
 
 
 
    }
 
 
 
}