ABSTRACT:Introducing the associated Bessel polynomials in terms of two nonnegative integers, we factorize their corresponding differential equation into a product of first-order differential operators by four different ways as shape invariance equations. Then, the radial part of the bound states of the Schrö dinger equation of a hydrogen-like atom is derived using one of the factorization methods in the framework of supersymmetric quantum mechanics. In this approach, we regenerate the radial bound states and their corresponding spectrum, which are consistent with the well-known facts. Based on the generalization of the supersymmetry idea, we shall show that two hydrogen-like atoms with the same energy of the electron possess three extra supersymmetric structures in addition to an ordinary one.