Como he oído, no es más que una forma de la piel de una integral. Y así, en este documento, se presenta un enfoque que se basa en la Polylogarithm Funciones. Vamos a participar en una serie de cuatro integración por partes que culminará en la respuesta $24\text{Li}_4(1)=24\zeta(4)$.
Para ello, comenzamos con la integral de interés $I$, el cual es dado por
$$I=\int_0^{\infty}\frac{x^4e^x}{(e^x-1)^2}\,dx$$
Hacer cumplir la sustitución de $x\to -\log (1-x)$ rendimientos
$$I=\int_0^1\frac{\log^4 (1-x)}{x^2}\,dx \tag 1$$
Ahora comenzamos una serie de cuatro sencilla integración por partes aplicaciones.
En primer lugar, dejando $u=\log^4(1-x)$ $v=-1/x$ $(1)$ rendimientos
$$I=-4\int_0^1\frac{\log^3 (1-x)}{x}\,dx \tag 2$$
donde nos tácitamente aplica $\lim_{x\to 1}\left(\log^4(1-x)(1-1/x)\right)=0$ para llegar a $(2)$.
Segundo, dejando $u=-4\log^3(1-x)$ $v=\log x$ $(2)$ rendimientos
$$I=-12\int_0^1 \log^2(1-x)\frac{\log x}{1-x}\,dx \tag 3$$
Tercero, dejando $u=-12\log^2(1-x)$ $v=\text{Li}_2(1-x)$ $(3)$ rendimientos
$$I=-24\int_0^1 \log(1-x)\frac{\text{Li}_2(1-x)}{1-x}\,dx \tag 4$$
Cuarto, dejando $u=-24\log(1-x)$ $v=-\text{Li}_3(1-x)$ $(4)$ rendimientos
$$\begin{align}
I&=-24\int_0^1 \frac{\text{Li}_3(1-x)}{1-x}\,dx\\\\ \tag 5
&=24\text{Li}_4(1)\\\\
&=24\zeta(4)\\\\
&=24\frac{\pi^4}{90}
\end{align}$$
Por lo tanto, hemos
$$\bbox[5px,border:2px solid #C0A000]{I=\frac{4\pi^4}{15}}$$
como era de esperar y tan fácil como uno, dos, tres, cuatro!