c++实现多元一次方程组

逆矩阵实现 实现三元及以上应该都没有问题。（之前写的了 想了想还是发上来）

#include <iostream>
using namespace std;
#define N 10
#define MAX 10
double X[MAX][MAX],Y[MAX][MAX],RESULT[MAX][MAX];
double des[N][N] ={};
//按第一章展开式计算|A|
double getA(double arcs[N][N],int n)
{
    if(n==1)
    {
        return arcs[0][0];
    }
    double ans = 0;
    double temp[N][N]={0.0};
    int i,j,k;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n-1;j++)
        {
            for(k=0;k<n-1;k++)
            {
                temp[j][k] = arcs[j+1][(k>=i)?k+1:k];
            }
         }
        double t = getA(temp,n-1);
        if(i%2==0)
        {
            ans += arcs[0][i]*t;
        }
        else
        {
            ans -=  arcs[0][i]*t;
        }
     }
     return ans;
}
//计算每一行每一列的每个元素所对应的余子式，组成A*
 void  getAStart(double arcs[N][N],int n,double ans[N][N])
 {
    if(n==1)
    {
        ans[0][0] = 1;
        return;
    }
    int i,j,k,t;
    double temp[N][N];
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            for(k=0;k<n-1;k++){
                for(t=0;t<n-1;t++)
                {
                    temp[k][t] = arcs[k>=i?k+1:k][t>=j?t+1:t];
                }
            }
            
            ans[j][i]  =  getA(temp,n-1);
            if((i+j)%2 == 1)
            {
                ans[j][i] = - ans[j][i];
            }
        }
     }
}
//得到给定矩阵src的逆矩阵保存到des中。
void GetMatrixInverse(double src[N][N],int n,double des[N][N])
 {
    double flag=getA(src,n);
    double t[N][N];
    if(flag==0)
        cout<<"error"<<endl;
    else
    {
        getAStart(src,n,t);
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                des[i][j]=t[i][j]/flag;
            }
        }
    }
}
void matrix()
{
    int a,b,c,i,j,k;
    //思路：先求x的逆矩阵des 然后des与Y相乘 得到result结果
    cout<<"请输入三个整数(a为矩阵A的行、b为矩阵A的列、c为矩阵B的列)："<<endl;
    cin>>a>>b>>c;
    cout<<"请输入X矩阵："<<endl;
    for(i = 0;i<a;i++){
        for(j=0;j<b;j++){
            cin>>X[i][j];
        }
    }
    cout<<"请输入Y矩阵："<<endl;
    for(i = 0;i<b;i++)
    {
        for(j=0;j<c;j++)
        {
            cin>>Y[i][j];
        }
    }
    cout<<endl;
    cout<<"A^-1 = "<<endl;
    GetMatrixInverse(X, a, des);
    //这里求到逆矩阵之后=>des
    for(int i=0;i<a;i++)
    {
        for(int j=0;j<b;j++)
        {
            cout<<des[i][j]<<"  ";
        }
        cout<<endl;
    }
//用逆矩阵与结果矩阵相乘
    for(i = 0;i<a;i++)
    {
        for(j = 0;j<c;j++)
        {
            RESULT[i][j] = 0;
            for(k=0;k<b;k++)
            {
                RESULT[i][j] += des[i][k]*Y[k][j];
            }
        }
    }
    cout<<"result:"<<endl;
    for(int i=0;i<a;i++)
    {
        for(int j=0;j<c;j++)
        {
            cout<<RESULT[i][j]<<"  ";
        }
        cout<<endl;
    }
}

int main()
{
    matrix();
}

运行结果： 

若有问题 欢迎提问👏