Ｓｋｌｌ - 経済統計練習帳

package mat;
 
 
import Jama.*; 
 
public class principal {
 
 double[] y=new double[100];
 double[][] a=new double[30][30];
 double[][] io=new double[3][3];
 double[] v=new double[100];
 double[] d=new double[100];
 double[] beta=new double[50];
  int number;
 
 
 
    void makedata(){
 
 
 
 
 
 double[][] alpha=new double[number][number];
 
 
 
 int s,sx;
 
for(s=0;s<number;s++){
for(sx=0;sx<number;sx++){
alpha[s][sx]=a[s+1][sx+1];
}    
}  
 
 
 
Matrix Ma = new Matrix(alpha);
 
 
EigenvalueDecomposition dd=new EigenvalueDecomposition(Ma);
 
 
 double[] x=dd.getRealEigenvalues();
 
 int maxpoint;
 maxpoint=0;
 double maxx;
 maxx=-999;
 for(s=0;s<number;s++){
 if(x[s]>maxx)maxpoint=s;    
 if(x[s]>maxx)maxx=x[s];  
 }
 
 Matrix Mb=dd.getV();
 
 
 for(sx=0;sx<number;sx++){
 beta[sx+1]=Mb.get(sx,maxpoint);
 }    
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 }
 
 
 
}
 
 
 
 
package mat;
 
import Jama.*; 
 
public class pro {
 
 
 double[][] a=new double[30][30];
 
 
 double[] beta=new double[50];
 
 public static void main(String[] args) {
       pro t=new pro();
    }
 
 
    pro(){
 
 
        a[1][1]=1;
        a[1][3]=-0.097;
        a[1][2]=0.376;
        a[2][2]=1;
   a[2][1]=0.376;
       a[2][3]=-0.359;
        a[3][3]=1;
   a[3][1]=-0.097;
       a[3][2]=-0.359;
 
 
        int number=3;
 
        principal sub=new principal();
        sub.number=number;
        sub.a=a;
        sub.makedata();
        beta=sub.beta;
 
        int s;
 
 
 for(s=1;s<number+1;s++){
 System.out.println(beta[s]);
 }   
 
 
 
 
 
 
 
 
 
 
 
 
 
 }
 
 
 
}