微分方程的計算

例如，我們要解以下微分方程：

\displaystyle \frac{dx}{dt}=kx

 In [1]: from sympy import *
 
 In [2]: from sympy.abc import x, t, k

 In [3]: x = Function('x')

 In [4]: dsolve(Eq(Derivative(x(t), t) - k * x(t), 0), x(t))
 Out[4]:
 \displaystyle x(t) = C_{1} e^{k t}

我們也可以設定初始條件，例如，在上例中，t = 0 時，x(0) = 0.5。這樣就可以求得常數。

 In [8]: dsolve(Eq(Derivative(x(t), t) - k * x(t), 0), x(t), ics = {x(0): 0.5})
 Out[9]:
 \displaystyle x(t) = 0.5 e^{k t}

參考閱讀：
https://docs.sympy.org/latest/modules/solvers/ode.html