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}.$$