Finite-difference time-domain (FDTD) schemes based on the bilinear z -transform are introduced for modelling time-domain wave propagation in dispersive Drude-critical points media. The accuracy and efficiency of the proposed technique are verified by comparisons with other FDTD algorithms and analytical solutions.