MULTIVARIABLE CALCULUS • MULTIPLE INTEGRALS

Triple Integrals: Cartesian — Set up and evaluate triple integrals in Cartesian coordinates

Extend integration to three dimensions and compute volumes, masses, and other quantities over solid regions.

Historical Context & Motivation

When you first learned about definite integrals, you discovered a way to add up infinitely many tiny quantities along a line — finding areas under curves, total distances, or accumulated change. Then double integrals extended that idea to two dimensions, letting you compute volumes under surfaces. The next natural question is: what happens when we need to accumulate a quantity throughout a three-dimensional solid region? That is exactly the problem that triple integrals were invented to solve.

Throughout the history of mathematics, the ability to handle higher-dimensional integrals grew alongside physics. Scientists needed to calculate the total mass of irregularly shaped objects, the gravitational pull of planets, and the charge distribution inside conductors. Each of these problems requires summing a function over a region of three-dimensional space.

1668
Newton & Leibniz Lay the Foundation
Isaac Newton and Gottfried Leibniz independently develop calculus, introducing the definite integral for functions of a single variable. This sets the stage for everything that follows.
1760s
Euler Explores Multiple Integration
Leonhard Euler extends integration to functions of two and three variables, applying iterated integrals to problems in mechanics and astronomy.
1811
Fourier's Heat Equation
Joseph Fourier uses triple integrals extensively while modeling heat flow through solid bodies, demonstrating that three-dimensional integration is essential for physics.
1900s
Rigorous Foundations
Henri Lebesgue and others formalize measure theory, giving triple (and higher) integrals a rigorous mathematical foundation. These tools now power modern engineering, physics, and data science.

The central question this lesson addresses is: How do we set up and evaluate an integral that accumulates a function f(x, y, z) throughout a solid region in Cartesian coordinates? By the end, you'll be able to identify the limits of integration, choose the right order of integration, and evaluate the resulting iterated integral step by step.

Core Principles & Definitions

Before diving into computation, let's pin down the key ideas. A triple integral extends the concept of single and double integrals into three-dimensional space. Just as a single integral sums thin slices and a double integral sums thin columns, a triple integral sums tiny volume elements (often written as dV) throughout a solid region E.

1

Volume Element dV

In Cartesian coordinates, the volume element is dV = dx dy dz. This represents a tiny rectangular box with side lengths dx, dy, and dz.
2

Region of Integration E

The solid region E is the 3D domain over which you integrate. It can be a box, a tetrahedron, or any bounded solid described by inequalities in x, y, and z.
3

Iterated Integration

A triple integral is evaluated as three nested single integrals. You integrate the innermost variable first, then the middle, then the outermost — just like peeling layers.
4

Order of Integration

There are six possible orders: dx dy dz, dx dz dy, dy dx dz, dy dz dx, dz dx dy, and dz dy dx. The choice depends on which order makes the limits simplest.
5

Fubini's Theorem

If f is continuous on E, you may evaluate the triple integral in any order of integration. This theorem guarantees all six orders yield the same answer.
KEY TAKEAWAY
Think of a triple integral like filling a swimming pool with tiny sugar cubes. Each cube has a volume dV = dx dy dz, and the function f(x, y, z) tells you how much "stuff" is at that location. The triple integral adds up the contribution of every cube throughout the entire pool. When f = 1, you're simply counting cubes — so the integral gives you the total volume of the region.

Visual Explanation

Understanding how to set up a triple integral becomes much easier once you can visualize the region of integration and see how the limits nest together. The diagram below shows a rectangular box region and how we slice it along each axis to build the iterated integral.

A box-shaped region E in 3D space, bounded by constants on each axis. The tiny highlighted cube represents the volume element dV = dx dy dz. The pink lines mark the x-limits, amber marks y-limits, and green marks z-limits.

In the diagram above, the box region E is the simplest type of integration domain. Each pair of limits is just a pair of constants, so the order of integration doesn't affect the difficulty. However, most real problems involve non-rectangular regions where the inner limits depend on the outer variables. In those cases, choosing the right order of integration is crucial for keeping the math manageable.

Mathematical Framework

Let's formalize the setup. Suppose you have a continuous function f(x, y, z) defined over a bounded solid region E in three-dimensional space. The triple integral of f over E is written as follows.

GENERAL TRIPLE INTEGRAL
∫∫∫_E f(x, y, z) dV
f(x, y, z) = the integrand (the function being integrated); E = the solid region of integration; dV = the volume element, equal to dx dy dz in Cartesian coordinates.

To actually evaluate this, we convert it into an iterated integral — three nested single-variable integrals. The most common setup describes E as a Type 1 region where z is sandwiched between two surfaces that depend on x and y, while x and y range over a flat region D in the xy-plane.

