在Sympy中找到微分方程的阶

Question

我有一个在Sympy（0.7.6.1）中定义的微分方程，例如：

>>> y = Function('y')
>>> x = Symbol('x')
>>> y0 = y(x)
>>> y1 = y0.diff(x)
>>> y2 = y1.diff(x)
>>> eq = y2 - 2*y1 + y0
>>> eq
                      2      
         d           d       
y(x) - 2*--(y(x)) + ---(y(x))
         dx           2          
                    dx   

我希望Sympy告诉我等式的顺序。 我在文档中什么都没找到。 我微弱地希望该degree可以工作，但事实并非如此：

>>> degree(eq)
1
>>> degree(eq, y)
0

上下文：出于教学目的，我编写了一个常数系数线性齐次ODE的求解器。 它找到特征方程并用Sympy的roots求解。 目前，我必须在ODE本身旁边手动传递ODE订单，这很不方便。

Answer 1

您可以使用ode_order 。 这为您提供了相对于函数的微分方程的阶数：

>>> ode_order(eq, y)
2