I observed that :
$$\int_0^\infty \frac{x}{e^x -1}\mathrm dx = \frac{π^2}{6} = \zeta(2)$$
$$\int_0^\infty \frac{x^2}{e^x -1}\mathrm dx = 2\zeta(3)$$
The first one can be evaluated by turning it into a geometric series and switching summation with integral. I derived second one by integration by parts and following same steps as in first one.
Question:
I want to know if it can be generalized as :
$$\int_0^\infty \frac{x^k}{e^x -1}\mathrm dx = k!\zeta(k+1)$$
a) for $k\in \mathbb N$
b) for (atleast for some) $k\in \mathbb R$
For $k\in\mathbb N$ , it suggests true due to pattern of integration by parts. But I am not sure. Can it be extended to some non-natural numbers ?