plot (exp(-r/(2*0.528E-10))/(r/0.528E-10))*((1-r/0.528E-10)) for r>0