Aquí hay otra solución usando el teorema de los residuos. Si tratamos de integrar la función $\frac{\sinh(az)}{\sinh(\pi z)}\frac{1}{1+n^2z^2}$ directamente, entonces obtenemos el resultado
$$\int_0^\infty \frac{\sinh(ax)}{\sinh(\pi x)}\frac{{\rm d}x}{1+n^2x^2} = \frac{\pi}{2n}\frac{\sin(a/n)}{\sin(\pi/n)} -\sum_{k=1}^\infty \frac{(-1)^{k+1}}{k^2n^2-1}\sin(ka)$$
que requiere el conocimiento de la identidad $\frac{\pi}{2n}\frac{\sin(a/n)}{\sin(\pi/n)} = \sum_{k=1}^\infty \frac{kn(-1)^{k+1}}{k^2n^2-1}\sin(ka)$ para obtener el deseado de identidad.
Podemos solucionar esto, en lugar de considerar el ligeramente modificada de la función
$$f(z) = \frac{\sinh(z)}{\sinh(\pi z)}\frac{1+untes}{1 + n^2z^2}$$ Bajo $z\a-z$ el plazo de dos en $f$ son pares e impares, respectivamente, que nos da
$$\int_{-\infty}^\infty f(x)\,{\rm d}x = 2\int_{0}^\infty \frac{\sinh(x)}{\sinh(\pi x)}\frac{{\rm d}x}{1 + n^2x^2}$$
y tanto de las integrales anteriores convergen mientras $|a|<\pi$.
Ahora podemos integrar $f$ a $z=-I$ a $z=R$ a lo largo del eje real y, a continuación, a lo largo de un semi-círculo de $C_R$ en la mitad superior del plano -. Tomamos $R=N+\frac{1}{2}$ $N\in\mathbb{N}$ para evitar la integración de más de uno de los polos de $f$. Para esta elección de $R$ tenemos $|f(z)| \leq \frac{1}{nR-1}$ junto $C_R$ entonces $\lim\limits_{R\to\infty}\int_{C_R}f(z)\,{\rm d}z = 0$.
Los polos de $f$ en la mitad superior del plano se encuentra en $z=ik$ para $k\in\mathbb{N}$ con
$$\text{Res}\left[\frac{\sinh(z)}{\sinh(\pi z)}\frac{1+untes}{1 + n^2z^2}; ik\right] = \frac{1}{\pi i}\frac{(-1)^{k+1}}{1+kn}\sin(ka)$$
y por el residuo teorema se sigue que
$$\int_{0}^\infty \frac{\sinh(x)}{\sinh(\pi x)}\frac{{\rm d}x}{1 + n^2x^2} = \sum_{k=1}^\infty \frac{(-1)^{k+1}}{1+kn}\sin(ka)$$