Python求解二阶微分方程的解析解

代码：

from sympy import symbols, Function, dsolve

# 定义自变量和因变量
x = symbols('x')
y = Function('y')(x)

# 定义微分方程
eq = y.diff(x, 2) + 4 * y.diff(x) + 3 * y - x
y = Function('y')
# 使用 dsolve 求解微分方程
solution = dsolve(eq, y(x))
print(solution)

结果：

Eq(y(x), C1exp(-3x) + C2*exp(-x) + x/3 - 4/9)