Foundations to FrontierPart VI: Calculus IIIBezhig

Module 39: Double integrals, changing order, and polar area elements

8 exercises, each with a complete solution. 1 figure.

On this page
  1. A. Prerequisites and observable goals
  2. B. The central idea: accumulate over an input region
  3. C. Iterated integrals and coordinates
  4. D. Why the polar factor is necessary
  5. E. A region-first workflow
  6. F. Worked examples
  7. G. Independent practice
  8. H. Complete solutions: Module 39
  9. I. Retention and readiness

A. Prerequisites and observable goals#

Before

Single-variable integration, planar regions, polar coordinates, and density.

After

Describe a region with correct variable bounds, evaluate iterated integrals, reverse order when useful, and include the polar area factor.

B. The central idea: accumulate over an input region#

A double integral ∬Df(x,y) dA\iint_Df(x,y)\dd A sums height times tiny area over a planar region. For nonnegative ff, it is the volume under a graph; for density, it is mass; for a general signed function, it is signed accumulation. The area element belongs to the input plane. It is not the curved surface area of the graph, which is a later topic.

C. Iterated integrals and coordinates#

For a continuous function on a rectangle, Fubini’s theorem permits either iterated order. On a vertically simple region D={(x,y):a≤x≤b, g1(x)≤y≤g2(x)}D=\{(x,y):a\le x\le b,\ g_1(x)\le y\le g_2(x)\},

∬Df dA=∫ab∫g1(x)g2(x)f(x,y) dy dx.\iint_D f\dd A=\int_a^b\int_{g_1(x)}^{g_2(x)}f(x,y)\dd y\dd x.

The inner variable changes first; the other is held constant. Horizontally simple regions use xx bounds depending on yy. A region may need splitting when changing order. For more general improper integrals, absolute integrability or other applicable theorems must be checked; swapping an order across singularities is not universally valid. [O3Double]

In polar coordinates, x=rcos⁡θx=r\cos\theta, y=rsin⁡θy=r\sin\theta and

 dA=r dr dθ,r≥0.\dd A=r\dd r\dd\theta,\qquad r\ge0.

Thus ∬Df dA=∬Ef(rcos⁡θ,rsin⁡θ)r dr dθ\iint_Df\dd A=\iint_E f(r\cos\theta,r\sin\theta)r\dd r\dd\theta on a suitable once-covered region. The factor rr changes area, not the function value. A nonnegative density ρ\rho gives mass m=∬Dρ dAm=\iint_D\rho\dd A, with centroid coordinates xˉ=∬xρ/m\bar x=\iint x\rho/m, yˉ=∬yρ/m\bar y=\iint y\rho/m. The average of ff over positive finite area A(D)A(D) is (1/A(D))∬Df dA(1/A(D))\iint_Df\dd A.

D. Why the polar factor is necessary#

A small polar cell has radial thickness Δr\Delta r and arc-length width approximately rΔθr\Delta\theta, so its area is approximately rΔrΔθr\Delta r\Delta\theta. Equal changes in angle sweep larger cells at larger radii. Ignoring rr would weight the disk’s center and edge incorrectly. The Jacobian module derives the factor algebraically.

Changing integration order does not mean mechanically swapping differentials and leaving the bounds unchanged. It means describing the same set of input points by a different family of slices. The region is the invariant object; the inequalities must be rebuilt from its geometry.

E. A region-first workflow#

Draw or describe the region with inequalities. Decide which slice direction makes the integrand and bounds easier. Integrate the inner variable while treating the outer one as constant, then integrate the resulting expression over its full outer interval. For order reversal, project the region onto the new outer axis and solve its boundary equations for the new inner variable. For a circular or radial region, test polar coordinates and explicitly include the area factor before integrating.

F. Worked examples#

Reverse order to make an integral elementary. Consider ∫01∫x1ey2 dy dx\int_0^1\int_x^1e^{y^2}\dd y\dd x. The region is 0≤x≤y≤10\le x\le y\le1. Reversing gives

