2019-05-19 线性方程组-高斯消元法

高斯消元法：（列主元高斯消元法）
1.正常的情况下，我们将矩阵化成上三角，然后化成行最简形式，这里算法略微不同；
2.这里处理的时候，每循环一次，消除一个未知数，没有意义的计算直接省略；

const double EPS=1e-8;
typedef vectorvec;
typedef vectormat;

//计算Ax=B
//当方程组无解或者无穷解时候,返回一个长度为0的数组 
vec gauss_jordan(const mat&A,const vec&b)
{
    int n=A.size();
    mat B(n,vec(n+1));//增广矩阵，只能对增广矩阵做初等行变换，得到的新的方程组同解
    for(int i=0;iabs(B[pivot][i])) pivot=j;
        }
        swap(B[i],B[pivot]);
        //无解或者无穷解 
        if(abs(B[i][i])