I'm having a little trouble calculating the double series given by the proposition below:
$$\int_0^1\frac{\left[\text{Li}_2(x)\right]^2}{x}\mathrm dx=\sum_{n=1}^\infty\sum_{k=1}^\infty\frac{1}{n^2k^2(n+k)}=2\zeta(2)\zeta(3)-3\zeta(5).\tag{1}$$
In fact, this integral has already been evaluated here and here, but i have no intention of evaluating it by integration or by using harmonic series. My interest is to attack the double sum in $(1)$ only using the manipulation of sums. Having clarified my interest in this topic, i naively proceeded as follows: $$\begin{align*}\sum_{n=1}^\infty\sum_{k=1}^\infty\frac{1}{n^2k^2(n+k)}&=\sum_{n,k\in\mathbb{N}_1}\frac{1}{n^2k^2(n+k)}\\ &=\sum_{\substack{{n,k\in\mathbb{N}_1}\\{n=k}}}\frac{1}{n^2k^2(n+k)}+\sum_{\substack{{n,k\in\mathbb{N}_1}\\{n>k}}}\frac{1}{n^2k^2(n+k)}+\sum_{\substack{{n,k\in\mathbb{N}_1}\\{n<k}}}\frac{1}{n^2k^2(n+k)}\\ &=\frac{\zeta(5)}{2}+\sum_{\substack{{n,k\in\mathbb{N}_1}\\{n>k}}}\frac{1}{n^2k^2(n+k)}+\sum_{\substack{{n,k\in\mathbb{N}_1}\\{n<k}}}\frac{1}{n^2k^2(n+k)}\tag{a}\end{align*}$$
Here's my difficulty with that approach! How can i solve the last two double sums in $(\text a)$ with the indexes $ n <k $ and $ n> k $ ?. Studying other posts (links at the end of the question) related to my question, i found here another possible approach using the fact that: $ \frac{1}{n^2k^2}=\frac{1}{(n+k)^2n^2}+\frac{1}{(n+k)^2k^2}+\frac{2}{(n+k)^3k}+\frac{2}{(n+k)^3n}$, then multiplying by $\frac{1}{n+k}$ on both sides and taking the double sum, we get: $$\sum_{n,k\in\mathbb{N}_1}\frac{1}{n^2k^2(n+k)}=2\sum_{n,k\in\mathbb{N}_1}\frac{1}{(n+k)^3n^2}+4\sum_{n,k\in\mathbb{N}_1}\frac{1}{(n+k)^4n}.\tag{b} $$
According to this article, $(\text b)$ is equal to $T(2,2,1)=2T(2,0,3)+4T(1,0,4)$, where $T(a,b,c)$ is Tornheim’s double series defined by $T(a,b,c)=\sum_{n,m=1}\frac{1}{n^am^b(n+m)^c}$. In this second approach, the expression in $(\text b)$ seems to me a little less aggressive, but if i apply the concepts of the mentioned article, i would feel a little difficult to understand, since it starts with already established formulas. In short, two questions hang in my mind: How to work $(\text a)$ and $(\text b)$ without using integrals or harmonic numbers?
Related links:
A. $\displaystyle \sum_{m,n=1\, m\neq n}^\infty{\frac{m^2+n^2}{mn(m^2-n^2)^2}}$
B. $\displaystyle \sum_{k=1}^{\infty} \sum_{n=1}^{\infty} \frac{1}{n^2k^2(n+k)^2}$