Метод решения СЛАУ

#include <iostream>
#include <cmath>
#include <fstream>
 
using namespace std;
 
double DetValue(double** A, int N)
{
    int p, h, k;
    double det = 0;
    double** B = new double*[N];
    for (int i = 0; i < N; i++)
    {
        B[i] = new double[N];
    }
    if (N == 1)
        return A[0][0];
    else if (N == 2)
    {
        det = (A[0][0] * A[1][1] - A[0][1] * A[1][0]);
        return det;
    }
    else
    {
        for (int p = 0; p < N; p++)
        {
            h = k = 0;
            for (int i = 1; i < N; i++)
            {
                for (int j = 0; j < N; j++)
                {
                    if (j == p)
                        continue;
                    B[h][k] = A[i][j];
                    k++;
                    if (k == N - 1)
                    {
                        h++;
                        k = 0;
                    }
                }
            }
            det = det + A[0][p] * pow(-1, p) * DetValue(B, N - 1);
        }
    }
    for (int i = 0; i < N; i++)
    {
        delete [] B[i];
    }
    delete [] B;
    return det;
}
 
int main()
{
    ifstream fin("input.txt");
    ofstream fout("output.txt");
    int N;
    double tmp, det, dTemp;
    fin >> N;
    double** A = new double*[N];
    for (int i = 0; i < N; i++)
    {
        A[i] = new double[N];
    }
    double*** B = new double**[N];
    for (int i = 0; i < N; i++)
    {
        B[i] = new double*[N];
        for (int j = 0; j < N; j++)
            B[i][j] = new double[N];
    }
    double* X = new double[N];
    double* root = new double[N];
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
        {
            fin >> A[i][j];
        }
        fin >> X[i];
    }
    for (int k = 0; k < N; k++)
    {
        for (int i = 0; i < N; i++)
        {
            for (int j = 0; j < N; j++)
            {
                B[k][i][j] = A[i][j];
            }
        }
    }
    for (int k = 0; k < N; k++)
    {
        for (int i = 0; i < N; i++)
        {
            B[k][i][k] = X[i];
        }
    }
    dTemp = DetValue(A, N);
    for (int i = 0; i < N; i++)
    {
        root[i] = DetValue(B[i], N) / dTemp;
        fout << root[i] << " ";
    }
    fout.close();
    fin.close();
    for (int i = 0; i < N; i++)
    {
        delete [] A[i];
    }
    delete [] A;
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
        {
            delete [] B[i][j];
        }
        delete [] B[i];
    }
    delete [] B;
    delete [] X;
    delete [] root;
    system("pause");
    return 0;
}