  Using Mathematica 3.0, the Schroedinger equation for bound states is solved.
The method of solution is based on a numerical integration procedure together
with convexity arguments and the nodal theorem for wave functions. The
interaction potential has to be spherically symmetric. The solving procedure is
simply defined as some Mathematica function. The output is the energy
eigenvalue and the reduced wave function, which is provided as an interpolated
function (and can thus be used for the calculation of, e.g., moments by using
any Mathematica built-in function) as well as plotted automatically.
