Question

我想在python中使用符号代数计算阶乘。我可以生成的最简单的问题版本如下：

from sympy.abc import j
from math import factorial
from sympy import summation
summation(factorial(j), (j, 1, 4))

我收到以下错误：

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "sympy/core/expr.py", line 194, in __int__
    r = self.round(2)
  File "sympy/core/expr.py", line 3042, in round
    raise TypeError('%s is not a number' % type(x))
TypeError: <class 'sympy.core.symbol.Symbol'> is not a number

从根本上说，我想要计算的是

summation(x**(j-1)/factorial(j-1), (j, 1, 3))

有什么建议吗？

Answer 1

使用Sympy自己的factorial函数（而不是math模块的阶乘函数）可能会返回您想要的内容。

按照初始设置但省略from math import factorial后，您可以写：

>>> from sympy import factorial, symbols
>>> x = symbols('x')
>>> summation(x**(j-1)/factorial(j-1), (j, 1, 3))
x**2/2 + x + 1

这将因子系列的总和减少为简单的二次方程。

我注意到你正在计算exp(x)幂级数展开的前几个项的总和：

1 + x + x**2/2! + x**3/3! + ...

Answer 2

我能想到避免循环的唯一解决方案是功能版本：S：

from math import factorial

reduce(int.__add__,map(factorial,xrange(1,5)))

Python中的求和因子

2 个答案: