VOLUME SPHERE/TORUS Maths21a, O. Knill 6) IN LAS VEGAS.

Maths21a, O. Knill 6) IN LAS VEGAS. The Monte Carlo Method is to shoot randomly onto the cube [L, L] [L, L] [L, L]
and see how many times we hit the sphere. Here an experiment with Mathematica:
We first calculate the volume of a sphere of radius R in different ways. Then we show how to calculate the R := (2Random[] 1); k = 0; Do[x = R; y = R; z = R; If[x2 + y2 + z2 < 1, k + +], {10000}]; k/10000
volume of the torus in three different ways. The page serves more as an illustration for the variety of tools
which are available. We did not cover all the material (yet) to understand all of the methods. Assume, we hit 5277 of 10000 the measured fraction of the volume of the sphere with the volume of the cube 8
is 0.5277. The volume of 1/8th of the sphere is /6 = 0.524

7) USING GAUSS THEOREM (see later) The vector field F (x, y, z) = (x, y, z) has divergence 3 Gauss theorem
tells that 3V is the flux of the vector field through the surface which is L times the surface area 4L2 . Therefore,
V = 4L3 /3

SEVEN WAYS TO COM-

PUTE THE VOLUME OF
THE SPHERE
FIVE WAYS TO COM-
PUTE THE VOLUME OF
THE TORUS

1) IN RECTANGULAR COORDINATES. The volume is

Z Z Z Z L Z
L2 x2 Z L2 x2 y2 1) WITH TORAL COORDINATES.
V = dxdydz = dzdydx
R L L2 x2 L2 x2 y 2
T (r, , ) = (x, y, z) = ((b + r cos()) cos(), (b + r cos()) sin(), r sin()) parameterizes the torus.
R L R L2 x2
After computing the most inner integral, we have V = 2(L2 x2 y 2 )1/2 dydx. Define a2 = L2 x2 The Jacobean is det(T ) = (x,y,z)
(r,,) = r(b + r cos(). The torus is the image of the cube [0, a] [0, 2] [0, 2]
L L2 x2
and use the integral in the box below:
under the map T . The change of variables formula gives
Z L
V = (L2 x2 ) dx = 2L3 2L3 /3 = 4L3 /3 .
L
Z a Z 2 Z 2 Z a
r(b + r cos()) dddr = (2)(2) br dr = 2 2 a2 b
x 0 0 0 0
Substitution = sin(u), dx = a cos(u)du gives:
a
Ra R /2 q R /2
a2 x2 dx = /2 a 1 sin2 (u)a cos(u) du = a2 /2 cos2 (u)du = a2 /2 .
a 2) USING CYLINDRICAL COORDINATES.
If we fix the z coordinate, we obtain
an annulus with
inner radius b a2 z 2 and outer radius 2 2
R a b + a z .
2 2
2 2 2 2
This annulus has the area (b + a z ) (b a z ) . Therefore, the volume is 4b a a z 2 dz =
2
2) IN CYLINDRICAL COORDINATES. At height z, we parameterize a disc of radius L2 z 2 , so that the
integral is 4b(a/2) = 2 2 a2 b.
Z L Z L2 z2 Z 2 Z L
rddrdz = (L2 z 2 )2dz = 4L3 2L3 /3 = 4L3 /3 . 3) USING PAPPUS CENTROID THEOREM. The volume of a solid of revolution generated by the revolution
L 0 0 L
of a region S in the x z plane around the z axis is equal to the product of the area of S and the arc length
2b of the circle on which the center of S moves.
3) IN SPHERICAL COORDINATES. In the case of the torus, the length of the curve is 2b. The area of the lamina is A = a2 . Therefore, the
volume is 2 2 a2 b.
Z R Z 2 Z Z R Z
2 sin() ddd = 2 2 d sin() d = 4L3 /3 . PROOF OF THE CENTROID THEOREM. We use a coordinate change transformation. In Polar coordinates,
0 0 0 0 0
the lamina S with center of mass (b, c) is parametrized by r and z. Introduce new coordinates T (u, v) = (u+b, v+
c) = (r, z) so that (0, 0) is the center of mass in theRnew
R coordinates. The
R R Jacobean of this coordinate
RR change is 1.
The volume of the solidR R of revolution is V = (2) S r drdz = (2) R
(u + b) Rdudv
R = 2b b dudv = 2bA,
4) WITH CAVALIERI. Cavalieri cuts the hemisphere at height z to obtain a disc of radius L2 z 2 with area where we used that
RR
u dudv = 0 because (u, v) = (0, 0) = ( R u dudv/A, R v dudv/A) is the center of
2 2 R
(L z ). He looked now at the complement of a cone of height L and radius L which when cut at height z mass of R.
gives a ring of outer radius L and inner radius z. The ring has area (L2 z 2 ). Cavalieri concludes (Cavalieri
principle) that the volume of that body is the same as the volume of the hemisphere. Since the difference of 4) MONTE CARLO AGAIN. Lets assume b = 2 and a = 1. If (x, y, z) is a random point in [3, 3] [3, 3]
the volume of the cylinder and the cone which is L3 L3 /3 the hemisphere has the volume 2L3 /3 and the [1, 1] then (r 2)2 + z 2 1 is the condition to be in the torus, where r2 = x2 + y 2 .
sphere has volume 4L3 /3. Assume, we hit 5484 of 100000 the measured fraction of the volume 72 of the box we estimate 725484/100000 =
3.94848 for the actual volume 4 2 = 39.4784.
5) CAS. Integrate[1, {x, L, L}, {y, Sqrt[L2 x2 ], Sqrt[L2 x2 ]}, {z, Sqrt[L2 x2 y 2 ], Sqrt[L2 x2 y 2 ]}]
5) CAS. Integrate[r, {r, 1, 3}, {theta, 0, 2Pi}, {z, Sqrt[1 (r 2)2 ], Sqrt[1 (r 2)2 ]}]