Relação entre a função de onda e o raio de Bohr

d1= D[r^2 Exp[-2r/a],{r}]

Solve[d1==0,r]

d2=D[r^2 (1-r/(2a))^2Exp[-r/a],{r}]

%==0

Solve[d2==0,r]

Reduce[d2==0,r]

3-Sqrt[5]//N

3+Sqrt[5]//N

Plot[ r^2(1-r/(2a))^2Exp[-r/a]/.a->1,{r,0,10},PlotRange->All]

Plot[ (1-r/(2a))^2Exp[-r/a]/.a->1,{r,0,10},PlotRange->All]

Plot[ (1-r/(2a))^2Exp[-r/a]/.a->1,{r,0,10},PlotRange->{0,0.04}]