模板 - 高斯约旦消元法

namespace Gauss_Jordan_Elimination{
    const int MAXN=105;
    double a[MAXN][MAXN],del;
    double ans[MAXN];

    const double eps=1e-10;

    bool gauss_jordan(int n){
        for(int i=1;i<=n;i++){
            int mx=i;
            for(int j=i+1;j<=n;j++)
                if(fabs(a[j][i])>fabs(a[mx][i]))
                    mx=j;
            for(int j=1;j<=n+1;j++)
                swap(a[i][j],a[mx][j]);
            if(fabs(a[i][i])<eps)
                return false;
            for(int j=1;j<=n+1;j++){
                if(j!=i){
                    double tmp=a[j][i]/a[i][i];
                    for(int k=i+1;k<=n+1;k++){
                        a[j][k]-=a[i][k]*tmp;
                    }
                }
            }
        }

        for(int i=1;i<=n;i++)
            ans[i]=a[i][n+1]/a[i][i];

        return true;
    }
}


using namespace Gauss_Jordan_Elimination;