Double Integrals: Rectangular & General Regions — Double Integrals Over Rectangular and General Regions - Calculus 3

SECTION 1

Historical Context & Motivation

The concept of integration in multiple variables arose naturally from problems in physics and geometry that single-variable calculus could not address. Calculating the mass of a plate with varying density, the volume beneath a surface, or the gravitational attraction of a lamina all require accumulating infinitesimal contributions across a two-dimensional region. The double integral was developed to formalize exactly this process, extending the Riemann sum paradigm from intervals on the real line to regions in the plane. The mathematical foundations were laid over two centuries of collaborative intellectual effort, beginning with the earliest attempts to compute volumes and areas by exhaustion, and culminating in the rigorous measure-theoretic framework of the twentieth century.

1680s • Newton & Leibniz Develop Calculus

The independent development of single-variable calculus by Isaac Newton and Gottfried Wilhelm Leibniz established the fundamental theorem of calculus, linking differentiation and integration on intervals. Their work set the stage for extending these ideas to higher dimensions.

1760s • Euler & Iterated Integrals

Leonhard Euler systematically computed double integrals by iterating single integrals—integrating first with respect to one variable and then the other. His computational approach became the standard technique and demonstrated the power of reducing multidimensional problems to sequences of one-dimensional ones.

1811 • Fourier's Analytical Theory of Heat

Joseph Fourier's work on heat conduction required integrating functions over two-dimensional domains, motivating rigorous treatment of double integrals in mathematical physics. His applications showcased the indispensability of multiple integration in modeling continuous phenomena.

1854 • Riemann's Integration Theory

Bernhard Riemann formalized the definition of the integral using partitions and Riemann sums. His framework extended naturally to functions of two variables by partitioning rectangular regions into subrectangles, providing the rigorous underpinning for double integrals.

1902 • Lebesgue Measure Theory

Henri Lebesgue introduced measure theory, generalizing integration to far more complex regions and functions. His framework unified double integrals over rectangular and non-rectangular domains under a single theoretical umbrella, resolving subtle convergence and integrability questions.

The central question this lesson addresses is both geometric and analytical: given a function f(x, y) defined over a region R in the xy-plane, how do we compute the total accumulated value of f across R? When R is a rectangle, the computation is relatively straightforward via iterated integrals. When R has curved or piecewise boundaries—a general region—we must carefully describe those boundaries as functions and adjust our limits of integration accordingly.

SECTION 2

Core Principles & Definitions

Understanding double integrals begins with generalizing the Riemann sum construction. In single-variable calculus, we partition an interval into subintervals, evaluate the function at sample points, multiply by subinterval width, and sum. For double integrals, we partition a planar region into small pieces, evaluate f(x, y) at sample points within each piece, multiply by the piece's area, and sum. The limit of this process, as the partition becomes infinitely fine, yields the double integral. The following core ideas form the foundation of this extension.

Riemann Sums in 2D

Partition the region R into n small subregions R_k with areas ΔA_k. Choose sample points (x_k^*, y_k^*) in each R_k. The Riemann sum Σ f(x_k^*, y_k^*) ΔA_k approximates the double integral, converging to the exact value as the mesh norm tends to zero.

Rectangular Regions

A rectangular region is defined by R = [a, b] × [c, d], where both variables have constant bounds. Iterated integrals over rectangles are computed by integrating with respect to one variable while treating the other as a constant, then integrating the result.

Fubini's Theorem

If f(x, y) is continuous on a rectangle R = [a, b] × [c, d], then the double integral equals either iterated integral: ∬_R f dA = ∫ₐᵇ ∫꜀ᵈ f(x,y) dy dx = ∫꜀ᵈ ∫ₐᵇ f(x,y) dx dy. The order of integration may be interchanged freely.

General (Non-Rectangular) Regions

A Type I region has variable y-bounds: a ≤ x ≤ b, g₁(x) ≤ y ≤ g₂(x). A Type II region has variable x-bounds: c ≤ y ≤ d, h₁(y) ≤ x ≤ h₂(y). These classifications determine how to set up the iterated integral's limits.

Geometric Interpretation

When f(x, y) ≥ 0, the double integral ∬_R f(x, y) dA represents the volume of the solid bounded above by z = f(x, y) and below by the xy-plane over the region R. When f = 1, the double integral simply gives the area of R.

✦

KEY TAKEAWAY

