Mi solución sólo vale para un número par de factores.
Se puede considerar un sistema de cuadrados anidados $\Gamma_n$ con centros en $0$ y vértices $z=\pm\pi(n+\frac{1}{2})\pm i\pi(n+\frac{1}{2})$ . Para ellos $d_n=\min\limits_{z\in\Gamma_n}|z|=\pi(n+\frac{1}{2})$ , longitud de $\Gamma_n$ es $S_n=8\pi(n+\frac{1}{2})$ y $\frac{S_n}{d_n}=8$ . Se puede demostrar que en este sistema de cuadrados el $\text{ctg} z$ está limitada por alguna constante $c$ independiente de $n$ .
Consideremos ahora la función $f(z)=\frac{\pi^2\text{ctg} \pi z}{(z+1)(z+2)(z+3)(z+4)}$ y por el teorema del residuo de Cauchy $\displaystyle\int\limits_{\Gamma_n^+}f(z)dz=2\pi i(\text{res}_0f(z)+\text{res}_{-1}f(z)+\text{res}_{-2}f(z)+\text{res}_{-3}f(z)+\text{res}_{-4}f(z)+\sum\limits_{k=-n}^n\text{res}_{ k}f(z))$ . En la última suma $k\ne0,-1,-2,-3,-4$ .
Calculamos los residuos $\text{res}_0f(z)=\frac{\pi}{24}$ , $\text{res}_{-1}f(z)=-\frac{11\pi}{36}$ , $\text{res}_{-2}f(z)=\frac{\pi}{4}$ , $\text{res}_{-3}f(z)=\frac{\pi}{4}$ , $\text{res}_{-4}f(z)=-\frac{11\pi}{36}$ , $\text{res}_{k}f(z)=\frac{\pi}{(k+1)(k+2)(k+3)(k+4)}$ .
Así que $\displaystyle\int\limits_{\Gamma_n^+}f(z)dz=2\pi i\left(-\frac{5\pi}{72}+\sum\limits_{k=-n}^{-5}\frac{\pi}{(k+1)(k+2)(k+3)(k+4)}+\sum\limits_{k=1}^{n}\frac{\pi}{(k+1)(k+2)(k+3)(k+4)}\right)$ .
Además, $\left|\displaystyle\int\limits_{\Gamma_n^+}f(z)dz\right|\leq\frac{c\pi^2S_n}{(d_n-1)(d_n-2)(d_n-3)}\to0$ Así que $\sum\limits_{k=1}^{\infty}\frac{1}{(k+1)(k+2)(k+3)(k+4)}+\sum\limits_{k=-\infty}^{-5}\frac{1}{(k+1)(k+2)(k+3)(k+4)}=\frac{5}{72}$ . Finalmente, la segunda suma se reduce a la primera.