#include "TMatrixD.h"
#include "fstream"
#include "iostream"
#include "math.h"
#include "TVector.h"
#include "TArrayD.h"
#include "iomanip"


using namespace std;


int main()
{
    TMatrixD matrixD(3,3);
    TMatrixD matrixD1(3,3);
    TMatrixD matrixD2(3,3);
    TMatrixD matrixD3(3,3);
//    TArrayD data1[3][3]{{1,1,1},{1,2,2},{2,1,2}};


    TArrayD data1(9);
    data1[0]=1;
    data1[1]=1;
    data1[2]=1;
    data1[3]=1;
    data1[4]=2;
    data1[5]=2;
    data1[6]=2;
    data1[7]=1;
    data1[8]=2;

    matrixD.SetMatrixArray(data1.GetArray());
    double DetD=matrixD.Determinant();


    matrixD1=matrixD;
    TMatrixDColumn(matrixD1,0)(0)=10;
    TMatrixDColumn(matrixD1,0)(1)=12;
    TMatrixDColumn(matrixD1,0)(2)=19;
    double DetD1=matrixD1.Determinant();


    matrixD2=matrixD;
    TMatrixDColumn(matrixD2,1)(0)=10;
    TMatrixDColumn(matrixD2,1)(1)=12;
    TMatrixDColumn(matrixD2,1)(2)=19;
    double DetD2=matrixD2.Determinant();



    matrixD3=matrixD;
    TMatrixDColumn(matrixD3,2)(0)=10;
    TMatrixDColumn(matrixD3,2)(1)=12;
    TMatrixDColumn(matrixD3,2)(2)=19;
    double DetD3=matrixD3.Determinant();
    
    double solution1=DetD1/DetD;
    double solution2=DetD2/DetD;
    double solution3=DetD3/DetD;


 
 
    for(int i=0;i<3;i++)
    {
        for(int j=0;j<3;j++)
        {
            cout<<setw(10)<<TMatrixDRow(matrixD,i)(j);
        }
        cout<<endl;
    }
    cout<<"ooooooooooooooooooooooooooooooooooooooooooooooooooooooo"<<endl;
    cout<<setw(15)<<solution1<<setw(15)<<solution2<<setw(15)<<solution3<<endl;
    cout<<"ooooooooooooooooooooooooooooooooooooooooooooooooooooooo"<<endl;



// ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
// Fazekas equations
// ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo

   ifstream ifile1("Aij86_matrix_Fazekas.txt");
   ifstream ifile2("Bi86_vector_Fazekas.txt");
   ofstream ofileFazekas("Solution_Fazekas.txt");
   TArrayD data2(81);
   TArrayD data3(9);



   TMatrixD matrixLarge(9,10);
   for(int i=0;i<9;i++)
   {
       for(int j=0;j<9;j++)
       {
           ifile1>>TMatrixDRow(matrixLarge,i)(j);
       }
   }
   double temp1;
   for(int i=0;i<9;i++)
   {
       ifile2>>TMatrixDRow(matrixLarge,i)(9);
       ifile2>>temp1;
   }
   ofstream ofile1("LargeMatrix.txt");
   for(int i=0;i<9;i++)
   {
       for(int j=0;j<10;j++)
       {
           ofile1<<setw(15)<<scientific<<TMatrixDRow(matrixLarge,i)(j);
       }
       ofile1<<endl;
   }

   

    TMatrixD matrixAFazekas(9,9);
    for(int i=0;i<9;i++)
    {
        for(int j=0;j<9;j++)
        {
            TMatrixDRow(matrixAFazekas,i)(j)=TMatrixDRow(matrixLarge,i)(j);
        }
    }


    
    TMatrixD matrixAFazekas1(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas1,i)(0)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=0)
        TMatrixDColumn(matrixAFazekas1,i)=TMatrixDColumn(matrixLarge,i);
    }



    TMatrixD matrixAFazekas2(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas2,i)(1)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=1)
        TMatrixDColumn(matrixAFazekas2,i)=TMatrixDColumn(matrixLarge,i);
    }




    TMatrixD matrixAFazekas3(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas3,i)(2)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=2)
        TMatrixDColumn(matrixAFazekas3,i)=TMatrixDColumn(matrixLarge,i);
    }




    TMatrixD matrixAFazekas4(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas4,i)(3)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=3)
        TMatrixDColumn(matrixAFazekas4,i)=TMatrixDColumn(matrixLarge,i);
    }




    TMatrixD matrixAFazekas5(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas5,i)(4)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=4)
        TMatrixDColumn(matrixAFazekas5,i)=TMatrixDColumn(matrixLarge,i);
    }




    TMatrixD matrixAFazekas6(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas6,i)(5)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=5)
        TMatrixDColumn(matrixAFazekas6,i)=TMatrixDColumn(matrixLarge,i);
    }




    TMatrixD matrixAFazekas7(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas7,i)(6)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=6)
        TMatrixDColumn(matrixAFazekas7,i)=TMatrixDColumn(matrixLarge,i);
    }





    TMatrixD matrixAFazekas8(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas8,i)(7)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=7)
        TMatrixDColumn(matrixAFazekas8,i)=TMatrixDColumn(matrixLarge,i);
    }




    TMatrixD matrixAFazekas9(9,9);
    for(int i=0;i<9;i++)
    {
        TMatrixDRow(matrixAFazekas9,i)(8)=TMatrixDRow(matrixLarge,i)(9);
    }
    for(int i=0;i<9;i++)
    {
        if(i!=8)
        TMatrixDColumn(matrixAFazekas9,i)=TMatrixDColumn(matrixLarge,i);
    }

    double a0=matrixAFazekas1.Determinant()/matrixAFazekas.Determinant();
    double a1=matrixAFazekas2.Determinant()/matrixAFazekas.Determinant();
    double a2=matrixAFazekas3.Determinant()/matrixAFazekas.Determinant();
    double a3=matrixAFazekas4.Determinant()/matrixAFazekas.Determinant();
    double a4=matrixAFazekas5.Determinant()/matrixAFazekas.Determinant();
    double a5=matrixAFazekas6.Determinant()/matrixAFazekas.Determinant();
    double a6=matrixAFazekas7.Determinant()/matrixAFazekas.Determinant();
    double a7=matrixAFazekas8.Determinant()/matrixAFazekas.Determinant();
    double a8=matrixAFazekas9.Determinant()/matrixAFazekas.Determinant();
    
    cout<<scientific<<a0<<endl
        <<scientific<<a1<<endl
        <<scientific<<a2<<endl
        <<scientific<<a3<<endl
        <<scientific<<a4<<endl
        <<scientific<<a5<<endl
        <<scientific<<a6<<endl
        <<scientific<<a7<<endl
        <<scientific<<a8<<endl;


    ofileFazekas<<scientific<<a0<<endl
        <<scientific<<a1<<endl
        <<scientific<<a2<<endl
        <<scientific<<a3<<endl
        <<scientific<<a4<<endl
        <<scientific<<a5<<endl
        <<scientific<<a6<<endl
        <<scientific<<a7<<endl
        <<scientific<<a8<<endl;


    return 100000000;

}