rw[n_,l_,r_]=(n r)^l LaguerreL[n-l-1,2 l+1,n r] Exp[-n r/2]

g1=Plot[{rw[1,0,r],rw[1,0,r]^2}, {r,0,5},PlotStyle->{Thickness[0.0075]},PlotRange->Full];
g2=PolarPlot[Abs[SphericalHarmonicY[0,0,θ,0]]^2,{θ,0,2π},PlotStyle->{Thickness[0.0075]}];
g3=DensityPlot[0.001rw[1,0,(Sqrt[x^2+y^2])]^2Abs[SphericalHarmonicY[0,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g4=DensityPlot[0.001rw[1,0,(Sqrt[x^2+y^2])]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g5=DensityPlot[Abs[SphericalHarmonicY[0,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
GraphicsArray[{{g1,g4},{g2,g5},{g3}}, ImageSize->Large]

g1=Plot[{rw[2,0,r],rw[2,0,r]^2}, {r,0,5},PlotStyle->{Thickness[0.0075]},PlotRange->Full];
g2=PolarPlot[Abs[SphericalHarmonicY[0,0,θ,0]]^2,{θ,0,2π},PlotStyle->{Thickness[0.0075]}];
g3=DensityPlot[0.001rw[2,0,(Sqrt[x^2+y^2])]^2Abs[SphericalHarmonicY[0,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g4=DensityPlot[0.001rw[2,0,(Sqrt[x^2+y^2])]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g5=DensityPlot[Abs[SphericalHarmonicY[0,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
GraphicsArray[{{g1,g4},{g2,g5},{g3}}, ImageSize->Large]

g1=Plot[{rw[2,1,r],rw[2,1,r]^2}, {r,0,5},PlotStyle->{Thickness[0.0075]},PlotRange->Full];
g2=PolarPlot[Abs[SphericalHarmonicY[1,0,θ,0]]^2,{θ,0,2π},PlotStyle->{Thickness[0.0075]}];
g3=DensityPlot[0.001rw[2,1,(Sqrt[x^2+y^2])]^2Abs[SphericalHarmonicY[1,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g4=DensityPlot[0.001rw[2,1,(Sqrt[x^2+y^2])]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g5=DensityPlot[Abs[SphericalHarmonicY[1,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
GraphicsArray[{{g1,g4},{g2,g5},{g3}}, ImageSize->Large]

g1=Plot[{rw[3,0,r],rw[3,0,r]^2}, {r,0,5},PlotStyle->{Thickness[0.0075]},PlotRange->Full];
g2=PolarPlot[Abs[SphericalHarmonicY[0,0,θ,0]]^2,{θ,0,2π},PlotStyle->{Thickness[0.0075]}];
g3=DensityPlot[0.001rw[3,0,(Sqrt[x^2+y^2])]^2Abs[SphericalHarmonicY[0,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g4=DensityPlot[0.001rw[3,0,(Sqrt[x^2+y^2])]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g5=DensityPlot[Abs[SphericalHarmonicY[0,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
GraphicsArray[{{g1,g4},{g2,g5},{g3}}, ImageSize->Large]

g1=Plot[{rw[3,1,r],rw[3,1,r]^2}, {r,0,5},PlotStyle->{Thickness[0.0075]},PlotRange->Full];
g2=PolarPlot[Abs[SphericalHarmonicY[1,0,θ,0]]^2,{θ,0,2π},PlotStyle->{Thickness[0.0075]}];
g3=DensityPlot[0.001rw[3,1,(Sqrt[x^2+y^2])]^2Abs[SphericalHarmonicY[1,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g4=DensityPlot[0.001rw[3,1,(Sqrt[x^2+y^2])]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g5=DensityPlot[Abs[SphericalHarmonicY[1,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
GraphicsArray[{{g1,g4},{g2,g5},{g3}}, ImageSize->Large]

g1=Plot[{rw[3,2,r],rw[3,2,r]^2}, {r,0,5},PlotStyle->{Thickness[0.0075]},PlotRange->Full];
g2=PolarPlot[Abs[SphericalHarmonicY[2,0,θ,0]]^2,{θ,0,2π},PlotStyle->{Thickness[0.0075]}];
g3=DensityPlot[0.001rw[3,2,(Sqrt[x^2+y^2])]^2Abs[SphericalHarmonicY[2,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g4=DensityPlot[0.001rw[3,2,(Sqrt[x^2+y^2])]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
g5=DensityPlot[Abs[SphericalHarmonicY[2,0,ArcTan[y/x],0]]^2,{x,-5,5},{y,-5,5},PlotPoints->100];
GraphicsArray[{{g1,g4},{g2,g5},{g3}}, ImageSize->Large]