The Schrodinger equation in 1-D is extended to 3-D in much the same way as the energy method in 1-D is extended to 2-D to deal with central-force motion. The kinetic-energy term associated with the angular variable(s) is written as a function of r and the angular momentum l, and is treated as a form of potential energy (the centrifugal potential energy). The resulting equation looks (somewhat in the case of the Schrodinger equation) like a 1-D problem with an "effective" potential energy. Although the program that finds energy levels and plots wavefunctions is itself very simple, its development involves a fair bit of detail (hint: skip to the program if you are not interested in detail).