Multivariable Calculus — Topic 05

Triple Integrals
Volume & Mass in Three Dimensions

Double integrals extend single-variable integration to flat regions. Triple integrals push that idea one dimension further, integrating over solid 3D regions to compute volumes, masses, centres of mass, and moments of inertia. The key challenge is always setting up the correct limits.

6Orderings of Integration
10Worked Examples
5SVG Diagrams
10Quiz Questions
Share this page
On This Page
  1. From Double to Triple Integrals
  2. Integration Over a Box, Fubini's Theorem
  3. Integration Over a General Region
  4. Setting Up Limits and Switching Order
  5. Applications: Volume, Mass, Centre of Mass
  6. Triple Integrals in Cylindrical Coordinates
  7. Triple Integrals in Spherical Coordinates
  8. Worked Examples
  9. Common Mistakes
  10. 10-Question Quiz

§ 01From Double to Triple Integrals

A double integral ∬ f(x,y) dA sums the values of a function over a 2D region. A triple integral ∭ f(x,y,z) dV does the same over a solid 3D region E.

The idea is the same: slice E into tiny volume elements dV, evaluate f at each point, multiply, and sum. In Cartesian coordinates, dV = dx dy dz. The resulting triple integral is evaluated as three nested single-variable integrals, working from the innermost outward.

Triple Integral — Definition
E f(x,y,z) dV = ∫∫∫ f(x,y,z) dx dy dz

When f = 1, the triple integral gives the volume of the region E. When f = ρ(x,y,z) (density), it gives the mass. The integrand determines what physical quantity is computed.

A solid region E subdivided into tiny boxes dV = dx dy dz
y z x E dV = dx dy dz E f dV = lim Σ f·ΔV
The triple integral is the limit of a Riemann sum: slice E into tiny boxes, evaluate f at a sample point in each box, multiply by the box's volume ΔV, and sum over all boxes inside E. As the boxes shrink, the sum converges to the triple integral.

§ 02Integration Over a Box — Fubini's Theorem

The simplest region for a triple integral is a rectangular box: B = [a,b] × [c,d] × [e,f], meaning a ≤ x ≤ b, c ≤ y ≤ d, e ≤ z ≤ f. All limits are constant numbers.

Fubini's Theorem for Triple Integrals
B f dV = ∫abcdef f(x,y,z) dz dy dx

On a rectangular box with continuous f, the triple integral equals the iterated integral in any of the six possible orderings. You are free to choose whichever order makes the arithmetic easiest.

All Six Integration Orders

For a box B = [a,b]×[c,d]×[e,f], the six orderings all give the same result:

Order 1
dz dy dx
Order 2
dz dx dy
Order 3
dy dz dx
Order 4
dy dx dz
Order 5
dx dz dy
Order 6
dx dy dz
How to Choose the Order on a Box If f(x,y,z) factors as g(x)·h(y)·k(z), then the triple integral factors into three independent single integrals — any order gives the same result trivially. If f does not factor, choose the order that produces the simplest inner integral first.

§ 03Integration Over a General Region

For a non-box region E, the limits of the inner integrals depend on the outer variables. There are three standard types, depending on which axis is "stacked".

Type 1 — z between surfaces, (x,y) in 2D region D
E f dV = ∬Du₁(x,y)u₂(x,y) f dz dA

E lies above the surface z = u₁(x,y) and below z = u₂(x,y). The projection of E onto the xy-plane is the 2D region D. Integrate z first (between the surfaces), then integrate the result over D using a double integral.

Type 2 — y between surfaces, (x,z) in 2D region D
E f dV = ∬Dv₁(x,z)v₂(x,z) f dy dA

Same idea but y is the "inner" variable, bounded by two surfaces y = v₁(x,z) and y = v₂(x,z). The 2D region D is the projection onto the xz-plane.

Type 3 — x between surfaces, (y,z) in 2D region D
E f dV = ∬Dw₁(y,z)w₂(y,z) f dx dA

