Почему я не могу это интегрировать (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()