Метод Гаусса

package operationofthesoftware;
import java.util.Scanner;
 
 
public class NewClass {
 
      public static float determinant_Gauss (int n, float[][] a, float det) {
          int k,i,j;
      for (k=0;k<n;k++){
     for (i=k+1;i<n;i++){
      for (j=k+1;j<n;j++){
        a[i][j]=a[i][j]-a[i][k] * a[k][j] / a[k][k];
        det=a[0][0];
        }
     }
      }
      for (i=1;i<n;i++){
       det=det*a[i][i];
      }
         return det;
    
}
 
 public static float[][] invers_Gauss(int n, float[][] a,float[][] add, float [][] invert)
    {
        int k,j;
        float det=0;
       determinant_Gauss (n,  a,  det) ;
        transp(a,n,add); // {транспонируем матрицу}
peresch(n,add,det,invert);  //{считаем дополнения}
 
for (k=0;k<n;k++)
for (j=0;j<n;j++)
invert[k][j]=invert[k][j]/det;// {создаем обратную матрицу}
          return invert;
 
}
 
 static void transp( float[][] a, int n, float[][] at){
int k,j;
for (k= 0;k<n;k++){
for (j=0;j< n;j++){
at[k][j]=a[j][k];
}
}
}
 
static void peresch(int n, float [][] b, float det, float[][] e){
int i,m,k,j;
   float z;
    float[][] d,c;
    d= new float [n][n];
    c= new float [n][n];
for (i=0;i<n;i++)
for (m=0;m<n;m++)
{
     for (j= 0;j<n;j++)//перестановка строк}
     {
         z=b[i][j];
         
         for (k=0;k>2;k--)
         d[k][j]=b[k-1][j];
         for (k=i+1;k<n;k++)
         d[k][j]=b[k][j];
         d[1][j]=z;
     }
     for (k=0;k<n;k++)// {перестановка столбцов}
     {
         z=d[k][m];
         for (j=m;j> 2;j--)
         c[k][j]=d[k][j-1];
         for (j=m+1;j<n;j++)
         c[k][j]=d[k][j];
         c[k][1]=z;
       ;
     determinant_Gauss(n,c,det);//{вычисление определителей}
     e[i][m]=det*znak1(i,m);//{вычисление дополнений}
     }
}
}
 
static int znak1(int i, int m)//{изменение знака при перестановке строк при нахождении дополнений}
{
    int  znak1;
if ((i+m)% 2==0){
    znak1 = 1;
}
else 
{znak1=-1;
}
          return znak1;
}
 
 
 
    public static float[] Gauss(int n, float [][] a, float [] b , float [] x ){
       int k,i,j; 
 
 
    for (i=0;i<n-1;i++)
    {
      for (k=i+1;k<n;k++) 
      {
          float del = a[k][i] / a[i][i]; 
         for (j=0;j<n;j++) 
         { a[k][j]= a[k][j] - a[i][j] * del;}
             b[k]= b[k] - b[i] * del;
      }
    }
        
  //ОБРАТНЫЙ ХОД
    for (i=n;i>0;i--)
    {
      for (j= i+1;j<n;j++){
          b[i]= b[i] - a[i][j] * x[j];}
          x[i]= b[i] / a[i][i];
    }
     return x;    
    }
 
    public static void main (String [] args){
    //Вычисление определителя матрицы методом Гаусса
        
        float det = 0;  // определитель матрицы
    int i,j;        // рабочие переменные
        Scanner in = new Scanner(System.in);
        System.out.println("Введите размер матрицы: ");
        //int  n=in.nextInt();
       int  n=3;
        float [][] a;
        a= new float[n][n];
        float [][] add;
        add= new float[n][n];
        float [][] invert;
        invert= new float[n][n];
        float [] b; //  матрица свободных членов
        b = new float[n];
        float [] x; //  вектор решения (вектор неизвестных)
        x= new float [n];
 
      /*  System.out.println("Введите исходную матрицу: ");
        for ( i = 0; i < n; i++) 
        {
         System.out.println("Строка "+ (i+1) +"("+n+" чисел)");
         for (j = 0; j < n; j++) 
         a[i][j]=in.nextFloat();
        }
             System.out.println();
             
             for (i = 0; i < n; i++) {
                    // цикл по второй размерности
                    for (j = 0; j < n; j++) {
                        //вывод элементов массива 
                    System.out.print(a[i][j] + "\t");
                    }
                    System.out.println();
                }
             
             System.out.println("Введите вектор свободных членов: ");
            for ( i = 0; i < n; i++) {
                System.out.println("b["+i+"]=");  
                b[i]=in.nextFloat();
            }*/
      a[0][0]=3;
      a[0][1]=3;
      a[0][2]=3;
      a[1][0]=2;
      a[1][1]=4;
      a[1][2]=4;
      a[2][0]=1;
      a[2][1]=3;
      a[2][2]=6;
      
      b[0]=1;
       b[1]=1;
        b[2]=1;
      
         float det1 = (float) determinant_Gauss(n, a, det); //вычисление определителя методом гаусса
         
        
      
           System.out.println("Значение определителя(по методу Гаусса)=" + det1);
           System.out.println("Обратная матрица методом Гаусса=" );
           invers_Gauss(n, a,add, invert);//вычисление обратной матрицы методом гаусса
            for (i = 0; i < n; i++) {
                    // цикл по второй размерности
                    for (j = 0; j < n; j++) {
                        //вывод элементов массива 
                    System.out.print(invert[i][j] + "\t");
                    }
                    System.out.println();
                }
    
 
 
        System.out.println("Решение СЛАУ методом Гаусса=" );
       x=Gauss(n, a, b,x);//Решение системы обыкновенных алгебраических уравнений методом Гаусса
          for (i = 0; i < n; i++) {
                        //вывод элементов массива 
                    System.out.print(x[i] + "\t");
                    }
                  
            
      
        
        //Вычисление обратной матрицы методом Гаусса
        //Решение системы обыкновенных алгебраических уравнений методом Гаусса
}
}