#include "stdafx.h"
#include <iostream>
#include <math.h>
using namespace std;
void main()
{
int i,j,k,m;
float max,sum,y,t;
float q;
const int n=3; // 声明一个常量变量，用来储存数组的长度
float a[n][n+1],x[n];
cout<<"请输入 a[i][j]: ";
for(i=0;i<=n-1;i++)
for(j=0;j<=n;j++)
cin>>a[i][j];
//输出a[i][j]
for(i=0;i<=n-1;i++)
{
for(j=0;j<=n;j++)
cout<<a[i][j]<<" ";
cout<<endl;
}
cout<<endl;
//把a[i][j]化为上三角矩阵
for(k=0;k<=n-2;k++)
{
for(max=0,i=k;i<=n-1;i++) // 求列主元
if(fabs(a[i][k])>max)
{
max=fabs(a[i][k]);
m=i;
}
if(m!=k) // 交换主元的位置
{
for(j=0;j<=n;j++)
{
t=a[k][j];
a[k][j]=a[m][j];
a[m][j]=t;
}
}
if(a[k][k]!=0) // 消元
{
for(i=k+1;i<=n-1;i++)
{
q=a[i][k];
for(j=k;j<=n;j++)
a[i][j]=a[i][j]-q*(a[k][j]/a[k][k]);
}
}
for(i=0;i<=n-1;i++) // 输出a[i][j]
{
for(j=0;j<=n;j++)
cout<<a[i][j]<<" ";
cout<<endl;
}
cout<<endl;
}
// 回代求x的解
for(i=n-1;i>=0;i--)
{
for(j=n,sum=a[i][j];j>i+1;j--)
sum=sum-a[i][j-1]*x[j-1];
x[i]=sum/a[i][i];
}
// 输出上三角矩阵
cout<<"变换后的矩阵为："<<endl;
for(i=0;i<=n-1;i++)
{
for(j=0;j<=n;j++)
cout<<a[i][j]<<" ";
cout<<endl;
}
cout<<endl;
// 输出方程的解
cout<<" 方程的解为："<<"x1 "<<"x2 "<<"......."<<endl;
for(i=0;i<=n-1;i++)
cout<<" "<<x[i];
cout<<endl;
}

不好意思，上传错了。这应该是求解 N 个未知数的线性方程组的程序。