We introduce the Dirac equation in four-dimensional gravity which is a generally covariant form. We choose the suitable variable and solve the corresponding equation. To solve such equation and to obtain the corresponding bispinor, we employ the factorization method which introduces the associated Laguerre polynomial. The associated Laguerre polynomials help us to write the Dirac equation of four-dimensional gravity in the form of the shape invariance equation. Thus we write the shape invariance condition with respect to the secondary quantum number. Finally, we obtain the spinor wave function and achieve the corresponding stability of condition for the four-dimensional gravity system.
