带有 sympy 的 n 阶导数

问题描述

我对 sympy 有点陌生
我想使用 sympy 计算表达式的 n 阶导数;但是,我不明白 diff 函数如何用于 n 阶导数:

from sympy import diff,symbols

x = symbols("x")
f = ((x**2-1)**5)

# for n = 2
# from the sympy docs,I do:
d_doc = diff(f,x,x)

# using the diff two times
d_2 = diff(diff(f,x),x)

我得到两个不同的结果:

>>> d_doc
10*(x**2 - 1)**3*(9*x**2 - 1)

>>> d_2
80*x**2*(x**2 - 1)**3 + 10*(x**2 - 1)**4

d_2 是这种情况下的正确答案。

这是为什么?
有没有办法制作一个接受 n 并返回 n 阶导数的函数

解决方法

一篇文章中的答案,(来自 Pranav Hosangadi 的评论):

同理,diff(f,x,x)简化了表达式

>>> simplify(diff(f,x))
(x**2 - 1)**3*(90*x**2 - 10)
>>> simplify(diff(diff(f,x),x))
(x**2 - 1)**3*(90*x**2 - 10)