行列式 - プログラム倉庫

package svm01;

class pro{

int m,n,mx,nx;
double a[][]=new double[4][4];
double b[]=new double[4];
double c[]=new double[4];
double e;
double z;

public static void main(String[] args) {
pro test=new pro();
}

pro(){

a[1][1]=1;
a[1][2]=1;
a[1][3]=1;
a[2][1]=0;
a[2][2]=1;
a[2][3]=1;
a[3][1]=0;
a[3][2]=0;
a[3][3]=1;
b[1]=5;
b[2]=2;
b[3]=1;
sol();

System.out.println(c[2]);

}

public void matrix(){
for(mx=1;mx<4;mx++){
for(m=1;m<4;m++){
z=(double)a[m][mx]/a[mx][mx];
if(m==mx)z=0;
for(n=1;n<4;n++){
a[m][n]=a[m][n]-z*a[mx][n];
}
b[m]=b[m]-z*b[mx];
}
}

}

public void sol(){

matrix();

for(m=1;m<4;m++){
c[m]=(double)b[m]/a[m][m];
}

}


    
}