MULTIVARIABLE CALCULUS • MULTIPLE INTEGRALS

Iterated Integrals — Set up and evaluate iterated integrals

Learn how to compute double integrals by breaking them into two successive single integrals.

Historical Context & Motivation

Long before calculators or computers existed, mathematicians faced a fundamental challenge: how do you find the volume under a curved surface, or the total mass of a flat plate whose density varies from point to point? Single-variable calculus gave us the definite integral to handle one-dimensional problems — areas under curves and accumulated quantities along a line. But the real world is at least two-dimensional, and extending integration to handle functions of two (or more) variables required a powerful new idea.

The breakthrough came from realizing that a two-dimensional integral can be broken into two iterated (repeated) one-dimensional integrals — one nested inside the other. This insight developed over centuries, beginning with early ideas about volumes and areas, and culminating in the rigorous theorems that make modern calculus work.

1635
Cavalieri's Principle
Italian mathematician Bonaventura Cavalieri proposed that volumes can be computed by slicing a solid into infinitely thin cross-sections and summing their areas — a precursor to iterated integration.
1680s
Newton & Leibniz Develop Calculus
Isaac Newton and Gottfried Wilhelm Leibniz independently created the foundations of single-variable calculus, including the definite integral and the Fundamental Theorem of Calculus.
1740s
Euler Extends Integration
Leonhard Euler applied repeated integration to problems in physics and mechanics, routinely computing double integrals by integrating first with respect to one variable, then the other.
1907
Fubini's Theorem
Italian mathematician Guido Fubini published the rigorous theorem guaranteeing that, under reasonable conditions, a double integral can always be evaluated as two iterated single integrals — and the order of integration can be swapped.

The central question that iterated integrals answer is straightforward: How do we actually compute a double integral? Defining a double integral as a limit of Riemann sums over rectangles is elegant in theory, but nearly impossible to compute directly. Iterated integrals give us a practical, step-by-step method — integrate one variable at a time — that converts a seemingly complex two-dimensional problem into a sequence of familiar single-variable integrations.

Core Principles & Definitions

Before diving into calculations, let's nail down the key ideas that make iterated integrals work. If you're comfortable evaluating a definite integral like ∫₀² x² dx, then you already have the core skill — iterated integrals simply ask you to do that process twice, one variable at a time.

1

Double Integral

The integral ∬R f(x, y) dA represents the total "accumulated value" of f over a 2-D region R. Geometrically, when f ≥ 0, it gives the volume under the surface z = f(x, y).
2

Iterated Integral

An iterated integral evaluates the double integral by performing two single integrals in sequence. You integrate with respect to one variable (the inner integral), then integrate the result with respect to the other variable (the outer integral).
3

Order of Integration

You can integrate with respect to y first (dy dx) or x first (dx dy). Fubini's Theorem guarantees both orders give the same answer when f is continuous on R.
4

Limits of Integration

The inner integral's limits may be constants (rectangular region) or functions of the outer variable (non-rectangular region). The outer integral's limits are always constants.
KEY TAKEAWAY
Think of an iterated integral like mowing a rectangular lawn in rows. First, you mow one complete strip from left to right (the inner integral). Then you step forward and mow the next strip, repeating until you've covered the whole lawn top to bottom (the outer integral). Each strip is a single-variable integral, and together they cover the entire two-dimensional region.

Visual Explanation

The diagram below shows how a double integral over a rectangular region R = [a, b] × [c, d] is broken into slices. We fix a value of x and integrate f(x, y) with respect to y from c to d — this produces a single number A(x) that represents the area of one vertical "slice." Then we integrate A(x) from a to b to sweep through all slices and accumulate the total volume.

The shaded cyan strip represents one "slice" at a fixed value x = x₀. The inner integral computes the area A(x₀) of that slice by integrating f(x₀, y) from c to d. The outer integral then sums all such slices from x = a to x = b.

Notice the key idea in the diagram: the inner integral treats x as a constant and produces a function of x alone. Once you have A(x), the outer integral is just an ordinary single-variable integral. This two-step process is why we call it an "iterated" integral — we iterate (repeat) the integration process.

Mathematical Framework

Let's formalize the ideas from the previous section. The mathematical machinery involves Fubini's Theorem and two standard forms for setting up iterated integrals.

