When we want to solve electromagnetic problems, the Finite Difference Time Domain (FDTD) method is a very useful numerical simulation technique to solve these problems. However, the traditional FDTD method is an explicit finite-difference scheme, so the method is limited by the Courant-Friedrich-Levy (CFL) stability condition. In other words, the minimum cell size will limit the maximum time-step size in a computational domain. Therefore, while simulating structures of fine scale dimensions, it will relatively result in a prohibitively high computation time generated by the maximum time-step size. The WLP-FDTD is based on the Weighted Laguerre Polynomials technique and the traditional FDTD algorithm. It is an implicit finite-difference equations. Therefore, it can completely avoid the stability constraint, and then improve calculation time by choosing relatively large time-step. In this thesis, we incorporate non-uniform grid method into the WLP-FDTD. By using them to simulate the structures of fine scale dimensions can reduce the computation time and memory usage. Further, we extend this method from two-dimensional to three-dimensional and add loss media into original formulations that will make the application of this method more widely.