Como se indicó anteriormente, esta suma puede ser evaluado utilizando los residuos de la teoría. Voy a decir el resultado: $f$ suficientemente "bien comportado" (lo que significa que se desvanece lo suficientemente rápido a lo largo de las secciones verticales de la típica rectangular de contorno se utiliza para derivar la siguiente relación):
$$\sum_{n=-\infty}^{\infty} f(n) = -\sum_k \operatorname*{Res}_{z=z_k} [\pi \, \cot{\pi z} \, f(z)]$$
donde $z_k$ es un no-reales(integral) de un polo de $f$ en el plano complejo.
En este caso, $f(z)=1/(z^4+a^4)$ y los polos están en $z = a\,e^{i (2 k-1) \pi/4}$, $k \in \{1,2,3,4\}$. Así, la evaluación de la suma se reduce a la generalización de los residuos en estas polos:
$$-\sum_k \text{Res}_{z=z_k} \pi \, \cot{\pi z} \, f(z) = -\sum_{k=1}^4 \frac{\pi \cot{(\pi e^{i (2 k-1) \pi/4})}}{4 a^3 e^{i 3 (2 k-1) \pi/4}} $$
Ahora, se puede deducir que
$$\cot{(b\, e^{i t})} = \frac{\sin(2 b \cos{t})}{\cosh(2 b \cos{t})-\cos(2 b \sin{t})}-i \frac{ \sinh(2 b \sin{t})}{\cosh(2 b \cos{t})-\cos(2 b \sin{t})}$$
Usted puede verificar esta fórmula en un programa como Wolfram Alpha o Mathematica, o puede derivar con el seno y coseno, además de teoremas.
El álgebra de los involucrados puede resultar tedioso y propenso a errores. Voy a esbozar aquí un par de consejos para obtener el resultado correcto. La reescritura de la suma de los residuos como (signo negativo):
$$\frac{\pi}{4 a^3} \left [ e^{-i 3 \pi/4} \frac{-\sin{(\sqrt{2} \pi a)} + i\sinh{(\sqrt{2} \pi a)} }{\cosh{(\sqrt{2} \pi a)} -\cos{(\sqrt{2} \pi a)}} + \\ e^{-i 9 \pi/4} \frac{\sin{(\sqrt{2} \pi a)} + i\sinh{(\sqrt{2} \pi a)} }{\cosh{(\sqrt{2} \pi a)} -\cos{(\sqrt{2} \pi a)}} + \\ e^{-i 15 \pi/4} \frac{\sin{(\sqrt{2} \pi a)} - i\sinh{(\sqrt{2} \pi a)} }{\cosh{(\sqrt{2} \pi a)} -\cos{(\sqrt{2} \pi a)}} + \\e^{-i 21 \pi/4} \frac{-\sin{(\sqrt{2} \pi a)} - i\sinh{(\sqrt{2} \pi a)} }{\cosh{(\sqrt{2} \pi a)} -\cos{(\sqrt{2} \pi a)}} \right ]$$
Debe ser claro que el exponenciales puede ser reducido, y se termina con dos pares de complejos conjugados si todo se hace correctamente. En este punto voy a dejar los detalles para que el lector y el estado del resultado final:
$$\sum_{n=-\infty}^{\infty} \frac{1}{n^4+a^4} = \frac{\pi}{\sqrt{2} \, a^3} \frac{\sinh{(\sqrt{2} \pi a)}+\sin{(\sqrt{2} \pi a)}}{\cosh{(\sqrt{2} \pi a)} -\cos{(\sqrt{2} \pi a)}}$$
BONO
Como comprobación, se puede derivar la conocida fórmula
$$\sum_{n=1}^{\infty} \frac{1}{n^4} = \frac{\pi^4}{90}$$
considerando el comportamiento del resultado anterior en el límite de $a \to 0$. Tenga en cuenta que
$$\sum_{n=1}^{\infty} \frac{1}{n^4+a^4} = \frac12 \left [ \sum_{n=-\infty}^{\infty} \frac{1}{n^4+a^4} - \frac{1}{a^4}\right]$$
Ahora Taylor ampliar el resultado lo suficiente como para conseguir un nonvanishing resultado:
$$ \frac{\sinh{(\sqrt{2} \pi a)}+\sin{(\sqrt{2} \pi a)}}{\cosh{(\sqrt{2} \pi a)} -\cos{(\sqrt{2} \pi a)}} \sim \frac{2 \sqrt{2} \pi a \left ( 1+ \frac{\pi^4 a^4}{30} \right )}{2 \pi^2 a^2 \left ( 1+ \frac{\pi^4 a^4}{90} \right )} \sim \frac{\sqrt{2}}{\pi a} \left ( 1+ \frac{\pi^4 a^4}{45} \right) $$
Poniendo todo esto, se obtiene la suma deseada.