1、LU 分解法求解线性方程：#includevoid solve(float l100,float u100,float b,float x,int n)int i,j;float t,s1,s2;float y100;for(i=1;i=1;i-) /* 第二次回代过程开始 */s2=0;for(j=n;ji;j-)t=-uij;s2=s2+t*xj;xi=(yi+s2)/uii; void main()float a100100,l100100,u100100,x100,b100;int i,j,n,r,k;float s1,s2;for(i=1;i=99;i+)/*将所有的数组置零，同时将
2、 L 矩阵的对角值设为 1*/for(j=1;j=99;j+)lij=0,uij=0;if(j=i) lij=1;printf (“input n:n“);/*输入方程组的个数*/scanf(“%d“,printf (“input array A:n“);/*读取原矩阵 A*/for(i=1;i=n;i+)for(j=1;j=n;j+)scanf(“%f“,printf (“input array B:n“);/*读取列矩阵 B*/for(i=1;i=n;i+)scanf(“%f“,for(r=1;r=n;r+)/*求解矩阵 L 和 U*/for(i=r;i=n;i+)s1=0;for(k=1
3、;k=r-1;k+)s1=s1+lrk*uki;uri=ari-s1;for(i=r+1;i=n;i+)s2=0;for(k=1;k=r-1;k+)s2=s2+lik*ukr;lir=(air-s2)/urr;printf(“array L:n“);/*输出矩阵 L*/for(i=1;i=n;i+)for(j=1;j=n;j+)printf(“%7.3f “,lij);printf(“n“);printf(“array U:n“);/*输出矩阵 U*/for(i=1;i=n;i+)for(j=1;j=n;j+)printf(“%7.3f “,uij);printf(“n“);solve(l,u
4、,b,x,n);printf(“解为:n“);for(i=1;i=n;i+)printf(“x%d=%fn“,i,xi);运行结果：input n:3input array A:2 2 34 7 7-2 4 5input array B:3 1 -7array L:1.000 0.000 0.0002.000 1.000 0.000-1.000 2.000 1.000array U:2.000 2.000 3.0000.000 3.000 1.0000.000 0.000 6.000解为:x1=2.000000x2=-2.000000x3=1.000000Press any key to continue
