Sympy中莱布尼茨法则衍生零积分未化简为零的问题求助

原因分析 #

当你对未求值的符号积分（即Integral(R*p(R), (R, b, r))这种未计算出原函数的积分对象）应用莱布尼茨求导法则时，SymPy会严格按照法则展开所有项——包括被积函数对求导变量的偏导积分项。哪怕这个偏导结果是0，SymPy也不会自动将这个Integral(0, ...)化简为0，因为这个积分对象处于“未求值”状态，默认不会主动触发这类特殊情况的化简。

而第一个例子中直接调用sym.integrate(0, ...)能得到0，是因为这是直接触发积分计算，SymPy会立即处理常数0的积分。

解决方法 #

你可以通过以下方式强制化简掉多余的积分项：

方法：调用doit()方法触发积分求值

对求导后的结果调用doit()，它会计算表达式中所有未求值的积分对象，包括Integral(0, ...)：

# 修改求导代码
i3s = sym.diff(i2s, r).doit()
print('i3s =', i3s)

运行后输出：

i3s = r*p(r)

补充说明 #

你提到的sym.simplify()调用多余是对的——默认的simplify()不会主动处理未求值的积分对象，必须显式触发求值操作。

内容的提问来源于stack exchange，提问作者sancho.s ReinstateMonicaCellio