#include "TMatrixD.h"
#include "fstream"
#include "iostream"
#include "math.h"
#include "TVectorD.h"
#include "TArrayD.h"
#include "iomanip"
#include "TDecompLU.h"


using namespace std;


int main()
{
    TMatrixD matrixD(3,3);
    

    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();
    
   
    TVectorD b(0,2,10.,12.,19.,"END");
    TDecompLU lu(matrixD);
    bool ok;
    lu.Solve(b,ok);
    cout<<"bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb"<<endl;
    b.Print();
    cout<<"bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb"<<endl;


    


    return 100000000;

}