Метод Гаусса (Java -> C++)

import java.util.Scanner;
import java.lang.Math.*;
public class Main {
    public static void printArr(float[][] arr,int n){
        for (int i=0;i<n;i++){
            for (int j=0;j<n;j++){
                System.out.printf("%f\t",arr[i][j]);
            }
            System.out.println();
        }
    }
    public static void printArrV(float[] arr, int n){
        for (int i=0;i<n;i++){
            System.out.printf("%f\t",arr[i]);
            if ((i%n==4)) System.out.print("\n");
        }
    }
    public static float[][] zapA(int n){
        float[][] a;
        a=new float[n][n];
        a[0][0]=a[1][1]=a[3][3]=a[4][4]=24;
        a[0][1]=a[0][2]=a[0][3]=a[0][4]=a[1][4]=a[2][0]=a[3][0]=a[3][1]=a[3][2]=2;
        a[1][0]=a[1][2]=a[2][1]=a[2][3]=a[4][1]=a[4][2]=a[4][3]=1;
        a[2][2]=32;
        a[1][3]=a[2][4]=a[3][4]=a[4][0]=3;
        return a;
    }
    public static float[] zapX(int n){
        float[] x;
        x=new float[n];
        for (int i=0;i<n;i++){
            x[i]=0;
        }
        return x;
    }
    public static float[] zapY(int n){
        float[] y;
        y=new float[n];
        y[0]=30;
        y[1]=55;
        y[2]=-28;
        y[3]=73;
        y[4]=-17;
        return y;
    }
    public static float[] zapY1(int n, int c){
        float[] y;
        y=new float[n];
        for (int i=0;i<n;i++){
            y[i]=0;
        }
        y[c]=1;
        return y;
    }
    public static float[] gauss(float[][]a, float[] x,float[] y, int n){
        float c;
        for (int k=0;k<n;k++){
            for (int i=k+1;i<n;i++){
                c=a[i][k]/a[k][k];
                for (int j=k;j<n;j++){
                    a[i][j]=a[i][j]-c*a[k][j];
                }
                y[i]=y[i]-c*y[k];
               
            }
        }
        x[n-1]=y[n-1]/a[n-1][n-1];
        for (int i=n-1;i>=0;i--){
            c=0;
            for (int j=i+1;j<n;j++){
                c=c+a[i][j]*x[j];
            }
            x[i]=(y[i]-c)/a[i][i];
        }
        return x;
    }
    public static float[] prIt(float[][]a, float[] x, float[] y, int n, float e){
        float aMax=a[0][0];
        float aMin=a[0][0];
        float b;
        for (int i=1;i<n;i++){
            if (a[i][i]>aMax){
                aMax=a[i][i];
            }
            if (a[i][i]<aMin){
                aMin=a[i][i];
            }
        }
        b=2/(aMin+aMax);
        float[] x1,x2;
        float sum1,sum2;
        int count=0;
        x2=new float[n];
        do {
            sum1=0;
            sum2=0;
            x1=zapX(n);
            for (int i=0;i<n;i++){
                x2[i]=x[i];
                for (int j=0;j<n;j++){
                    x1[i]=x1[i]+a[i][j]*x[j];
                }
                x1[i]=x1[i]-y[i];
                x[i]=x[i]-b*x1[i];
            }
            count++;
            for (int i=0;i<n;i++){
                sum1=sum1+x[i];
                sum2=sum2+x2[i];
            }
        }while (Math.abs(sum1-sum2)/5>e);
        System.out.printf("Метод простой итерации завершился за %d итераций\n",count);
        return x;
    }
    public static float[] zeid(float[][]a, float[] x, float[] y, int n, float e){
        float[] x1,x2;
        int count=0;
        float sum1,sum2,s;
        x2=new float[n];
        x1=zapX(n);
       do{
            sum1=0;
            sum2=0;
            for (int i=0;i<n;i++){
                s=0;
                for (int j=0;j<n;j++){
                    if ((j==i)&&(j<n-1)) j++;
                    if ((j==i)&&(j==n-1)) break;
                    s=s+a[i][j]/a[i][i]*x[j];
                }
                x[i]=y[i]/a[i][i]-s;
                x2[i]=x1[i];
                x1[i]=x[i];
            }
            count++;
            for (int i=0;i<n;i++){
                sum1=sum1+x[i];
                sum2=sum2+x2[i];
            }
        } while (Math.abs(sum1-sum2)/5>e);
        System.out.printf("Метод Зейделя завершился за %d итераций\n",count);
        return x;
    }
    public static void main(String[] args) {
        float[] y,x;
        float[][] a,x1;
        int n=5;
        x1=new float[n][n];
        a=zapA(n);
        y=zapY(n);
        x=zapX(n);
        System.out.print("Исходный массив А:\n");
        printArr(a,n);
        System.out.print("Исходный вектор y:\n");
        printArrV(y,n);
        gauss(a, x, y, n);
        System.out.print("Вектор x(по gauss):\n");
        printArrV(x,n);
        System.out.printf("Проверка корней:\n");
        y=zapY(n);
        a=zapA(n);
        float[] x2;
        x2=new float[n];
        for (int i=0;i<n;i++){
            x2[i]=0;
            for (int j=0;j<n;j++){
                x2[i]=x2[i]+a[i][j]*x[j];
            }
        }
        printArrV(x2,n);
        printArrV(y,n);
        for (int i=0;i<n;i++){
            y=zapY1(n,i);
            a=zapA(n);
            x=zapX(n);
            x1[i]=gauss(a,x,y,n);
        }
        System.out.print("Обратная матрица(транспонированная)(по gauss):\n");
        printArr(x1,n);
        a=zapA(n);
    
        float[][] c;
        c=new float[n][n];
        for (int i=0;i<n;i++){
            for (int j=0;j<n;j++){
                c[i][j]=0;
            }
        }
        for (int i=0;i<n;i++){
            for (int j=0;j<n;j++){
                for (int k=0;k<n;k++){
                    c[i][j]=c[i][j]+a[i][k]*x1[j][k];
                }
            }
        }
        System.out.printf("Произведение матриц:\n");
        printArr(c,n);
        System.out.print("Введите точность для метода простой итерации и Зейделя:\n");
        float e;
        Scanner scan = new Scanner(System.in);
        e=scan.nextFloat();
        a=zapA(n);
        y=zapY(n);
        x=zapX(n);
        x=prIt(a,x,y,n,e);
        System.out.print("Вектор х(по простой итерации):\n");
        printArrV(x,n);
        a=zapA(n);
        y=zapY(n);
        x=zapX(n);
        x=zeid(a,x,y,n,e);
        System.out.print("Вектор х(по Зейделю):\n");
        printArrV(x,n);
    }
}