where the complex amplitudes \(E_n\)  are to be determined. This is feasible by employing the Jacobi-Anger expansion [1]}:
\(\mathrm {e}^{i x \sin {y}} = \sum _{n=-\infty }^{+\infty } J_n(x) ~\mathrm {e}^{iny},\) 
