# Evaluating $\sum_{i=1}^{\infty} \sum_{j=1}^{\infty} \frac{1}{ij(i+j)^2}$

I was playing around with double sums and encountered this problem: Evaluate $$\sum_{i=1}^{\infty} \sum_{j=1}^{\infty} \frac{1}{ij(i+j)^2}$$ It looks so simple I thought it must have been seen before, but after using ApproachZero I only found extremely similar variants of the problem:

$$\sum_{j=1}^{\infty} \sum_{k=1}^{\infty} \frac{1}{j^2(j+k)^2}$$ here

$$\sum_{k=1}^{\infty} \sum_{n=1}^{\infty} \frac{1}{n^2k^2(n+k)^2}$$ here

$$\sum_{n,k \in \mathbb{N_1}} \frac{1}{n^2k^2(n+k)}$$ here

Some general facts about $$T(r,s,t)=\sum_{u,v=1}^{\infty} \frac{1}{u^rv^s(u+v)^t}$$ here

While the last link looks promising, I don't believe that it answers my specific question(but I haven't read the paper fully, only skimmed).

I have an answer to this question(shared Q&A style), but I am curious as to whether there is another way to evaluate the sum. I have been experimenting with double sums for a few weeks now, so many techniques to attack "simple" sums such as this would be much appreciated! In terms of level, I have officially finished a Calculus 2 course but I have self-taught myself Calculus 3 and Complex Analysis, and worked through a few "integration problems" books so any further integral techniques would be fine.

• Here's one way to solve it, assuming you know how to calculate $\zeta(4)$ (it's in the calculation of the quantity named $J$ in ysharifi's post): Integral of squared infinite sum (AoPS thread) - found using Approach$0$. Sep 1, 2023 at 17:29
• Interesting, thank you! I'll check it out Sep 1, 2023 at 19:07

My own answer to the above question.

We can do this by analyzing a particular $$n=i+j$$. Instead of summing in the row/column way(when the double sum is written out) we sum by the diagonals from the top row to the left column as shown below

