We study the stability properties of, and the phase error present in, a finite element scheme for Maxwell’s
equations coupled with a Debye or Lorentz polarization model. In one dimension we consider a second order
formulation for the electric field with an ordinary differential equation for the electric polarization added...