PH 314 Fall 2002 MJM      Ellipse and ellipsoid integrals.                               10/13/02

The equation of ellipse is (x/a)2 + (y/b)2 =1. Now let x' = x/a, and y'=y/b, so that the ellipse equation becomes
    x'2 + y'2 =1.
This is a circle of radius 1.

When we calculate the area of an ellipse we write A = ò dx dy over the bounding curve. The bounding curve starts out as (x/a)2 + (y/b)2 =1 , which makes for messy substitutions and a somewhat complicated calculation. But! When we go to the new coordinates we have A = ab ò dx' dy' over the new curve.  Since the bounding curve in the new coordinates is a unit circle, the area is obviously p which makes the area of the ellipse A = p ab .

Ellipsoid calculations go the same way, starting with (x/a)2 + (y/b)2 + (z/c)2 = 1 .
When we change to x' = x/a, y' = y/b, and z' = z/c, the bounding surface is a unit sphere  Then the volume of an ellipsoid is easily worked out to be 4/3 p abc

Calculating integrals like I1 = ellipsoid ò x2 dx dy dz or I2 = ellipsoid ò y2 dx dy dz now becomes easy when we change variables to x', y', and z'. In the new coordinates

I1 = a3bc unit sphereò x'2 dx' dy' dz' and I2 = ab3cunit sphereò y'2 dx' dy' dz'.

For integrals over the unit sphere, unit sphereò x'2 dx' dy' dz' = unit sphereò y'2 dx' dy' dz' because x', y' and z' are all equivalent directions when integrating over the volume of a sphere of radius 1.

(Back to the ellipse for a moment.) For an ellipse, consider an integral like I3 = ò x2 dx dy . Over the ellipse when we change variables this becomes I3 = a3b unit circleò x'2 dx' dy' .

Because ò x'2 dx' dy' = ò y'2 dx' dy' over the unit circle, and r'2 = x'2 + y'2, within the unit circle, we can conclude that I3 = a3b unit circleò1/2 r'2 dx' dy' . But when we integrate, we can use strips of area of length 2p r' and width dr', so then I3 becomes

    I3 = 1/2 a3b 0ò1 r'2 2p r' dr' = p /4 a3b.

Higher powers of x in an integral are a little more difficult, but still much easier over the unit circle.
Take I4 = ellipseò x4 dx dy. This integral becomes I4 = a5 b unit circleò(r' cos q ')4 r' dr' dq '. The r' part is easy and we have I4 = 1/6 a5 b 0ò2p(cos q ')4 dq ' . The angle integral is easily handled with double angle formulas or something similar ( to give 3p /4) .

For integrals over the ellipsoid, things are perhaps simpler than for the circle. When we transform to the unit sphere the integral dx' dy' dz' may be written as r' dq ' r' sin q ' dj ' dr'. If there is no j dependence in the integral, then the volume element is 2 p r'2 dr' sin q ' dq ' . The substitution w' = cos q ' is quite convenient because dw' = - sin q ' dq ' . Then the w integral goes from 1 to -1, and when we flip the limits the negative sign disappears inside the integral.

Exercises:
1) Show that IA = ellipsoid ò x2 dx dy dz = a3bc unit sphereò x'2 dx' dy' dz' = 4p /15 a3bc .
2) Evaluate IB = ellipsoid òx4 dx dy dz .
3) Evaluate IB = ellipsoid òx2 y2 dx dy dz .