The following python program is an example of how one can obtain the eigenvalues for a single-nucleon moving in a harmonic oscillator potential. It is rather easy to change the onebody-potential with ones like a Woods-Saxon potential.

• The c++ and Fortran versions of this program can be found at https://github.com/NuclearStructure/PHY981/tree/master/doc/pub/spdata/programs.
• The c++ program uses the c++ library armadillo, see http://arma.sourceforge.net/.