Multivariable Calculus · Undergraduate

Double Integrals over Rectangles

Quick answer

The double integral of f(x, y) over a rectangle adds up f times area, the way a single integral adds up f times length. Cutting the rectangle into small pieces and standing a box on each gives a Riemann sum, and the limit is the volume under the surface. Fubini's theorem turns the double integral into two single integrals, one inside the other, taken in either order. Dividing by the area gives the average value of f.

What you'll learn

  • Estimate a double integral with a Riemann sum of boxes
  • Interpret a double integral as a signed volume
  • Evaluate a double integral as an iterated integral in either order
  • Find the average value of a function over a rectangle

Volume under a surface

A single integral ∫abf(x) dx\int_a^b f(x)\,dx gives the area under a curve. One dimension up, a function f(x,y)≥0f(x, y) \ge 0 defined on a rectangle R=[a,b]×[c,d]R = [a, b] \times [c, d] has a graph z=f(x,y)z = f(x, y), a surface. The solid between RR and the surface has a volume, and the double integral ∬Rf(x,y) dA\iint_R f(x, y)\,dA is that volume.

The solid under z = 16 − x² − 2y² over a square The surface z = 16 − x² − 2y², drawn above the square from 0 to 2 in both x and y. It is highest, at 16, above the origin corner and lowest, at 4, above the corner (2, 2). The solid between the square and the surface has volume 48. 2 4 6 1 8 1 10 12 2 14 16 2 x y z z = 16 − x² − 2y²
The solid under z = 16 − x² − 2y² over a square

Riemann sums with boxes

To estimate the volume, cut RR into small rectangles, pick a sample point in each, and stand a box on each rectangle, as tall as ff at the sample point. The volume of a box is f(xi∗,yj∗) ΔAf(x_i^*, y_j^*)\,\Delta A, so the solid’s volume is about

∑i∑jf(xi∗,yj∗) ΔA\sum_{i}\sum_{j} f(x_i^*, y_j^*)\,\Delta A

a double Riemann sum. As the rectangles shrink, the sums approach the double integral:

∬Rf(x,y) dA=lim⁡m,n→∞∑i=1m∑j=1nf(xi∗,yj∗) ΔA\iint_R f(x, y)\,dA = \lim_{m, n \to \infty} \sum_{i=1}^{m}\sum_{j=1}^{n} f(x_i^*, y_j^*)\,\Delta A

The limit exists for every continuous ff. When ff takes negative values, boxes below the xyxy-plane count as negative, and the integral is a signed volume.

Four boxes estimating the volume The square from 0 to 2 in x and y, cut into four unit squares. On each stands a box as tall as 16 − x² − 2y² at the square's upper right corner: 13 on the square nearest the origin, 10 and 7 on the two side squares, and 4 on the far square. Their total volume, 34, estimates the volume under the surface. 2 4 6 1 8 1 10 12 2 14 16 2 x y z
Four boxes estimating the volume

Why two single integrals are enough

Slice the solid with a plane x=constantx = \text{constant}. The slice is a region under the curve z=f(x,y)z = f(x, y) for c≤y≤dc \le y \le d, and its area is a single integral in yy, with xx held fixed:

A(x)=∫cdf(x,y) dyA(x) = \int_c^d f(x, y)\,dy

A thin slab of thickness Δx\Delta x has volume about A(x) ΔxA(x)\,\Delta x, and adding the slabs is the slicing idea behind the disk and washer methods: the volume is ∫abA(x) dx\int_a^b A(x)\,dx. Slicing the other way, with planes y=constanty = \text{constant}, gives the same volume. Fubini’s theorem says this holds for every continuous ff on a rectangle:

∬Rf(x,y) dA=∫ab∫cdf(x,y) dy dx=∫cd∫abf(x,y) dx dy\iint_R f(x, y)\,dA = \int_a^b \int_c^d f(x, y)\,dy\,dx = \int_c^d \int_a^b f(x, y)\,dx\,dy

The right-hand sides are iterated integrals. Work from the inside out: the inner integral treats the outer variable as a constant, the same way a partial derivative does. The order of dy dxdy\,dx tells which variable goes first, and each pair of limits belongs to its own variable.

Worked examples

Average value

The average value of a function of one variable is its integral divided by the length of the interval. With two variables, divide by area:

favg=1area(R)∬Rf(x,y) dAf_{\text{avg}} = \frac{1}{\text{area}(R)}\iint_R f(x, y)\,dA

A box on RR as tall as favgf_{\text{avg}} has the same volume as the solid under the surface.

Common mistakes

