#include <iostream>
#include <fstream>

using namespace std;

int main()
{
    unsigned rozmiarmacierzy;
    double** A;
    double* b;
    double* p;

    while(true)
    {
        int wybor;
        cout << "\nMENU:" <<endl;
        cout <<"1 - wczytanie danych" << endl;
        cout <<"2 - metoda eliminacji Gaussa bez wyboru elementu"<< endl;
        cout <<"3 - metoda eliminacji Gaussa z wyborem elementu w wierszach"<<endl;
        cout <<"4 - metoda eliminacji Gaussa z wyborem elementu w kolumnach"<<endl;
        cout <<"5 - metoda eliminacji Gaussa z pelnym wyborem elementu"<< endl;
        cout <<"6 - wyjscie z programu"<<endl;

        cin>>wybor;

        switch(wybor)
        {
        case 1:

        {

            ifstream plik_zrodlowy("test.csv");
            if (!plik_zrodlowy.is_open())
            {
                cout <<"Plik nie zosstal otwarty!"<<endl;
                return -1;
            }

            plik_zrodlowy >> rozmiarmacierzy;

            A = new double*[rozmiarmacierzy];

            A[0] = new double[rozmiarmacierzy*rozmiarmacierzy];

            for(unsigned i = 1; i< rozmiarmacierzy; i++)
                A[i] = A[i-1] + rozmiarmacierzy;

            b = new double[rozmiarmacierzy];

            char pol_kolumny;

            for (unsigned i = 0; i < rozmiarmacierzy+1; i++)
                plik_zrodlowy >> pol_kolumny;

            for (unsigned i = 0; i < rozmiarmacierzy; i++)
            {
                for (unsigned j = 0; j < rozmiarmacierzy; j++)
                {
                    plik_zrodlowy >> A[i][j];
                    plik_zrodlowy >> pol_kolumny;
                }
                plik_zrodlowy >> pol_kolumny;
                plik_zrodlowy >> b[i];
            }
            plik_zrodlowy.close();

            cout << "\nPobrano dane"<<endl;
            for(int i=0; i<rozmiarmacierzy; i++)
            {
                for(int j=0; j<rozmiarmacierzy; j++)
                {
                    cout <<"\t"<< A[i][j];
                }
                cout<<endl;
            }
            cout << rozmiarmacierzy<<endl;

            break;
        }


        case 2:

            int i,j;
            p = new double[rozmiarmacierzy];

            for(int zmienna=0; zmienna < rozmiarmacierzy; zmienna++)
            {
                for(i=zmienna; i<rozmiarmacierzy; i++)
                {
                    for(j=zmienna; j<rozmiarmacierzy; j++)
                    {
                          p[i]=A[i][i];
                        if(i>j)
                        {
                            p[i]=A[i][j];
                            p[i]=p[i]/p[j];
                            cout<<p[i]<<endl;
                            A[i][j]=A[i][j]-A[i-1][j]*p[i];

                            A[i][j]=0;
                            j++;

                        }

                    }
                }
            }


        for(int i=0; i<rozmiarmacierzy; i++)
            {
                for(int j=0; j<rozmiarmacierzy; j++)
                {
                    cout <<"\t"<< A[i][j];
                }
                cout<<endl;
            }
            break;
        case 3:
                break;
        case 4:
                break;
        case 5:
                break;
        case 6:
                return 0;
        }
    }
}