gauss.rar_gauss column c++_列主元 高斯消元_高斯消元法

#include<stdio.h> #include<math.h> void main() { float a[10][10],b[10],d = 0,t,sum = 0; int k = 1,l = 0,i,j = 0,m,n; printf("请输入方程组组数："); scanf_s("%d",&n); printf("请输入方程组未知数个数："); scanf_s("%d",&m); printf("请输入方程组各项系数：\n"); for(i=1;i<=n;i++) { printf("第%d组方程系数：",i); for(j=1;j<=m;j++) scanf_s("%f",&a[i][j]); scanf_s("%f",&b[i]); } while(k<=n) {d=a[k][k];l=k; for(i=k+1;i<=n;i++) if(fabs(a[i][k])>fabs(d)) { d=a[i][k];l=i; } if(l!=k) { for(j=k;j<=n;j++) { t=a[l][j];a[l][j]=a[k][j];a[k][j]=t;} t=b[k];b[k]=b[l];b[l]=t; } for(j=k+1;j<=n;j++) a[k][j]=a[k][j]/a[k][k]; b[k]=b[k]/a[k][k]; for(i=k+1;i<=n;i++) { for(j=k+1;j<=n;j++) a[i][j]=a[i][j]-a[i][k]*a[k][j]; j=1; b[i]=b[i]-a[i][k]*b[k]; } k++; } for(i=n-1;i>=1;i--) { sum=0; for(j=i+1;j<=n;j++) sum=sum+a[i][j]*b[j]; b[i]=b[i]-sum; } printf("结果：\n"); for(i=1;i<=n;i++) printf("x%d=%4.2f\n",i,b[i]); system("pause"); }