四阶龙格库塔法

四阶龙格-库塔方法的原理

对于微分方程

采用四阶龙格-库塔法的计算公式为

其中，h为求解步长。

matlab编程

一、

用这个来举例：

结果

二、用ode45

Matlab内部函数ode45直接进行常微分方程求解，函数格式为：
[X,Y]=ode(odefun,tspan,Y0)
Odefun设置为微分方程中需要积分的函数，
tspan为微分方程积分的范围，
Y0为Y的初始值。

结果