Think of a double integral like measuring total rainfall over a city. A rain gauge at each point records the rate of rainfall f(x, y), and you want the total volume of water collected over the entire area. You mentally slice the city into thin north–south strips (integrating in one direction), compute the rainfall along each strip (the inner integral), and then sum all strips together (the outer integral). For a rectangular park, each strip spans the same width; for an irregularly shaped district, the strip endpoints vary from position to position—that is the difference between rectangular and general regions.

SECTION 3

Visualizing Double Integrals Over Rectangles

The following diagram illustrates the fundamental construction of a double integral over a rectangular region. The rectangle R = [a, b] × [c, d] in the xy-plane is partitioned into a grid of subrectangles. Over each subrectangle, a rectangular prism (box) is erected whose height equals the function value f(x_k^*, y_k^*) at a chosen sample point. The sum of all box volumes approximates the double integral; the approximation improves as the grid becomes finer. In the limit, the sum converges to the exact volume beneath the surface z = f(x, y).

Each colored box represents a term f(x_k^*, y_k^*) × ΔA in the Riemann sum. The dashed green curve represents the surface z = f(x, y). As the partition is refined, the box tops approach the surface, and their total volume converges to ∬_R f(x, y) dA.

The diagram above captures the essential geometric meaning of the double integral. When we write ∬_R f(x, y) dA, we are computing the signed volume between the surface and the xy-plane. Over a rectangle, the partition is a simple grid whose subrectangles each have area ΔA = Δx × Δy. The iterated integral structure—integrating first in y then in x, or vice versa—corresponds to summing along columns (or rows) of this grid and then accumulating those column (or row) sums.

SECTION 4

Mathematical Framework

We now establish the formal definitions and key theorems that govern double integrals over both rectangular and general regions. The theoretical backbone is Fubini's Theorem, which allows us to evaluate a double integral as two successive single integrals, and the characterization of general regions as Type I or Type II domains.

DOUBLE INTEGRAL DEFINITION

∬_R f(x, y) dA = lim_{‖P‖→0} Σₖ f(xₖ*, yₖ*) ΔAₖ

P denotes a partition of R into subrectangles, ‖P‖ is the mesh norm (diameter of the largest subrectangle), (xₖ*, yₖ*) is a sample point in the k-th subrectangle, and ΔAₖ is its area. If f is continuous (or bounded with discontinuities on a set of measure zero), this limit exists and is independent of the sample-point choice.

FUBINI'S THEOREM — RECTANGULAR REGIONS

∬_R f(x, y) dA = ∫_a^b ∫_c^d f(x, y) dy dx = ∫_c^d ∫_a^b f(x, y) dx dy

Here R = [a, b] × [c, d] is a closed rectangle and f is continuous on R. The theorem guarantees that both orders of iteration yield the same result. When one order is computationally simpler than the other, Fubini's Theorem permits choosing the more convenient order.

TYPE I REGION (ITERATED INTEGRAL)

∬_D f(x, y) dA = ∫_a^b ∫_{g₁(x)}^{g₂(x)} f(x, y) dy dx

D is a Type I region: D = {(x, y) | a ≤ x ≤ b, g₁(x) ≤ y ≤ g₂(x)}, where g₁ and g₂ are continuous functions. The inner integral's limits depend on x, reflecting the curved boundaries of D.

TYPE II REGION (ITERATED INTEGRAL)

∬_D f(x, y) dA = ∫_c^d ∫_{h₁(y)}^{h₂(y)} f(x, y) dx dy

D is a Type II region: D = {(x, y) | c ≤ y ≤ d, h₁(y) ≤ x ≤ h₂(y)}, where h₁ and h₂ are continuous functions. The inner integral runs over x with limits that depend on y. Many regions can be described as both Type I and Type II, and the choice affects the difficulty of the computation.

💡

When to Switch the Order of Integration

If the inner integral in one order leads to an antiderivative that is difficult or impossible to express in closed form, try reversing the order. Sketch the region D, re-describe it as the other type, and set up new limits. This strategy is especially effective when the integrand contains expressions like

e^{y²}

or sin(x²) that lack elementary antiderivatives in one variable but simplify in the other.

SECTION 5

Type I vs. Type II Regions — A Detailed Breakdown

