Here $\lambda \in \mathbb{C}$ and $$a_n = \frac{1}{\pi}\int_0^\pi e^{\lambda \cos(t)} \cos(nt)\,dt$$
I'm clueless about how to reach that expression for $a_n$. The things that came to my mind were to develop the Laurent series for some of the expressions, but didn't find myself nearer. I know the following:
$$e^{\frac{1}{2}\lambda(z+\frac{1}{z})}=\sum_{n=0}^\infty \frac{(\frac{1}{2} \lambda(z+\frac{1}{z}))^n}{n!}$$
$$e^{\frac{1}{2}\lambda(z+\frac{1}{z})}=e^{\frac{1}{2}\lambda z}e^{\frac{1}{2}\lambda\frac{1}{z}}=\sum_{n=0}^\infty\frac{(\frac{1}{2}\lambda z)^n}{n!}\sum_{n=0}^\infty\frac{(\frac{1}{2}\lambda \frac{1}{z})^n}{n!}=\sum_{n=0}^\infty\sum_{k=0}^n\frac{\lambda^nz^{2k-n}}{2^nk!(n-k)!}$$
where the last equality is because of the Cauchy Product.
I think it's easier to use the fact that if $f$ is holomorphic in a punctured neighbourhood $U$ around $z = 0$ and $f(z) = \sum_{n = -\infty}^\infty a_n z^n$ is the Laurent series at $z = 0$, then $$ a_n = \frac{1}{2\pi i}\oint_{\gamma} \frac{f(z)}{z^{n+1}} dz$$ for any closed contour $\gamma$ in $U$ that winds around $z = 0$ once anticlockwise.
For $f(z) = e^{\tfrac 1 2 \lambda (z + 1/z)}$, you probably want to take $\gamma$ to the unit circle, parameterised as $\gamma(t) = e^{it}$ for $ t \in [0, 2\pi].$ I hope this works out!