The bound state solutions of the Dirac equation with equal scalar and vector Makarov potentials are obtained. It is shown that angular component of the Dirac equation can be solved with the gactorization method， which enables us to find immediately the eigenvalues and at the same time the manipulation process for the normalized eigenfunctions. The radial bound state solutions are expressed in terms of the confluent hypergeometric functions and the energy equation is derived from the boundary condition satisfied by the radial wavefunctions.