使用列主元消元法解方程组c语言代码.docx

使用列主元消元法解方程组 c语言代码 # includestdio.h # includemath.h # define N 3 main(){ int i,j,k,h,s,m,n,z; float max(float *y); float A[N][N+1],B[N],ma,t,M[N]={0},X; float x[N]; /*输入系数矩阵 */ for(i=0;iN;i++) for(j=0;jN+1;j++) scanf(%f,A[i][j]); printf(向命令窗口输入的曾广矩阵是：\n); for(i=0;iN;i++) { for(j=0;jN+1;j++) printf(%2.1f ,A[i][j]); printf(\n); } /*正消过程*/ for (i=0;iN-1;i++) { /*找出每行最大的值 */ for (j=i;jN;j++) B[j-i]=A[j][i]; ma=max(B); for (s=0;sN;s++) B[s]=0; /* 交换最大值行和当前行 */ if(ma!=0){ for (k=i;kN;k++) if (A[i][i]!=ma A[k][i]==ma ma!=0) { for (h=0;hN+1;h++) { t=A[i][h]; A[i][h]=A[k][h]; A[k][h]=t; } } } else { printf(A是非奇异矩阵！\n); break; } /* 消去过程 */ for (m=i+1;mN;m++) M[m]=A[m][i]/A[i][i]; for (n=i+1;nN;n++) { for (z=0;zN+1;z++) A[n][z]=A[n][z]-M[n]*A[i][z]; } for (m=0;mN;m++) M[m]=0; } /*输出三角矩阵 */ printf(通过高斯列主元消去法化成的上三角曾广矩阵是：\n); for(i=0;iN;i++) { for(j=0;jN+1;j++) printf(%2.1f ,A[i][j]); printf(\n); } /* 回带过程 */ x[N-1]=A[N-1][N]/A[N-1][N-1]; for(i=N-2;i0 || i==0;i--) { for(j=N-1;ji;j--) X=X+A[i][j]*x[j]; x[i]=(A[i][N]-X)/A[i][i]; X=0; } /* 将方程的解输出 */ printf(该方程组的解依次为：\n); for(i=0;iN;i++){ printf(x[%d]=%2.1f ,i+1,x[i]); } } /* 计算最大值函数 */ float max(float *y){ int w; float max=y[0]; for (w=0;wN;w++) if (fabs(y[w])fabs(max)) max=y[w]; return max; }