In this investigation a simple method developed by introducing spin to Schrodinger equation to study the relativistic hydrogen atom. By separating Schrodinger equation to radial and angular parts, we modify these parts to the associated Laguerre and Jacobi differential equations, respectively. Bound state Energy levels and wave functions of relativistic Schrodinger equation for Hydrogen atom have been obtained. Calculated results well matched to the results of Dirac’s relativistic theory. Finally the factorization method and supersymmetry approaches in quantum mechanics, give us some first order raising and lowering operators, which help us to obtain all quantum states and energy levels for different values of the quantum numbers n and m.