I read somewhere that
$$ \int_{0}^{\infty} \frac{x^3}{e^{x} - 1} = \frac{\pi^4}{15}$$
Does anyone see a way to prove this? My first idea was doing a contour integration and use the residue theorem, but that seems to be a lot of work so if anyone has a better idea I'd be glad to hear it.
For any $n\in\mathbb{N}^+$, $$ \int_{0}^{+\infty}\frac{x^n}{e^x-1}\,dx = \sum_{m\geq 1}\int_{0}^{+\infty}x^n e^{-mx}\,dx = n!\cdot\sum_{m\geq 1}\frac{1}{m^{n+1}} = n!\cdot\zeta(n+1).$$ In our case, $6\cdot\zeta(4)=\frac{\pi^4}{15}$.