Solving the Schroedinger equation for bound states with Mathematica 3.0
  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. 
Wolfgang Lucha, F. F. Schoberl