TYPE 1 ITERATED INTEGRAL (dz dy dx)
∫ₐᵇ ∫_{g₁(x)}^{g₂(x)} ∫_{u₁(x,y)}^{u₂(x,y)} f(x, y, z) dz dy dx
a, b = constant limits on x (outermost); g₁(x), g₂(x) = lower and upper bounds on y (may depend on x); u₁(x, y), u₂(x, y) = lower and upper bounds on z (may depend on x and y).
VOLUME OF A REGION
Volume(E) = ∫∫∫_E 1 dV
When f(x, y, z) = 1, the triple integral gives the volume of the region E. This is one of the most common applications.
TRIPLE INTEGRAL OVER A BOX
∫ₐᵇ ∫꜀ᵈ ∫ₑᶠ f(x, y, z) dz dy dx
When E is a box [a, b] × [c, d] × [e, f], all six limits are constants and the integration order doesn't matter. This is the simplest case and a good starting point for practice.
💡 Integration Strategy
Always start from the innermost integral and work outward. When evaluating the inner integral, treat the other variables as constants — exactly like partial differentiation in reverse. After integrating with respect to z, you'll have a function of x and y, which you then integrate with respect to y, and so on.

Setting Up the Limits of Integration

The hardest part of triple integrals isn't the integration itself — it's setting up the correct limits. Here's a systematic approach you can follow for any region.

  1. Step 1 — Sketch the region. Draw or visualize the solid E. Identify its bounding surfaces.
  2. Step 2 — Choose the order of integration. Decide which variable to integrate first (innermost). Usually pick the variable that is "sandwiched" between two surfaces.
  3. Step 3 — Find the inner limits. For the innermost variable, express its lower and upper bounds in terms of the remaining variables.
  4. Step 4 — Project onto a plane. Collapse the region onto the plane of the two remaining variables. This gives you the domain D for the middle and outer integrals.
  5. Step 5 — Find the middle and outer limits. Set up limits for D just like you would for a double integral.
A tetrahedron bounded by the plane x + y + z = 1 and the three coordinate planes. The right panel shows how the limits are determined from outer to inner. Notice that inner limits depend on outer variables, forming a nesting structure.

The tetrahedron example above is one of the most common non-rectangular regions you'll encounter. The key insight is that projecting the solid onto the xy-plane gives a triangle, which determines the middle and outer limits. The inner variable z then ranges from the bottom surface (z = 0) up to the top surface (z = 1 − x − y).

⚠️ Common Mistake
Don't write all limits as constants unless the region is a box. If the region has curved or slanted boundaries, at least some limits must be functions of the outer variables. Always check that your limits correctly describe every point (x, y, z) inside E.

Worked Example

Let's work through a complete triple integral from start to finish. We'll evaluate the integral of f(x, y, z) = 6xy over the tetrahedron bounded by x + y + z = 1, x = 0, y = 0, and z = 0.

Evaluate ∫∫∫_E 6xy dV over the tetrahedron with x + y + z ≤ 1, x ≥ 0, y ≥ 0, z ≥ 0
1
Step 1 — Identify the region and choose an orderThe region E is the tetrahedron with vertices at (0,0,0), (1,0,0), (0,1,0), and (0,0,1). We choose the order dz dy dx. The limits are: x from 0 to 1, y from 0 to 1 − x, z from 0 to 1 − x − y.
∫₀¹ ∫₀^(1−x) ∫₀^(1−x−y) 6xy dz dy dx
2
Step 2 — Integrate with respect to z (innermost)Since 6xy contains no z, the z-integration is straightforward. Treat 6xy as a constant with respect to z: ∫₀^(1−x−y) 6xy dz = 6xy × z |₀^(1−x−y) = 6xy(1 − x − y).
Inner integral = 6xy(1 − x − y)
3
Step 3 — Integrate with respect to y (middle)Now integrate 6xy(1 − x − y) with respect to y from 0 to 1 − x. First expand: 6xy(1 − x − y) = 6x(y − xy − y²). Integrate term by term: ∫₀^(1−x) 6x(y − xy − y²) dy = 6x [y²/2 − xy²/2 − y³/3] from 0 to 1 − x. Substitute y = 1 − x: y² = (1−x)², y³ = (1−x)³. = 6x [(1−x)²/2 − x(1−x)²/2 − (1−x)³/3] = 6x × (1−x)² [1/2 − x/2 − (1−x)/3] = 6x × (1−x)² [(3 − 3x − 2 + 2x)/6] = 6x × (1−x)² [(1 − x)/6] = x(1−x)³
Middle integral = x(1 − x)³
4
Step 4 — Integrate with respect to x (outermost)Expand x(1 − x)³ = x(1 − 3x + 3x² − x³) = x − 3x² + 3x³ − x⁴. Now integrate: ∫₀¹ (x − 3x² + 3x³ − x⁴) dx = [x²/2 − x³ + 3x⁴/4 − x⁵/5] from 0 to 1 = 1/2 − 1 + 3/4 − 1/5 = 10/20 − 20/20 + 15/20 − 4/20 = 1/20
∫∫∫_E 6xy dV = 1/20
CHECKING YOUR WORK
Always verify that your answer makes sense. Since 6xy is non-negative throughout the tetrahedron (where x ≥ 0 and y ≥ 0), the integral should be positive — and 1/20 is indeed positive. You can also check units or dimensions: the tetrahedron has volume 1/6, and the average value of 6xy inside it is (1/20) ÷ (1/6) = 6/20 = 3/10, which is reasonable since x and y range from 0 to 1.

