  The supersymmetry in quantum mechanics and shape invariance condition are
applied as an algebraic method to solve the Dirac-Coulomb problem. The ground
state and the excited states are investigated using new generalized ladder
operators.