x is the inner variable. In practice this ordering is less common but occasionally the most natural choice.

Type 1 region: z sandwiched between two surfaces
y z x D z = u₁(x,y) z = u₂(x,y) z
For a Type 1 region, integrate z first between the two bounding surfaces, then integrate over the projected region D in the xy-plane.
Reading off limits: from innermost to outermost
x: a → b (constants) y: g₁(x) → g₂(x) (fn of x) z: u₁(x,y)→u₂(x,y) (fn of x,y) f(x, y, z) dz dy dx Innermost limits may depend on all outer vars.
The key nesting rule: the innermost limits can depend on all outer variables, the middle limits can depend only on the outermost variable, and the outermost limits must be constants.

§ 04Setting Up Limits and Switching Order

Setting up the correct limits is the hardest skill in triple integration. A reliable 4-step strategy works for most problems.

The 4-Step Strategy for Setting Up Limits

1
Sketch the region E in 3D (or at least its bounding surfaces). Identify the top and bottom surfaces, and any side constraints.
2
Choose the order of integration. Usually z is integrated first (between the top and bottom surfaces), then the remaining double integral is over the projection D in the xy-plane.
3
Write the z-limits: z runs from the bottom surface z = u₁(x,y) to the top surface z = u₂(x,y).
4
Set up the double integral over D: find the projection of E onto the xy-plane, determine whether D is Type I (x fixed, y from curve to curve) or Type II (y fixed, x from curve to curve), and write the remaining limits accordingly.

Switching the Order of Integration

Sometimes an integral in one order is impossible to evaluate in closed form, but in another order it becomes elementary. The strategy for switching order is:

1
From the original limits, describe the region E in plain English: what are the geometric constraints on x, y, and z?
2
Sketch E (or its projection onto each coordinate plane) to confirm your description.
3
Re-derive the limits in the new desired order by asking: for a fixed value of the outermost variable, what range does the next variable take? Then for those two fixed, what range does the innermost variable take?
The Limits Must Always Form a Consistent Nest Outermost limits: constants only. Middle limits: may depend on the outermost variable only. Innermost limits: may depend on both outer variables. Breaking this rule — for example, putting an x in the outermost limits when it has not yet been integrated — is the most common structural error.

§ 05Applications: Volume, Mass, and Centre of Mass

Triple integrals are the natural tool for any quantity that is "spread through" a 3D region.

Application 1
Volume
V = ∭E 1 dV

Setting f = 1 integrates the constant function 1 over all of E. The result is the total volume. This is almost always the simplest possible triple integral.

Application 2
Mass
M = ∭E ρ(x,y,z) dV

If ρ is the mass density (mass per unit volume), integrating ρ over the solid gives the total mass. Uniform density ρ = constant gives M = ρ·V.

Application 3
Centre of Mass
x̄ = (1/M) ∭ x ρ dV

The x-coordinate of the centre of mass is the weighted average of x over E. Similarly for ȳ and z̄. The centre of mass (x̄, ȳ, z̄) is the balance point of the solid.

Application 4
Moments of Inertia
Iz = ∭E (x²+y²) ρ dV

The moment of inertia about the z-axis measures resistance to rotation. Similar formulas with (y²+z²) and (x²+z²) give Ix and Iy.

§ 06Triple Integrals in Cylindrical Coordinates

When the region E has cylindrical symmetry (a cylinder, cone, or anything involving x² + y²) cylindrical coordinates (r, θ, z) transform the integral into a much simpler form.

Cylindrical Coordinates Formula
E f dV = ∫∫∫ f(r cosθ, r sinθ, z) · r dz dr dθ

The key factor is the Jacobian r. Every cylindrical-coordinate triple integral must include this factor. Also: x² + y² becomes r², which dramatically simplifies most integrands.

