用SymPy求第二个隐式导数

Question

我们在Python核Python3中有以下细胞：

from sympy.interactive.printing import init_printing
from sympy.core.symbol import symbols
from sympy.geometry.util import idiff
init_printing(use_unicode=True)
x = symbols('x')
y = symbols('y')

然后就跑：

^{pr2}$

我们得到正确答案：2x-y/x-2y

但如果运行代码：

eqd2 = idiff(eqd1, y, x)
eqd2

我们得到：x/y，但答案是：18/（x-2y）²或6（x²-xy+y²）/（x-2y）³

x/y和18/（x-2y）²看起来并不是一个相等的表达式或一个是另一个的阶乘

Answer 1

跑步

eqd2 = idiff(eq, y, x, n=2)
simplify(expand(eqd2))

生产

6(x²-xy+y²)/(x³-6x²y-12xy²-8y³)

这是正确的结果（经过一点简化）。在