==> geometry/points.in.sphere.s <==
Use spherical polar coordinates, and w.l.o.g. choose the polar axis
through one of the points. Now the distance between the two points is
sqrt ( r1^2 + r2^2 - 2 r1 r2 cos(theta))
and cos(theta) is (conveniently) uniformly distributed between -1 and
+1, while r1 and r2 have densities 3 r1^2 d(r1) and 3 r2^2 d(r2). Split
the total integral into two (equal) parts with r1 < r2 and r1 > r2, and
it all comes down to integrating polynomials.
More generally, the expectation of the n'th power of the distance
between the two points is
2^n . 72 / ((n+3)(n+4)(n+6))
So the various means are:
the (arithmetic) mean distance is 36/35 = 1.028571...
the root mean square distance is sqrt(6/5) = 1.095445...
the geometric mean distance is 2exp(-3/4) = 0.944733...
the harmonic mean distance is 5/6 = 0.833333...
the inverse root mean inverse square distance is
2/3 = 0.666666...