Алгоритм Краута

package javaapplication4;
 
public class Main {
 
   public static void main(String[] args) {
       float[][] A = new float[][]{{4,1,1,1},{1,4,1,1},{1,1,4,1},{1,1,1,4}};
            float[] X=new float[4];
            float[] Z=new float[4];
            float[] B=new float[]{1,2,3,4};
            float[][] L= new float[4][4];
            float[][] U = new float[4][4];
 
            for (int i = 0; i < 4; i++) {
                for (int j = 0; j < 4; j++){
                  L[i][j]=0;
                  U[i][j]=0;
                  if(i==j){
                      U[i][i]=1;
                  }
 
                }
 
            }
             for (int i = 0; i < 4; i++) {
                for (int j = 0; j < 4; j++){
                    L[i][0]=A[i][0];
                }
             }
 
                for (int j = 1; j < 4; j++){
                    U[0][j]=A[0][j]/L[0][0];
                }
        
            for(int k=1;k<3;k++){
                for(int i=k;i<4;i++){
                    float sum1=0;
                    for(int m=0;m<k-1;m++){
                        sum1+=L[i][m]*U[m][k];
                    }
                    L[i][k]=A[i][k]-sum1;
                }
                for(int j=k+1;j<4;j++){
                    float sum1=0;
                    for(int m=0;m<k-1;m++){
                        sum1+=L[k][m]*U[m][j];
                    }
                    U[k][j]=(A[k][j]-sum1)/L[k][k];
                }
            }
 
float sum=0;
for(int m=0;m<3;m++)
{
    sum+=L[3][m]*U[m][3];
}
L[3][3]=A[3][3]-sum;
 
 
 
Z[0]=B[0]/L[0][0];
for(int k=1;k<4;k++){
    float sum1=0;
    for(int j=0;j<k-1;j++){
        sum1+=L[k][j]*Z[j];
 
    }
    Z[k]=(B[k]-sum)/L[k][k];
}
 
 
X[3]=Z[3];
for(int k=2;k>0;k--){
    float sum1=0;
    for(int j=k+1;j<3;j++)
    {
       sum1+=U[k][j]*X[j];
    }
    X[k]=Z[k]-sum1;
}
 
    }
 
}

...            
 
int k1 = 0;
        while (k1<3) {
                k1++;
                for (int i=k1; i<4; i++) {
                    float sum1 = 0;
                    for(int m=0; m<k1; m++) {
                        sum1 += L[i][m] * U[m][k1];
                    }
                    L[i][k1] = A[i][k1] - sum1;
                }
                for (int j=k1+1; j<4; j++) {
                    float sum1 = 0;
                    for (int m=0; m<k1; m++) {
                        sum1 += L[k1][m] * U[m][j];
                    }
                    U[k1][j] = (A[k1][j] - sum1) / L[k1][k1];
                }
            }
 
...

Z[0]=B[0]/L[0][0];
for(int k=1;k<4;k++){
    float sum1=0;
    for(int j=0;j<k-1;j++){
        sum1+=L[k][j]*Z[j];
 
    }
    Z[k]=(B[k]-[B]sum1[/B])/L[k][k];
}

Z[0]=B[0]/L[0][0];
for(int k=1;k<4;k++){
    float sum1=0;
    for(int j=0;j<k-1;j++){
        sum1+=L[k][j]*Z[j];
 
    }
    Z[k]=(B[k]-[B]sum1[/B])/L[k][k];
}