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 ?

I think it is reasonable to start with:
$$\sum_{k=1}^{+\infty}\frac{H_k^{(2)}H_k}{k^2}=\sum_{k=1}^{+\infty}\frac{H_k}{k^4}+\sum_{k=1}^{+\infty}\frac{H_k}{k^2}\sum_{1\leq j< k}\frac{1}{j^2},\tag{1}$$
that leads to:
$$\sum_{k=1}^{+\infty}\frac{H_k^{(2)}H_k}{k^2}=\left(\sum_{k=1}^{+\infty}\frac{H_k}{k^2}\right)\left(\sum_{j=1}^{+\infty}\frac{1}{j^2}\right)-\sum_{k=1}^{+\infty}\frac{1}{k^2}\sum_{1\leq j< k}\frac{H_j}{j^2},\tag{2}$$
Now since:
$$\operatorname{Li}_2(x)+\frac{\log^2(1-x)}{2}=\sum_{k=1}^{+\infty}\frac{H_k}{k}x^k,\tag{3}$$
$$\frac{\log^2(1-x)}{2}=\sum_{k=1}^{+\infty}\frac{H_{k-1}}{k}x^k,\tag{4}$$
follows. By dividing by $x$ and integrating between $0$ and $1$ we get:
$$\sum_{k=1}^{+\infty}\frac{H_{k-1}}{k^2}=\frac{1}{2}\int_{0}^{1}\frac{\log^2(x)}{1-x}dx=\frac{1}{2}\int_{0}^{+\infty}\frac{u^2}{e^u-1}du=\zeta(3),\tag{5}$$
so:
$$\sum_{k=1}^{+\infty}\frac{H_k^{(2)}H_k}{k^2}=2\zeta(2)\zeta(3)-\sum_{k=1}^{+\infty}\frac{1}{k^2}\sum_{1\leq j< k}\frac{H_j}{j^2}.\tag{6}$$
For the last term consider:
$$-\frac{\log(1-xy)}{y(1-xy)}=\sum_{k=1}^{+\infty}H_k x^k y^{k-1}, \tag{7}$$
multiply both terms by $-\log(y)$ and integrate between $0$ and $1$ with respect to $y$:
$$\int_{0}^{1}\frac{\log(y)\log(1-xy)}{y(1-xy)}dy = \sum_{k=1}^{+\infty}\frac{H_k}{k^2}x^k.\tag{8}$$ 
Multiplying both sides by $-\frac{\log x}{1-x}$ and integrating between $0$ and $1$ with respect to $x$ should do the trick. For the last part it is only required to find an appropriate birational diffeomorphism of the unity square that puts the integral in a nicer form - a sort of "reverse Viola-Rhin method".
