plot exp(-r/(0.528E-10)), r=0..infinity