He aquí otro enfoque, ligeramente más largo. Utilizará una integral doble para evaluar la suma.
Observando que \int_0^1 x^{n - 1} \, dx = \frac{1}{n} \quad \text{and} \quad \int_0^1 y^n \, dy = \frac{1}{n + 1}, la suma puede escribirse como \sum_{n = 1}^\infty \frac{H_n}{n(n + 1)} = \int_0^1 \int_0^1 \frac{1}{x} \sum_{n = 1}^\infty H_n (xy)^n \, dx dy.\tag1 Aquí se ha cambiado el orden de la suma con la integral doble y se puede hacer debido al teorema de convergencia dominada.
Utilizando el función generadora para el número armónico, a saber -\frac{\log(1-x)}{1 - x}~=~\sum_{n=1}^\infty H_n x^n, la suma infinita que aparece dentro de la integral en (1) puede escribirse como \sum_{n = 1}^\infty \frac{H_n}{n(n + 1)} = -\int_0^1 \int_0^1 \frac{\ln (1 - xy)}{x(1 - xy)} \, dx dy. Ejecución de una sustitución de x \mapsto x/y da \begin{align} \sum_{n = 1}^\infty \frac{H_n}{n(n + 1)} &= -\int_0^1 \int_0^y \frac{\ln (1 - x)}{x (1 - x)} \, dx dy\\ &= -\int_0^1 \int_0^y \left [\frac{\ln (1 - x)}{x} + \frac{\ln (1 - x)}{1 - x} \right ] \, dx \, dy\\ &= \int_0^1 \operatorname{Li}_2 (y) \, dy + \frac{1}{2} \int_0^1 \ln^2 (1 - y) \, dy\\ &= I_1 + \frac{1}{2} I_2. \end{align} Aquí \operatorname{Li}_2 (y) es el función dilogaritmo .
Para la primera de estas integrales, escribiendo la función dilogaritmo en términos de su representación en serie, tenemos \begin{align} I_1 &= \int_0^1 \sum_{n = 1}^\infty \frac{y^n}{n^2} \, dy\\ &= \sum_{n = 1}^\infty \frac{1}{n^2} \int_0^1 y^n \, dy\\ &= \sum_{n = 1}^\infty \frac{1}{n^2(n + 1)}\\ &= \sum_{n = 1}^\infty \left [\frac{1}{n^2} + \frac{1}{n + 1} - \frac{1}{n} \right ]\\ &= \sum_{n = 1}^\infty \frac{1}{n^2} - \sum_{n = 1}^\infty \left (\frac{1}{n} - \frac{1}{n + 1} \right )\\ &= \frac{\pi^2}{6} - 1, \end{align} donde el primero de la serie es el conocido problema de Basilea, mientras que el segundo es telescópico.
Para la segunda de las integrales, defina I(a) = \int_0^1 (1 - y)^a \, dy, \quad a > - 1. Observe que I''(0) = \int_0^1 \ln^2 ( 1 - y) \, dy, donde la derivada es con respecto a a . Ahora, como I(a) = \operatorname{B} (1, a + 1) = \frac{\Gamma (a + 1)}{\Gamma (a + 2)} = \frac{1}{a + 1}, donde \operatorname{B}(x,y) es el Función beta se tiene fácilmente I''(a) = \frac{2}{(a + 1)^3}, que se puede encontrar en I''(0) = I_2 = 2 .
Así, \sum_{n = 1}^\infty \frac{H_n}{n (n + 1)} = \frac{\pi^2}{6} - 1 + \frac{1}{2} \cdot 2 = \frac{\pi^2}{6}.