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

loading...


2

我们在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²-xyy²)/(x-2y)³

x / y和18 /(x-2y)²看起来不是一个相同的表达式或另一个是因子的

1回答

  • 3

    运行

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

    产生

    6(x2-xy y2)/(x3-6x2y-12xy2-8y3)

    这是正确的结果(稍微简化之后) .

loading...

评论

暂时没有评论!