Todos sabemos que las sumas de la forma $\displaystyle{\sum\frac1{(n+a)(n+b)}}$ o $\displaystyle{\sum\frac1{(n+a)(n+b)(n+c)}}$, etc., son telescópicas en la naturaleza, y su solución es trivial: para el entero de los valores de un, b, y c, que es ! Pero, ¿y si un, b, y c son no enteros? Entonces, ¿qué? De repente, las cosas no son tan simples y triviales ya, y el anteriormente banal y engaño tame de problemas toma en nuevo giro, ganando totalmente inesperado dimensiones de profundidad, el significado y la comprensión. La clave, cuando se trata con estas sorprendentes giros de los acontecimientos, es solo dar un paso atrás, y tratar de analizar y reformular el viejo y gastado de solución en términos de lo que podría ser útil o relevante para eludir los obstáculos planteados por el recién encontrado situación.
Por ejemplo, no es difícil demostrar que $\displaystyle{\sum_{n=1}^\infty\frac1{(n+a)(n+b)}}=\frac{H_a-H_b}{a-b}$ o que $$\sum_{n=1}^\infty\frac1{(n+a)(n+b)(n+c)}=\frac{(a-b)\cdot H_c\ +\ (b-c)\cdot H_a\ +\ (c-a)\cdot H_b}{(a-b)(b-c)(c-a)}$$ por los valores naturales de la un, b, c. De hecho, el usuario de la Variable Aleatoria ya ha demostrado en este hilo, aunque ahora me doy cuenta de que mi fórmula inicial es ligeramente equivocado, en el sentido de que la suma debe comenzar en $0$, $\gamma$ presente, o en $1$, $\gamma$ ausente. Ahora, la cuestión es: ¿cómo redefinir $H_m$ , así como ser capaz de extender su significado no-natural argumentos. Natural de argumentos, tenemos la conocida fórmula de $H_m=\displaystyle\sum_{k=1}^m\frac1k$ . Ahora, tomemos la función simple de $f_k(x)=\displaystyle\frac{x^k}k$ , y el aviso de que $f_k'(x)=x^{k-1}$ . Y ya sabemos que $\displaystyle{\sum_{k=1}^mx^{k-1}=\frac{1-x^m}{1-x}}$ , por fin, este nos permite concluir que $\displaystyle{H_m=\int_0^1\frac{1-x^m}{1-x}dx}$ , que, a diferencia de la fórmula anterior, se puede extender fácilmente a los no-natural de los argumentos. Esto ya ha sido hecho por Euler dos siglos y medio, así que es apenas un nuevo territorio. Ahora, en nuestro caso, $\{a,b,c\}=\left\{\frac13,\frac23,1\right\}$ , y el valor de $H_1$ es $1$, por lo que todo lo que queda por hacer es calcular $H_\frac13$ y $H_\frac23$ utilizando la fórmula anterior, ya que $$\sum_{n=0}^\infty\frac1{(3n+1)(3n+2)(3n+3)}=\frac16+\frac1{27}\cdot\sum_{n=1}^\infty\frac1{\left(n+\frac13\right)\left(n+\frac23\right)(n+1)}$$
$$H_\frac13=\int_0^1\frac{1-\sqrt[3]x}{1-x}dx=\int_0^1\frac{1-t}{1-t^3}d\left(t^3\right)=\int_0^1\frac1{1+t+t^2}d\left(t^3\right)=\int_{\sqrt[3]0}^{\sqrt[3]1}\frac{3\,t^2}{1+t+t^2}dt=$$
$$=3\int_0^1\left(1-\frac{1+t}{1+t+t^2}\right)dt=3\,\bigg[\int_0^11\cdot dt-\tfrac12\int_0^1\frac{2+2t}{1+t+t^2}dt\bigg]=$$
$$=3\,\bigg[t|_0^1-\tfrac12\bigg(\int_0^1\frac1{1+t+t^2}dt+\int_0^1\frac{1+2t}{1+t+t^2}dt\bigg)\bigg]=$$
$$=3\,\bigg[1-\tfrac12\int_0^1\frac1{\left(t+\frac12\right)^2+\frac34}dt-\tfrac12\cdot\ln\left(1+t+t^2\right)_0^1\bigg]=$$
$$=3\,\bigg[1-\tfrac12\cdot\frac1{\sqrt3/2}\cdot\arctan\bigg(\frac{t+\frac12}{\sqrt3\big/2}\bigg)_0^1-\frac{\ln3}2\bigg]=3-\frac\pi{2\sqrt3}-\tfrac32\ln3.$$
Lo mismo para $\displaystyle{H_{\frac23}=\tfrac32+\frac\pi{2\sqrt3}-\tfrac32\ln3}$ . Entonces, sustituyendo estos valores en la fórmula original de arriba, por fin llegamos al resultado deseado $\displaystyle{I=\frac{\pi\sqrt3-3\ln3}{12}}$ .