Abstract
An alternative formulation for a frequency dependent (FD) ADI-FDTD method is proposed, based on an electric flux field D tridiagonal matrix rather than an electric field matrix. The procedure for D source excitation is explained. Results show that the accuracy of the method is unaffected by the inclusion of frequency dependency in the model, compared with FD-FDTD.