上面公式为赶法计算公式 实验目的 通过编写三介 Runge-Kuta法、中点法解常微分方程(初值问题)计算程序,通过编写 追赶法解三对角程组程序进一步熟悉了解递推算法基本程序模块结构以达到能套用该 基本程序模块编同类程序;通过三阶 Runge-Kuta法、中点法计算结果的比较,对算法的 稳定性有较深刻感性认识通过追赶法程序和 Gauss消元法程序比较两算法的效率 实验内容 1.套用递推鲜法基本程序模块编写追赶法程序; 2.套用递推鲜法基本程序模块编写三阶 Runge-Kutta法解常傲分方程程序; 3.编写中点法解常微分方程程序 4给出不同步长及不同的递推次数比较2和3的结果 实验步骤 1.调试并运行你所编写的三阶 Runge-Kutta法程序,尽量将td取大一些h取小一些, 即n值取大一些记下y(t)值; 2调试你所纸写的中点法程序和三阶 Runge-Kutta法用同一程序用例,同-t值,h 值计算并记下各:t)值; 3.调试你所组写的追赶法程序.正确录人Gass列主元消元法程序.用同一程序用例 比较其计算结果·長好在两程序外加一重循环(须修改说明),让循环次数相同用手表或跑 表测量两程序运有谢时间 实验结论 思考题 1.为什么h轼小时,n较大时三阶 Runge-Kutta法程序和中点法程序结果不一样? 2.追赶法算法是否也可列入广义递推算法? 实验四广义递推算法 预备知识 Gas列主元消元法可计算行列式值.其计算过程及计算公式为: 1·按列选主元 算法:尝试法 选第i列主元,尝试范围为一重循环初值取j=i,终值为j=n;定解条件是找|a 最大者 ·260·