Estoy tratando de averiguar cómo evaluar las siguientes: $$ J=\int_{0}^{\infty}\frac{x^3}{e^x-1}\ln(e^x - 1)\,dx $$ Estoy intentado considerando $I(s) = \int_{0}^{\infty}\frac{x^3}{(e^x-1)^s}\,dx\implies J=-I'(1)$, pero yo no podía entender lo $I(s)$ fue. Mi otra idea era contorno de integración, pero no estoy seguro de cómo tratar con el logaritmo. Mathematica dice que $J\approx24.307$.
Le he pedido a una pregunta similar y la respuesta involucrados $\zeta(s)$, por lo que sospecho que esto va así.