在学习《数学实验与建模》这本书的时候,第8.3节微分方程建模实例中
Malthus模型和Logistic模型中,都提高了符号解,但是并没有给出符号解的计算code,
在这里将相关代码给出,作为学习笔记,期间查阅了该博主的博客,表示感谢。
https://blog.csdn.net/weixin_45870904/article/details/113151534
Malthus模型
import sympy as sp
from sympy import diff,dsolve,simplify
x=sp.symbols('x',cls=sp.Function)
t=sp.symbols('t')
r=sp.symbols('r')
eq=diff(x(t),t)-(r*x(t))
x=dsolve(eq)
print(simplify(x))
Logistic模型
import sympy as sp
import math
from sympy import diff,dsolve,simplify
x=sp.symbols('x',cls=sp.Function)
t=sp.symbols('t')
r=sp.symbols('r')
xm=sp.symbols('xm')
c1=sp.symbols('c1')
x0=sp.symbols('x0')
eq=diff(x(t),t)-r*(1-(x(t)/xm))*x(t)
x=dsolve(eq)
print(simplify(x))