如何计算 Python 中的 10 个积分

Question

我需要使用数值方法确定性地计算 python 中的以下总和。

我有以下功能：

有没有类似scipy.integrate.tplquad的可以算出这么多积分的？

Answer 1

使用 sympy

import sympy as sym

x = sym.symbols('x1:11')
expr = sym.prod(x) * sym.exp(-sum(x))
integral = sym.integrate(expr, *((xi, 0, 1) for xi in x))

print(integral) 将显示明确的表达式，尽管 integral.factor() 将显示更简洁的形式。

(-2 + E)**10*exp(-10)    # E is the mathematical constant

您还可以评估任何您想要的精度。

>>> integral.evalf()  # this can take an argument `n` for precision digits
1.65960210179703e-6