Can you please help me in the step-by-step calculation of this improper intergal?
$$\int_{0}^{\infty}\frac{x^2 \ln x}{e^x-1} \ \mathrm d x$$
I can solve a similar integral where I replaced $\ln(x)$ with $x$. So the integral becomes, $$\int_{0}^{\infty}\frac{x^3}{e^x-1} \ \mathrm d x$$ $$=\int_{0}^{\infty}\frac{x^3e^{-x} }{1-e^{-x} } \ \mathrm d x$$ so that $$=\sum_{n=0}^{\infty} \int_{0}^{\infty}x^3e^{-x(n+1) }\mathrm d x$$ $$=Γ(3+1)\sum_{n=0}^{\infty} \frac{1}{(n+1)^{3+1}}$$ $$=Γ(4)ζ(4)$$ $$=6\frac{π^4}{90}$$ $$\approx 6.493$$
Thank you in advance for your detailled explanation and derivation.