龙格库塔法求微分方程(共13页).doc

精选优质文档-倾情为你奉上MATLAB程序设计实践课程考核一、编程实现“四阶龙格库塔（R-K）方法求常微分方程”，并举一例应用之。【实例】采用龙格-库塔法求微分方程：1、算法说明：在龙格-库塔法中，四阶龙格-库塔法的局部截断误差约为o(h5)，被广泛应用于解微分方程的初值问题。其算法公式为：其中：2、流程图：2.1、四阶龙格库塔（R-K）方法流程图：输入待求微分方程、求解的自变量范围、初值以及求解范围内的取点数等。确定求解范围内的步长k = 取点数？否求解：求解并输出：是结束算法2.2、实例求解流程图：输入求解的自变量范围求出待求简单微分方程的真值解用MATLAB自带函数ode23求解待求微分方程结束用自编函数四阶龙格库塔（R-K）方法求解待求微分方程开始3、源程序代码3.1、四阶龙格库塔（R-K）方法源程序：function x,y = MyRunge_Kutta(fun,x0,xt,y0,PointNum,varargi