计算物理作业(龙格库塔四阶求解初值).docVIP

EAST CHINA INSTITUTE OF TECHNOLOGY 计算物理作业 学生学号： 2014203701 学生姓名： 专 业： 物理学 班 级: 1423701 2017年 1月6 日 -一编写程序求解微分方程 用四阶龙格一库塔法求解微分方程的初值问题 dy dy dx y p(2.0) = 1 2.0 x 2.6 问题描述： 已知的一个微分方程及初值条件如下所示: y=x/y y(2.o)= 1 2.0 x 2.6 编写一个程序求解该微分方程，并分别以步长0.06计算y的值，然 后以图像的形式将通过程序运算得到的结果与精确解r=x2-3进行 比较。 主要算法与公式: 在这里，我们求解微分方程的方法是四阶Runge-Kutta法 y(k +1) = y(k) + h * (k 1 + 2 * k2 + 2 * k3 + k4)/6; kl = f(x(k),y(k)); k2 = f(x(k) + 0.5 * h, y(k) + 0.5 * h * k 1); k3 = f(x(k) + 0.5 * h, y(k) + 0.5 * h * k2); k4 = f(x(k) + h,y(k) + h*k3); 采用这种方法计算的结果一般比较精确，并且编写的程序具有占用内 存小的优点。 源代码:function SDeqRunKu4(dfun, a, b, yO