Practice problems

  1. Estimate ∬R(x+y) dA\iint_R (x + y)\,dA over R=[0,2]×[0,2]R = [0, 2] \times [0, 2] with four unit squares and midpoints. Then find the exact value.

    Answer

    Estimate 88; exact value 88

    Full solution

    At the midpoints (0.5,0.5)(0.5, 0.5), (1.5,0.5)(1.5, 0.5), (0.5,1.5)(0.5, 1.5) and (1.5,1.5)(1.5, 1.5), x+yx + y equals 11, 22, 22 and 33, for a total of 88. Exactly, ∫02(2x+2) dx=4+4=8\int_0^2 (2x + 2)\,dx = 4 + 4 = 8. The midpoint rule is exact for linear functions: on each square, the surface rises above the box on one side as much as it dips below on the other.

  2. Evaluate ∫01∫02(x+y) dy dx\displaystyle\int_0^1 \int_0^2 (x + y)\,dy\,dx.

    Answer

    33

    Full solution

    The inner integral is [xy+y22]02=2x+2\Big[xy + \tfrac{y^2}{2}\Big]_0^2 = 2x + 2, and ∫01(2x+2) dx=1+2=3\int_0^1 (2x + 2)\,dx = 1 + 2 = 3.

  3. Evaluate ∫03∫016xy2 dy dx\displaystyle\int_0^3 \int_0^1 6xy^2\,dy\,dx.

    Answer

    99

    Full solution

    The inner integral is 6x⋅13=2x6x \cdot \tfrac{1}{3} = 2x, and ∫032x dx=9\int_0^3 2x\,dx = 9.

  4. Evaluate ∬Rex+y dA\iint_R e^{x + y}\,dA over R=[0,1]×[0,1]R = [0, 1] \times [0, 1]. Give the exact value and a decimal.

    Hint

    ex+y=exeye^{x + y} = e^x e^y.

    Answer

    (e−1)2≈2.9525(e - 1)^2 \approx 2.9525

    Full solution

    The integrand is a product, so the integral is (∫01ex dx)(∫01ey dy)=(e−1)(e−1)\left(\int_0^1 e^x\,dx\right)\left(\int_0^1 e^y\,dy\right) = (e - 1)(e - 1).

  5. Find the volume of the solid under the plane z=4−x−yz = 4 - x - y and over R=[0,1]×[0,2]R = [0, 1] \times [0, 2].

    Answer

    55

    Full solution

    The plane stays above z=1z = 1 on RR, so the integral is a true volume. The inner integral is ∫02(4−x−y) dy=8−2x−2=6−2x\int_0^2 (4 - x - y)\,dy = 8 - 2x - 2 = 6 - 2x, and ∫01(6−2x) dx=6−1=5\int_0^1 (6 - 2x)\,dx = 6 - 1 = 5.

  6. Find the average value of f(x,y)=xyf(x, y) = xy on R=[0,2]×[0,4]R = [0, 2] \times [0, 4].

    Answer

    22

    Full solution

    The integral is (∫02x dx)(∫04y dy)=2⋅8=16\left(\int_0^2 x\,dx\right)\left(\int_0^4 y\,dy\right) = 2 \cdot 8 = 16, and the area is 88, so the average is 22.

  7. Evaluate ∫12∫01(2x+y) dx dy\displaystyle\int_1^2 \int_0^1 (2x + y)\,dx\,dy.

    Answer

    52\tfrac{5}{2}

    Full solution

    The inner integral, in xx, is [x2+xy]01=1+y\Big[x^2 + xy\Big]_0^1 = 1 + y, and ∫12(1+y) dy=1+32=52\int_1^2 (1 + y)\,dy = 1 + \tfrac{3}{2} = \tfrac{5}{2}.

  8. Without integrating, decide whether ∬R(y−1) dA\iint_R (y - 1)\,dA over R=[0,3]×[0,2]R = [0, 3] \times [0, 2] is positive, negative or zero.

    Answer

    Zero

    Full solution

    On the half of RR with y>1y > 1 the surface lies above the xyxy-plane, and on the half with y<1y < 1 it lies below by the same amounts, mirrored across y=1y = 1. The signed volumes cancel. Directly, ∫02(y−1) dy=0\int_0^2 (y - 1)\,dy = 0.

  9. A student evaluates ∫02∫12x2y dx dy\displaystyle\int_0^2 \int_1^2 x^2 y\,dx\,dy and expects the answer 44 from Example 3. They get 143\tfrac{14}{3}. Who is right?

    Hint

    Which variable runs from 11 to 22 here?

    Answer

    The student computed correctly, but this integral is over a different rectangle, so 44 does not apply.

    Full solution

    With dxdx inside, the limits 11 to 22 belong to xx and 00 to 22 belong to yy: the rectangle is [1,2]×[0,2][1, 2] \times [0, 2]. The value is (∫12x2 dx)(∫02y dy)=73⋅2=143\left(\int_1^2 x^2\,dx\right)\left(\int_0^2 y\,dy\right) = \tfrac{7}{3} \cdot 2 = \tfrac{14}{3}. Switching the order of integration means moving the limits along with their differentials.

Frequently asked questions

What does a double integral measure?

When f(x, y) ≥ 0, the double integral of f over a region is the volume of the solid between the region and the surface z = f(x, y). In general it is a signed volume: parts below the xy-plane count as negative.

What is an iterated integral?

Two single integrals, one inside the other. The inner integral is taken with the other variable held constant; its result is a function of the outer variable, which the outer integral then integrates.

What does Fubini's theorem say?

For a continuous function on a rectangle, the double integral equals the iterated integral in either order: integrate in y first and then x, or in x first and then y.

How do you find the average value of f over a rectangle?

Divide the double integral of f by the area of the rectangle.

When does a double integral split into a product?

When the region is a rectangle and f(x, y) = g(x)h(y). Then the double integral is the integral of g times the integral of h.

What to learn next