Radial wavefunction at origin for excited states $|{R_{nl}}^2| (0)$ for quarkonium

Question

I am solving the radial Schrodinger wave equation for quarkonium containing a quark and antiquark. The system is non-relativistic as quarks are charm and bottom which are heavy. The wave equation is numerically solved using the R K method in mathematica. I have wavefunction in numerical form. I want to find the radial wavefunctions for excited orbital angular momentums $|{R_{nl}}^2| (0)$ at the origin. I know the form given in eq 4.22 by Quigg and Rosner, but I am not fully confident on how to use this equation. I have calculated the masses already.