Система линейных уравнений. Крамер C++.

#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <tchar.h>
using namespace std;
 
 
//Функция возведения числа в нужную степень, для вычисления минора. [url]http://www.mathelp.net/MA1.htm[/url]
int mypow(int x)
{ 
    int y=1;
    for(int z=1;z<=x;z++)
    {
        y*=-1;
    }
    return y;
}
 
//Функия получения новой матрицы из старой, путём вычёркивания первой (математически) нулевой(Си-шно) строки и х`овой колонки.
int GetNew(double m[100][100],int msize,double result[100][100],int st)
{
    int rsize;
    int z=0;
    rsize=msize-1;
    for(int x=1;x<=msize;x++)
        {
         for(int y=0;y<=msize;y++)
         {
             if(y!=st)
             {
                 result[x-1][z]=m[x][y]; // z здесь затем, чтобы не попаст в у=st x-номер строки
                 z++;
             }
         }
         z=0;
    }
    return rsize;
}
 
//Рекурсивная функция вычисления определителя:
double GetDet(double m[100][100], int msize)
{
    double result=0;
    if(msize==1)
        {
            return result=(m[0][0]*m[1][1]-m[0][1]*m[1][0]);
        }
        else //Вычисляем определитель путём вычёркивания строк и столбцов. Берём элементы из первой строки.
        {
            double a[100][100];
            int asize;
            for(int x=0;x<=msize;x++)
                {
                    asize=GetNew(m,msize,a,x);
                    result+= mypow(x+1)*a[0][x]*(GetDet(a,asize));
                }
            return result;
        }
}
 
 
 
//Функция присвоения матрицы
int EqlMatrix(double m[100][100], int msize, double s[100][100])
{
    for(int x=0;x<msize;x++)
    {
        for(int y=0;y<msize;y++)
        {
            m[x][y]=s[x][y];
        }
    }
    return msize;
}
 
//Функция замены колонки на матрицу-столбец
void SetColumn(double m[100][100], int msize, int column , double source[100][100])
{   
    for(int x=0;x<=msize;x++)
    {
        m[x][column]=source[x][0];
    }
}
 
int main(int argc, char* argv[])
{   
    int flag, size,flag2;
    double  first[100][100],second[100][100],temp[100][100];
    int firstsize,secondsize,tempsize;
    double h, element1, element2, maindet;
 
    bool c=1;
    cout<<"Metod Kramera dlya resheniya sistem lineynih uravneniy.\n\n";
    while(c)
    {
        //Ввод размера
        cout<<"Vvedite razmer matrici (Otchet razmera vedetcya s ""0"", ot 1 do 7):";
        cin>> size;
 
        //Отсев размера
        do
        {
            if((size>=1)&&(size<=7))
            {
                cout<< "\n";
                break;
                flag=1;
            }  
            else
            {
                cout<< "Vvedite x ot 1 do 7, skazali ghe! \n";
                flag=0;
                cin>> size;
            }
        }
        while (flag=1);
        
        //Заполнение квадратной матрицы из коэффициентов при переменных
        firstsize=size;
                for(int x=0;x<=size;x++)
                        for(int y=0;y<=size;y++)
                        {
                                cout<<"\nVvedite ["<<x<<";"<<y<<"] element pervoy matrici:";
                                cin>>element1;
                                first[x][y]=element1;
                        }
 
        //Заполнение столбца свободных элементов
        secondsize=size;
        for(int x=0; x<=size ; x++)
        {
            cout<<"\nVvedite "<<x<<" element vtotoy matrici:";
            cin>>element2;
            second[x][0]=element2;
        }
 
        //Вывод системы уравнений
        for(int x=0;x<=size;x++)
        {
            for(int y=0;y<=size;y++)
            {
                cout<<first[x][y]<<" ";
            }
            cout<<" |"<<second[x][0]<<"\n";
        }
 
        //вычисляем основной определитель
        maindet = GetDet(first, firstsize);
        if(maindet!=0)
        {
            for(int x=0;x<=size;x++)
            {
                //Заменяем колонку, вычисляем определитель , делим на главный определитель и выводим 
                tempsize = EqlMatrix(first, firstsize,temp);
                SetColumn(temp/*куда*/, tempsize,x/*какую*/,second/*откуда*/);
                h=GetDet(temp,tempsize)/maindet;
                cout<<"\nElement "<<x<<" = "<<GetDet(temp,tempsize)<<"/"<<maindet<<"="<<h<<"";
            }
        }
        else
        {
            cout<<"\n \n Glavniy determenant osnovnoy matrici raven 0,\n Metod Kramera ne radotaet i mi moghem otvetit' lish na vopros: skolko resheniy.\n";
            for(int x=0;x<=size;x++)
            {
                tempsize = EqlMatrix(first, firstsize,temp);
                SetColumn(temp,tempsize,x,second);
                GetDet(temp,tempsize);
                if ((GetDet(temp,tempsize))!=0) 
                {
                    flag2=0;
                    break;
                }
                else
                {
                    flag2=1;
                }
            }
            if(flag2==1)
            {
                cout<<" Beskonechno mnogo resheniy";
            }
            else
            {
                cout<<" Resheniy net";
            }
        }
        cout<<"\n \n Hotite reshit' eshe 1 sistemu(1-yes. 0- no)?:";
        cin>>c;
    }
    return 0;
    cin.get();
    cin.get();
}