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; }