Conversions
x = r cosθ, y = r sinθ, z = z
The xy part is standard 2D polar. z is unchanged.
Key identity
x² + y² = r²
Replace immediately whenever this expression appears in the integrand or limits.
Volume element
dV = r dz dr dθ
The Jacobian r is essential and must never be omitted.
Standard limits
0 ≤ r, 0 ≤ θ ≤ 2π, z free
For a full cylinder of radius R: 0 ≤ r ≤ R, 0 ≤ θ ≤ 2π, h₁ ≤ z ≤ h₂.

§ 07Triple Integrals in Spherical Coordinates

When E is a sphere, hemisphere, or any region where x² + y² + z² appears, spherical coordinates (ρ, θ, φ) give the simplest form. See the Cylindrical & Spherical page for full derivations.

Spherical Coordinates Formula
E f dV = ∫∫∫ f(ρ sinφ cosθ, ρ sinφ sinθ, ρ cosφ) · ρ² sinφ dρ dφ dθ

The Jacobian is ρ² sinφ. Both factors are required. For a full sphere of radius R: 0 ≤ ρ ≤ R, 0 ≤ θ ≤ 2π, 0 ≤ φ ≤ π.

Region TypeBest SystemKey Reason
Box / planesCartesianConstant limits; no simplification from other systems
Cylinder, coneCylindricalx²+y² = r²; Jacobian r; limits become constant
Sphere, hemisphereSphericalx²+y²+z² = ρ²; Jacobian ρ² sinφ; limits constant
Mixed (cone + sphere)SphericalBoth surfaces simplify in spherical — cone is φ=const, sphere is ρ=const
ParaboloidCylindricalz = r² describes the paraboloid; limits become simple

§ 08Worked Examples

Example 01 Triple integral over a box — ∭B xyz dV, B = [0,2]×[0,1]×[0,3]
Setup

All limits are constants (a box), so Fubini applies immediately. Choose order dz dy dx.

020103 xyz dz dy dx
Factorise

The integrand xyz = x · y · z factors completely. Apply Fubini's product rule:

= [∫02 x dx] · [∫01 y dy] · [∫03 z dz]
Each factor
02 x dx = [x²/2]02 = 2 01 y dy = [y²/2]01 = 1/2 03 z dz = [z²/2]03 = 9/2
Result
∭ xyz dV = 2 · (1/2) · (9/2) = 9/2
Example 02 Volume of the tetrahedron with vertices (0,0,0), (1,0,0), (0,1,0), (0,0,1)
Region

The tetrahedron is bounded by x ≥ 0, y ≥ 0, z ≥ 0, and the plane x + y + z = 1 (equivalently z ≤ 1 − x − y).

Set up

Integrate z first from 0 to 1−x−y. Then the projection D onto the xy-plane is the triangle x ≥ 0, y ≥ 0, x + y ≤ 1.

V = ∫0101−x01−x−y dz dy dx
z-integral
01−x−y dz = 1 − x − y
y-integral
01−x (1−x−y) dy = [(1−x)y − y²/2]01−x = (1−x)²/2
x-integral
01 (1−x)²/2 dx = [−(1−x)³/6]01 = 1/6
Result
V = 1/6

Confirmed by the formula V = (1/6)|det[edges]| = 1/6 for this standard tetrahedron ✓

Example 03 E z dV where E is under z = 4 − x² − y² and above the xy-plane
Region

The paraboloid z = 4 − x² − y² meets z = 0 when x² + y² = 4 (circle of radius 2). So E: 0 ≤ z ≤ 4−x²−y², (x,y) in the disk x²+y² ≤ 4.

z-integral
04−x²−y² z dz = [z²/2]04−x²−y² = (4−x²−y²)²/2
Switch to polar

The remaining double integral is over the disk x²+y² ≤ 4. Switch to polar: x²+y² = r², dA = r dr dθ, r ∈ [0,2], θ ∈ [0,2π].

∬ (4−r²)²/2 · r dr dθ = ∫002 r(4−r²)²/2 dr dθ
r-integral

Let u = 4−r², du = −2r dr:

