Así, con el fin de obtener la respuesta, intenté aplicar algunas expansiones de Taylor, lo que me llevó a nada en realidad. Después de un tiempo traté de usar el teorema de la suma
$\sum_{n=-\infty}^{+\infty}{f\left(n\right)}=-\sum_{i=1}^{m}{Res_{z=z_i}{\pi\cot\left(\pi z\right)f\left(z\right)}}$ $f\left(z\right)$'s de los polos.
El residuo de a $z=-\frac{1}{2}$ es igual al residuo de a $z=\frac{1}{2}$,$-\frac{\pi ^2}{16}$. Lo que tengo parece ser no la respuesta correcta después de todo:
$-\sum_{i=1}^{m}{Res_{z=z_i}{\pi\cot\left(\pi z\right)f\left(z\right)}} = -\frac{\pi ^2}{8}$, por lo que el $\sum_{n=-\infty}^{+\infty}{f\left(n\right)}=\frac{\pi ^2}{8}$ Puesto que el $f\left(z\right)$ es incluso, llego $\sum_{n=0}^{+\infty}{f\left(n\right)}=\frac{\pi ^2}{16}$
Sin embargo, la tarea inicial fue encontrar la suma de $\sum_{n=1}^{+\infty}{f\left(n\right)}$. Supongo que $\sum_{n=1}^{+\infty}{f\left(n\right)}=\sum_{n=0}^{+\infty}{f\left(n\right)}-f\left(0\right)=\frac{\pi ^2}{16}-1$, lo cual es incorrecto, ya que de acuerdo a Mathematica, que debería obtener $\frac{\pi ^2}{16}-\frac{1}{2}$. No sé donde está mi error en este punto. Tal vez, todo el enfoque no está bien ejecutado.