  A general form for ladder operators is used to construct a method to solve
bound-state Schr\"odinger equations. The characteristics of supersymmetry and
shape invariance of the system are the start point of the approach. To show the
elegance and the utility of the method we use it to obtain energy spectra and
eigenfunctions for the one-dimensional harmonic oscillator and Morse potentials
and for the radial harmonic oscillator and Coulomb potentials.