02 r(4−r²)²/2 dr = (1/2)∫40 u²(−du/2) = (1/4)∫04 u² du = (1/4)(64/3) = 16/3
θ-integral
0 (16/3) dθ = 32π/3
Result
∭ z dV = 32π/3
Example 04 Cylindrical: volume inside x² + y² = 9 between z = 0 and z = x + 4
Convert

x² + y² = 9 → r = 3. The plane z = x + 4 becomes z = r cosθ + 4. Limits: 0 ≤ r ≤ 3, 0 ≤ θ ≤ 2π, 0 ≤ z ≤ r cosθ + 4.

z-integral
0r cosθ+4 dz = r cosθ + 4
r-integral
03 (r cosθ + 4) · r dr = cosθ [r³/3]03 + 4[r²/2]03 = 9 cosθ + 18
θ-integral
0 (9 cosθ + 18) dθ = [9 sinθ + 18θ]0 = 0 + 36π = 36π
Result
V = 36π

The ∫cosθ term vanishes over a full revolution, the tilted plane contributes zero net volume above the symmetric cylinder.

Example 05 Cylindrical: ∭ (x²+y²) dV over the solid cone z = √(x²+y²), z ≤ 2
Convert

Cone: z = r. Region: r ≤ z ≤ 2, which requires r ≤ 2. So 0 ≤ r ≤ 2, r ≤ z ≤ 2, 0 ≤ θ ≤ 2π. Integrand x²+y² = r².

002r2 r² · r dz dr dθ
z-integral
r2 dz = 2 − r
r-integral
02 r³(2−r) dr = ∫02 (2r³−r⁴) dr = [r⁴/2 − r⁵/5]02 = 8 − 32/5 = 8/5
θ-integral
0 (8/5) dθ = 16π/5
Result
∭ (x²+y²) dV = 16π/5
Example 06 Spherical: volume of the unit ball ∭B dV
Setup

Unit ball: 0 ≤ ρ ≤ 1, 0 ≤ θ ≤ 2π, 0 ≤ φ ≤ π. Integrand f = 1.

V = ∫00π01 ρ² sinφ dρ dφ dθ
Factor

The integrand ρ² sinφ separates completely:

= [∫01 ρ² dρ] · [∫0π sinφ dφ] · [∫0 dθ] = (1/3) · 2 · 2π = 4π/3
Result
V = 4π/3 ✓

The classical formula for the volume of a sphere of radius 1, derived in three lines.

Example 07 Spherical: ∭E e^(x²+y²+z²)^(3/2) dV, E is the ball of radius 2
Recognise

x²+y²+z² = ρ², so (x²+y²+z²)^(3/2) = ρ³. Integrand becomes e^(ρ³). This is impossible in Cartesian.

00π02 e^(ρ³) · ρ² sinφ dρ dφ dθ
Separate
= [∫02 ρ²e^(ρ³)dρ] · [∫0πsinφ dφ] · [∫0dθ]
ρ-integral

Let u = ρ³, du = 3ρ² dρ → ρ² dρ = du/3:

02 ρ²e^(ρ³)dρ = (1/3)[e^(ρ³)]02 = (e⁸ − 1)/3
Combine
= (e⁸−1)/3 · 2 · 2π = 4π(e⁸−1)/3
Result
∭ e^(ρ³) dV = 4π(e⁸−1)/3
Example 08 Switching order: ∫0101y1 e^(x²) dx dy dz
Problem

The inner integral ∫y1 e^(x²) dx has no closed form (e^(x²) has no elementary antiderivative). We must switch order.

Describe E

From the limits: 0 ≤ z ≤ 1 (z-variable is free, constant limits), 0 ≤ y ≤ 1, y ≤ x ≤ 1. The region in the xy-plane: 0 ≤ y ≤ x, 0 ≤ x ≤ 1 (triangle). So E = {(x,y,z) : 0 ≤ z ≤ 1, 0 ≤ y ≤ x, 0 ≤ x ≤ 1}.

Switch to dx last

New order: dz dy dx. Since z is independent, keep 0 ≤ z ≤ 1. For the xy-triangle, switch to: 0 ≤ x ≤ 1, 0 ≤ y ≤ x.