FUBINI'S THEOREM (RECTANGULAR REGION)
∬_R f(x, y) dA = ∫_a^b ∫_c^d f(x, y) dy dx = ∫_c^d ∫_a^b f(x, y) dx dy
If f(x, y) is continuous on the rectangle R = [a, b] × [c, d], then the double integral equals either iterated integral. Here dA represents the area element, and the two forms correspond to integrating y-first or x-first.
ITERATED INTEGRAL — dy dx ORDER
∫_a^b [ ∫_c^d f(x, y) dy ] dx
Step 1 (inner): Hold x constant, integrate f(x, y) with respect to y from c to d. This yields a function of x alone. Step 2 (outer): Integrate that result with respect to x from a to b.
GENERAL REGION — TYPE I (dy dx)
∫_a^b ∫_{g₁(x)}^{g₂(x)} f(x, y) dy dx
For a Type I region, y ranges between two curves y = g₁(x) and y = g₂(x), while x ranges between constants a and b. The inner limits depend on x.
GENERAL REGION — TYPE II (dx dy)
∫_c^d ∫_{h₁(y)}^{h₂(y)} f(x, y) dx dy
For a Type II region, x ranges between two curves x = h₁(y) and x = h₂(y), while y ranges between constants c and d. The inner limits depend on y.
💡 Reading the Notation
Always read an iterated integral from the inside out. The innermost integral sign pairs with the innermost differential (dy or dx). When you see ∫₀¹ ∫₀ˣ f(x, y) dy dx, you first integrate with respect to y (from 0 to x), then integrate the result with respect to x (from 0 to 1).

Type I vs. Type II Regions

Understanding the shape of your region of integration is the most important step in setting up an iterated integral. Regions are classified as Type I or Type II based on which variable has the "moving" limits. Many regions can be described as either type, and choosing wisely can make the integral much easier to evaluate.

Left: a Type I region where, for each fixed x, y ranges between two curves g₁(x) and g₂(x). Right: a Type II region where, for each fixed y, x ranges between h₁(y) and h₂(y). The pink arrows show the direction of the inner integration.
Comparison of Type I and Type II iterated integral setups
FeatureType I (dy dx)Type II (dx dy)
Inner variableyx
Inner limitsg₁(x) to g₂(x)h₁(y) to h₂(y)
Outer variablexy
Outer limitsConstants a to bConstants c to d
Best when...Region is easily described by top/bottom curves in yRegion is easily described by left/right curves in x

Worked Example

Let's walk through a complete example to see the process in action. We will evaluate the double integral of f(x, y) = 2x + 3y over the rectangular region R = [1, 3] × [0, 2].

Evaluate ∬_R (2x + 3y) dA over R = [1, 3] × [0, 2]
1
Step 1 — Set Up the Iterated IntegralSince R is a rectangle, we can use either order. Let's choose dy dx. The integral becomes ∫₁³ ∫₀² (2x + 3y) dy dx. The inner integral integrates y from 0 to 2, and the outer integral integrates x from 1 to 3.
∫₁³ ∫₀² (2x + 3y) dy dx
2
Step 2 — Evaluate the Inner Integral (with respect to y)Treat x as a constant and integrate with respect to y. The antiderivative of 2x with respect to y is 2xy, and the antiderivative of 3y is (3/2)y². Evaluate from y = 0 to y = 2: [2x(2) + (3/2)(2²)] − [2x(0) + (3/2)(0²)] = 4x + 6 − 0 = 4x + 6.
Inner integral = 4x + 6
3
Step 3 — Evaluate the Outer Integral (with respect to x)Now integrate 4x + 6 with respect to x from 1 to 3. The antiderivative is 2x² + 6x. Evaluate: [2(3²) + 6(3)] − [2(1²) + 6(1)] = [18 + 18] − [2 + 6] = 36 − 8 = 28.
∬_R (2x + 3y) dA = 28
4
Step 4 — Verify by Switching OrderFor a rectangle, we can switch to dx dy: ∫₀² ∫₁³ (2x + 3y) dx dy. Inner integral: integrate with respect to x → [x² + 3xy] from x = 1 to x = 3 gives (9 + 9y) − (1 + 3y) = 8 + 6y. Outer integral: ∫₀² (8 + 6y) dy = [8y + 3y²] from 0 to 2 = 16 + 12 = 28. ✓ Same answer!
Both orders give 28 — Fubini's Theorem confirmed

Non-Rectangular Region Example

