<s>[INST] Find the closed form for $\int_{0}^{1}\left(x^k+{\ln{x}\over 1-x}\right)^2dx=f(k)$ Find the closed form for

$$\int_{0}^{1}\left(x^k+{\ln{x}\over 1-x}\right)^2dx=f(k)\tag0$$

Setting $k=0$ then $f(0)=1$.
$f(k)$ seems to be rational numbers for all values of $k\ge0.$
We are going to prove that 

$$I=\int_{0}^{1}\left(1+{\ln{x}\over 1-x}\right)^2dx=\color{red}{1}?\tag1$$

Expand $(1)$
$$I=\int_{0}^{1}\left(1+{2\ln{x}\over 1-x}+{\ln^2{x}\over (1-x)^2}\right)dx=1$$
Standard integral $$\int_{0}^{1}{\ln{x}\over 1-x}dx=-{\pi^2\over 6}$$
Then
$$I=1-{\pi^2\over 3}+\int_{0}^{1}{\ln^2{x}\over (1-x)^2}dx=1\tag2$$
$$J=\int_{0}^{1}{\ln^2{x}\over (1-x)^2}dx=\sum_{n=0}^{\infty}(1+n)\int_{0}^{1}x^n\ln^2{x}dx\tag3$$
Standard integral $$\int_{0}^{1}x^n\ln^2{x}dx={2\over (1+n)^3}$$
$$J=\sum_{n=0}^{\infty}(1+n)\cdot{2\over (1+n)^3}={\pi^2\over 3}\tag4$$
Substitute $(4)$ into $(2)$
$$I=1\tag5$$
 [/INST] \begin{align*}
\int_{0}^{1}\left ( x^k + \frac{\ln x}{1-x} \right )^2  \, {\rm d}x&= \int_{0}^{1}x^{2k} \, {\rm d}x + 2\int_{0}^{1}\frac{x^k \ln x} {1-x} \, {\rm d}x + \int_{0}^{1}\frac{\ln^2 x}{(1-x)^2} \, {\rm d}x \\ 
 &= \frac{1}{2k+1} +  2\int_{0}^{1}\frac{x^k \ln x} {1-x} \, {\rm d}x + \int_{0}^{1}\frac{\ln^2 x}{(1-x)^2} \, {\rm d}x
\end{align*}
Now the second integral is easy. 
\begin{align*}
\int_{0}^{1}\frac{x^k \ln x}{1-x} \, {\rm d}x &=\int_{0}^{1} x^k \ln x \sum_{n=0}^{\infty} x^n \, {\rm d}x \\ 
 &= \sum_{n=0}^{\infty} \int_{0}^{1}x^{n+k} \ln x \, {\rm d}x\\ 
 &= -\sum_{n=0}^{\infty} \frac{1}{\left ( n+k+1 \right )^2}\\ 
 &= -\psi^{(0)} (k+1) \\
 &= -\zeta(2) + \mathcal{H}^{(2)}_k
\end{align*}
And the last one using the $\displaystyle \sum_{n=0}^{\infty}x^n =\frac{1}{1-x}$ and differentiating once we have that:
\begin{align*}
\int_{0}^{1}\frac{\ln^2 x}{\left ( 1-x \right )^2} &=\int_{0}^{1} \ln^2 x \sum_{n=1}^{\infty} n x^{n-1} \, {\rm d}x \\ 
 &=\sum_{n=1}^{\infty} n \int_{0}^{1} x^{n-1} \ln^2 x \, {\rm d}x \\ 
 &= 2\sum_{n=1}^{\infty} \frac{n}{n^3}\\ 
 &=\frac{\pi^2}{3}
\end{align*}
Thus:
$$\int_{0}^{1}\left ( x^k + \frac{\ln x}{1-x} \right )^2 \, {\rm d}x = \frac{1}{2k+1} +2 \mathcal{H}_k^{(2)}$$
Here $\displaystyle \mathcal{H}_k^{(2)}=\sum_{n=1}^{k} \frac{1}{n^2}$
Woops. I came third. Never mind, since I wrote all this I leave the message.
Edit:  Fixed error.
</s>