= ∫010x01 e^(x²) dz dy dx
z-integral
01 dz = 1
y-integral
0x e^(x²) dy = x e^(x²)

e^(x²) is treated as a constant when integrating in y.

x-integral
01 x e^(x²) dx = [e^(x²)/2]01 = (e−1)/2
Result
∭ e^(x²) dV = (e−1)/2 ≈ 0.859
Example 09 Mass of a solid hemisphere ρ ≤ 3, z ≥ 0, with density δ = z
Setup

Upper hemisphere: 0 ≤ ρ ≤ 3, 0 ≤ θ ≤ 2π, 0 ≤ φ ≤ π/2. Density δ = z = ρ cosφ.

M = ∫00π/203 ρ cosφ · ρ² sinφ dρ dφ dθ
Separate
= [∫03ρ³ dρ][∫0π/2sinφ cosφ dφ][∫0dθ]
Each factor
03ρ³ dρ = 81/4 0π/2sinφ cosφ dφ = [sin²φ/2]0π/2 = 1/2 0dθ = 2π
Result
M = (81/4)(1/2)(2π) = 81π/4
Example 10 Centre of mass of a solid pyramid: base [0,1]×[0,1] in xy-plane, apex at (0,0,1), uniform density
Region

The pyramid has square base 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, and apex at z = 1. At height z, the cross-section is a square of side (1−z). So: 0 ≤ z ≤ 1, 0 ≤ x ≤ 1−z, 0 ≤ y ≤ 1−z.

Volume V
V = ∫0101−z01−z dy dx dz = ∫01(1−z)² dz = [−(1−z)³/3]01 = 1/3
By symmetry

The pyramid is symmetric about x = (1−z)/2 and y = (1−z)/2. So x̄ = ȳ = 1/4 by the symmetry of a uniform square section.

∭ z dV = ∫01 z(1−z)² dz = ∫01(z − 2z² + z³)dz = 1/2 − 2/3 + 1/4 = 1/12 z̄ = (1/12) / (1/3) = 1/4
Centre of Mass
(x̄, ȳ, z̄) = (1/4, 1/4, 1/4)

The centre of mass lies one-quarter of the way from the base to the apex, consistent with the general rule that a solid pyramid's centroid is at 1/4 of the height.

§ 09Common Mistakes

Mistake 1 — Wrong Nesting of Limits The outermost limits must be constants. The middle limits may depend only on the outermost variable. The innermost limits may depend on both outer variables. Putting a function of y in the outermost integral limits, for example, is a fundamental structural error that makes the integral meaningless.
Mistake 2 — Forgetting the Jacobian in Cylindrical or Spherical In cylindrical: dV = r dr dθ dz — the factor r must appear. In spherical: dV = ρ² sinφ dρ dφ dθ — both factors ρ² and sinφ must appear. Omitting either Jacobian gives an answer that is wrong by a multiplicative factor that grows with the domain.
Mistake 3 — Treating All Orders as Interchangeable for Non-Box Regions Fubini's theorem lets you swap order freely on a box with constant limits. On a general region E, the limits themselves change when you change order — you must re-derive the bounds from the geometric description of E, not just swap dz and dy naively.
Mistake 4 — Incorrect φ Limits in Spherical For the full sphere, φ runs from 0 to π. For the upper hemisphere only, φ runs from 0 to π/2. Using π where π/2 is needed doubles the answer. Always draw the region and check which polar angles are included.
Mistake 5 — Using Cartesian When a Better System Exists Computing the volume of a sphere in Cartesian coordinates requires evaluating three nested square-root expressions — it is technically possible but extremely tedious. Always check: does the region involve x²+y² (use cylindrical) or x²+y²+z² (use spherical)? Choosing the right coordinate system is not optional — it is the key skill.

§ 1010-Question Quiz

Test your understanding of iterated triple integrals, limits, applications, and coordinate systems.

Quiz — Triple Integrals  |  10 Questions

Cookie Settings