为什么不能同意这一点（Mathematica可以）

Question

为什么sympy不能执行此积分：

import sympy as sp
x = sp.symbols('x', real=True, nonzero=True)
sp.integrate(x**3/(sp.exp(x)-1), (x, 0, sp.oo))

这应该返回pi ^ 4 / 15.

Answer 1

使用Integral，然后使用`evalf'

评估函数

import sympy as sp
from sympy import Integral, oo, exp, pi, integrate
from sympy.abc import a, x, y
x = Integral(x**3/(exp(x)-1), (x, 0, oo))
x.evalf()