高斯消元求逆矩阵板子。
#include#include #include #include using namespace std;const double EPS=0.00000001;#define N 105int n;double B[N][N*2],A[N][N*2],C[N][N];bool guass_jordan(){ memcpy(B,A,sizeof(A)); for(int i=1;i<=n;++i){ for(int j=n+1;j<=n*2;++j){ if(i==j-n){ B[i][j]=1.0; } else{ B[i][j]=0.0; } } } for(int i=1;i<=n;++i){ int pivot=i; for(int j=i+1;j<=n;++j){ if(fabs(B[j][i])>fabs(B[pivot][i])){ pivot=j; } } swap(B[i],B[pivot]); if(fabs(B[i][i])