问题描述
我对 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)