∫01∫0yey2 dx dy=∫01yey2 dy=e−12.\int_0^1\int_0^y e^{y^2}\dd x\dd y =\int_0^1y e^{y^2}\dd y=\frac{e-1}{2}.

The original inner antiderivative is not elementary, but integration across horizontal widths introduces the factor needed for substitution.

Exact triangle 0≤x≤y≤10\le x\le y\le1. A vertical slice runs from y=xy=x to y=1y=1; a horizontal slice runs from x=0x=0 to x=yx=y. Bounds change with the order, while the region does not.

A disk integral. Over x2+y2≤4x^2+y^2\le4, the integrand x2+y2x^2+y^2 becomes r2r^2. Therefore ∬D(x2+y2) dA=∫02π∫02r3 dr dθ=8π\iint_D(x^2+y^2)\dd A=\int_0^{2\pi}\int_0^2r^3\dd r\dd\theta=8\pi. The r3r^3 combines the transformed height r2r^2 with the area factor rr.

G. Independent practice#

F foundation C standard course work T transfer S optional stretch L runnable lab or project. Attempt the complete set before opening its solutions.

M39.01FSolution

Evaluate ∫02∫03(x+2y) dy dx\int_0^2\int_0^3(x+2y)\dd y\dd x.

M39.02CSolution

Evaluate ∬Dxy dA\iint_Dxy\dd A where D={x≥0,y≥0,x+y≤2}D=\{x\ge0,y\ge0,x+y\le2\}.

M39.03TSolution

Reverse the order and evaluate ∫01∫x1yey3 dy dx\int_0^1\int_x^1 y e^{y^3}\dd y\dd x. Explain why the new order makes substitution useful.

M39.04CSolution

Evaluate ∬x2+y2≤1(x2+y2) dA\iint_{x^2+y^2\le1}(x^2+y^2)\dd A using polar coordinates.

M39.05TSolution

Find the area of the annular sector 1≤r≤31\le r\le3, 0≤θ≤π/30\le\theta\le\pi/3, using a double integral.

M39.06TSolution

On the triangle x≥0,y≥0,x+y≤1x\ge0,y\ge0,x+y\le1, density is ρ=1+x\rho=1+x. Find mass and both centroid coordinates.

M39.07CSolution

Find the average value of f=x2+y2f=x^2+y^2 on [0,2]×[0,1][0,2]\times[0,1].

M39.08SSolution

The region is described vertically by 0≤x≤20\le x\le2, 0≤y≤min⁡(x,2−x)0\le y\le\min(x,2-x). Give an iterated integral for its area in each order and explain why one description needs a split.

Optional small hints

For M39.02, a vertical slice has 0≤y≤2−x0\le y\le2-x. For M39.06, moments require multiplying by the coordinate as well as density. For M39.08, the roof changes formula at x=1x=1.

H. Complete solutions: Module 39#

Compare the first step where your work differs. A corrected answer is not yet an independent success.

M39.01Show solution

Hold the outer variable constant in the inner integral

Integrating in yy gives [xy+y2]03=3x+9[xy+y^2]_0^3=3x+9. Then ∫02(3x+9) dx=[3x2/2+9x]02=24\int_0^2(3x+9)\dd x=[3x^2/2+9x]_0^2=\boxed{24}. The rectangular area is 66; an average height of 44 is consistent with the affine integrand’s center value.

Back to M39.01
M39.02Show solution

Translate a triangular inequality into bounds

Use 0≤x≤20\le x\le2 and 0≤y≤2−x0\le y\le2-x. Then

I=∫02x(2−x)22 dx=12∫02(4x−4x2+x3) dx=12(8−32/3+4)=2/3.I=\int_0^2\frac{x(2-x)^2}{2}\dd x =\frac12\int_0^2(4x-4x^2+x^3)\dd x =\frac12(8-32/3+4)=\boxed{2/3}.

The integrand is nonnegative, so the final sign fits the geometry.

Back to M39.02
M39.03Show solution

Use slice width to create a substitution factor

