I want to calculate the series of the Basel problem $\displaystyle{\sum_{n=1}^{\infty}\frac{1}{n^2}}$ by applying the Euler series transformation. With some effort I got that
$$\displaystyle{\frac{\zeta (2)}{2}=\sum_{n=1}^{\infty}\frac{H_n}{n2^n}}.$$
I know that series like the $\displaystyle{\sum_{n=1}^{\infty}\frac{H_n}{n2^n}}$ are evaluated here, but the evaluations end up with some values of the $\zeta$ function, like $\zeta (2),\zeta(3).$
First approach: Using the generating function of the harmonic numbers and integrating term by term, I concluded that
$$\displaystyle{\sum_{n=1}^{\infty}\frac{H_n}{n2^n}=\int_{0}^{\frac{1}{2}}\frac{\ln (1-x)}{x(x-1)}dx},$$
but I can't evaluate this integral with any real-analytic way.
First question: Do you have any hints or ideas to evaluate it with real-analytic methods?
Second approach: I used the fact that $\displaystyle{\frac{H_n}{n}=\sum_{k=1}^{n}\frac{1}{k(n+k)}}$ and then, I changed the order of summation to obtain
$$\displaystyle{\sum_{n=1}^{\infty}\frac{H_n}{n2^n}=\sum_{k=1}^{\infty}\frac{2^k}{k}\left(\sum_{m=2k}^{\infty}\frac{1}{m2^m}\right)}.$$
To proceed I need to evaluate the
$$\int_{0}^{\frac{1}{2}}\frac{x^{2k-1}}{1-x}dx,$$
since $\displaystyle{\sum_{m=2k}^{\infty}\frac{1}{m2^m}=\int_{0}^{\frac{1}{2}}\frac{x^{2k-1}}{1-x}dx}.$
Second question: How can I calculate this integral?
Thanks in advance for your help.
$$ \sum_{n\geq 1}\frac{H_{n-1}}{n}x^n = \frac{1}{2}\log(1-x)^2 \tag{1} $$ follows from the termwise integration of $\sum_{n\geq 1}H_n x^n = \frac{-\log(1-x)}{1-x}.$ It leads to $$ \sum_{n\geq 1}\frac{H_{n-1}}{n 2^n} = \frac{1}{2}\log^2(2).\tag{2}$$ On the other hand $$ \sum_{n\geq 1}\frac{1}{n^2 2^n} = \text{Li}_2\left(\frac{1}{2}\right)=\frac{\pi^2}{12}-\frac{\log^2(2)}{2}\tag{3} $$ follows from the dilogarithm reflection formula, and by summing $(2)$ and $(3)$ the identity we usually derive from Euler's acceleration method $$ \sum_{n\geq 1}\frac{H_n}{n 2^n} = \frac{\pi^2}{12}\tag{4} $$ easily follows.
Late addendum (2020). It is interesting to notice that $$ H_n = \gamma-\lim_{m\to 0^+} \frac{d}{dm}\left(\frac{1}{(n+1)_m}\right) \tag{5}$$ where $(n+1)_m$ is the rising Pochhammer symbol and $\gamma=-\Gamma'(1)$ is the Euler-Mascheroni constant.
We have $$\begin{eqnarray*} \sum_{n\geq 1}\frac{1}{n 2^n (n+1)_m}&=&\sum_{n\geq 1}\frac{1}{2^n(n)_{m+1}}=\sum_{n\geq 1}\frac{\Gamma(n)}{2^n \Gamma(n+m+1)}\\&=&\frac{1}{\Gamma(m+1)}\sum_{n\geq 1}\frac{B(n,m+1)}{2^n}\\&=&\frac{1}{\Gamma(m+1)}\int_{0}^{1}\sum_{n\geq 1}\frac{(1-x)^m x^{n-1}}{2^n}\,dx\\&=&\frac{1}{\Gamma(m+1)}\int_{0}^{1}\frac{(1-x)^m}{2-x}\,dx\\&=&\frac{1}{\Gamma(m+1)}\int_{0}^{1}\frac{x^m}{1+x}\,dx\end{eqnarray*} \tag{6}$$ hence by the chain rule and differentiation under the integral sign $$\begin{eqnarray*} \sum_{n\geq 1}\frac{H_n}{n 2^n}&=&\gamma\log(2)-\gamma\int_{0}^{1}\frac{dx}{1+x}-\int_{0}^{1}\frac{\log(x)}{1+x}\,dx\\&=&\int_{0}^{1}-\log(x)\sum_{n\geq 0}(-1)^{n} x^n\,dx=\sum_{n\geq 0}\frac{(-1)^n}{(n+1)^2}\\&=&\sum_{n\geq 1}\frac{(-1)^{n+1}}{n^2}=\zeta(2)-2\sum_{n\geq 1}\frac{1}{(2n)^2}=\left(1-\frac{1}{2}\right)\zeta(2)=\color{red}{\frac{\pi^2}{12}}\end{eqnarray*} \tag{7}$$ "undoing" Euler acceleration method.