Numerical solver for the 3 D stationary nuclear Schrödinger equation