python求微分方程的解

在 Python 中求微分方程的解可以使用 sympy 模块。下面是一个简单的例子：

假设有一个一阶常微分方程 dy/dx = x^2，初始条件为 y(0) = 1。我们可以使用 sympy 模块中的 dsolve 函数来求解：

from sympy import symbols, Function, diff, dsolve

x = symbols("x")
y = Function("y")(x)

eq = diff(y, x) - x**2
ic = {y.subs(x, 0): 1}

sol = dsolve(eq, ics=ic)
print(sol)

输出结果为：

Eq(y(x), x**3/3 + 1)

即该微分方程的解为 y = x^3/3 + 1。

需要注意的是，sympy 模块的 dsolve 函数只能求解一些简单的微分方程，对于复杂的微分方程可能无法求解。