Archimedean volume ratio

And although he made many excellent discoveries, he is said to have asked his kinsmen and friends to place over the grave where he should be buried a cylinder enclosing a sphere, with an inscription giving the proportion by which the containing solid exceeds the contained.

Sphere and cylinder

(Quote from Plutarch; see references on the Archimedes site by Chris Rorres. Image from Wikipedia)

The proportion is 3:2, as Archimedes discovered. More generally, let A_n be the ratio of the volumes of n-dimensional unit ball B^n and of the circumscribed cylinder B^{n-1}\times [-1,1]. This sequence begins with A_1=1, A_2=\pi/4, A_3=2/3, A_4=3\pi/16, A_5=8/15, A_6=5\pi/12\approx 0.49\dots. Note that beginning with 6 dimensions, the ball occupies less than half of the cylinder. It is not hard to see that the sequence monotonically decreases to zero. Indeed, the ratio of (n-1)-dimensional cross-sections by the hyperplane x_{n}=t, for -1\le t\le 1, is (1-t^2)^{(n-1)/2}. Hence A_n = \frac{1}{2} I_{(n-1)/2} where I_p:=\int_{-1}^1 (1-t^2)^p\,dt. The integrand monotonically converges to 0 a.e. as p\to\infty, and therefore I_p does the same.

The Archimedean ratio A_n is rational when n is odd, and is a rational multiple of \pi when n is even. Indeed, integrating I_p by parts, we get
\displaystyle I_p = \int_{-1}^1 (1-t^2)^p\,dt = 2p \int_{-1}^1 t^2(1-t^2)^{p-1}\,dt = 2p(I_{p-1}-I_{p})
hence I_p = \frac{2p}{2p+1}I_{p-1}. In terms of A_n = \frac{1}{2} I_{(n-1)/2}, we have A_{n}= \frac{n-1}{n}A_{n-2}. Since A_1=1 and A_2=\pi/4 are known (the former being trivial and the latter being one of definitions of \pi), the claim follows.

Combining the monotonicity with the recurrence relation, we find that \displaystyle \frac{n-1}{n} =\frac{A_{n}}{A_{n-2}} \le \frac{A_{n}}{A_{n-1}}\le 1. Hence \displaystyle \frac{A_{n}}{A_{n-1}} \to 1, which yields the Wallis product formula for \pi. Indeed, the recurrence shows that for n even, \displaystyle  \frac{A_{n}}{A_{n-1}} is a particular rational multiple of \displaystyle \frac{A_{2}}{A_{1}} = \frac{\pi}{4}.

The above is essentially the inductive proof of Wallis’s formula I remember from my calculus (/real analysis) class. Unfortunately, between integration by parts and the squeeze lemma we lose geometry. In geometric terms, Wallis’s formula says that the volumes of B^{n+1}\times B^{n-1} are B^n\times B^n are asymptotically equivalent: their ratio tends to 1. Is there a way to see that these solids have a large overlap? Or, to begin with, that B^{n+1}\times B^{n-1} has the smaller volume? I don’t see this even with B^3\times B^1 vs B^2\times B^2. By the way, the latter is the bidisc from complex analysis.

Aside: Wallis’s formula for statistics students by Byron Schmuland presents a nice 2-page proof (which uses the relation between \pi and the Gaussian).

Added: Another way to generalize the Archimedean volume ratio to higher dimensions is to divide the volume of B^{2n+1} by the volume of B^{n}\times B^{n+1}. The quotient is always rational, namely \displaystyle \frac{(n+1)!}{(2n+1)!!}. The sequence begins with \displaystyle 1, \frac23, \frac25, \frac{8}{35}, \frac{8}{63}, \frac{16}{231}, \frac{16}{429} and obviously tends to zero. It does not appear to be related to \pi.

One thought on “Archimedean volume ratio”

  1. Note sure if this is rigorous, but we may think B^3 as an imaginary number w_1 (its projection to the imaginary plane) with another real number (the real part of imaginary number w_2). Then
    B^3 \times B^1 = \{(w_1,w_2) | Re(w_1)^2+Im(w_1)^2+Re(w_2)^2<1, Im(w_2)^2<1\}
    bidisc = \{(w_1,w_2) | Re(w_1)^2+Im(w_1)^2<1, Re(w_2)^2+Im(w_2)^2<1\}
    That is, for a fixed w_1 on the unit disc, w_2 is confined to lie on a rectangle with width 2(\sqrt{1-|w_1|^2}) and height 2 for B^3 \times B^1, and on the unit disc for bidisc. Intuitively the rectangle shrinks as w_1 moves away from the origin, so B^3 \times B^1 should span a smaller volume than bidisc.
    V(B^3 \times B^1)=\int_0^{2 \pi} \int_0^1 4 \sqrt{1-r^2} r dr d\theta = 8/3 \pi
    V(B^2 \times B^2)=\pi^2
    And V(B^3 \times B^1)/V(B^2 \times B^2)=8/(3\pi)=A_3/A_2

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

This site uses Akismet to reduce spam. Learn how your comment data is processed.