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 form ...