用c++编高斯消元法解m行n列的线性方程组

好的，下面是用C++编写的高斯消元法求解m行n列的线性方程组的程序：

#include <iostream>
#include <cmath>

using namespace std;

const double eps = 1e-6;

int Gauss(double a[][N], int m, int n, double ans[]) {
    int i, j, k, r;
    for (k = 0; k < min(m, n); k++) {
        r = k;
        for (i = k + 1; i < m; i++) {
            if (fabs(a[i][k]) > fabs(a[r][k]))
                r = i;
        }
        if (fabs(a[r][k]) < eps)
            return 0;
        if (r != k)
            for (j = k; j <= n; j++)
                swap(a[k][j], a[r][j]);
        for (j = k + 1; j <= n; j++)
            a[k][j] /= a[k][k];
        a[k][k] = 1;
        for (i = 0; i < m; i++) {
            if (i != k)
                for (j = k + 1; j <= n; j++)
                    a[i][j] -= a[i][k] * a[k][j];
            a[i][k] = 0;
        }
    }
    for (i = 0; i < min(m, n); i++)
        ans[i] = a[i][n];
    return 1;
}

其中，a为系数矩阵，m为方程组的行数，n为方程组的列数，ans为解向量。

需要注意的是，代码中的N需要根据实际情况进行修改。此外，代码中的eps表示精度，可以根据需要进行调整。

使用时，只需要将系数矩阵a和解向量ans传入函数Gauss中即可，返回值为1表示方程组有解，ans中存储的即为解向量；返回值为0表示方程组无解。