How to write and solve system of differential equations?

Question

Kon Dig 2017-1-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/321579-how-to-write-and-solve-system-of-differential-equations

回答： Torsten 2017-1-23

thesis equations.JPG

I have the following two equations that describe how drug "d" is released from a sphere of radius "r" made of a polymer that diodegrades into "m" monomers: (you can see the equations in the jpg file).

I want to feed matlab with the equations so that I can extract a release time result. Need help with how to present the equations in matlab, which solver to use and any feedback that can make the system clear to my understanding.

FYI: (1.) a: most probably will be 1, (2.) Dm / Dd / k_DH: constants, (3.) the second term of both equations (v_r * (∂C_d)/∂r ) could be omitted at first

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Lateef Adewale Kareem 2017-1-23

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/321579-how-to-write-and-solve-system-of-differential-equations#answer_251638

This is a pde. To solve this problem, I suggest you use method of lines. Discretize the equation in space. let the solution at each point in space be represented by a member of a vector. then write the equation for each member as a function of such member, its neighbors and time (use boundary condition to account for neighbours of most extreme members. then pass the problem to ODE45.

Although you are calling an ode solver, but the solution to the system is actually the solution to the pde and you can animate the Concentration versus r plot over time.