Python解微分方程

1.求解常微分方程的步驟：

from sympy import *
init_printing()
#定義符號常量x 與 f(x) g(x)。這里的f g還可以用其他字母替換，用於表示函數
x = Symbol('x')
f, g = symbols('f g', cls=Function)

#用diffeq代表微分方程： f''(x) − 2f'(x) + f(x) = sin(x) 
diffeq = Eq(f(x).diff(x, x) - 2*f(x).diff(x) + f(x), sin(x))
#調用dsolve函數,返回一個Eq對象，hint控制精度
print(dsolve(diffeq, f(x),hint='1st_linear'))