The classification of a region as Type I or Type II dictates the structure of the iterated integral. A Type I region is bounded on the left and right by vertical lines x = a and x = b, while the bottom and top boundaries are curves y = g₁(x) and y = g₂(x). For each fixed x in [a, b], the variable y ranges between these two curves. A Type II region is bounded above and below by horizontal lines y = c and y = d, with left and right boundaries given by curves x = h₁(y) and x = h₂(y). The diagram below contrasts these two perspectives on the same region.

Left: a Type I region with vertical slices — for each fixed x, y ranges from g₁(x) to g₂(x). Right: a Type II region with horizontal slices — for each fixed y, x ranges from h₁(y) to h₂(y). The pink arrows indicate the direction of the inner integral in each case.

Feature	Type I Region	Type II Region
Outer variable	x ranges from a to b (constants)	y ranges from c to d (constants)
Inner variable	y ranges from g₁(x) to g₂(x)	x ranges from h₁(y) to h₂(y)
Slice direction	Vertical slices (parallel to y-axis)	Horizontal slices (parallel to x-axis)
Best used when	Boundaries are naturally expressed as y = f(x)	Boundaries are naturally expressed as x = g(y)
Example region	Region between y = x² and y = x, 0 ≤ x ≤ 1	Region between x = y and x = √y, 0 ≤ y ≤ 1

Comparison of Type I and Type II region descriptions

ℹ

Regions That Are Both Types

Many regions can be expressed as either Type I or Type II. For instance, the region bounded by y = x² and y = 1 can be described as Type I with 0 ≤ x ≤ 1 and x² ≤ y ≤ 1, or as Type II with 0 ≤ y ≤ 1 and −√y ≤ x ≤ √y. Some regions require decomposition into subregions, each of which is Type I or Type II, before the integral can be set up as a sum of iterated integrals.

SECTION 6

Worked Examples

Example 1: Double Integral Over a Rectangle

Evaluate ∬_R (2xy + 3) dA, where R = [1, 3] × [0, 2].

Evaluating a Double Integral Over a Rectangle

1
Step 1 — Identify the Region and Set Up

R is the rectangle [1, 3] × [0, 2], so x ranges from 1 to 3 and y ranges from 0 to 2. Since the region is rectangular, both iteration orders are equally straightforward. We choose dy dx: ∫₁³ ∫₀² (2xy + 3) dy dx.
2
Step 2 — Evaluate the Inner Integral (with respect to y)

Treating x as a constant, integrate (2xy + 3) with respect to y from 0 to 2: ∫₀² (2xy + 3) dy = [xy² + 3y]₀² = x(2)² + 3(2) − 0 = 4x + 6.

Inner integral = 4x + 6
3
Step 3 — Evaluate the Outer Integral (with respect to x)

Now integrate the result from Step 2 with respect to x from 1 to 3: ∫₁³ (4x + 6) dx = [2x² + 6x]₁³ = (2(9) + 18) − (2(1) + 6) = (18 + 18) − (2 + 6) = 36 − 8 = 28.

∬_R (2xy + 3) dA = 28

Example 2: Double Integral Over a General Region (Type I)

Evaluate ∬_D x dA, where D is the region bounded by y = x² and y = x.

Evaluating a Double Integral Over a Type I Region

1
Step 1 — Sketch and Identify the Region

The curves y = x² and y = x intersect where x² = x, giving x(x − 1) = 0, so x = 0 and x = 1. For 0 ≤ x ≤ 1, the parabola y = x² lies below the line y = x. Thus D is a Type I region: 0 ≤ x ≤ 1, x² ≤ y ≤ x.
2
Step 2 — Set Up the Iterated Integral

∬_D x dA = ∫₀¹ ∫_{x²}^{x} x dy dx.
3
Step 3 — Evaluate the Inner Integral

Since x is constant with respect to y: ∫_{x²}^{x} x dy = x · (x − x²) = x² − x³.

Inner integral = x² − x³
4
Step 4 — Evaluate the Outer Integral

∫₀¹ (x² − x³) dx = [x³/3 − x⁴/4]₀¹ = 1/3 − 1/4 = 4/12 − 3/12 = 1/12.

∬_D x dA = 1/12

SECTION 7

Strengths, Limitations & Strategy Selection

Double integrals via iterated integrals are a powerful and general technique, but their practical effectiveness depends on the choice of integration order and the description of the region. Understanding the strengths and limitations of each approach helps you select the optimal strategy for a given problem.

