Compile and run your routines with the main program ``getEdebug.c'' in
Appendix A. This program accepts as input the
number of integration intervals N and the number of states you want,
and then returns all of the energies. At N=1000 intervals, you should
be able to recover the first five states to within 5 decimal places.
You should also be able to verify the scaling of your error in going
from N=1000 to N=10000.