Write a main program which makes the call
with the potential Veff[k]=-1./r[k] for a total number of intervals N=200 and our standard change of variables. Verify that for n=0, 1 and 2 you get values in the ranges (-0.9,-0.5], (-0.5,-0.125] and (-0.125,-0.055556], respectively.E=rtbisp480(func_SchrodingerNodes,-0.9,0.,1e-12,n,Veff,r,dr,N);