python数学建模导论2.2 常微分方程符号解与数值解的计算（未放入章节为微分方程介绍章节，自行百度）

from sympy import *
y = symbols('y', cls=Function)
x = symbols('x')
eq = Eq(y(x).diff(x,2)+2*y(x).diff(x,1)+y(x),x*x)
print(dsolve(eq, y(x)))

结果

from sympy import *
y = symbols('y', cls=Function)
x = symbols('x')
eq = Eq(y(x).diff(x,2)+4*y(x).diff(x,1)+29*y(x), 0)
print(dsolve(eq, y(x)))

C1=symbols('C1')
C2=symbols('C2')
f=(C1*sin(5*x) + C2*cos(5*x))*exp(-2*x)    #第五行eq的解
print(f.diff(x,1)) #eq的一阶导的形式

'''
(0*C1+1*C2)*1=0
-2*(C1*0+C2*1)*1+(5*C1*1-5*C2*0)*1=15
'''