高斯消去法.docVIP

测试数据2:方程无解的情况 三、程序清单 新的高斯消去法.cpp #includeiostream.h #includemath.h const int s=100; void shuru(double a[s][s],double b[s],int n); //函数输入原型声明 void xiaoyuan(double a[s][s],double b[s],int n); //函数消元原型声明 void huidai(double a[s][s],double b[s],double x[s],int n); //函数回带原型声明 void shuchu(double x[s],int n,double a[s][s]); //函数输出原型声明 void main() //主函数的定义 { double a[s][s],b[s],x[s]; //变量数组的定义 int n; cout输入方程组的阶数:\n; cinn; shuru(a,b,n); //调用函数输入 xiaoyuan(a,b,n); //调用函数消元 huidai(a,b,x,n); //调用函数回带 shuchu(x,n,a); //调用函数输出 } void shuru(double a[s][s],double b[s],int n) //函数输入的定义 { int i,j; cout输入系数矩阵:\n; i=0; while(in) { j=0; while(jn) { cina[i][j]; //输入矩阵到系数矩阵 j++; } i++; } cout输入常数矩阵:\n; i=0; while(in) { cinb[i]; //输入矩阵到常数矩阵 i++; } } void xiaoyuan(double a[s][s],double b[s],int n) //函数消元的定义 { int i,j,l; double p,q,max; for(l=0;ln-1;l++) { max=fabs(a[l][l]); //调用绝对值函数，选出绝对值最大者 j=l; for(i=l;in;i++) if(maxfabs(a[i][l])) { max=fabs(a[i][l]); j=i; } for(i=l;in;i++) //交换系数行 { p=a[l][i]; a[l][i]=a[j][i]; a[j][i]=p; } p=b[l]; //交换常数行 b[l]=b[j]; b[j]=p; for(i=l+1;in;i++) //消元计算 { q=a[i][l]/a[l][l]; for(j=l;jn;j++) a[i][j]=double(a[i][j]-q*a[l][j]); //消元后的系数矩阵 b[i]-=q*b[l]; //消元后的常数矩阵 } } } void huidai(double a[s][s],double b[s],double x[s],int n) //函数回带的定义 { double p; int i,j; x[n-1]=b[n-1]/a[n-1][n-1]; i=n-2; while(i=0) { p=0; j=i+1;