龙格库塔法在matlab中能不能实现一阶常微分方程的计算?

来源:百度知道 编辑:UC知道 时间:2024/07/04 23:59:52
比如说:“dp/dt=-3*p/(100+t)",p(0)=7;怎么用龙格库塔法在matlab中实现。

ode45和ode23都是用龙哥库塔法求解常微分方程的命令,直接用即可,具体用法可参看其帮助文档。
[t,y]=ode45(@(t,p)-3*p./(100+t),[0 200],7);
plot(t,y);xlabel('t');