高斯消元法

学过至少一元二次方程解法的同学们肯定知道，最常用的求根方法是加减消元法。

所谓加减消元法，就是通过方程的等式性质，将某一未知数的系数同一化，再通过两个等式的加减消去一个元，从而将多元方程不断降元，从而解出一个未知数，不断代回，从而求出所有根。

此加减消元法就是高斯消元法。

所以程序思路也尽可知。

以下直接附上代码

代码转自https://blog.csdn.net/qq_33929112/article/details/52233470

#include <iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const double EPS=1e-8;
#define maxn 105
double A[maxn][maxn];
double B[maxn][maxn+1];
double x[maxn];
double b[maxn];
int n;
//求解Ax=b，其中A是方阵
void gauss_jordan()
{
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            B[i][j]=A[i][j];
        }
    }
    for(int i=0;i<n;i++){
        B[i][n]=b[i];
    }
    for(int i=0;i<n;i++){
        //把正在处理的未知数系数的绝对值最大的式子换到第i行
        int pivot=i;
        for(int j=i;j<n;j++){
            if(abs(B[j][i]>abs(B[pivot][i]))) pivot=j;
        }
        for(int j=0;j<=n;j++){
            swap(B[i][j],B[pivot][j]);
        }
        //无解或者无穷多解
        if(abs(B[i][i])<EPS) return;
        //把正在处理的未知数的系数变成1
        for(int j=i+1;j<=n;j++){
            B[i][j]/=B[i][i];
        }
        for(int j=0;j<n;j++){
            if(i!=j){
                //从第j个式子中消去第i个未知数
                for(int k=i+1;k<=n;k++){
                    B[j][k]-=B[j][i]*B[i][k];
                }
            }
        }
    }
    for(int i=0;i<n;i++){
        x[i]=B[i][n];
    }
    return;
}
int main()
{
    cin>>n;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            cin>>A[i][j];
        }
    }
    for(int i=0;i<n;i++){
        cin>>b[i];
    }
    gauss_jordan();
    for(int i=0;i<n;i++){
        cout<<x[i]<<" ";
    }
    printf("\n");
    return 0;
}