Plot[(4*r^2)/(0.53*10^(-11))^3*Exp[(-2*r)/(0.53*10^(-11))], {r, 0, 5*10^(-11)}]