Написать решение систем линейных уравнений методом Крамера

#include <tchar.h>
#include <iostream>
 
using namespace std;
 
/*функция возведения в определённую степень числа -1*/
int mypow(int x){ 
    int y=1;
    for(int z=1;z<=x;z++){
  y*=-1;
    }
    return y;
}
 
/*класс матрица*/
class matrix{
public:
 
/*Конструктор: все элементы 0, размерность 0 */
 matrix(){
  size=0;
  
  for(int x=0;x<=100;x++)
    for(int y=0;y<=100;y++)
    m[x][y];
 }
/*Рекурсивная функция вычисления определителя */
 int GetDet(){
    int result=0;
  if(size==1){ //Если у матрицы размерность 2(математическая, С`шная 1), то вычисляем определитель
return m[0][0]*m[1][1]-m[0][1]*m[1][0];
  }
  else{
/*Вычисляем определитель путём вычёркивания строк и столбцов. Берём элементы из первой строки. */   
   for(int x=0;x<=size;x++){
result+= mypow(x+2)*m[0][x]*(GetNew(x)).GetDet(); //рекурсия...
   }
  return result;
  }
    }
 
 /*Функия получения новой матрицы из старой, путём вычёркивания первой(нулевой) строки и х`овой колонки*/
 matrix GetNew(int st){
  matrix result;
  int z=0;
  result.size=size-1;
 
  for(int x=1;x<=size;x++){
   for(int y=0;y<=size;y++){
  if(y!=st){
   result.m[x-1][z]=m[x][y];
   z++;
  }
  
   }
   z=0;
  }
  return result;
 }
 
/*Функция вывода на экран квадратной матрицы*/
  void echo(void){
   for(int x=0;x<=size;x++){
    for(int y=0; y<=size;y++){
     cout<<m[x][y]<<" ";
    }
    cout<<"\r\n";
   }
  }
/*Функция вывода на экран матрицы-столбца */
  void echoEx(void){
    for(int x=0;x<=size;x++){
     cout<<m[x][0]<<"\r\n";
    }
   }
 
/*Функция установки размерности матрицы*/
   inline void SetSize(int r){
    size=r;
   }
 
/*Функция установки значения элемента*/
   inline void SetValue(int value, int x , int y){
    m[x][y]=value;
   }
 
   /*Функция замены колонки на матриуц-столбец*/
   inline void SetColumn(int column , matrix source){    
    for(int x=0;x<=size;x++){
     m[x][column]=source.m[x][0];
    }
   }
 
private:
    int m[100][100]; //сама матрица
    int size; //размерность
    
};
int _tmain(int argc, _TCHAR* argv[])
{
matrix first,second,temp; //наши матрицы
int element,size,maindet; //необходимые переменные
bool c=1; //для цикла
cout<<"The Kramer method for the system of linear equations. By Ky3bMu4 (c) 2007.\r\n\r\n";
while(c){
cout<<"\r\nInsert size of matrix:"; 
/*тут небольшая неровность: в математике матрица начинается с 1 столбцы, а в С - с 0-го   */
cin>>size;
first.SetSize(size);
second.SetSize(size);
 
 
/* вводим основную матрицу */
for(int x =0; x<=size; x++){
    for( int y=0; y<=size;y++){
  cout<<"\r\nInsert ["<<x<<";"<<y<<"] element of first matrix:";
  cin>>element;
  first.SetValue(element,x,y);
  }
}
 
/* вводим матрицу-столбец с числами*/
for(int x=0; x<=size; x++){
    cout<<"\r\nInsert "<<x<<" element of second matrix:";
    cin>>element;
    second.SetValue(element,x,0);
}
maindet = first.GetDet(); //вычисляем основной определитель(на него будем делить)
if(maindet!=0){
for(int x=0;x<=size;x++){ //кол-во элементов(переменных в линейном уравнении) == размерности главной матрицы.
//Заменяем колонку, вычисляем определитель , делим на главный определитель, пробуем сократить и выводим 
temp = first;
temp.SetColumn(x,second); //заменянем колонку
cout<<"\r\nElement "<<x<<" = "<<temp.GetDet()<<"/"<<maindet; // выводим на экран  
}
}else{ //облом, главный определитель == 0, Крамер не прокатит :(((
    cout<<"\r\nSorry, but the main determinat of first matrix is 0.Stop.";
}
 
// The End
cout<<"\r\n\r\nFinish!";
cout<<"\r\nDo you want to again calculate new matrix(1/0)?:";
cin>>c;
}
 
    return 0;
}

#include <tchar.h>

#include <iostream>
#include <cmath>
using namespace std;
double det(int n,double **a);
void Kramer(int n,double **aMain,double *freeVar,double *answer);
int main()
{
    int n;
    cout<<"Enter a number of variables ";cin>>n;
    double **aMain;
    aMain= new double* [n];
    double *freeVar;
    freeVar = new double [n];
 
    double *answer = new double [n];
 
 
 
    for(int i=0;i<n;i++)
        aMain[i]=new double [n];
 
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
            {cout<<"A["<<i<<"]["<<j<<"]=";cin>>aMain[i][j];}
 
    cout<<endl<<"Please enter free variables"<<endl;
    for(int i=0;i<n;i++)
    {
        cout<<"b["<<i<<"]=";cin>>freeVar[i];
    }
 
    Kramer(n,aMain,freeVar,answer);
 
 
 
 
    delete [] aMain;
    delete [] freeVar;
    delete [] answer;
    return 0;
 
 
}
 
 
double det(int n, double **mat)
{
    double d=0;
    int c, subi, i, j, subj;
    double **submat;
    submat= new double*[n];
    for(int i=0;i<n;i++)
        submat[i]= new double [n];
 
    if (n == 2)
    {
        return( (mat[0][0] * mat[1][1]) - (mat[1][0] * mat[0][1]));
    }
    else
    {
        for(c = 0; c < n; c++)
        {
            subi = 0;
            for(i = 1; i < n; i++)
            {
                subj = 0;
                for(j = 0; j < n; j++)
                {
                    if (j == c)
                    {
                        continue;
                    }
                    submat[subi][subj] = mat[i][j];
                    subj++;
                }
                subi++;
            }
        d = d + (pow(-1 ,c) * mat[0][c] * det(n - 1 ,submat));
        }
    }
    delete [] submat;
    return d;
}
 
void Kramer(int n,double **aMain,double *freeVar,double *answer)
{
    double *detArray = new double [n];
    double deterMain;
    double *temp= new double [n];
     cout<<"Finding determinant of main matrix"<<endl;
    cout<<"deterMain = ";deterMain=det(n,aMain);cout<<deterMain<<endl;
 
    for(int i=0;i<n;i++)
        {for(int j=0;j<n;j++)
    {
        temp[j]=aMain[i][j];
        aMain[i][j]=freeVar[j];
 
 
    }
    detArray[i]=det(n,aMain);
    for(int k=0;k<n;k++)
    {
        aMain[i][k]=temp[k];
    }
        }
 
        for(int i=0;i<n;i++)
        {answer[i]=detArray[i]/deterMain;
        cout<<answer[i]<<endl;
        }
        delete [] detArray;
        delete [] temp;
 
}