Mathematica suggests that integrating the Logarithamic Integral, $\text{li}(x)=\int_0^x \frac{dt}{\log t}$, multiplied by $x^n$, between the limits $0$ and $1$ leads to the following result
$$\int_0^1 x^n\, \text{li}(x) \, dx= - \frac{\log(n + 2)}{(n + 1)}$$
(See this question for how this improper integral was found)
Does anyone recognise this result from the mathematical literature or know of a straight forward proof?
I don't fully understand how Mathematica arrives at this result, particularly in respect of the improper limit of the integral at $1$.
We can integrate by parts to solve it, however pay attention that we run into divergence issues if we take $\left(\frac{x^{n+1}}{n+1}\right)'=x^n$.
To avoid that we will go with: $\left(\frac{x^{n+1}-1}{n+1}\right)'=x^n$. $$I(n)=\int_0^1 \left(\frac{x^{n+1}-1}{n+1}\right)'\text{li}(x) dx\overset{}=\underbrace{\left(\frac{x^{n+1}-1}{n+1}\right)\text{li}(x)\bigg|_0^1}_{=0} -\frac{1}{n+1}\int_0^1 \frac{x^{n+1}-1}{\ln x}dx$$
Here's an elementary approach for the last integral in order to get $\ln(n+2)$ and the result follows.