The Finite-Diierence Time Domain (FDTD) method has received considerable attention recently. The popularity of the FDTD method stems from the fact that it is not limited to a speciic geometry and it does not restrict the constitutive parameters of a scatterer. Furthermore, it provides a direct solution to problems with transient illumination, but can also be used for harmonic analysis. However,...