Python求解常微分方程

sympy、numpy、scipy、matplotlib是强大的处理数学问题的库，可以执行积分、求解常微分方程、绘图等功能，其开源免费的优势可以与MATLAB媲美。

• 一阶常微分方程

from sympy import *
f = symbols('f', cls=Function)#定义函数标识符
x = symbols('x')#定义变量
eq = Eq(diff(f(x),x,1),f(x))#构造等式，即dy/dx=y
#diff(函数,自变量,求导次数)
print(dsolve(eq, f(x)))
pprint(dsolve(eq, f(x)))#以"pretty"形式打印方程的解

输出结果：

大家可以看出两种打印方式的区别。

• 二阶常微分方程

from sympy import *
f = symbols('f', cls=Function)#定义函数标识符
x,p,q =