N[{BesselJZero[1/2,1],BesselJZero[3/2,1],BesselJZero[5/2,1],BesselJZero[1/2,2],BesselJZero[7/2,1],BesselJZero[3/2,2]}]
Out[1]= {3.14159, 4.49341, 5.76346, 6.28319, 6.98793, 7.72525}

sBesselj[l_,x_]=Sqrt[Pi/(2 x)] BesselJ[l+(1/2),x]
Plot[sBesselj[0,x],{x,0,10}]
Plot[sBesselj[1,x],{x,0,10}]

now easier:
Plot[SphericalBesselJ[0,x],{x,0,10}]
Plot[SphericalBesselJ[1,x],{x,0,10}]

The Figure 2.11 energies need to be multipled by Pi^2 to compare to our results
or our results need to be divided by Pi^2 to compare to Figure 2.11

%1^2/Pi^2
Out[5]= {1., 2.04575, 3.36563, 4., 4.94763, 6.0468}

R[n_,l_,r_]=Sqrt[2 Factorial[n]/Factorial[l+1/2+n]] r^l LaguerreL[n, l +1/2, r^2] Exp[-r^2/2]
Table[Integrate[R[n,0,r] r^2 R[n,0,r] r^2, {r,0,Infinity}],{n,0,10}]
Table[Integrate[R[n,1,r] r^2 R[n,1,r] r^2, {r,0,Infinity}],{n,0,10}]
Table[Integrate[R[n,2,r] r^2 R[n,2,r] r^2, {r,0,Infinity}],{n,0,10}]
Table[Integrate[R[n,3,r] r^2 R[n,3,r] r^2, {r,0,Infinity}],{n,0,10}]

         3  7  11  15  19  23  27  31  35  39  43
Out[2]= {-, -, --, --, --, --, --, --, --, --, --}
         2  2  2   2   2   2   2   2   2   2   2

In[3]:= 
         5  9  13  17  21  25  29  33  37  41  45
Out[3]= {-, -, --, --, --, --, --, --, --, --, --}
         2  2  2   2   2   2   2   2   2   2   2

In[4]:= 
         7  11  15  19  23  27  31  35  39  43  47
Out[4]= {-, --, --, --, --, --, --, --, --, --, --}
         2  2   2   2   2   2   2   2   2   2   2

In[5]:= 
         9  13  17  21  25  29  33  37  41  45  49
Out[5]= {-, --, --, --, --, --, --, --, --, --, --}
         2  2   2   2   2   2   2   2   2   2   2
	 
(2 n + l + 3/2) Note this is the same as E...see that <r^2>=2<V>=<E> by Virial Thm

Table[Sqrt[Pi] 4^n Integrate[R[n,0,r] r R[n,0,r] r^2, {r,0,Infinity}],{n,0,10}]
http://www.research.att.com/~njas/sequences
reports this is
In[11]:=  Table[n Binomial[2 n,n],{n,1,11}]

n=0 is  Gamma[2 + l]/Gamma[l+3/2] 
n=1 is (l+9/4) Gamma[2 + l]/Gamma[l+5/2]
n=2 is (32 l^2+176 l +225)/32 Gamma[2 + l]/Gamma[l+7/2]
n=3 is (128 l^3+1248 l^2+3820 l +3675)/128 Gamma[2 + l]/Gamma[l+9/2]
n=4 is (2048 l^4+30720 l^3+164224 l^2+370752 l+297675)/2048 Gamma[2 + l]/Gamma[l+11/2]
n=5 is (8192 l^5+174080 l^4+1415680 l^3+5503360 l^2+10207068 l+7203735)/8192 Gamma[2 + l]/Gamma[l+13/2]

clearly (l^n+lower terms) Gamma[2 + l]/Gamma[l+2n+3/2]
                      ^^we could produce a formula for the constant term using the l=0 result above