This routine will ultimately have two uses: (1) to compute the
solution 
to Schrödinger's equation for a given energy E,
(2) to return information to allow us to determine the energies E. In
this problem set, we focus on (2), determining E, and thus one of the
arguments of the function (Psiout[]) will not be used until later.