Evaluate ∫₀¹ ∫₀ˣ (x + y) dy dx
1
Step 1 — Identify the RegionThis is a Type I setup. The outer variable x goes from 0 to 1. For each fixed x, y goes from 0 to x. This describes the triangular region below the line y = x and above the x-axis, from x = 0 to x = 1.
Region: triangle with vertices (0,0), (1,0), (1,1)
2
Step 2 — Inner Integral (with respect to y)Integrate (x + y) with respect to y from 0 to x. Antiderivative: xy + y²/2. Evaluate: [x·x + x²/2] − [0] = x² + x²/2 = 3x²/2.
Inner integral = 3x²/2
3
Step 3 — Outer Integral (with respect to x)Integrate 3x²/2 from 0 to 1. Antiderivative: (3/2)(x³/3) = x³/2. Evaluate: 1³/2 − 0 = 1/2.
∫₀¹ ∫₀ˣ (x + y) dy dx = 1/2

Common Strategies & Pitfalls

Even though the idea is simple — integrate twice — students often run into trouble with setup errors and algebraic mistakes. The table below highlights the most common issues alongside strategies to avoid them.

Strategies and common pitfalls when working with iterated integrals
Strategy / StrengthCommon Pitfall
Sketch the region before writing limits. A quick drawing reveals whether to use Type I or Type II.Setting up limits without a sketch, leading to incorrect bounds that don't match the actual region.
The inner limits can be functions of the outer variable (e.g., y goes from 0 to x).Putting functions of the inner variable in the outer limits — outer limits must always be constants.
Treat the outer variable as a constant during the inner integration, just like a parameter.Accidentally integrating the outer variable during the inner step (e.g., integrating x when you should treat it as constant).
Switching the order of integration can simplify a difficult integral dramatically.Switching the order without correctly re-deriving the new limits from a sketch of the region.
Check your answer by trying both orders on a rectangular region.Assuming both orders always have the same limits — they only do for rectangles.
KEY TAKEAWAY
Think of setting up the limits like giving GPS directions. The outer integral says "drive from Street A to Street B" (constant limits). The inner integral says "at each street, walk from one side of the block to the other" — and the block width might change from street to street (variable limits). If you mix up which directions are constant and which change, you'll end up in the wrong neighborhood.

Connection to Advanced Topics

Iterated integrals are not just a technique — they are a gateway to a wide range of advanced calculus ideas. Understanding how to set up and evaluate them prepares you for triple integrals, alternative coordinate systems, and real-world applications in physics and engineering.

How iterated integrals connect to more advanced topics
This LessonAdvanced Extension
Double integrals over rectangular regions (dx dy or dy dx)Triple integrals (∫∫∫ f(x,y,z) dV) for computing volumes and masses of 3-D objects
Cartesian coordinates (x, y)Polar coordinates (r, θ) with dA = r dr dθ — ideal for circular regions
Switching order of integration (Fubini's Theorem)Change of variables with Jacobians — generalizes coordinate transformations
Volume under z = f(x, y) when f ≥ 0Surface area, center of mass, moments of inertia, probability distributions

In physics, for example, double integrals let you compute the total charge on a plate with non-uniform charge density, or the gravitational force exerted by a flat mass distribution. The iterated integral is the practical tool that makes these computations feasible. Mastering it now gives you a solid foundation for all of these extensions.

Practice Problems

PROBLEM 1CONCEPTUAL
In the iterated integral ∫₂⁵ ∫₁⁴ f(x, y) dy dx, which variable is integrated first, and what are its limits? Which variable is integrated second?
PROBLEM 2BASIC CALCULATION
Evaluate the iterated integral ∫₀² ∫₀³ (4xy) dy dx.
PROBLEM 3INTERMEDIATE
Set up and evaluate the double integral of f(x, y) = x² over the triangular region with vertices at (0, 0), (2, 0), and (2, 2). Use the dy dx order.
PROBLEM 4APPLIED
A metal plate occupies the rectangular region [0, 3] × [0, 2] (in meters). Its density at point (x, y) is ρ(x, y) = 1 + x + y kg/m². Set up and evaluate an iterated integral to find the total mass of the plate.
PROBLEM 5CRITICAL THINKING
Consider the iterated integral ∫₀¹ ∫ₓ¹ eʸ² dy dx. Explain why this integral is difficult to evaluate in its current order, then switch the order of integration and evaluate.

Lesson Summary

An iterated integral converts a double integral into two successive single-variable integrals. You evaluate the inner integral first (treating the other variable as a constant), then feed the result into the outer integral. Fubini's Theorem guarantees that both orders of integration (dy dx and dx dy) yield the same result when the integrand is continuous.

For rectangular regions, all limits are constants and either order works equally well. For non-rectangular regions, the inner limits become functions of the outer variable. Always sketch the region, classify it as Type I or Type II, and choose the order of integration that makes the calculation simplest. Sometimes switching the order is the key to evaluating an otherwise impossible integral.

Varsity Tutors • Multivariable Calculus • Iterated Integrals — Set up and evaluate iterated integrals