An Efficient Numerical Method for Time Domain Computational Electromagnetic Simulation