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

void Invers_Array ( double ** & Matrix, const short int & size )
{
    double a, b, sum = 0;
    double ** Invers_Matrix = new double * [ size ];     /* Единичная матрица */
    for ( int i = 0; i < size; i ++ )
    {
        Invers_Matrix [i] = new double [ size ];
        for ( int j = 0; j < size; j ++ )
            Invers_Matrix [i][j] = 0;
            Invers_Matrix [i][i] = 1;
    }
    for ( int i = 0; i < size; i ++ )     /* Приведём матрицу к треугольной при помощи алгоритма Гаусса */
    {
        a = Matrix [i][i];
        for ( int j = i + 1; j < size; j ++ )
        {
            b = Matrix [j][i];
            for ( int k = 0; k < size; k ++ )
            {
                Matrix [j][k] = Matrix [i][k] * b - Matrix [j][k] * a;
                Invers_Matrix [j][k]= Invers_Matrix [i][k] * b - Invers_Matrix [j][k] * a;
            }
        }
    }
    for ( int i = 0; i < size; i ++ )     /* Обратный ход вычисления элементов обратной матрицы */
    {
        for ( int j = size - 1; j >= 0; j -- )
        {
            for ( int k = size - 1; k > j; k -- )
                sum += Matrix [j][k] * Invers_Matrix [k][i];
                if ( Matrix [j][j] == 0 )
                {
                    for ( int i = 0; i < size; i ++ )
                        delete [] Invers_Matrix [i];
                        delete [] Invers_Matrix;
                }
                Invers_Matrix [j][i] = ( Invers_Matrix [j][i] - sum ) / Matrix [j][j];
        }
    }
    printf ( "\n%s", " Обратная матрица \n\n" );
    Output_Array ( Invers_Matrix, size );
}