Difficulty: Intermediate | Prerequisites: Single-variable integration, partial derivatives, basic vector operations (dot product, cross product).
This material sits in the second half of a standard Calculus III course, where you move from integrating over simple rectangular regions to integrating over regions whose boundaries are curves or lines. You should already be comfortable setting up and evaluating definite integrals, and you need a working knowledge of partial derivatives. If you can take a cross product of two vectors and sketch a region in the xy-plane from a pair of inequalities, you are ready.
Double integrals let you add up a quantity (population, mass, probability) spread across a two-dimensional region. Surface area extends the idea to measuring how much "sheet" a surface occupies in 3D. Both require careful attention to the shape of the integration region and the correct setup of bounds.
Double integral
The integral ∬_D f(x, y) dA, which sums f over a region D in the plane. Think of it as stacking up thin slices of area, each weighted by the value of f at that point.
Density function ρ(x, y)
A function giving the "amount per unit area" at each point, whether that is population, mass, charge, or probability. In simple terms, it tells you how concentrated something is at each location.
Region of integration (D)
The two-dimensional domain over which you integrate. It can be rectangular, triangular, or bounded by curves. Getting the bounds right is usually the hardest part of the problem.
Iterated integral
A double integral written as two nested single integrals with specific bounds. The inner integral is evaluated first, and its limits may depend on the outer variable.
Surface area (of a graph z = f(x, y))
The area of a surface sitting above a region R in the xy-plane, computed as ∬_R √(1 + (∂z/∂x)² + (∂z/∂y)²) dA. In simple terms, it measures how much the tilting of the surface stretches the flat area beneath it.
Cross product
The vector operation a × b that produces a vector perpendicular to both a and b, whose magnitude equals the area of the parallelogram they span. Used here to find the area element on a surface defined by two edge vectors.
Normal vector to a surface
A vector perpendicular to the surface at a given point. For a plane z = ax + by + c, the normal direction is constant everywhere, which simplifies the surface-area integral.
Splitting a region by sign of x
When the region D is symmetric about the y-axis but the boundary equations change form depending on whether x ≥ 0 or x < 0, split the integral into two pieces. For the triangular region bounded by y = 5x, y = -5x, and |x| ≤ 1:
When x ≥ 0: the bounds on y run from 5x up to 5, with x from 0 to 1.
When x < 0: the bounds on y run from -5x up to 5, with x from -1 to 0.
This is a common pattern for V-shaped or diamond-shaped regions.
Setting up the population integral
Given a density function ρ(x, y) = 10^(3-y), the total population over D is P = ∬_D ρ(x, y) dA. The integral becomes two iterated integrals:
P = ∫₋₁⁰ ∫₋₅ₓ⁵ 10^(3-y) dy dx + ∫₀¹ ∫₅ₓ⁵ 10^(3-y) dy dx
Substitution for exponential integrands
To integrate 10^(3-y) with respect to y, let t = 3 - y, so dt/dy = -1. The antiderivative is -10^t / ln(10) + C. This substitution turns an awkward exponential into a standard form.
After evaluating the inner integral and substituting the bounds, you get expressions involving 10^(3+5x) and 10^(3-5x). These require a second substitution (u = 3 + 5x or v = 3 - 5x) for the outer integral.
Final result
The total population works out to approximately 75 people. The exact expression involves terms with 1/ln²(10), but the numerical answer is what matters on most exams.
The setup
Given the plane 2x + y - z + 4 = 0 (equivalently z = 2x + y + 4), with 1 ≤ y ≤ 4 and 2 ≤ x ≤ 6, find the surface area.
Identifying the corner points
Plug the boundary values of x and y into z = 2x + y + 4 to get the four corners of the surface:
P = (6, 1, 17)
Q = (2, 1, 9)
R = (6, 4, 20)
S = (2, 4, 12)
These define a parallelogram in 3D that sits above the rectangle [2, 6] × [1, 4] in the xy-plane.
Using the cross product to find area
Form two edge vectors from a common vertex (P):
a = PQ = ⟨-4, 0, -8⟩
b = PR = ⟨0, 3, 3⟩
The cross product a × b = ⟨24, 12, -12⟩, and its magnitude |a × b| = √(24² + 12² + 12²) = 12√6.
This gives the area of the parallelogram directly: 12√6.
Using the surface-area integral
Alternatively, apply the surface-area formula. Since z = 2x + y + 4:
∂z/∂x = 2
∂z/∂y = 1
Area = ∬_R √(1 + 4 + 1) dA = √6 · ∬_R 1 dA = √6 · (area of rectangle) = √6 · (4 × 3) = 12√6.
Both methods agree, as they must. The integral method is more general (it works for curved surfaces too), while the cross-product method is quicker for flat surfaces.
Why the answers match
For a plane, the partial derivatives are constants, so √(1 + f_x² + f_y²) is constant and factors out of the integral. The surface-area integral then reduces to that constant times the area of the base region. The cross-product method computes the same quantity from the other direction.
Population / mass from a density function
P = ∬_D ρ(x, y) dA
Antiderivative of 10^t
∫ 10^t dt = 10^t / ln(10) + C
Surface area of z = f(x, y) over a region R
Area = ∬_R √(1 + (∂z/∂x)² + (∂z/∂y)²) dA
Magnitude of a cross product
|a × b| = √(a₁² + a₂² + a₃²), where a × b = ⟨a₂b₃ - a₃b₂, a₃b₁ - a₁b₃, a₁b₂ - a₂b₁⟩
Cross product for area
The area of the parallelogram spanned by vectors a and b equals |a × b|.
Double integrals with density functions are how engineers and scientists compute total quantities from spatially varying data. Estimating the population of a city from census-tract densities, finding the total charge on a semiconductor wafer, or calculating the mass of a plate with non-uniform thickness all use exactly this setup.
Surface area calculations appear in heat transfer (how much surface is exposed to airflow), material science (coating a curved part with a uniform layer), and architecture (estimating cladding for a curved roof).
Students often set up the bounds as though y always runs from a number to another number. When the region is not rectangular, the inner bounds depend on the outer variable. Sketch the region first.
Forgetting the ln(10) denominator when integrating 10^t. The antiderivative of a^t is a^t / ln(a), not just a^t.
Confusing the surface-area formula with the arc-length formula. Surface area has 1 + f_x² + f_y² under the root (two partial derivatives), while arc length has 1 + (dy/dx)² (one ordinary derivative).
Assuming the cross-product method and the integral method are different techniques for different problems. For a plane, they always give the same answer; the integral method generalises to curved surfaces.
⚠️ Setting up bounds for non-rectangular regions is tested heavily. Expect at least one problem where you must sketch the region and determine which variable to integrate first.
⚠️ The surface-area formula √(1 + f_x² + f_y²) is a standard exam question. Know it cold.
⚠️ Substitution inside a double integral (as with the 10^(3-y) integrand) is a common source of lost marks. Practise the chain of substitutions until it feels mechanical.
True or false: the surface area of z = f(x, y) over a flat region R is always greater than or equal to the area of R itself. (True, because √(1 + f_x² + f_y²) ≥ 1.)
Fill in the blank: the antiderivative of 10^t with respect to t is ____. (10^t / ln(10) + C)
True or false: if you swap the order of integration in a double integral, the answer changes. (False, provided you adjust the bounds correctly.)
Fill in the blank: for the plane z = 2x + y + 4, the value of √(1 + (∂z/∂x)² + (∂z/∂y)²) is ____. (√6)
True or false: the magnitude of a × b gives the area of the triangle formed by a and b. (False, it gives the area of the parallelogram. The triangle is half that.)
Q: A region D in the xy-plane is bounded by y = 3x, y = -3x, and y = 6. Set up (but do not evaluate) the double integral of f(x, y) over D using dy dx order.
A: Split by sign of x. For x ≥ 0: ∫₀² ∫₃ₓ⁶ f(x, y) dy dx. For x < 0: ∫₋₂⁰ ∫₋₃ₓ⁶ f(x, y) dy dx.
Q: Find the surface area of the plane z = 3x + 4y + 1 over the rectangle [0, 2] × [0, 5].
A: ∂z/∂x = 3, ∂z/∂y = 4. Area = ∬ √(1 + 9 + 16) dA = √26 · (2 × 5) = 10√26.
Q: Two edge vectors of a flat surface are a = ⟨1, 0, 2⟩ and b = ⟨0, 3, 1⟩. What is the area of the parallelogram they define?
A: a × b = ⟨0·1 - 2·3, 2·0 - 1·1, 1·3 - 0·0⟩ = ⟨-6, -1, 3⟩. |a × b| = √(36 + 1 + 9) = √46.
Q: Why does the integrand of the surface-area formula equal 1 when the surface is the xy-plane itself (z = 0)?
A: Because ∂z/∂x = 0 and ∂z/∂y = 0, so √(1 + 0 + 0) = 1. A flat horizontal surface has the same area as its shadow below, which makes geometric sense.
Double integrals are the foundation for triple integrals (the next major topic), where you add a third variable and integrate over a volume. The surface-area formula reappears when you study surface integrals and flux in vector calculus. The cross-product technique used here for area carries directly into computing normal vectors for surface integrals.
The substitution skills practised in the population-density problem (chained u-substitutions inside an iterated integral) are exactly what you need for cylindrical and spherical coordinate transformations.
double integral, iterated integral, surface area formula, cross product area, population density integral, non-rectangular region, integration bounds, 10^x antiderivative, exponential integrand substitution, parallelogram area, normal vector, MATH 215, Calculus III, multivariable calculus, surface area of a plane