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 .