While the second-order in time and in space Yee algorithm has been the primary FDTD algorithm used to date, higher-order FDTD methods have been proposed in order to reduce the numerical phase error. In 1989 Fang [10] proposed both a second-order accurate in time, fourth-order accurate in space FDTD algorithm and a fourth-order accurate in time, fourth-order accurate in space FDTD algorithm. Independently, Deveze et al. [163,164] published a similar higher-order FDTD method and also developed an absorbing boundary condition for the method.