Recordar la representación integral de la Polylogarithm dada por
$$\operatorname{Li}_s(z)~=~\frac1{\Gamma(s)}\int_0^\infty \frac{t^{s-1}}{e^t/z-1}\mathrm dt$$
Además tenga en cuenta que
$$\frac{\partial}{\partial z}\left[\frac{t^{s-1}}{e^t/z-1}\right]=\frac{t^{s-1}e^t}{(e^t-z)^2}$$
que es precisamente la estructura de su integrando. Por lo tanto, mediante la aplicación de una variación de Feynman, el Truco de la que obtenemos
\begin{align*}
\frac{\mathrm d}{\mathrm dz}\Gamma(s)\operatorname{Li}_s(z)&=\frac{\mathrm d}{\mathrm dz}\int_0^\infty \frac{t^{s-1}}{e^t/z-1}\mathrm dt\\
\Gamma(s)\frac{\mathrm d}{\mathrm dz}\operatorname{Li}_s(z)&=\int_0^\infty \frac{\partial}{\partial z}\frac{t^{s-1}}{e^t/z-1}\mathrm dt\\
\Gamma(s)\frac{\operatorname{Li}_{s-1}(z)}z&=\int_0^\infty\frac{t^{s-1}e^t}{(e^t-z)^2}\mathrm dt
\end{align*}
Ahora conectar $s=5$ e $z=1$ rendimientos
\begin{align*}
\Gamma(5)\frac{\operatorname{Li}_{5-1}(1)}1&=\int_0^\infty\frac{t^{5-1}e^t}{(e^t-1)^2}\mathrm dt\\
4!\operatorname{Li}_4(1)&=\int_0^\infty\frac{t^{4}e^t}{(e^t-1)^2}\mathrm dt
\end{align*}
Tenga en cuenta que el lado derecho es la integral nos interesa, mientras que el lado izquierdo puede ser evaluado directamente en términos de Riemann Zeta Función. En general esto se reduce a utilizar el valor de la $\zeta(4)$, como J. G.'s respuesta hizo así, desde la $\operatorname{Li}_s(1)=\zeta(s)$ como se puede comprobar mediante el uso de la serie de la representación de la Polylogarithm.
$$\therefore~\int_0^\infty\frac{t^{4}e^t}{(e^t-1)^2}\mathrm dt~=~4!\operatorname{Li}_4(1)~=~\frac{4\pi^4}{15}$$