The original inequalities are 0≤x≤10\le x\le1 and x≤y≤1x\le y\le1. Equivalently 0≤y≤10\le y\le1 and 0≤x≤y0\le x\le y. At fixed yy, the integrand is constant in xx, so

I=∫01∫0yyey3 dx dy=∫01y2ey3 dy.I=\int_0^1\int_0^y y e^{y^3}\dd x\dd y =\int_0^1 y^2e^{y^3}\dd y.

Now u=y3u=y^3 gives  du=3y2 dy\dd u=3y^2\dd y, with bounds 00 and 11. Hence I=13∫01eu du=(e−1)/3I=\tfrac13\int_0^1e^u\dd u=\boxed{(e-1)/3}. The horizontal slice contributes an additional factor yy, exactly the factor that makes the derivative of y3y^3 available. Merely changing the differentials without changing the bounds would not describe this region.

Back to M39.03
M39.04Show solution

Transform both the integrand and area element

The unit disk corresponds to 0≤r≤10\le r\le1, 0≤θ≤2π0\le\theta\le2\pi. Thus I=∫02π∫01r2 r dr dθ=2π(1/4)=π/2I=\int_0^{2\pi}\int_0^1r^2\,r\dd r\dd\theta=2\pi(1/4)=\boxed{\pi/2}. Omitting the area factor would produce a different and unjustified radial weighting.

Back to M39.04
M39.05Show solution

Area is an integral of one with the correct area element

The area is ∫0π/3∫13r dr dθ=(π/3)(9−1)/2=4π/3\int_0^{\pi/3}\int_1^3r\dd r\dd\theta=(\pi/3)(9-1)/2=\boxed{4\pi/3}. Independently, this is one-sixth of an annulus of area π(9−1)=8π\pi(9-1)=8\pi, giving the same value.

Back to M39.05
M39.06Show solution

Density-weighted area and moments

Using 0≤x≤10\le x\le1, 0≤y≤1−x0\le y\le1-x, m=∫01(1+x)(1−x) dx=2/3m=\int_0^1(1+x)(1-x)\dd x=2/3. The numerator for xˉ\bar x is ∫01x(1−x2) dx=1/4\int_0^1x(1-x^2)\dd x=1/4. The numerator for yˉ\bar y is 12∫01(1+x)(1−x)2 dx=5/24\tfrac12\int_0^1(1+x)(1-x)^2\dd x=5/24. Therefore xˉ=3/8, yˉ=5/16\boxed{\bar x=3/8,\ \bar y=5/16} and m=2/3\boxed{m=2/3} in consistent mass units. Relative to the uniform triangle’s centroid, mass shifts toward larger xx.

Back to M39.06
M39.07Show solution

Divide the accumulation by input area

The double integral is ∫02[x2y+y3/3]01 dx=∫02(x2+1/3) dx=10/3\int_0^2[x^2y+y^3/3]_0^1\dd x=\int_0^2(x^2+1/3)\dd x=10/3. The rectangle’s area is 22, so the average is 5/3\boxed{5/3}. The average of a nonconstant quadratic is not necessarily its value at the rectangle’s center.

Back to M39.07
M39.08Show solution

A piecewise roof becomes one horizontal description

Vertical slices require A=∫01∫0x1 dy dx+∫12∫02−x1 dy dxA=\int_0^1\int_0^x1\dd y\dd x+\int_1^2\int_0^{2-x}1\dd y\dd x. Horizontal slices use 0≤y≤10\le y\le1 and y≤x≤2−yy\le x\le2-y, so A=∫01∫y2−y1 dx dy=∫01(2−2y) dy=1A=\int_0^1\int_y^{2-y}1\dd x\dd y=\int_0^1(2-2y)\dd y=\boxed1. The split arises because the upper vertical boundary changes from one line to the other at x=1x=1, not because Fubini’s theorem changes the region.

Back to M39.08

I. Retention and readiness#

Retrieve M17.05 and M32.08. Explain why the inner bounds may depend on the outer variable but not on the variable currently being integrated after its endpoint evaluation. Advance when you can draw the region implied by someone else’s bounds and detect an omitted polar factor.