多元一次方程解法 C++

#include<iostream>
#include<math.h>
#include<fstream>
#include<stdlib.h>
using namespace std;
#define MaxElement  4

void print(double (*pArray)[MaxElement + ], int iWidth,int iHigh);
void main(){
    int n,m;
    double a[MaxElement][MaxElement + ] = {
        //3
//         {100, 10, 1, 10},
//         {400, 20, 1, 20},
//         {900, 30, 1, 10}, 

        //
        {,,,,},
        {,,,,-},
        {,,,,-},
        {,,-,,-},
    };//第四列是增广矩阵
    int i,j;
    n = MaxElement;
    cout<<"输入方程组介数：";
    cout<<n<<endl;
    cout<<"输入增广矩阵："<<endl;
    for(i = ; i < n; i++){
        for(j = ; j < n + ;j++){
            cout<<a[i][j]<<"  ";
        }
        cout<<endl;
    }
    for(j = ; j < n; j++)
    {
        double max = ;
        double imax = ;
        for(i = j; i < n; i++)
        {
            if(imax < fabs(a[i][j])){
                imax = fabs(a[i][j]);
                max = a[i][j];//得到各行中所在列最大元素
                m = i;
            }
        }
        if(fabs(a[j][j]) != max)
        {
            double b = ;
            for(int k = j;k < n + ; k++){
                b = a[j][k];
                a[j][k] = a[m][k];
                a[m][k] = b;
            }
        }
        print(a, MaxElement, MaxElement + );
        for(int r = j;r < n + ;r++)
        {
            a[j][r] = a[j][r] / max;//让该行的所在列除以所在列的第一个元素，目的是让首元素为1
        }
        print(a, MaxElement, MaxElement + );
        for(i = j + ;i < n; i++)
        {
            double c = a[i][j];
            if(c == ) continue;
            for(int s = j;s < n + ;s++){
                double tempdata = a[i][s];
                a[i][s] = a[i][s] - a[j][s] * c;//前后行数相减，使下一行或者上一行的首元素为0
                print(a, MaxElement, MaxElement + );
            }
            print(a, MaxElement, MaxElement + );
        }
        print(a, MaxElement, MaxElement + );
    }
    for(i = n - ; i >= ; i--)
    {
        for(j = i + ;j < n; j++)
        {
            double tempData = a[i][j];
            double data1 = a[i][n];
            double data2 = a[j][n];
            a[i][n] = a[i][n] - a[j][n] * a[i][j];
            print(a, MaxElement, MaxElement + );
        }
    }
    print(a, MaxElement, MaxElement + );
    cout<<"方程组的解是："<<endl;
    for(int k = ; k < n; k++){
        cout<<"x"<<k<<" = "<<a[k][n]<<endl;
    }
}
void print(double (*pArray)[MaxElement + ], int iWidth,int iHigh) {
    std::cout<<"Array: "<<"\n";
    for(int i = ; i < iWidth; i++)
    {
        for(int j = ; j < iHigh;j++)
        {
            cout<<pArray[i][j]<<"  ";
        }
        cout<<endl;
    }
}