I ran into this integral today
$$\int_0^1 \left( x \ln x \right)^{2020} \,\mathrm{d}x \overset{?}{=} \frac{\Gamma(2021)}{2021^{2021}}$$
My solution goes along these lines. Recalling the identity
$$\ln x = \lim_{n \rightarrow +\infty} n \left( x^{1/n} - 1 \right)$$
one can write the integral as follows:
\begin{align*} \int_{0}^{1} \left ( x \ln x \right )^{2020}\, \mathrm{d}x &= \lim_{n \rightarrow +\infty} n^{2020} \int_{0}^{1} \left ( x^{2020} \left ( x^{1/n} -1 \right )^{2020} \right )\, \mathrm{d}x \\ &\!\!\!\!\!\!\overset{u=x^{1/n}}{=\! =\! =\! =\!} \lim_{n \rightarrow +\infty} n^{2021} \int_{0}^{1} u^{2021n-1} \left ( 1-u \right )^{2020} \, \mathrm{d}u \\ &=\lim_{n \rightarrow +\infty} n^{2021} \int_{0}^{1} u^{2021n-1} \left ( 1-u \right )^{2021-1} \, \mathrm{d}u \\ &=\lim_{n \rightarrow +\infty} n^{2021} \mathrm{B} \left ( 2021n, 2021 \right ) \\ &=\lim_{n \rightarrow +\infty} n^{2021} \; \frac{\Gamma \left ( 2021 n \right ) \Gamma \left ( 2021 \right )}{\Gamma \left ( 2021 n + 2021 \right )} \\ &= \Gamma \left ( 2021 \right ) \lim_{n \rightarrow +\infty} n^{2021} \frac{\Gamma \left ( 2021 n \right )}{\Gamma \left ( 2021 n + 2021 \right )} \end{align*}
Using Gautschi's inequality one has that
$$n^{2021}\left ( 2021n -1 \right )^{1-2022}<\frac{ n^{2021} \Gamma\left ( 2021n -1 + 1 \right )}{\Gamma\left ( 2021n-1 + 2022 \right )}< n^{2021}\left ( 2021n\right )^{1-2022}$$
and thus by the squeeze theorem one has the result I stated above.
My doubt however is at the application of the inequality. Is it correct? Can you suggest other ways of solving the problem?