A novel finite-difference time-domain (FDTD) method, called shift operator FDTD (SO-FDTD) method is developed for anisotropic magnetized dispersive media. The recursive relation between operators is used. In this paper, some expressions containing the dielectric constants of magnetized dispersive media are written as rational polynomial function. The SO-FDTD formulation for anisotropic magnetized plasma is derived. The high efficiency and effectiveness of the method are confirmed by computing the reflection and transmission through a magnetized plasma layer, with the direction of the propagation parallel to the direction of the biasing field. A comparison with frequency domain analytic results is included. The CPU time was several times shorter than that of the JEC method.