Обратная матрица методом Гаусса-Жордана

#include <iostream>
 
void inversion(double **A, int N)
{
    double temp;
 
    double **E = new double *[N];
 
    for (int i = 0; i < N; i++)
        E[i] = new double [N];
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
        {
            E[i][j] = 0.0;
 
            if (i == j)
                E[i][j] = 1.0;
        }
 
    for (int k = 0; k < N; k++)
    {
        temp = A[k][k];
 
        for (int j = 0; j < N; j++)
        {
            A[k][j] /= temp;
            E[k][j] /= temp;
        }
 
        for (int i = k + 1; i < N; i++)
        {
            temp = A[i][k];
 
            for (int j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (int k = N - 1; k > 0; k--)
    {
        for (int i = k - 1; i >= 0; i--)
        {
            temp = A[i][k];
 
            for (int j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            A[i][j] = E[i][j];
 
    for (int i = 0; i < N; i++)
        delete [] E[i];
 
    delete [] E;
}
 
int main()
{
    int N;
 
    std::cout << "Enter N: ";
    std::cin >> N;
 
    double **matrix = new double *[N];
 
    for (int i = 0; i < N; i++)
        matrix[i] = new double [N];
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
        {
            std::cout << "Enter matrix[" << i << "][" << j << "] = ";
            std::cin >> matrix[i][j];
        }
 
    inversion(matrix, N);
 
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
            std::cout << matrix[i][j] << "  ";
 
        std::cout << std::endl;
    }
 
    for (int i = 0; i < N; i++)
        delete [] matrix[i];
 
    delete [] matrix;
 
    std::cin.get();
    return 0;
}

for(int i=0;i<n;i++){                     //заполняется единичная матрица
        for(int j=0;j<n;j++){
            if (i==j){Ed[i][j]=std::complex<double>(1,0);}
            else {Ed[i][j]=std::complex<double>(0,0);}
        }
    }       
        for (int k = 0; k < n; k++)  //обратная матрица
    {
        temp = std::complex<double>(AYT[k][k]);
        for (int j=0;j<n;j++)
        {
            AYT[k][j]=std::complex<double>(AYT[k][j]/temp);
            Ed[k][j]=std::complex<double>(Ed[k][j]/temp);
        }
 
        for (int i=k+1;i<n;i++)
        {
            temp = std::complex<double>(AYT[i][k]);
            for (int j=0;j<n;j++)
            {
                AYT[i][j]=std::complex<double>(AYT[i][j]-(AYT[k][j]*temp));
                Ed[i][j]=std::complex<double>(Ed[i][j]-(Ed[k][j]*temp));
            }
        }
    }
 
    for (int k=n-1;k>0;k--)
    {
        for (int i=k-1;i>=0;i--)
        {
            temp = std::complex<double>(AYT[i][k]);
            for (int j=0;j<n;j++)
            {
                AYT[i][j]=std::complex<double>(AYT[i][j]-(AYT[k][j]*temp));
                Ed[i][j]=std::complex<double>(Ed[i][j]-(Ed[k][j]*temp));
            }
        }
    }
 
    for (int i=0;i<n;i++){
        for (int j=0;j<n;j++){
            AYT[i][j] = std::complex<double>(Ed[i][j]);}
    }
 
       
    cout << "OBRATNAYA!!!!!!!!!!!!!!!!!!!!!!\n\n\n\n";
    for (int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            cout<<AYT[i][j]<<"\t";
        }
        cout<<"\n";
    }