Нахождение обратной матрицы на основе LU- разложения

public static Matrix Reverse(Matrix A){
        Matrix L, U, A1 = new Matrix(A.rows);
        int rows = A.rows;
        L = new Matrix(A.L);
        U = new Matrix(A.U);
        LU(A);
        L.set(A.L);
        U.set(A.U);
        for (int i = rows-1; i >= 0; i--) {
            for (int j = rows-1; j >= 0; j--) {
                if(i == j) {
                    double sum = 0.0;
                    for (int k = j - 1; k >= 0; k--) {
                        sum += U.matrix[j][k]*A1.matrix[k][j];                     
                    }
                    A1.matrix[j][j] = (1 - sum) / U.matrix[j][j];
                }
                else if (i < j){
                    double sum = 0.0;
                    for (int k = i - 1; k >= 0; k--) {
                        sum += U.matrix[i][k]*A1.matrix[k][i];
                    }
                    A1.matrix[i][j] = sum / (-U.matrix[i][i]);
                }
                else{
                    double sum = 0.0;
                    for (int k = j - 1; k >= 0; k--) {
                        sum += L.matrix[k][j]*A1.matrix[i][k];
                    }
                    A1.matrix[i][j] = -sum;
                }
            }
        }
        return A1;
    }