1. 程式人生 > >Python解微分方程

Python解微分方程

ini diff als lns sympy dso mbo color symbol

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))

Python解微分方程