用追赶法解三对角线方程组.docVIP

云南大学数学与统计学实验教学中心 实验报告 课程：—2012学年第一学期 成绩： 指导教师： 学生姓名：名称：：2011/10/29 实验学时： 3 学院： 数学与统计学院 专业： 数学与应用数学 年级： 一、目的 内容 三、环境四.实验方法 追赶法；其算法如下： 计算： （i=2,3,4……n-1） 解 (i=2,3,4……,n) 解 （i=n-1,n-2,……,2,1） 五、实验过程 1实验步骤 (1）编程: 根据所用算法及选用语言编出源程序 2）. 开机, 打开所用语言系统输入所编源程序. 3）. 调试程序, 修改错误直至能正确运行. 4）. 运行程序并输出计算结果. #include stdio.h #include conio.h #include math.h /*预处理指令*/ #define MAX_n 100 /*宏定义，即最多可输入100阶方阵*/ #define PRECISION 0.000001 /*精度*/ void SulutionOutput(float x[],int n) /*输出方程组的解*/ /*带分号的为声明*/ void TriDiagonalMatrixInput(float a[],float b[],float c[],float f[],int n) /*三对角元素的输入*/ int Z_G_method(float a[],float b[],float c[],float f[],int n) /*三对角方程组求解 — 追赶法 */ void main() /*主函数*/ { int n; float a[MAX_n],b[MAX_n],c[MAX_n], d[MAX_n],f[MAX_n]; /*定义数组*/ printf(\nInput n=); scanf(%d,n); TriDiagonalMatrixInput(a,b,c,,d,f,n); /*以下三行均为调用函数声明*/ Z_G_method(a,b,c,,d,f,n); SulutionOutput(f,n); } 3 调试过程 1.本打算用递归调用，但是对那个不太熟悉，就改为用循环了； 2. 首先用do……while循环，出现语法错误太多；无法调试； 后改为了do……while循环了； 5.逐步调试修改程序，达到预期结果； 六、总结 Floating point error: Domain. 原因分析： 编写程序时出现了语法错误 （2）编写程序时没去分大小写，造成多处语法错误；原因分析：都怪自己粗心马虎 3．体会和收获。 无论干什么事都要认真，踏实！努力做好自己该做的。 上课要认真听讲，专心学习，做作业才省力！ 各门课程之间相互联系，应搞好各科的学习。 七：程序源代码： #include stdio.h #include conio.h #include math.h #define MAX_n 100 #define PRECISION 0.000001 //解输出 void SulutionOutput(float x[],int n) { int i; for(i=1;i=n;++i) printf(\nx[%d]=%f,i,x[i]); } //三对角方程组元素输入 void TriDiagonalMatrixInput(float a[],float b[],float c[],float f[],int n) { int i; printf(Input b[1],c[1],f[1]:); scanf(%f %f %f,b[1],c[1],f[1]); for(i=2;in;++i) { printf(Input a[%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(Input a[%d],b[%d],f[%d],f[%d]:,n,n,n,n); scanf(%f %f