Gauss Elimination for NxM matrix
问题 /* Program to demonstrate gaussian <strong class="highlight">elimination</strong> on a set of linear simultaneous equations */ #include <iostream> #include <cmath> #include <vector> using namespace std; const double eps = 1.e-15; /*Preliminary pivoting strategy Pivoting function */ double pivot(vector<vector<double> > &a, vector<double> &b, int i) { int n = a.size(); int j=i; double t=0; for(int k=i; k<n; k+=1) { double aki = fabs(a[k][i]); if(aki>t) { t=aki; j=k; } } if(j>i) { double dummy;