FDTD is playing a key role in electromagnetic field computation. Moreover, it has been demonstrated that FDTD is special case of MRTD (multiresolution in time domain) if the base function and weight function is Haar scaling function in the course of solving Maxell equations by Wavelet-Galerkin scheme. Higher-order FDTD algorithm corresponds to Haar wavelet packet function.