我试图使用sympys diff计算f（y（x））w.r.t到x的更高导数。
from sympy import *
from IPython.display import display
init_printing(use_latex=True)
x = symbols('x')
f, y = symbols('f, y', cls=Function)
d2 = diff(f(y(x)),x,2)
print(d2)
print(d2.doit())
Sympy回归：
Derivative(y(x), x)**2*Derivative(f(y(x)), y(x), y(x)) + Derivative(y(x), x, x)*Subs(Derivative(f(_xi_1), _xi_1), (_xi_1,), (y(x),))
Derivative(f(y(x)), y(x))*Derivative(y(x), x, x) + 2*Derivative(y(x), x)**2*Derivative(f(y(x)), y(x), y(x))
胶乳图片：Sympy result。
虽然第一个结果似乎是正确的，但我不理解doit（）操作后第二个表达式中的因子2。
答案 0 :(得分：0)
看起来你偶然发现了一个was just fixed a few weeks ago。
的错误您可以通过将[0-9]{2} [0-9]{2} [0-9]{4}，f和y分别替换为某些函数或值来对此进行测试（附加到您的代码中）：
x打印值应该相同。
这个问题可以进一步分离到：
f_ex = Lambda(x, x**2)
y_ex = Lambda(x, sin(x))
x_ex = 2
substitutions = [ (f,f_ex), (y,y_ex), (x,x_ex) ]
print( d2.subs(substitutions).doit().n() ) #-1.30728724172722
print( d2.doit().subs(substitutions).doit().n() ) #-0.960930862590836
在这里，普通print((Derivative(f(y(x)), x, x)))
print((Derivative(f(y(x)), y(x), y(x))).doit())
只添加因子2，这显然是错误的。