where each box corresponds to a particular choice of i,j. Now, because $$n=i+j$$ will hit all integers from $$2$$ to $$\infty$$, we can rewrite the summation as $$S = \sum_{n=2}^{\infty} \frac{1}{n^2}f(n)$$ where $$f(n)=\sum_{i=1}^{n-1} \frac{1}{i(n-i)}$$ where we set $$j$$ to be $$n-i$$ because $$i+j=n$$. This is the "sum by diagonals" representation. Now, we can find an explicit formula for $$f(n)$$. Note that $$\frac{1}{i(n-i)}=\frac{1}{n}(\frac{1}{n-i}+\frac{1}{i})$$ so, we have that $$f(n)=\frac{1}{n} \sum_{i=1}^{n-1} \frac{1}{i}+\frac{1}{n-i}=$$ $$\frac{1}{n} (\sum_{i=1}^{n-1} \frac{1}{i}+ \sum_{i=1}^{n-1} \frac{1}{n-i})=$$ $$\frac{1}{n}(H_{n-1}+\sum_{i=1}^{n-1} \frac{1}{n-i})=$$ $$\frac{1}{n}(H_{n-1}+\sum_{i=1}^{n-1} \frac{1}{i})=$$ $$\frac{2 H_{n-1}}{n}$$ where the second-to-last equality follows from the reflection formula $$\sum_{i=a}^{b} f(i)=\sum_{i=a}^{b} f(a+b-i)$$. From this, we have that $$S=\sum_{n=2}^{\infty} \frac{1}{n^2} \cdot \frac{2 H_{n-1}}{n}=$$ $$2 \sum_{n=2}^{\infty} \frac{H_{n-1}}{n^3}=$$ $$2 \sum_{n=2}^{\infty} \frac{H_n-\frac{1}{n}}{n^3}=$$ $$2(\sum_{n=2}^{\infty} \frac{H_n}{n^3} - \sum_{n=2}^{\infty} \frac{1}{n^4})$$ Here, by noting that $$H_n=\frac{1}{n}$$ for $$n=1$$, we also have that $$S=2(\sum_{n=1}^{\infty} \frac{H_n}{n^3} - \sum_{n=1}^{\infty} \frac{1}{n^4})$$ as the $$n=1$$ terms cancel each other out. Now, the second summation is clearly $$\zeta(4)=\frac{\pi^4}{90}$$, so we have that $$S=2(\sum_{n=1}^{\infty} \frac{H_n}{n^3} - \frac{\pi^4}{90})=$$ $$2 \sum_{n=1}^{\infty} \frac{H_n}{n^3} - \frac{\pi^4}{45}$$ Now, all that is left is to evaluate $$\sum_{n=1}^{\infty} \frac{H_n}{n^3}$$. To do this, we use the method outlined in Problem 3.58 of "Limits, Series, and Fractional Part Integrals". Recall the Dilogarithm function $$\mathrm{Li_2}(z)=\sum_{n=1}^{\infty} \frac{z^n}{n^2} = - \int_{0}^{z} \frac{\ln(1-t)}{t} dt$$ For example, note that $$\mathrm{Li_2}(1)=\zeta(2)=\frac{\pi^2}{6}$$ and $$\mathrm{Li_2}(0)=0$$. It is a known result(proved in an earlier problem in the book) that $$\sum_{n=1}^{\infty} \frac{H_n}{n^3} = - \sum_{n=1}^{\infty} \frac{1}{n^2} \int_{0}^{1} x^{n-1} \ln(1-x) dx$$ and we know that $$- \sum_{n=1}^{\infty} \frac{1}{n^2} \int_{0}^{1} x^{n-1} \ln(1-x) dx=- \int_{0}^{1} \frac{\ln(1-x)}{x} (\sum_{n=1}^{\infty} \frac{x^n}{n^2}) dx=$$ $$- \int_{0}^{1} \frac{\ln(1-x)}{x} \mathrm{Li_2}(x) dx$$ Now, sub $$u=\mathrm{Li_2}(x), du = \frac{-\ln(1-x)}{x} dx$$, and using the special values of $$\mathrm{Li_2}(x)$$ above for the bounds, we see that $$\sum_{n=1}^{\infty} \frac{H_n}{n^3} = \int_{0}^{\frac{\pi^2}{6}} u du=$$ $$\frac{u^2}{2} \Big|_0^{\frac{\pi^2}{6}}=$$ $$\frac{\pi^4}{72}$$ So, we finish with $$S = 2 \cdot \frac{\pi^4}{72} - \frac{\pi^4}{45}=$$ $$\boxed{\frac{\pi^4}{180}}=$$ $$\boxed{\frac{1}{2} \zeta(4)}$$

• Very elegant approach for sure ! (+1) Sep 2, 2023 at 8:35

$$\sum_{j=1}^{\infty} \frac{1}{i\,j\,(i+j)^2}=\frac{\gamma+\psi ^{(0)}(i+1)-i \psi ^{(1)}(i+1) }{i^3}$$

$$S_p=\sum_{i=1}^{p}\frac{\gamma}{i^3}=\gamma\, H_p^{(3)}$$

$$\sum_{i=1}^{p}\frac{\psi ^{(1)}(i+1)}{i^2}=\frac{1}{120} \left(\pi^4-60 \psi ^{(1)}(p+1)^2+10 \psi^{(3)}(p+1)\right)$$

$$\sum_{i=1}^{\infty}\frac{\psi ^{(0)}(i+1)}{i^3}=\frac{\pi ^4}{72}-\gamma \zeta (3)$$

Expanding as series

$$S_p=\frac{\pi ^4}{180}+\frac{1-\gamma }{2 p^2}+\frac{3 \gamma -4}{6 p^3}+\frac{13-6 \gamma }{24 p^4} +O\left(\frac{1}{p^5}\right)$$ which is in a relative error smaller than $$1.0$$% as soon as $$p>17$$ and smaller than $$0.1$$% as soon as $$p>65$$.

• Very interesting technique! When trying to evaluate the sum, I got to the first line of your solution before I gave up; splitting it up and using the series expansion is super smart! Sep 2, 2023 at 15:43