Aspect	Strength	Limitation
Rectangular regions	Constant limits make both orders equally easy; no region-sketching required.	Only applicable when the domain has no curved boundaries — a special case.
Type I setup	Natural when boundaries are given as y = g(x); aligns with vertical cross-sections.	Fails when the region cannot be described by a single pair of y-bounds for each x; may require decomposition.
Type II setup	Natural when boundaries are given as x = h(y); may simplify the integrand.	Requires solving for x as a function of y, which may be algebraically difficult.
Switching order of integration	Can convert an impossible inner integral into a tractable one (e.g., ∫ $e^{y²}$ dy becomes ∫ $e^{y²}$ dx).	Requires re-describing the region in the other type, which demands careful geometric reasoning.
Cartesian coordinates overall	Conceptually straightforward; directly applicable to most polynomial and algebraic boundaries.	Circular and polar regions are awkward; polar coordinates (r, θ) are often superior for circular symmetry.

Strengths and limitations of double integral strategies in Cartesian coordinates

🧭

STRATEGIC INSIGHT

Choosing the right integration order is like choosing the right direction to mow a oddly shaped lawn. If the lawn narrows along the north–south axis but has clean straight edges east to west, it makes sense to mow east–west strips and stack them north to south. In double integrals, the 'natural' direction is the one where the limits of integration are simplest—ideally functions you can antidifferentiate easily. When one order leads to an integrand you cannot antidifferentiate, sketch the region, reverse the order, and try again.

SECTION 8

Connections to Advanced Theory

Double integrals over rectangular and general regions form the foundation upon which several more advanced topics in multivariable calculus and analysis are built. Understanding the Cartesian framework thoroughly prepares you for coordinate transformations, triple integrals, and surface integrals. The table below maps these connections explicitly.

This Lesson	Advanced Extension	Key New Idea
Double integral over rectangle [a,b]×[c,d]	Triple integral over box [a,b]×[c,d]×[e,f]	Add a third iterated integral; dA becomes dV
Cartesian coordinates (x, y)	Polar coordinates (r, θ), dA = r dr dθ	Jacobian factor r arises from the coordinate transformation
General regions described by bounding curves	Change of variables via Jacobians	Any smooth coordinate change T: (u,v) → (x,y) introduces \|∂(x,y)/∂(u,v)\| in the integrand
∬ f dA for volume under a surface	Surface integrals ∬_S f dS	Integration over curved surfaces in 3D rather than flat regions in the plane
Fubini's Theorem (switching order)	Green's Theorem, Stokes' Theorem	Relate double integrals over regions to line integrals along boundaries

How double integrals connect to advanced topics in multivariable calculus

Perhaps the most immediately relevant extension is the transition to polar coordinates. When the region D is a disk, annulus, or sector, expressing it in Cartesian coordinates leads to square-root expressions in the limits that complicate integration. In polar coordinates, these regions become rectangles in (r, θ)-space, and the area element becomes dA = r dr dθ, incorporating the Jacobian of the coordinate transformation. The Jacobian concept generalizes far beyond polar coordinates—it is the key tool for any substitution in multiple integrals, analogous to u-substitution in single-variable calculus but involving a determinant of partial derivatives.

ℹ

Looking Ahead

Green's Theorem (covered later in Calculus 3) states that under appropriate smoothness conditions, ∬_D (∂Q/∂x − ∂P/∂y) dA = ∮_C P dx + Q dy. This remarkable result converts a double integral over a region into a line integral around its boundary—directly leveraging the region descriptions (Type I/II) you are learning now.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why Fubini's Theorem allows us to compute a double integral as two successive single integrals. Under what conditions does the theorem apply, and what does it tell us about the order of integration?

Show answer

Fubini's Theorem states that if f(x, y) is continuous (or more generally, integrable and the iterated integrals exist) on a region R, then the double integral ∬_R f dA can be evaluated by integrating first with respect to one variable and then the other. The theorem holds for continuous functions on closed, bounded regions—including both rectangles and general regions described by continuous boundary curves. The crucial consequence is that the order of integration can be reversed: ∫∫ dy dx = ∫∫ dx dy, giving computational flexibility. This works because continuity ensures the function is well-behaved enough that slicing the volume in either direction produces the same total.

Question 2: BASIC CALCULATION

Evaluate ∬_R (x² + y) dA, where R = [0, 2] × [−1, 1].

Show answer

