解线性方程组.docVIP

//线性方程组列主元高斯消元法求解 #include stdafx.h #includememory #if _DEBUG #includeiostream //测试时使用 #endif using namespace std; //绝对值函数 __inline double _abs(double v) { return v 0 ? -v : v; } #if _DEBUG void test(int n, const double *pd) { int i, j, k = 0; for (i = 0; i n; i++) { for (j = 0; j = n; j++) { cout pd[k] ; k++; } cout endl; } } #endif //线性方程组列主元高斯消元法 //n 方程元数；pCoef 系数，必须以行主序方式存放的二维数组； //pOut 长度为 n 的一维数组(调用者负责维护)，用于输出数据 //返回值：0 成功，-1 无解，1 申请内存失败, 2 不定解。 int GaussJordanElimination(int n, const double *pCoef, double *pOut) { double *pcf; int rows = n, columns = n +