<s>[INST] Triple Euler sum result $\sum_{k\geq 1}\frac{H_k^{(2)}H_k }{k^2}=\zeta(2)\zeta(3)+\zeta(5)$ In the following thread
I arrived at the following result 
$$\sum_{k\geq 1}\frac{H_k^{(2)}H_k }{k^2}=\zeta(2)\zeta(3)+\zeta(5)$$
Defining 
$$H_k^{(p)}=\sum_{n=1}^k \frac{1}{n^p},\,\,\, H_k^{(1)}\equiv H_k $$
But, it was after long evaluations and considering many variations of product of polylogarithm integrals.
I think there is an easier approach to get the solution, any ideas ?
 [/INST] Consider the integral $$I= - \int_0^1 \frac{\ln(1-x^2)}{\sqrt{1-x^2}} (\sin^{-1} x)^4 \,dx.$$
Since $$(\sin^{-1} x)^4 = \frac32 \sum_{n=1}^{\infty} \cfrac{2^{2n} H_{n-1}^{(2)}}{n^2 \binom{2n}{n}} \,x^{2 n} \tag{1}$$
and 
$$-\int_0^1 \frac{\ln(1-x^2)}{\sqrt{1-x^2}} x^{2n}\,dx= \frac{\pi}{2} \binom{2n}{n} \frac{(H_n + 2\ln2)}{2^{2n}}, \tag{2}$$
we have 
$$\begin{align*}  
&I= - \frac32 \sum_{n=1}^{\infty} \cfrac{2^{2n} H_{n-1}^{(2)}}{n^2 \binom{2n}{n}} \int_0^1 \frac{\ln(1-x^2)}{\sqrt{1-x^2}} x^{2n}\,dx
\\&= \frac{3 \pi}{4} \sum_{n=1}^{\infty} \frac{H_{n-1}^{(2)}}{n^2} ( H_n +2 \ln2 )
\\& = \frac{\pi^5}{80} \ln2 + \frac{3 \pi}{4} \sum_{n=1}^{\infty} \frac{H_{n-1}^{(2)}\,H_n}{n^2}.
\end{align*}$$
However, substituting $x\mapsto \sin x$ and employing the fourier expansion of $\ln \cos x$:
$$\begin{align*}
& I= -2 \int_0^{\pi/2} x^4 \, \ln\cos x\, dx
\\&= 2 \int_0^{\pi/2} x^4 \left(\ln2 + \sum_{n=1}^{\infty} \frac{(-1)^n \cos(2 x n)}{n} \right)dx
\\&=  \frac{\pi^5}{80}\ln2  + 2 \sum_{n=1}^{\infty} \frac{(-1)^n}{n} \int_0^{\pi/2} x^4 \cos(2 x n) dx
\\&= \frac{\pi^5}{80}\ln2 + 2 \sum_{n=1}^{\infty} \frac{(-1)^n}{n} \frac{(-1)^n}{n^2}\left(\frac{\pi^3}{8}-\frac{3 \pi}{4 n^2}\right)
\\&= \frac{\pi^5}{80}\ln2 + \frac{\pi^3}{4}\zeta(3) - \frac{3 \pi}{2} \zeta(5). \end{align*}$$
Therefore,

$$\sum_{n=1}^{\infty} \frac{H_{n-1}^{(2)}\,H_n}{n^2} = 2\zeta(2)\,\zeta(3)-2\zeta(5).$$

Finish off using Euler's formula for $\sum H_n/n^q $.

Notes.
You may find a proof of $(1)$ here, and $(2)$ is just the derivative of a beta function. The swap of the sum and the integral should be justified.
I found this proof while exploring series involving $H_n^{(2)}$. Using the same method, I also obtained the following
related results:
$$\sum_{n=1}^{\infty} \frac{H_{n-1}^{(2)}\,H_{2n}}{n^2} =\frac{11}{4}\zeta(2)\,\zeta(3)-\frac{47}{16}\zeta(5) \tag{3}$$
$$\sum_{n=1}^{\infty} \frac{H_{n-1}^{(2)}\,H_{n}^2}{n^2} = 4 \zeta(3)^2 - \frac{5}{8} \zeta(6) \tag{4}$$
$$\sum_{n=1}^{\infty} \frac{H_n \left(H_{n-1}^{(2)2}-H_{n-1}^{(4)}\right)}{n^2} = 3\,\zeta(3)\,\zeta(4)-4\,\zeta(2)\,\zeta(5)+4\,\zeta(7) \tag{5}$$
and others.
</s>