Mira la respuesta final. La respuesta final te dice que esta serie no se puede calcular de forma sencilla. Aquí propongo una solución aunque no es una buena solución:
$$\Sigma_{n=1}^\infty \frac1{(n+a)(n+b)}=\frac{H_a-H_b}{a-b}$$
Dónde $H_x$ es número armónico de $x$ .
(También se puede sustituir por la función digamma. Aunque $\psi_0(x)\neq H_x$ el resultado sigue siendo el mismo).
Informática $H_x$ tampoco es fácil. Es tan difícil como el cálculo de integrales. Pero si tienes su tabla, puedes calcular la serie.
En este caso, si usted sabe $H_{\frac16}$ y $H_{\frac23}$ puede resolver esta serie (¡condición impar!).
$$\Sigma_{n=0}^\infty \frac{(-1)^n}{(3n+1)}=\Sigma_{n=0}^\infty (\frac1{(6n+1)}-\frac1{(6n+4)})=\Sigma_{n=0}^\infty \frac3{(6n+1)(6n+4)}=1-\frac14+\frac3{36}\Sigma_{n=1}^\infty \frac1{(n+\frac16)(6n+\frac46)}=1-\frac14+\frac1{12}\times\frac{H_{\frac46}-H_{\frac16}}{\frac46-\frac16}=\frac34+\frac16(H_{\frac23}-H_{\frac16})$$
Necesito conocer estos dos números armónicos (la gran desventaja de este método):
$$H_{\frac23}=\frac32+\frac\pi{2 \sqrt3}-\frac{3 \log 3}2$$
$$H_{\frac16}=6-\frac{\sqrt3 \pi}2-\frac{3 \log 3}2-\log4$$
$$H_{\frac23}-H_{\frac16}=-\frac92+\frac{2\pi \sqrt3}3+2\log 2$$
$$\Sigma_{n=0}^\infty \frac{(-1)^n}{(3n+1)}=\frac34+\frac16 (-\frac92+\frac{2\pi \sqrt3}3+2\log 2)$$
$$\therefore ~~~\Sigma_{n=0}^\infty \frac{(-1)^n}{(3n+1)}=\frac{\pi \sqrt3}9+\frac{\log2}3$$