首页 > 其他 > 详细

高斯消元

时间:2018-07-08 17:33:29      阅读:124      评论:0      收藏:0      [点我收藏+]
//高斯消元 
//可解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\n");
            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\n",num[i]);
    return 0;
}

 

高斯消元

原文:https://www.cnblogs.com/water-radish/p/9280606.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!