For all, $n\geq 1$, prove that
$$\int_0^1 x^{2n-1}\ln(1+x)\,dx =\frac{H_{2n}-H_n}{2n}\tag1$$ This identity I came across to know from here, YouTube which is proved in elementary way.
Trying to make alternate efforts to prove it, we can observe that $$\ln(1+x)=\ln\left(\frac{1-x^2}{1-x}\right)=\ln(1-x^2)-\ln(1-x)$$On multiplying by $x^{2n-1}$ and integrating from $0$ to $1$ , we have first order partial derivatives of beta function.
How can we prove that identity $(1)$ without the use of derivatives of beta function?
It is not necessary to use digamma function and multiplication formula of harmonic numbers, indeed
$$\int_0^1 x^{2n-1}\ln(1+x)dx=\\=\sum_{k=1}^{\infty}\frac{(-1)^{k+1}}{k}\int_0^1 x^{2n+k-1}dx =\sum_{k=1}^{\infty}\frac{(-1)^{k+1}}{k(k+2n)}=\\=\sum_{k=1}^{\infty}\left(\frac{1}{(2k-1)(2n+2k-1)}-\frac{1}{2k(2n+2k)}\right)=\\=\frac{1}{2n}\sum_{k=1}^{\infty}\left(\frac{1}{2k-1}-\frac{1}{2(n+k)-1}\right)-\frac{1}{4n}\sum_{k=1}^{\infty}\frac{n}{k(n+k)}=\\=\frac{1}{2n}\sum_{k=1}^n\frac{1}{2k-1}-\frac{1}{4n}\sum_{k=1}^{\infty}\left(\frac{1}{k}-\frac{1}{n+k}\right)=\\=\frac{1}{2n}\sum_{k=1}^n\frac{1}{2k-1}-\frac{1}{4n}\sum_{k=1}^n\frac{1}{k}=\\=\frac{1}{2n}\left(\sum_{k=1}^n\frac{1}{2k-1}-\sum_{k=1}^n\frac{1}{2k}\right)=\\=\frac{1}{2n}\left(\sum_{k=1}^n\frac{1}{2k-1}+\sum_{k=1}^n\frac{1}{2k}-\sum_{k=1}^n\frac{1}{k}\right)=\\=\frac{1}{2n}\left(\sum_{h=1}^{2n}\frac{1}{h}-\sum_{k=1}^n\frac{1}{k}\right)=\frac{H_{2n}-H_n}{2n}\;.$$