高斯消元法c++程序

//这是一个高斯消元法的函数，不足之处是只能对方阵阵进行操作
//操作规则：先确定方阵的行数，再在宏定义中写明N的值
//能有效进行高斯消元。

#include
#define N 5
using namespace std;

//////////////////////定义一个输入函数/////////////////////////////

int input(double array[N][N])
{
// double array[N][N];
int i,j;
for (i=0;ifor (j=0;jcin>>array[i][j];
return 0;
}

///////////////////////一个变幻函数////////////////////////////////

void vary(double array[N][N]) //变换成上三角
{
int row1,row2,col;
//row1表示行，row2表示在row1的基础上的下一行，col表示相应的列
double m;
for(row1=0;row1{
for(row2=row1+1;row2{
m=array[row2][row1]/array[row1][row1];
for(col=row1;col{
array[row2][col]=array[row2][col]-array[row1][col]*m;
}
}
}
}
//////////////////////////主函数///////////////////////////////////
int main()
{
double a[N][N];
cout<<"\n\n";
cout<<"\t\t\t\t请输入一个 "<input(a);
vary(a);
for(int i=0;i{
for(int j=0;j{
cout<// printf("%d",a[i][j]);
}
cout<}

}