#include#include #include #include using namespace std; double a[110][110]; int n; const double esp=1e-6; int gauss() { int c,r; for(r=0,c=0;c fabs(a[t][c])) t=i; if(fabs(a[t][c]) =c;i--) a[r][i]/=a[r][c]; //1 for(int i=r+1;i =c;j--) //倒叙 不然a[i][c]会变成0 a[i][j]-=a[r][j]*a[i][c]; // 进行运算 } r+=1; } if(r esp) return 2;//无解 return 1;//无穷解 } //解出唯一解 for(int i=n-1;i>=0;i--) { for(int j=i+1;j >n; for(int i=0;i >a[i][j]; } int t=gauss(); if(t==0) { for(int i=0;i



