使用Sympy进行代数递归

Question

我想获得可以通过使用SymPy递归定义的多项式的代数项。

作为使用Maple的示例，Hermite多项式定义为

H(0,x):=1;
H(1,x):=2*x;
for n from 1 to 6 do
    H(n_1,x):= 2*x*H(n,x) - 2*n*H(n-1,x);
    print( simplify(%) );
end do:

产生所需的代数项，但是，我无法在SymPy中找到相同的方法。请您帮忙？

Answer 1

就像在Maple中一样，您可以循环，但是在这里使用一个函数似乎是合理的。您可以使用标准的python函数：

clean_number

当然，如果需要，您可以在函数中进行简化。但是，由于这可能导致表达式无法取消，因此我只在脚本底部简化了一次。