$TextStyle = {FontFamily -> Times, FontSize -> 14}

j[n_, x_] := 
      BesselJ[n + 1/2, x]
Plot[{j[0, x], j[1, x], j[2, x]}, {x, 0, 10}]

<< NumericalMath`BesselZeros`
Clear[t]
z[l_] := BesselJZeros[l + 1/2, 5]
En[l_, n_] := z[l][[n]]^2
ztable = 
      Table[
Table[{l, z[l][[n]]}, {n, 1, 
            5}], 
{l, 0, 4}]
Etable = {}
Do[
        Etable = 
          Join[Etable, 
Table[{l, En[l, n]}, {n, 1, 
                5}]], 
{l, 0, 4}];

TableForm[ztable]
\(Print[
      StyleForm["Compare table on page 180.", 
FontSize -> 18, FontColor -> RGBColor[1, 0, 0]]]\)

(* Energy level diagram *)
Clear[
        spd]
\(spd = {"s", "p", "d", "f", "g"}\)
spdl = {RGBColor[0, 0, 
          1], 
Table[
          Text[spd[[i + 1]], {0.1 + i, 450}], {i, 0, 
            4}]}
\(ListPlot[
      Etable, 
PlotStyle -> {PointSize[
            0.03]}, 
AxesLabel -> {"l", "E"}, 
Epilog -> {spdl}]\)

3.77*69

(* Radial wave functions *)
l = 
        0
Do[
zero = z[
            l][[n]]; 
u[n_, l_, x_] := 
        j[l, x*zero]; 
plotrf[n, l] = 
        Plot[u[n, l, x], {x, 0, 1}], 
{n, 1, 
        5}]
Show[
      Table[plotrf[n, l], {n, 1, 5}]]