Решение СЛАУ методом Крамера

#include <time.h>
#include <iostream>
#include <math.h>
#include <stdlib.h>
 
using namespace std;
//extern int X;
 
int **FIND_Minor( int **matrix, int size, int str, int stb )
 
{
 
    
    int **minor = new int *[ size - 1 ];
 
    //m_str, m_stb - коэффициенты, с помощью которых строится сам минор
    int m_str = 0;
    int m_stb;
    for( int i = 0; i < size; i++ )
    {
        //Элементы из строки с номером str не записываются в минор
        if( i != str )
        {
            m_stb = 0;
            minor[ m_str ] = new int[ size - 1 ];
            for( int j = 0; j < size; j++ )
            {
                //Элементы из столбца с номером stb не записываются в минор
                if( j != stb )
                {
                    minor[ m_str ][ m_stb ] = matrix[ i ][ j ];
                    m_stb++;
                }
            }
            m_str++;
        }
    }
 
    return minor;
}
 
long long FIND_Det( int **matrix, int size )
{
 
    //Проверка условия выхода из рекурсии
    if( size == 2 )
    {
        return matrix[ 0 ][ 0 ] * matrix[ 1 ][ 1 ] - matrix[ 0 ][ 1 ] * matrix[ 1 ][ 0 ];
    }
    else if( size == 1 )
    {
        return matrix[ 0 ][ 0 ];
    }
 
   //Момент оптимизации - ищем строку, в которой больше всего нулей, чтобы по ней находить определитель
    //Работает в связи с условием разложения по строке, которое пояснено ниже
    int max_0, str_0 = 0, max_0_b = 0;
    for( int i = 0; i < size; i++ )
    {
        max_0 = 0;
        for( int j = 0; j < size; j++ )
        {
            if( matrix[ i ][ j ] == 0 )
            {
                max_0++;
            }
        }
        if( max_0 > max_0_b )
        {
            str_0 = i;
            max_0_b = max_0;
        }
    }
 
    //Само рекурсивное разложение по строке
    int result = 0;
 
    for( int j = 0; j < size; j++ )
    {
        if( matrix[ str_0 ][ j ] )
        {
            result += matrix[ str_0 ][ j ] * ( unsigned )pow( -1.f, ( unsigned )j ) * FIND_Det( FIND_Minor( matrix, size, str_0, j ), size - 1 );
 
        }
    }
 
    return result;
}
 
 
 //  int column_of_coefficients[ size ];
 //   for( int i = 0; i < size; i++ )
   // {
     //    column_of_coe for( int j = 0; j < size; j++ )w
 
 
 
 
 
 
 
int main()
{
    srand( ( unsigned )time(0) );
 
    float begin,end;
 
 
    int size;
    int det[ size ];
 
    cout << "Input size: "<< endl;
    cin >> size;
    int **matrix = new int*[ size ];
 
    char sw;
    cout << "\n1Avtomaticly plenum - \n2Manual plenum - \n";
    cin >> sw;
 
    switch( sw )
    {
    case '1':
        cout << "\nMatrix:\n";
        for( int i = 0; i < size; i++ )
        {
            matrix[ i ] = new int[ size ];
            for( int j = 0; j < size; j++ )
            {
                matrix[ i ][ j ] = rand()%10;
                cout << matrix[ i ][ j ] << " ";
            }
            cout << "\n";
        }
        break;
 
    case '2':
        cout << "\nYour Matrix:\n";
        for( int i = 0; i < size; i++ )
        {
            matrix[ i ] = new int[ size ];
            for( int j = 0; j < size; j++ )
            {
                cin >> matrix[ i ][ j ];
            }
        }
        break;
 
    default:
        cout << "\n";
        return 0;
    }
 
// Создаю двумерный массив и оттождествляю его с данной матрицей, чтобы потом использовать в замене столбцов.
    int col_of_coeffic[ size ][ size ];
    for ( int i = 0; i < size; i ++)
    {
      for ( int j = 0; j < size; j ++)
      {
          col_of_coeffic[ i ][ j ] = matrix[ i ][ j ] ;
        }
    }
// Создаю рандомный столбец решений СЛАУ, его я буду переставлять, дабы находить определители по методу Крамера.
  int column_of_coefficients[ size ];
  for( int i = 0; i < size; i++ )
  {
       column_of_coefficients[ i ] = rand()%10;
       cout <<  column_of_coefficients[ i ] << endl;
  }
        for ( int j = 0; j < size; j ++)
        {
 
         for ( int i = 0; i < size; i++ )
          {
             for ( int j = 0; j < size; j ++)
             {
// Собственно, вот и сама перестановка.
                 matrix[ i ][ j ] = column_of_coefficients[ i ] ;
             }
 
// Нахожу и-тый определитель, затем все повторяю, заменив матрицу на прежнюю.
         det [ i ] = FIND_Det( matrix, size ) ;
             matrix[ i ][ j ] = col_of_coeffic[ i ][ j ];
        }
      }
 
 // cout << "column solutions : " << endl;
     int X[ size ];
// Самая последняя часть - ввожу пустой одномерный вертикальный массив, заполняю его по формуле - Δi / Δ 
 
       for ( int i = 0; i < size; i++ )
       {
          X[ i ] = det [ i ] / FIND_Det( matrix, size );
       }
    for ( int i = 0; i < size; i++ )
       {
          cout << X[ i ] << endl;
        }
// Вывожу его и определитель исходной рандомной матрицы.
    cout << "\nDeterminate= ";
    begin = clock();
    cout << FIND_Det( matrix, size ) << endl;
    end = clock();
 
    cout << "\nTime  - ";
    cout << ( end - begin ) / 1000 << " seconds" << endl;
 
    return 0 ;
}