En la introducción sobre la pregunta de cómo se puede ir sobre la computación $$\sum_{n=1}^\infty \frac{1}{3n(3n-1)}.$$
Puedo contribuir con un asintótica de expansión y una serie infinita de esta discusión, mostrado y demostrado a continuación. Además de publicar esta respuesta también estoy planteando dos preguntas I. se Puede verificar esta expansión, utilizando una prueba técnica? II. Podría incluso ser una referencia a la siguiente fórmula en algún lugar?
Esta es la prueba, que es de procedimiento, con el resultado final. Supongamos que buscamos para evaluar
$$T(q) = \sum_{n\ge 1} \frac{1}{qn(qn-1)} =
\sum_{n=1}^p \frac{1}{qn(qn-1)} + \sum_{n\ge 1} \frac{1}{(qn+pq)(qn+pq-1)} \\=
\sum_{n=1}^p \frac{1}{qn(qn-1)} + \frac{1}{p^2}\sum_{n\ge 1} \frac{1}{(n+p)(n+(pq-1)/q)}$$
con $p, q\ge 2$ dos números enteros.
La suma plazo, la llame $S(x)$, es armónica y puede ser evaluado por la inversión de sus Mellin transformar.
Recordar el Mellin transformar la identidad para armónica sumas con la base de la función $g(x)$, que es
$$\mathfrak{M}\left(\sum_{k\ge 1}\lambda_k g(\mu_k x); s\right) =
\left(\sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s}\right) g^*(s)$$
donde $g^*(s)$ es la transformada de Mellin $g(x).$
En el presente caso tenemos
$$\lambda_k = 1, \quad
\mu_k = k \quad \text{y} \quad
g(x) = \frac{1}{(x+p)(x+(pq-1)/q)}.$$
La transformada de Mellin $g(x)$ es
$$\int_0^\infty \frac{1}{(x+p)(x+(pq-1)/q)} x^{m-1} dx.$$
Elegimos $p\ge 2$ de modo que la base de la función no tiene un polo en $x=-1$, lo que provocaría la expansión sobre el infinito que hemos llegado a obtener no convergen. (Como se evalúa la armónica suma de $x=1$ esto sería situado en el límite entre los dos discos alrededor de cero y sobre el infinito, causando una evaluación en $x=1$ a fallo).
Evaluamos la Mellin transformar siguiente. Utilizar un contorno circular para obtener
$$ g^*(s) (1 - e^{2\pi i (s-1)}) \\=
2\pi i \left(
\operatorname{Res}\left(\frac{x^{m-1}}{(x+p)(x+\frac{pq-1}{p})}; s=-p\right) +
\operatorname{Res}\left(\frac{x^{m-1}}{(x+p)(x+\frac{pq-1}{p})}; s=-\frac{pq-1}{q}\right)
\right) \\=
2\pi i \left(-q(-p)^{m-1} + p(-(pq-1)/q)^{m-1}\right)
= 2\pi i e^{\pi i (s-1)} p \left(((pq-1)/q)^{m-1}-p^{m-1}\right).$$
Esto implica que
$$g^*(s) = 2\pi i q \frac{-e^{\pi i s}}{1 - e^{2\pi i s}}
\left(\left(\frac{pq-1}{q}\right)^{m-1} - p^{m-1}\right)\\=
-2\pi i q \frac{1}{e^{-\pi i s} - e^{\pi i s}}
\left(\left(\frac{pq-1}{q}\right)^{m-1} - p^{m-1}\right) \\=
q\pi\frac{2}{e^{\pi i s} - e^{-\pi i s}}
\left(\left(\frac{pq-1}{q}\right)^{m-1} - p^{m-1}\right) =
q\frac{\pi}{\sin(\pi s)} \left(\left(\frac{pq-1}{q}\right)^{m-1} - p^{m-1}\right).$$
De ello se desprende que la transformada de Mellin la suma plazo para $S(x)$ es dada por
$$ Q(s) =
q \frac{\pi}{\sin(\pi s)}
\left(\left(\frac{pq-1}{q}\right)^{m-1} - p^{m-1}\right) \zeta(s)
\quad\text{porque}\quad \sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s} = \zeta(s).$$
Finalmente podemos invertir la Mellin transformar con la inversión integral
$$\frac{1}{2\pi i}\int_{3/2-i\infty}^{3/2+i\infty} P(s)/x^s ds$$
para obtener una expansión sobre infinidad de $S(x)$ comenzando con la pole en $s=2$ y conseguir
$$ S(x) \sim
- q \sum_{m\ge 2} (-1)^m
\left(\left(\frac{pq-1}{q}\right)^{m-1} - p^{m-1}\right) \frac{\zeta(m)}{x^m}.$$
La configuración de $x=1$ y sustituyendo en la ecuación original por fin tenemos
$$ T(q) = \sum_{n=1}^p \frac{1}{qn(qn-1)}
- \frac{1}{q} \sum_{m\ge 2} (-1)^m
\left(\left(\frac{pq-1}{q}\right)^{m-1} - p^{m-1}\right) \zeta(m).$$
Anexo Wed Apr 2 20:14:50 CEST 2014. El siguiente MSE enlace se muestra cómo trabajar con distintos Mellin transforma donde hay un poste en el eje real positivo (incluso).