We derive the spectral decomposition of the Lippmann-Schwinger equation for electrodynamics, obtaining the fields as a sum of eigenmodes. The method is applied to cylindrical geometries.