高斯消元

//高斯消元 
//可解n元一次方程 
//消成左上到右下的对角线 精度相对于上三角损失更大 
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define eps 1e-7
using namespace std;
int n;
double coe[101][101],num[101];
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
            scanf("%lf",&coe[i][j]);
        scanf("%lf",&num[i]);
    }
    for(int i=1;i<=n;i++)
    {
        int ma=i;
        for(int j=i;j<=n;j++)
            if(coe[ma][i]-coe[j][i]<eps)
                ma=j;
        swap(coe[ma],coe[i]);
        swap(num[ma],num[i]);
        if(fabs(coe[i][i])<eps)
        {
            printf("No Solution
");
            return 0;
        }
        for(int j=i+1;j<=n;j++)
            coe[i][j]/=coe[i][i];
        num[i]/=coe[i][i];
        for(int j=1;j<=n;j++)
        {
            if(j==i) continue;
            for(int k=i+1;k<=n;k++)
                coe[j][k]-=coe[i][k]*coe[j][i];
            num[j]-=num[i]*coe[j][i];
        }
    }
    for(int i=1;i<=n;i++)
        printf("%.2lf
",num[i]);
    return 0;
}