《解线性方程组.doc》由会员分享，可在线阅读，更多相关《解线性方程组.doc（18页珍藏版）》请在金锄头文库上搜索。
1、#include#include#define n 3/2-1 顺序解线性方程组void Cgauss(float ann,float bn,float xn)int i,j,k;float sum=0.0;for(k=0;k=0;i-)sum=0.0;for(j=n-1;ji;j-)sum=sum+aij*xj;xi=bi-sum;/2-2 顺序解线性方程组void Cgauss1(float ann,float bn,float xn)int i,j,k;float temp,sum=0.0;for(k=0;kfabs(akk)&akk!=0)for(j=k;j=0;i-)sum=0.0;
2、for(j=n-1;ji;j-)sum=sum+aij*xj;xi=bi-sum;/2-3 顺序解线性方程组void Cgauss2(float ann,float bn,float xn)int i,j,is,jsn,k;float temp,max,t,sum=0.0;for(k=0;k max)max=t,jsk=j,is=i; if(max+1.0=1.0) printf(矩阵异常!);elseif(jsk!=k)/列交换, 每次列交换都是全部交换，因此 i 是从 0 开始for(i=0;i=0;i-)sum=0.0;for(j=n-1;ji;j-)sum=sum+aij*xj;xi=
3、bi-sum;void main()int i,j;float ann=1,2,-2,2,1,2,3,0,4,bn=2,3,1,xn;printf(方程各元素如下 :n);for(i=0;i#include#include#define n 3void Crout(float ann,float bn,float xn)int r,i,k;float Lnn, Unn,yn,sum;for(i=0;i=0;r-)sum=0.0;for(i=r+1;i#include#include#define n 3void Doolittle(float ann,float bn,float xn)int
4、 r,i,k,j;float Lnn=0, Unn=0,yn=0,sum;/*初始化矩阵 l*/ for(i=0; i=0;r-)sum=0.0;for(i=r+1;i#include#include#define n 3void Ssqrt(float ann,float bn,float xn)int j,i,k;float Lnn=0,yn,sum;/ L00=(float)sqrt(a00);for(j=0;j=0;i-)sum=0.0;for(k=i+1;k#include#include#define n 4void Chase(float Ann,float fn,float x
5、n)int i,j;float an, bn,cn,a1n, b1n,c1n,yn;/系数矩阵用三个一维数组表示for(i=1,j=0;i=0;i-)xi=yi-c1i*xi+1;printf(方程解为:n);for(i=0;in;i+)printf(x%d=%.4f ,i,xi);printf(n);void main()int i,j;float ann=3,1,0,0,2,3,1,0,0,2,3,1,0,0,1,3,bn=1,0,1,0,xn;printf(-追赶法分解法解线性方程组-n);printf(方程各元素如下:n);for(i=0;in;i+)for(j=0;jn;j+)printf(%.2f ,aij);printf(n);printf(方程常数依次为:n);for(i=0;in;i+)printf(%.2f ,bi);printf(n);Chase(a,b,x);
