La suma infinita en la respuesta de Chen Wang, es decir, $ \displaystyle \sum_{n=1}^{\infty} \frac{H_{4n}}{n^{2}}$ puede evaluarse mediante integración de contornos considerando la función $$f(z) = \frac{\pi \cot(\pi z) [\gamma + \psi(-4z)]}{z^{2}}, $$
donde $\psi(z)$ es el función digamma y $\gamma$ es el Constante de Euler-Mascheroni .
La función $f(z)$ tiene polos de orden $2$ en los enteros positivos, polos simples en los enteros negativos, polos simples en los cuartos de enteros positivos y un polo de orden $4$ en el origen.
La función $\psi(-4z)$ tiene polos simples en los semienteros positivos, pero se anulan con los ceros de $\cot( \pi z)$ .
Consideremos ahora un cuadrado en el plano complejo (llamémoslo $C_{N}$ ) con vértices en $\pm (N + \frac{1}{2}) \pm i (N +\frac{1}{2})$ .
En los lados del cuadrado, $\cot (\pi z)$ está uniformemente acotada.
Y cuando $z$ es de gran magnitud y no se encuentra en el eje real positivo, $\psi(-4z) \sim \ln(-4z)$ .
Así que $ \displaystyle \int_{C_{N}} f(z) \ dz $ desaparece cuando $N \to \infty$ a través de los enteros positivos.
Por lo tanto,
$$\sum_{n=1}^{\infty} \text{Res} [f(z), n] + \sum_{n=1}^{\infty} \text{Res}[f(z),-n] + \text{Res}[f(z),0] + \sum_{n=0}^{\infty} \text{Res}\Big[f(z), \frac{2n+1}{4} \Big] =0 .$$
Para determinar los residuos, necesitamos las siguientes expansiones de Laurent.
En los enteros positivos,
$$ \gamma + \psi (-4z) = \frac{1}{4} \frac{1}{z-n} + H_{4n} + \mathcal{O}(z-n) $$
y
$$ \pi \cot (\pi z) = \frac{1}{z-n} + \mathcal{O}(z-n) .$$
En el origen,
$$ \gamma+ \psi(-4z) = \frac{1}{4z} -4 \zeta(2) z -16 \zeta(3) z^{2} + \mathcal{O}(z^{3})$$
y $$ \pi \cot (\pi z) = \frac{1}{z} - 2 \zeta(2) z + \mathcal{O}(z^{3}) .$$
Y en los cuartos enteros positivos,
$$ \gamma + \psi(-4z) = \frac{1}{4} \frac{1}{z-\frac{2n+1}{4}} + \mathcal{O}(1)$$
y
$$ \pi \cot (\pi z) = (-1)^{n} \pi + \mathcal{O}\Big(z- \frac{2n+1}{4} \Big) .$$
Entonces en los enteros positivos,
$$f(z) = \frac{1}{z^{2}} \Big( \frac{1}{4} \frac{1}{(z-n)^{2}} + \frac{H_{4n}}{z-n} + \mathcal{O}(1) \Big), $$
lo que implica
$$\begin{align} \text{Res} [f(z),n] &= \text{Res} \Big[ \frac{1}{4z^{2}} \frac{1}{(z-n)^{2}} , n \Big] + \text{Res} \Big[ \frac{1}{z^{2}} \frac{H_{4n}}{z-n}, n \Big] \\ &= - \frac{1}{2n^{3}} + \frac{H_{4n}}{n^{2}} .\end{align}$$
En los enteros negativos,
$$ \text{Res}[f(z),-n] = \frac{\gamma + \psi(4n)}{n^{2}} = \frac{H_{4n-1}}{n^{2}} = \frac{H_{4n}}{n^{2}} - \frac{1}{4n^{3}} . $$
En el origen,
$$ f(z) = \frac{1}{z^{2}} \Big( \frac{1}{4z^{2}} - \frac{\zeta(2)}{2} - 4 \zeta(2) - 16 \zeta(3) z + \mathcal{O}(z^{2}) \Big),$$
lo que implica
$$\text{Res}[f(z),0] = -16 \zeta(3) .$$
Y en los cuartos enteros positivos,
$$ f(z) = \frac{\pi}{4z^{2}} \frac{(-1)^{n}}{z- \frac{2n+1}{4}} + \mathcal{O}(1),$$
lo que implica
$$ \begin{align} \text{Res} \Big[ f(z),\frac{2n+1}{4} \Big] &= \text{Res} \Big[\frac{\pi}{4z^{2}} \frac{(-1)^n}{z- \frac{2n+1}{4}}, \frac{2n+1}{4} \Big] \\ &= 4 \pi \ \frac{(-1)^{n}}{(2n+1)^{2}} . \end{align} $$
Juntando todo, tenemos
$$ - \frac{1}{2} \sum_{n=1}^{\infty} \frac{1}{n^{3}} + 2 \sum_{n=1}^{\infty} \frac{H_{4n}}{n^{2}} - \frac{1}{4} \sum_{n=1}^{\infty} \frac{1}{n^{3}} - 16 \zeta(3) + 4 \pi \sum_{n=0}^{\infty} \frac{(-1)^{n}}{(2n+1)^{2}} $$
$$ = - \frac{1}{2} \zeta(3) + 2 \sum_{n=1}^{\infty} \frac{H_{4n}}{n^{2}} - \frac{1}{4} \zeta(3) - 16 \zeta(3) + 4 \pi G = 0 .$$
Por lo tanto,
$$ \sum_{n=1}^{\infty} \frac{H_{4n}}{n^{2}} = \frac{67}{8} \zeta(3) - 2 \pi G .$$
EDITAR :
Encontré la expansión Laurent de $\psi(-4z)$ en los enteros positivos utilizando la ecuación funcional de la función digamma para expresar $\psi(4z)$ como
$$ \psi(4z) = \psi(4z+4n+1) - \frac{1}{4z+4n} - \frac{1}{4z+4n-1} - \ldots - \frac{1}{4z} .$$
Entonces evalué el límite $$\lim_{z \to -n} (z+n) \psi(4z) = - \frac{1}{4}$$ y el límite $$\lim_{z \to -n} \Big(\psi(4z) + \frac{1}{4} \frac{1}{z+n} \Big) = - \gamma +H_{4n} .$$
Esto conduce a la expansión $$\gamma + \psi (-4z) = \frac{1}{4} \frac{1}{z-n} + H_{4n} + \mathcal{O}(z-n) .$$
Hice algo similar para encontrar la expansión en los cuartos de enteros positivos.