I need to take the residue of the following function at infinity: $$ f(z) = \oint_\infty \left(\frac{e^{-\alpha/z}e^{-\alpha z}}{z}\right)dz $$ Which, up to a sign, is invariant under inversions $z\rightarrow 1/w$, since $dz \rightarrow -\frac{1}{w^2} dw$.
How can I proceed to take a Laurent expansion of this function at $z = \infty$ (or $z = 0$)?
Around $z = 0$, the exponential is the generating function for the Modified Bessel function, e.g. $$ e^{\frac{t}{2}(z + 1/z)} = \sum_{n=-\infty}^\infty I_n(t)z^n. $$
Thus, we find
$$ f(z) = \oint_\infty \left(\frac{e^{-\alpha/z}e^{-\alpha z}}{z}\right)dz = \oint_\infty \frac{dz}{z}\sum_{n=-\infty}^\infty I_n(2\alpha)z^n = I_0(2\alpha). $$ Which is the correct answer.