Strengths, Limitations & Common Pitfalls

Cartesian coordinates are the default choice for triple integrals, but they aren't always the best choice. Understanding when they shine and when they struggle will save you time and effort.

Cartesian triple integrals: when to use and when to consider alternatives
AspectStrengthLimitation
Region shapeExcellent for boxes, tetrahedra, and regions bounded by planes or polynomials.Awkward for spheres, cones, and cylinders — limits become messy square roots.
IntegrandWorks well when f(x,y,z) is a polynomial or separable function of x, y, z.Difficult when f involves x² + y² or x² + y² + z², which simplify in cylindrical or spherical.
Setup complexityLimits are straightforward when boundaries are planes or simple surfaces.Changing the order of integration may require re-deriving limits from scratch.
Volume elementdV = dx dy dz — no extra Jacobian factor needed.Other coordinate systems require a Jacobian (r for cylindrical, ρ² sin φ for spherical).
🔄 WHEN TO SWITCH COORDINATES
If the region or integrand involves circular symmetry (like x² + y² ≤ r²), switch to cylindrical coordinates. If the region has spherical symmetry (like x² + y² + z² ≤ R²), switch to spherical coordinates. Think of it like choosing the right tool: you can hammer in a screw, but a screwdriver is faster.

Connection to Other Coordinate Systems

The Cartesian triple integral you learned in this lesson is the foundation for all triple integrals, regardless of coordinate system. When you move to cylindrical or spherical coordinates, the same idea applies — you still set up three nested integrals — but the volume element and limits change. Here's a quick comparison to show where your Cartesian skills fit in the bigger picture.

Comparison of coordinate systems for triple integrals
FeatureCartesianCylindricalSpherical
Variablesx, y, zr, θ, zρ, φ, θ
Volume elementdx dy dzr dr dθ dzρ² sin φ dρ dφ dθ
Best forBoxes, tetrahedra, planar boundariesCylinders, cones, circular symmetrySpheres, hemispheres, radial symmetry
Jacobian1 (no extra factor)rρ² sin φ

As you continue in multivariable calculus, you'll learn the change of variables theorem (also called the substitution theorem for multiple integrals), which generalizes all coordinate transformations. The Jacobian factor — like the r in cylindrical or ρ² sin φ in spherical — accounts for how volumes distort under the transformation. Mastering the Cartesian case first gives you a solid reference point, because in Cartesian coordinates the Jacobian is simply 1.

🚀 Looking Ahead
Triple integrals also connect to powerful results like the Divergence Theorem (also known as Gauss's Theorem), which relates a triple integral over a solid region to a surface integral over its boundary. This is a cornerstone of physics and engineering.

Practice Problems

Test your understanding with the following five problems, arranged from conceptual to challenging. Work through each one carefully, and check your reasoning against the provided answers.

PROBLEM 1CONCEPTUAL
Explain in your own words why the innermost limits of a triple integral can depend on the outer variables, but the outermost limits must always be constants.
PROBLEM 2BASIC CALCULATION
Evaluate ∫₀² ∫₀³ ∫₀¹ (x + 2y + 3z) dz dy dx.
PROBLEM 3INTERMEDIATE
Find the volume of the solid region E bounded by z = 0, z = x, y = 0, and y = 4 − x², where x ≥ 0. Set up and evaluate the integral using the order dz dy dx.
PROBLEM 4APPLIED
A solid block occupies the region 0 ≤ x ≤ 2, 0 ≤ y ≤ 1, 0 ≤ z ≤ 3. The density at any point is given by δ(x, y, z) = 4xz kg/m³. Find the total mass of the block by evaluating the appropriate triple integral.
PROBLEM 5CRITICAL THINKING
Consider the integral ∫₀¹ ∫₀ˣ ∫₀ʸ f(x, y, z) dz dy dx. Rewrite this integral with the order of integration changed to dx dy dz. (Hint: first describe the region E as a set of inequalities, then re-derive limits for the new order.)

Lesson Summary

A triple integral in Cartesian coordinates evaluates ∫∫∫_E f(x, y, z) dV by converting the expression into three nested single integrals — an iterated integral. The volume element is dV = dx dy dz, requiring no extra Jacobian factor. By Fubini's Theorem, any of the six possible integration orders yields the same result, provided f is continuous on the bounded region E.

To set up the integral, sketch the region, choose an order, and determine limits from outer (constants) to inner (functions of outer variables). The innermost integral is evaluated first, treating outer variables as constants. When f = 1, the triple integral gives the volume of E. For problems with circular or spherical symmetry, consider switching to cylindrical or spherical coordinates for simpler limits and integrands.

Varsity Tutors • Multivariable Calculus • Triple Integrals: Cartesian