The formula you give is not specifically the distribution of one kind of orbital.
Rather, this formula transforms Ψ, whose square is a probability
per volume unit, into a probability
per radius unit. That is, find the electron in this particular cubic picometre for instance, versus find the electron in this chosen picometre interval of distance to the nucleus, in any direction. The measurement volume is then a picometre thickness at a sphere around the nucleus.
The 4πr
2 just tells that the volume of the "picometre times sphere area" increases as the sphere's area for a given thickness.
Worth investing some time to grasp, because it tells you why the S orbitals have the maximum probability density right at the nucleus, but when drawn versus the distance, probability is zero there. This is because of the r
2 only.
If you've read about density of states for phonons or electrons in 3D crystals (semiconductors can also be produced 2D or 1D), it's the same r
2 story here.
It's also the same r
2 story for the kinetic energy distribution (hi Boltzman) in gas. In each direction it decreases exponentially, but the "non-oriented" energy distribution is zero around zero energy, and goes first to a maximum before decreasing.
Only s orbitals have nonzero density at the nucleus and contribute to radioactive decay by electron capture, as an illustration.
As for p, d, f... orbitals, ψ depends also on the angle, not just the distance as for S orbitals. Fortunately, ψ is just a product of the angle's function and the distance's function, but the distance's function
depends on the shell: 2p, 3p... - just as 1s, 2s... are not the same. Very nice pictures there:
http://winter.group.shef.ac.uk/orbitron/click a shell in the list at left.