Abstract-We describe the implementations of Drude-critical point model for describing dispersive media into finite difference time domain algorithm using piecewise-linear recursive-convolution and auxiliary differential equation methods. The advantages, accuracy and stability of both implementations are analyzed in detail. Both implementations were applied in studying the transmittance and reflectance of thin metal films, and excellent agreement is observed between analytical and numerical results.