Question

龙格库塔法matlab编程

里面的Sin=250（mg/L）。我现在需要的是得到S和X的一组数据。怎么编程求解？

Answer 1

用matlab的四、五阶龙格库塔方法（ode45）可以求得其数值解。设定其微分方程组的初始条件X(0)=0，S(0)=0。解决方法：

1、常微分方程组的自定义函数odefun（t，y），其主要内容

dy(1)=3*y(2)/(60+y(2))*y(1)-0.16*y(1);

dy(2)=4*250-4*y(2)-5*y(2)/(60+y(2))*y(1)-0.06*y(1);

2、t、y1、y2的求解

[t,y]=ode45(@odefun,[0 10],y0)

3、计算结果及t与S、X的关系图，S与X的关系图

追问

能不能把完整程序发一下