三对角方阵的追赶法

《三对角方阵的追赶法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、班级：数学114学号：201100144407姓名：徐美荣3一：实验目的加强"追赶法"在解方程组中的应用，掌握多种不同的解线性方程组的解法以及编程的应用二：实验内容用追赶法解下列严格对角优势的三对角线方程组2-100x10-13-20x2=10-12-1x3000-35x41三：程序源代码：#include#include#include#defineMAX_n100#definePRECISION0.000001//解输出voidSuluti

2、onOutput(floatx[],intn){inti;for(i=1;i<=n;++i)printf("x[%d]=%f",i,x[i]);}//三对角方程组元素输入voidTriDiagonalMatrixInput(floata[],floatb[],floatc[],floatf[],intn){inti;printf("Inputb[1],c[1],f[1]:");scanf("%f%f%f",&b[1],&c[1],&f[1]);for(i=2;i

3、"Inputa[%d],b[%d],c[%d],f[%d]:",i,i,i,i);scanf("%f%f%f%f",&a[i],&b[i],&c[i],&f[i]);}printf("Inputa[%d],b[%d],f[%d],f[%d]:",n,n,n,n);scanf("%f%f%f%f",&a[n],&b[n],&c[n],&f[n]);}//三对角方程组求解—追赶法3intZ_G_method(floata[],floatb[],floatc[],floatf[],intn){inti

4、;c[1]/=b[1];for(i=2;i0;--i)f[i]-=c[i]*f[i+1];}voidmain(){intn;floata[MAX_n],b[MAX_n],c[MAX_n],d[MAX_n],f[MAX_n];printf("Inputn=");sca

5、nf("%d",&n);TriDiagonalMatrixInput(a,b,c,,d,f,n);Z_G_method(a,b,c,,d,f,n);SulutionOutput(f,n);}3