Abstract
Efforts have been made to solve the Dirac equation with axially deformed scalar and vector Woods–Saxon potentials in the coordinate space with the imaginary time step method. The results of the single-particle energies thus obtained are consistent with those calculated with the basis expansion method, which demonstrates the feasibility of the imaginary time step method for the relativistic static problems.