Simple treatment of multi-term dispersion in FDTD

Three new simple and efficient algorithms are proposed for the numerical treatment of the multi-term Debye or Lorentz dispersion in the FDTD method. The formulation is based on the auxiliary differential equation, but requires much fewer operations than the published schemes. The approach is equivalent to the best higher order recursive schemes in terms of memory and computational expense, but without the linearity assumptions.