《C语言写高斯消元法》由会员分享,可在线阅读,更多相关《C语言写高斯消元法(2页珍藏版)》请在金锄头文库上搜索。
1、高斯消元法C语言程序及其输出结果杨学玉( 33)C语言程序:#include #include #include #define N 10double aN+1N+1,bN+1,xN+1;/* 用在此处以便可以让下面的程序也调用 */ main() int n,i,j,k; double temp;printf( 请输入方程组的维数 :n); scanf(%d,&n);if(nN)printf( 错误:元素超过初设定的值 %dn,N); exit(0);printf( 开始输入各元素的值: ); for(i=1;i=n;i+) for(j=1;j=n;j+)printf(请输入元素%d%d,i
2、,j); scanf(%lf,&aij);printf( 请输入常数 b%d,j); scanf(%lf,&bi);for(j=1;j=n-1;j+) for(i=j+1;i=n;i+) temp=aij/ajj;for(k=1;k=1;j-)xj=bj;for(k=n;k=j+1;k-)xj=xj-xk*ajk;xj=xj/ajj;for(j=1;j=n;j+)printf(x%d=%fn,j,xj);/* 列出高斯消元后得到的数组 */fun(int m,int n,double aN+1N+1,double bN+1)int i,j;printf( 消元后的矩阵是: n);for(i=1;i=m;i+)for(j=1;j=n;j+)printf(%f,aij);printf(n) ;C语言输出结果:Mat-lab 程序及输出结果: