int SolveSys(CMatrix *A, CVector *b){
const int N = A->GetnEq();
const int LDA = A->GetnVar();
CMatrix* Mat;
Mat = new CMatrix(N,LDA,A->GetMat());
int NRHS(1);
int *IPIV;
int LDB = b->GetSize();
int INFO(-1);
IPIV = new int[N];
INFO = LAPACKE_dgesv(LAPACK_COL_MAJOR,N,NRHS,Mat->GetMat(),LDA,IPIV,b->GetVec(),LDB);
delete Mat;
Mat = NULL;
}
