Double Integrals and Fubini's Theorem

The double integral as the Riemann construction with area elements in place of widths, Fubini's theorem as what reduces that two-dimensional limit to two ordinary integrals, and why the freedom to choose an order is a property of the rectangle rather than of integration.

Definition

The double integral. Over a rectangle R = [ a , b ] × [ c , d ] ,

∬ R f d A = lim ∑ f ( x i ∗ , y j ∗ ) Δ A ,

the same Riemann construction as the single-variable integral, with area elements Δ A in place of widths Δ x . The region is chopped into small rectangles, each contributes its area times a sample height, and the limit is taken as the pieces shrink.

Fubini's theorem evaluates that two-dimensional limit as two one-dimensional ones, in either order, when f is continuous on R :

∬ R f d A = ∫ c d ∫ a b f ( x , y ) d x d y = ∫ a b ∫ c d f ( x , y ) d y d x .

The inner integration holds the outer variable fixed, so its symbols pass through as constants and the result is an expression in the outer variable alone.

Over a non-rectangular region the inner limits are functions of the outer variable. Reversing the order then means redescribing the region, not exchanging the written limits.

Assumptions and scope

  • Fubini's theorem requires continuity on the rectangle, or absolute integrability more generally. Without it the two iterated integrals can differ.

  • Over a non-rectangular region the inner limits depend on the outer variable, and reversing the order requires redescribing the region rather than swapping the limits.

Worked material

Example

Integrands and regions worth recognising

Four cases, chosen because each one settles a different question about what the order and the region cost.

---

1. A separable integrand, where the double integral is a product. Over R = [ 0 , 2 ] × [ 0 , 3 ] with f ( x , y ) = x y 2 , the inner integration carries the x through as a constant:

∫ 0 3 x y 2 d y = x [ y 3 3 ] 0 3 = 9 x , ∫ 0 2 9 x d x = [ 9 x 2 2 ] 0 2 = 18 .

But because the integrand factors as g ( x ) h ( y ) and the limits are constants, the whole thing factors too:

∬ R x y 2 d A = ( ∫ 0 2 x d x ) ( ∫ 0 3 y 2 d y ) = ( 2 ) ( 9 ) = 18   ✓

The factorisation needs both conditions. ∫ 0 2 ∫ 0 3 ( x + y 2 ) d y d x = 24 , and the corresponding product of separate integrals is ( 2 ) ( 9 ) = 18 , which is not it: a sum does not factor. And over a non-rectangular region the limits carry the outer variable, so even a product integrand does not separate.

---

2. A region where the order decides how much work it is. Consider ∫ 0 1 ∫ y 1 e x 2 d x d y .

As written, the inner integral ∫ e x 2 d x has no elementary antiderivative, and the calculation stops. Reversing the order rescues it. The region is { ( x , y ) : 0 ≤ y ≤ 1 ,   y ≤ x ≤ 1 } , which is the triangle below the line y = x ; described the other way, x runs over [ 0 , 1 ] and for each x , y runs from 0 to x :

∫ 0 1 ∫ 0 x e x 2 d y d x = ∫ 0 1 x e x 2 d x = [ e x 2 2 ] 0 1 = e − 1 2 ≈ 0.859141 .

The inner integration produced the factor x that the substitution needed. Both orders are valid and give the same number; only one is computable in closed form. Check: a 3000 × 3000 midpoint sum over the triangle returns 0.85914 … .

---

3. A region needing two pieces in one order and one in the other. Let S be bounded by y = x , y = 2 − x and y = 0 : a triangle with vertices ( 0 , 0 ) , ( 1 , 1 ) , ( 2 , 0 ) .

Sweeping in vertical strips, the upper boundary changes at x = 1 , so d y d x needs two integrals:

∫ 0 1 ∫ 0 x 1 d y d x + ∫ 1 2 ∫ 0 2 − x 1 d y d x = 1 2 + 1 2 = 1 .

Sweeping in horizontal strips, each strip at height y runs from x = y to x = 2 − y throughout, so d x d y needs one:

∫ 0 1 ∫ y 2 − y 1 d x d y = ∫ 0 1 ( 2 − 2 y ) d y = [ 2 y − y 2 ] 0 1 = 1   ✓

Both give the triangle's area of 1. The choice of order is not about correctness; it is about how many pieces the region breaks into.

---

4. Where Fubini has nothing to say. Over [ 0 , 1 ] × [ 0 , 1 ] take

f ( x , y ) = x 2 − y 2 ( x 2 + y 2 ) 2 .

The two iterated integrals both exist and come to π / 4 and − π / 4 . They disagree, and the double integral does not exist at all, because ∬ | f | d A diverges near the origin where f is unbounded.

Nothing in the arithmetic of either iterated integral signals this. Each is a perfectly ordinary calculation returning a finite number. Only checking the hypothesis does, and the hypothesis fails because f is not continuous on the closed rectangle and not absolutely integrable over it. This is the case that makes Fubini a theorem rather than a convention.

---

What the four have in common. The region and the integrand together decide what is possible: whether the integral factors, whether either order is computable, how many pieces are needed, and whether the orders agree at all. None of that is visible from the notation ∬ R f d A , which is why the first step is always to describe the region and check the integrand rather than to start integrating.

Non-example

Errors an iterated integral invites

Swapping the limits of a non-rectangular double integral. Over a rectangle the limits are constants and either order works. Over a region where the inner limits depend on the outer variable, exchanging the order requires describing the region the other way round, reversing the written limits produces a different region and a different number.

Treating Fubini as unconditional. Over a rectangle with a continuous integrand the two orders agree, and disagreement means an arithmetic slip. That is a consequence of the hypothesis, not a general fact about integration: without continuity or absolute integrability the two iterated integrals can genuinely differ, and the two-dimensional integral may not exist at all. Quoting "either order works" without the rectangle and the continuity is quoting a conclusion without its premise.

Reading the inner limits as fixed when they are not. Over a region bounded by y = x 2 and y = x , the inner limits in ∫ ∫ d y d x run from x 2 to x , both functions of the outer variable. Substituting them leaves an expression in x alone, which is the point. Writing constant inner limits there integrates over a rectangle the region does not fill.

Related units

Requires

Learn this topic

Used in

Sources

Results update as you type. Use the up and down arrow keys to move between results, Enter to open one, and Escape to close.

Type to search.

Settings

Appearance

Interface density

Your record

Your progress is stored in this browser and nowhere else: an identifier, the answers you have given, the mastery states and review schedule derived from them, and the lesson you last opened. Clearing it makes you a new learner on this device. It cannot be undone, and it will not affect your appearance or density settings.

Focus timer

Focus--minutes remaining

Phase

Kept in this browser only, and used to label the session in your own history.

Today

Nothing recorded yet. Finish a focus session and it will appear here.

Settings

Focus sessions between long breaks.

Sessions you are aiming for in a day.

Notifications

Your history

Sessions are stored in this browser and nowhere else. They are not evidence and never reach your mastery record.