I will discuss the linear stability of weakly charged and slowly rotating Kerr-Newman black holes under coupled gravitational and electromagnetic perturbations. We show that the solutions to the linearized Einstein-Maxwell equations decay at an inverse polynomial rate to a linearized Kerr-Newman solution plus a pure gauge term. The proof uses tools from microlocal analysis and a detailed description of the resolvent of the Fourier transformed linearized Einstein-Maxwell operator at low frequencies.