Triple integrals: iterated limits, cylindrical and spherical coordinates
Ten triple integrals worked by hand: a product over a box, the volume and centroid of a tetrahedron, limits read off a slanted plane, the solid under a paraboloid and its centroid in cylindrical coordinates, a moment of inertia, the volume of a ball in spherical coordinates, an ice-cream cone, the integral of ρ² over a ball, a change of order that makes an impossible integral easy, and the probability that a 3-D Gaussian lands in the unit ball, with worked solutions and the mistakes of dropping the r, leaving out sin φ, running φ to 2π, putting numbers where the inner limits need variables, and taking the cone's radius from the sphere.
Before you start
A triple integral adds up a quantity over a solid: its volume when the integrand is ,1, its mass when the integrand is a density, a probability when the integrand is a joint density of three variables. It is computed as three single integrals inside each other, and everything that goes wrong goes wrong in the limits or in the volume element. The inner limits may depend on the outer variables and the outer limits must be numbers. A round solid becomes a box in cylindrical or spherical coordinates, but only if the volume element changes with it: dV=rdzdrdθ in cylindrical and dV=ρ2sinφdρdφdθ in spherical. These ten problems start with boxes and tetrahedra, move to paraboloids, cylinders, balls and cones, change the order of integration to rescue an integral that has no antiderivative, and end with the probability mass of a three-dimensional Gaussian inside the unit ball.
Iterated integrals. Over a region ,E,∭EfdV is computed one variable at a time, innermost first. The innermost limits may depend on the two outer variables, the middle limits on the outer variable only, and the outer limits are numbers. To read off the limits, find where a line in the innermost direction enters and leaves ,E, then describe the shadow of E on the remaining plane as in a double integral.
Volume and averages.∭E1dV is the volume of .E. The centroid has ,xˉ=V1∭ExdV, and the same for ,yˉ,.zˉ. A product f(x)g(y)h(z) over a box splits into three single integrals.
Cylindrical coordinates.,x=rcosθ,,y=rsinθ,,z=z, with r≥0 and .0≤θ<2π. Then x2+y2=r2 and :dV=rdzdrdθ: polar coordinates in the plane with z unchanged.
Spherical coordinates.,x=ρsinφcosθ,,y=ρsinφsinθ,,z=ρcosφ, where ρ≥0 is the distance from the origin, φ∈[0,π] is the angle down from the positive z-axis and θ∈[0,2π) is the polar angle in the xy-plane. Then x2+y2+z2=ρ2 and .dV=ρ2sinφdρdφdθ.
The cone z=x2+y2 is φ=4π in spherical coordinates and z=r in cylindrical. The standard normal CDF is ,Φ, with .Φ(1)=0.8413.
Evaluate ∭BxyzdV over the box .B=[0,1]×[0,2]×[0,3].
·
T is the tetrahedron ,x,y,z≥0,.x+y+z≤1. Find its volume and ,xˉ, the x-coordinate of its centroid.
··
Find the volume of the solid in the first octant under the plane ,2x+3y+z=6, writing the limits in the order .dzdydx.
··
E is the solid between the paraboloid z=x2+y2 and the plane .z=4. Find its volume and the height zˉ of its centroid, in cylindrical coordinates.
·
A solid cylinder ,x2+y2≤4,0≤z≤3 has density .1. Find its moment of inertia about the z-axis, .Iz=∭(x2+y2)dV.
·
Use spherical coordinates to show that a ball of radius a has volume .34πa3.
··
Find the volume of the "ice-cream cone" inside the sphere x2+y2+z2=4 and above the cone ,z=x2+y2, first in spherical and then in cylindrical coordinates.
··
For the unit ball ,B, evaluate ,∭B(x2+y2+z2)dV, and use symmetry to get .∭Bz2dV.
···
Evaluate .∫01∫z1∫y1ex3dxdydz.
···
X=(X1,X2,X3) has independent standard normal coordinates, with density .(2π)−3/2e−∥x∥2/2. Find the probability that X lies in the unit ball.
Answers
29=4.5
V=61 and xˉ=41
V=6
V=8π≈25.13 and zˉ=8π64π/3=38
Iz=24π≈75.40
V=34πa3
V=38π(2−2)≈4.907
∭Bρ2dV=54π and ∭Bz2dV=154π
6e−1≈0.2864
P(∥X∥≤1)≈0.1987
Worked solutions
Problem 1
Evaluate ∭BxyzdV over the box .B=[0,1]×[0,2]×[0,3].
∭BxyzdV=∫01xdx∫02ydy∫03zdzThe integrand is a product of one function of each variable and the limits are all numbers, so the integral splits.
,∫01xdx=21,,∫02ydy=2,∫03zdz=29Three single integrals.
21⋅2⋅29=29Multiply.
Iterated directly: ,∫03zdz=29, then ,∫0229ydy=9, then ∫019xdx=29The same number, innermost first, as a check.
29=4.5Splitting works only when both the integrand and the region are products; over a tetrahedron the limits tie the variables together.
Problem 2
T is the tetrahedron ,x,y,z≥0,.x+y+z≤1. Find its volume and ,xˉ, the x-coordinate of its centroid.
z runs from 0 up to ;1−x−y; the shadow on the xy-plane is the triangle ,x,y≥0,x+y≤1A vertical line through (x,y) enters T at the floor and leaves through the slanted face.
V=∫01∫01−x(1−x−y)dydx=∫012(1−x)2dx=61The inner y-integral of c−y from 0 to ,c, with ,c=1−x, is .2c2.
∭TxdV=∫01x2(1−x)2dx=21(21−32+41)=241The same inner integrals, with the extra factor x carried along.
xˉ=1/61/24=41Divide by the volume.
V=61 and xˉ=41By symmetry the centroid is :(41,41,41): a quarter of the way from the base to the opposite vertex, as for any tetrahedron.
Problem 3
Find the volume of the solid in the first octant under the plane ,2x+3y+z=6, writing the limits in the order .dzdydx.
0≤z≤6−2x−3yFloor ,z=0, roof the plane solved for .z.
0≤y≤36−2xThe shadow is where the roof is above the floor: 2x+3y≤6 with .x,y≥0. A line in the y-direction leaves it at .3y=6−2x.
0≤x≤3The shadow's extent in ,x, where 2x=6 on the x-axis. Only this outer range is a pair of numbers.
∫03∫0(6−2x)/3(6−2x−3y)dydx=∫036(6−2x)2dx=6With ,c=6−2x, the inner integral of c−3y from 0 to 3c is ;3c2−6c2=6c2; then .∫036(6−2x)2dx=61⋅6216=6.
V=6Check with the formula for a tetrahedron, 61 of the product of the intercepts: .61⋅3⋅2⋅6=6. Taking the limits as 0 to ,6,0 to 2 and 0 to 3 gives the box around it, 36 (Mistake 4).
Problem 4
E is the solid between the paraboloid z=x2+y2 and the plane .z=4. Find its volume and the height zˉ of its centroid, in cylindrical coordinates.
,r2≤z≤4,,0≤r≤2,0≤θ≤2πThe paraboloid is ;z=r2; it meets z=4 at .r=2.
V=∫02π∫02∫r24rdzdrdθ=2π∫02(4−r2)rdrThe factor r is part of ;dV; nothing in the integrand depends on .θ.
2π[2r2−4r4]02=2π(8−4)=8πLeaving out the r gives 332π (Mistake 1).
∭EzdV=2π∫02216−r4rdr=π(32−332)=364πThe inner integral is .∫r24zdz=216−r4.
V=8π≈25.13 and zˉ=8π64π/3=38The centroid is above the middle height 2 because the solid is wide at the top and narrows to a point at the bottom.
Problem 5
A solid cylinder ,x2+y2≤4,0≤z≤3 has density .1. Find its moment of inertia about the z-axis, .Iz=∭(x2+y2)dV.
x2+y2=r2 and ,dV=rdzdrdθ, with ,0≤z≤3,0≤r≤2In cylindrical coordinates the cylinder is a box.
Iz=∫02π∫02∫03r2⋅rdzdrdθThe integrand r2 times the r from dV gives .r3.
=2π⋅424⋅3=24πThe integral splits: ,∫02πdθ=2π,,∫02r3dr=4,.∫03dz=3.
The mass is ,M=π⋅22⋅3=12π, so Iz=21MR2 with R=2The textbook formula for a solid cylinder, as a check.
Iz=24π≈75.40Mass far from the axis counts with the square of its distance, which is why a hollow cylinder of the same mass has twice this moment.
Problem 6
Use spherical coordinates to show that a ball of radius a has volume .34πa3.
,0≤ρ≤a,,0≤φ≤π,0≤θ≤2πφ only needs to reach ,π, the negative z-axis; θ supplies the full turn.
V=∫02π∫0π∫0aρ2sinφdρdφdθ.dV=ρ2sinφdρdφdθ.
,∫0aρ2dρ=3a3,,∫0πsinφdφ=2,∫02πdθ=2πThe integrand is a product and the limits are numbers, so it splits.
V=2π⋅2⋅3a3Multiply.
V=34πa3Without the sinφ the result is 32π2a3≈6.58a3 (Mistake 2), and running φ to 2π gives 0 (Mistake 3).
Problem 7
Find the volume of the "ice-cream cone" inside the sphere x2+y2+z2=4 and above the cone ,z=x2+y2, first in spherical and then in cylindrical coordinates.
Spherical: ,0≤ρ≤2,,0≤φ≤4π,0≤θ≤2πThe cone makes the angle 4π with the z-axis, and the solid is the part of the ball within that angle.
V=2π[−cosφ]0π/4⋅323=2π(1−22)38The integral splits as in Problem 6.
Cylindrical: ,r≤z≤4−r2, and the cone meets the sphere where ,r2+r2=4, at r=2So ,0≤r≤2, not .2.
2π∫02(4−r2−r)rdr=2π(38−22−322)∫02r4−r2dr=[−31(4−r2)3/2]02=38−22 and .∫02r2dr=322.
V=38π(2−2)≈4.907Both routes agree. Spherical coordinates turn the solid into a box, which is why the first route is one line. Running r to 2 in the second route gives 0 (Mistake 5).
Problem 8
For the unit ball ,B, evaluate ,∭B(x2+y2+z2)dV, and use symmetry to get .∭Bz2dV.
x2+y2+z2=ρ2The integrand depends only on the distance from the origin.
∭Bx2dV=∭By2dV=∭Bz2dVThe ball looks the same after relabeling the axes, so the three integrals are equal and add to .54π.
∭Bρ2dV=54π and ∭Bz2dV=154πDirectly: .∫02π∫0π∫01ρ4cos2φsinφdρdφdθ=2π⋅32⋅51=154π. Dividing by the volume ,34π, the average of ρ2 over the ball is .53.
Problem 9
Evaluate .∫01∫z1∫y1ex3dxdydz.
As written the inner integral ∫y1ex3dx has no elementary antiderivativeSo change the order.
The region is 0≤z≤y≤x≤1Read the inequalities off the limits: z≤y from the middle, y≤x from the inner, and all of them in .[0,1].
In the order :dzdydx:,0≤z≤y,,0≤y≤x,0≤x≤1For fixed x and ,y,z runs up to ;y; for fixed ,x,y runs up to .x.
∫01∫0x∫0yex3dzdydx=∫01ex3∫0xydydx=∫012x2ex3dxex3 is constant in the two inner integrals.
6e−1≈0.2864With ,u=x3,.∫012x2ex3dx=61∫01eudu. The factor 2x2 is the area of the triangle 0≤z≤y≤x in the yz-plane, and it is exactly what the substitution needed.
Problem 10
X=(X1,X2,X3) has independent standard normal coordinates, with density .(2π)−3/2e−∥x∥2/2. Find the probability that X lies in the unit ball.
P(∥X∥≤1)=(2π)−3/2∫02π∫0π∫01e−ρ2/2ρ2sinφdρdφdθThe density depends only on ,ρ, so spherical coordinates turn the ball into a box.
=(2π)−3/2⋅4π∫01ρ2e−ρ2/2dρ=π2∫01ρ2e−ρ2/2dρThe angles give ,4π, the area of the unit sphere, and .4π(2π)−3/2=2/π.
∫01ρ2e−ρ2/2dρ=[−ρe−ρ2/2]01+∫01e−ρ2/2dρ=−e−1/2+2π(Φ(1)−21)Parts with u=ρ and ;dv=ρe−ρ2/2dρ; the last integral is a normal probability.
P=2Φ(1)−1−π2e−1/2≈0.6827−0.4839Multiply through by .2/π.
P(∥X∥≤1)≈0.1987Each coordinate is within 1 of 0 with probability ,0.6827, yet the vector is within distance 1 less than 20% of the time. In d dimensions ∥X∥ concentrates near ,d, so a high-dimensional Gaussian has almost no mass near its mode.
Where this goes wrong
1. Dropping the r in cylindrical coordinates
The integrand is rewritten in r and ,θ, and dV is taken as .dzdrdθ.
,r2≤z≤4,,0≤r≤2,0≤θ≤2πRight so far: Problem 4.
“.V=∫02π∫02∫r24dzdrdθ.”The step that causes the mistake: a box dzdrdθ has volume ,rdzdrdθ, because the arc it spans grows with .r.
V=2π∫02(4−r2)dr=332π≈33.51Compare the cylinder ,r≤2,,0≤z≤4, of volume :16π≈50.27: the paraboloid cuts away exactly half of it, so the answer is .8π≈25.13. Without the ,r, a thin ring near the axis counts as much as a ring of the same width at the rim.
2. Spherical coordinates without sin φ
The volume element is taken as .ρ2dρdφdθ.
,0≤ρ≤a,,0≤φ≤π,0≤θ≤2πRight so far: Problem 6.
“.dV=ρ2dρdφdθ.”The step that causes the mistake: a circle of latitude at angle φ has radius ,ρsinφ, which shrinks to 0 at the poles, and the volume element has to shrink with it.
V=2π⋅π⋅3a3=32π2a3≈6.58a3The true volume is .34πa3≈4.19a3. Without ,sinφ, the thin polar regions count as if they were as wide as the equator.
3. φ running from 0 to 2π
Both angles are given the full turn.
dV=ρ2sinφdρdφdθ and 0≤θ≤2πRight so far: Problem 6.
“.0≤φ≤2π.”φ is measured from the positive z-axis and reaches the negative z-axis at ;π;θ already sweeps all the way around.
,∫02πsinφdφ=0, so V=0A volume of 0 for a ball is a sure sign. For ,π<φ<2π,sinφ<0 and the ball is counted a second time with a negative sign. With ,0≤φ≤π,.V=34πa3.
4. Numbers where the inner limits need variables
The limits are read off the intercepts, as if the solid were a box.
The plane 2x+3y+z=6 meets the axes at ,x=3,,y=2,z=6Right so far: Problem 3.
“,0≤z≤6,,0≤y≤2,.0≤x≤3.”Only the outer limits can be numbers; the roof z=6−2x−3y comes down as x and y grow.
V=3⋅2⋅6=36That is the box around the tetrahedron, six times its volume. With z up to 6−2x−3y and y up to ,36−2x,.V=6.
5. The cone's radius taken from the sphere
In cylindrical coordinates the outer radius is read off the sphere.
r≤z≤4−r2Right so far: Problem 7.
“The sphere has radius ,2, so .0≤r≤2.”The solid only reaches ,r=2, where the cone meets the sphere. Beyond it the "floor" z=r is above the "roof" .4−r2.
V=2π∫02(4−r2−r)rdr=2π(38−38)=0The negative part from 2 to 2 exactly cancels the solid. Set the floor equal to the roof to find the edge of the shadow: r=4−r2 gives r=2 and .V=38π(2−2)≈4.907.