In this paper, we transform quasi-state frequency domain diffusion equations to wave equations in the fictitious wave domain. A complex frequency shifted perfectly matched layer (CFS-PML) boundary condition is adopted to the high-order finite difference time domain (FDTD) algorithm for reducing storage requirements and improving computational efficiency. We divide the dipole source by the pseudo-delta function for simulating electromagnetic response of arbitrary orientation electric dipole. By comparing the numerical solution with 1D analytical solution and 3D frequency controlled source electromagnetic (CSEM) results, we verify the validity, the accuracy and the efficiency of the algorithm, and then investigate and analyze the influence of different grid parameters and boundary conditions on numerical result at different frequency.