使用SymPy获得第二个隐式衍生物

Question

我们在jupyter内核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')

然后运行：

eq = x ** 2 - x * y + y ** 2 - 3
eqd1 = idiff(eq, y, x)
eqd1

我们得到了正确答案：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 ³ ）

这是正确的结果（稍微简化之后）。