Figure for ps3, p2.

f[n_,r_] := (2/(n!)) (r^(2n+1)) E^(-r^2) ;

f[#, r] &/@ Range[2,6,2]
Integrate[ f[#, r], {r,0,Infinity} ] &/@ {1,3,5}

p = With[
{nn=Range[2,6,2]},
 Plot[  Evaluate[f[#, r] &/@ nn], {r, 0, 5}, PlotRange -> Full, PlotLegends->Placed[(Row[{"n = ", #}] &/@ nn),{Right, Top}], AxesLabel->{r, f[n,r]}]

]
peeters`exportForLatex["ps3Fig2",p]