How to solve a system of three ODE
4 次查看(过去 30 天)
显示 更早的评论
I was looking into epidemiology, and I came to a system of three ODEs to solve the spread. For the life of me I can't remember how to do this in Matlab. Here is my system, with the initial conditions.
S(0) = 19,465,197
E(0) = 0
I(0) = 0
dS/dt=-147199.4991*S+10*S*I
dE/dt=10*S*I-147203.5121*E
dI/dt=4*E-147199.5197*I
0 个评论
回答(2 个)
Benjamin Schwabe
2012-6-15
Hi James,
start writing it in a vectorform with x=(S,E,I); Then you have x(0)=x0 as defined and you can write dx/dt=f(x); create a function doin this. the you might use ode45 oder another ode solver for you problem. The help is quite detailed here.
Hope that helps, Benjamin
0 个评论
Walter Roberson
2012-6-15
WIth the symbolic toolbox, it would look somewhat like
syms S(t) E(t) I(t)
dsolve(S(0) == 19465197, E(0) == 0, I(0) == 0, diff(S(t), t) == -1.471994991*10^5*S(t)+10*S(t)*I(t), diff(E(t), t) == 10*S(t)*I(t)-1.472035121*10^5*E(t), diff(I(t), t) = 4*E(t)-1.471995197*10^5*I(t))
with solution
E(t) = 0, I(t) = 0, S(t) = 19465197*exp(-147199.4991*t)
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!