Question

求个用高斯消元法解方程的C语言程序。（要在TC理能运行）谢谢了，大神帮忙啊

求个用高斯消元法解方程的C语言程序。一定要在TC里能运行才可以。

Answer 1

#include "Stdio.h" #include "Conio.h" #include "Stdlib.h" /*L是矩阵的行减1,从程序上看是最外层循环的次数 N 对应矩阵的行数，M对应矩阵的列数 可以通过改变L、N、M来控制矩的阶数 */ #define L 3 #define N 4 #define M 5 void gauss(double a[N][M],double x[N]) { int i,j,l,n,m,k=0; double temp[N]; /*第一个do-while是将增广矩阵消成上三角形式*/ do{n=0; for(l=k;l<L;l++)temp[n++]=a[l+1][k]/a[k][k]; for(m=0,i=k;i<N;i++,m++) for(j=k;j<M;j++)a[i+1][j]-=temp[m]*a[k][j]; k++; }while(k<N) ; /*第二个do-while是将矩阵消成对角形式，并且重新给k赋值*/ k=L-1; do {n=0; for(l=k;l>=0;l--)temp[n++]=a[k-l][k+1]/a[k+1][k+1]; for(m=0,i=k;i>=0;i--,m++) for(j=k;j<M;j++)a[k-i][j]-=temp[m]*a[k+1][j]; k--; }while(k>=0) ; /*下一个for是解方程组*/ for(i=0;i<N;i++) x[i]=a[i][N]/a[i][i]; } main() { int i,j; double a[N][M]={0},answer[N]; printf("!!The size of Maxrix is %d * %d,each line enter %d element:\n ",N,M,M); for(i=0;i<N;i++) { printf("Enter the Matrix's %d line:\n",i); for(j=0;j<N+1;j++) scanf("%lf",&a[i][j]); } printf("\nthe corss matrix is:\n_ _ _ _ _\n"); gauss(a,answer); for(i=0;i<N;i++) { for(j=0;j<M;j++) printf("%-2lf ",a[i][j]); putchar('\n'); } printf("_ _ _ _ _\nthe solve is:\n"); for(i=0;i<N;i++) printf("x%d=%lf\n",i+1,answer[i]); getch(); }

希望采纳