SymPy:KroneckerDelta和指数的衍生品

Question

请参考以下MWE：

import sympy as s
x = s.IndexedBase('x')
y = s.IndexedBase('y')
i,j,k = map(s.Idx,['i','j','k'])
a = s.exp(x[i]*y[j]*s.KroneckerDelta(i,j))
b = a.diff(x[j])

b的值是 b 的 LaTeX 渲染。由于我还不允许嵌入图像，这里是文本形式

((Derivative(KroneckerDelta(i, j), i)*Derivative(i, x[j]) + Derivative(KroneckerDelta(i, j), j)*Derivative(j, x[j]))*x[i]*y[j] + KroneckerDelta(i, j)*y[j])*exp(KroneckerDelta(i, j)*x[i]*y[j])

关键的一点是，KroneckerDelta的衍生品i和j指数存在未经评估的衍生品，i和j指数的衍生品x[i]。为什么这些不0？

Answer 1

我认为，我所期望的索引变量符号区分的行为尚未在 Sympy 中完全实现。我继续使用其他工具，如Maxima或Mathematica（甚至笔和纸）进行这些计算。