Thank you again for your help. Another question, though. In the link you gave me, both functions rR(r) and (rR(r))^2 are plotted, but with no explanation or definition.
So, I assumed that rR(r) was the radial probability, and I took the formula of rR(r) and differentiated it with respect to r. Then I set the resulting derivative equal to zero, and solved for r, which would be the r for the maximum of the function rR(r). This should give me what I want--the maximum radial probability as a function of Z. What I got was rmax = a0/Z. (I'm assuming a0 is the Bohr radius of H1 atom, 5.29E-11 meters.) If all this is right, then for U238, where Z=92, then this maximum radial probability of the s1 inner radius is:
(5.29E-11)/(92)=(2.62E-13) meters.
Which, of course, is a quite a bit smaller than the Bohr radius for the H1 atom.
Does this make sense? Am I doing this right, or do I need to work with the function (rR(r))^2 instead?
Thanks again for your help.