Set up as ∫₀² ∫₋₁¹ (x² + y) dy dx. Inner integral: ∫₋₁¹ (x² + y) dy = [x²y + y²/2]₋₁¹ = (x² + 1/2) − (−x² + 1/2) = 2x². Outer integral: ∫₀² 2x² dx = [2x³/3]₀² = 16/3. The answer is 16/3. Note that the y-term contributes zero because the region is symmetric about y = 0 and y is an odd function on [−1, 1].

Question 3: INTERMEDIATE

Evaluate ∬_D y dA, where D is the region bounded by y = √x and y = x/2. Set up the integral as both Type I and Type II, then evaluate using whichever is simpler.

Show answer

First, find the intersection points: √x = x/2 gives x = x²/4, so x² − 4x = 0 and x = 0 or x = 4. As Type I: D = {(x, y) | 0 ≤ x ≤ 4, x/2 ≤ y ≤ √x}. The integral is ∫₀⁴ ∫_{x/2}^{√x} y dy dx. Inner: [y²/2]_{x/2}^{√x} = x/2 − x²/8. Outer: ∫₀⁴ (x/2 − x²/8) dx = [x²/4 − x³/24]₀⁴ = (4 − 64/24) = (4 − 8/3) = 4/3. As Type II: D = {(x, y) | 0 ≤ y ≤ 2, y² ≤ x ≤ 2y}. Both setups yield the answer 4/3.

Question 4: APPLIED

A thin plate occupies the triangular region with vertices at (0, 0), (2, 0), and (0, 4). Its density at point (x, y) is ρ(x, y) = 1 + x (in units of kg/m²). Find the total mass of the plate by evaluating an appropriate double integral.

Show answer

The triangle's hypotenuse connects (2, 0) and (0, 4), which is the line y = −2x + 4, or equivalently y = 4 − 2x. As a Type I region: 0 ≤ x ≤ 2, 0 ≤ y ≤ 4 − 2x. The mass is M = ∬_D (1 + x) dA = ∫₀² ∫₀^{4−2x} (1 + x) dy dx. Inner integral: (1 + x)(4 − 2x) = 4 + 2x − 2x². Outer integral: ∫₀² (4 + 2x − 2x²) dx = [4x + x² − 2x³/3]₀² = 8 + 4 − 16/3 = 12 − 16/3 = 20/3 kg. The total mass is 20/3 ≈ 6.67 kg.

Question 5: CRITICAL THINKING

The integral ∫₀¹ ∫_y¹ sin(x²) dx dy cannot be evaluated directly because sin(x²) has no elementary antiderivative with respect to x. Reverse the order of integration and evaluate the resulting integral.

Show answer

The original limits describe a Type II region: 0 ≤ y ≤ 1, y ≤ x ≤ 1. Sketch this: it is the triangle with vertices (0,0), (1,0), (1,1) above the line y = x. Rewriting as Type I: 0 ≤ x ≤ 1, 0 ≤ y ≤ x. The reversed integral is ∫₀¹ ∫₀ˣ sin(x²) dy dx. Inner integral: sin(x²) × x (since ∫₀ˣ dy = x). So we compute ∫₀¹ x sin(x²) dx. Let u = x², du = 2x dx: (1/2)∫₀¹ sin(u) du = (1/2)[−cos(u)]₀¹ = (1/2)(−cos 1 + 1) = (1 − cos 1)/2 ≈ 0.2298. This classic problem demonstrates the power of reversing the order of integration to make an otherwise impossible integral tractable.

SUMMARY

Lesson Summary

Summary

A double integral ∬_R f(x, y) dA extends single-variable integration to two dimensions, computing the accumulated value of f over a planar region R. Over a rectangular region R = [a, b] × [c, d], Fubini's Theorem guarantees that the double integral equals either iterated integral, and the order of integration can be freely interchanged. When f(x, y) ≥ 0, the integral represents the volume of the solid beneath the surface z = f(x, y) above R.

For general (non-rectangular) regions, the domain D is classified as Type I (y-bounds depend on x) or Type II (x-bounds depend on y). The choice determines the limits of the inner integral and often affects computational difficulty. When one order leads to an intractable inner integral, reversing the order of integration — by re-describing the region as the other type — is a powerful problem-solving strategy. These techniques form the foundation for polar-coordinate double integrals, triple integrals, change-of-variable formulas via the Jacobian, and the deep theorems of vector calculus including Green's Theorem.