Code for exactly solving neuteron's binding states in different ways.
A fortran compiler (e.g., gfortran) is installed.
This code directly solves the equation.
gfortran direct_solution.f90
./a.out
This code solves the equation by expanding the wave function with eigen functions of harmonic oscillator.
gfortran ho_basis.f90
./a.out