Multiple Integrals

Learning Objectives

  • Interpret a double integral as signed accumulation over a planar region.
  • Evaluate double integrals by iterated integration and Fubini's Theorem.
  • Describe Type I and Type II regions and reverse the order of integration when advantageous.
  • Convert circular and sector-shaped regions to polar coordinates using dA=r dr dθdA=r\,dr\,d\theta.
  • Apply double integrals to area, volume under a surface, mass, first moments, centroids, and planar moments of inertia.
  • Set up introductory surface-area integrals for smooth graphs and parametrized surfaces.
  • Distinguish an ordinary region integral from an iterated integral written with reversed (oriented) limits.

A double integral accumulates a scalar quantity over a two-dimensional region. Depending on the integrand, that accumulation can represent signed volume, area, mass, a first moment, or another engineering quantity. This topic develops the planar geometry and coordinate skills needed before triple integrals and general changes of variables are introduced later.

Double Integrals and Signed Accumulation

Double Integral

For an integrable function ff on a planar region RR, the double integral ∬Rf(x,y) dA\iint_R f(x,y)\,dA is the limit of signed Riemann sums ∑f(xi∗,yj∗) ΔA\sum f(x_i^*,y_j^*)\,\Delta A.

When f(x,y)≥0f(x,y)\ge 0, ∬Rf dA\iint_R f\,dA can be interpreted as the geometric volume between z=f(x,y)z=f(x,y) and the xyxy-plane. When ff changes sign, contributions above the plane are positive and contributions below the plane are negative. The integral is therefore a signed quantity, not the sum of absolute box volumes.

Double Integral

Signed accumulation of f over a planar region R.

∬Rf(x,y) dA\iint_R f(x,y)\,dA

Variables

SymbolDescriptionUnit
RRPlanar region of integration-
f(x,y)f(x,y)Scalar integrand-
dAdADifferential area element-

Interactive Simulation

The simulation below uses midpoint Riemann sums on the square [−2,2]×[−2,2][-2,2]\times[-2,2]. Boxes below the xyxy-plane contribute negatively, so the displayed approximation represents the signed double integral over the entire square.

Signed Double Integral Visualizer

Midpoint cells above the xy-plane contribute positively; cells below it contribute negatively. Increasing the resolution improves the Riemann-sum approximation.

Signed integral over [-2, 2] × [-2, 2]

21.7600
Exact value
21.3333
Absolute error
4.27e-1
∬R4−x2−y2 dA≈∑f(xi∗,yj∗) ΔA\iint_R 4-x^2-y^2\,dA \approx \sum f(x_i^*,y_j^*)\,\Delta A
Positive contributionNegative contribution

Geometric volume versus signed integral

If the goal is the geometric volume between a surface and the plane and the surface crosses the plane, split the region where ff changes sign and integrate ∣f∣|f|, or otherwise account for the positive and negative portions separately. By contrast, ∬Rf dA\iint_R f\,dA itself keeps the sign of ff.

Iterated Integrals and Fubini's Theorem

Iterated Integral

An iterated integral evaluates a multiple integral through successive single-variable integrations, treating variables not currently integrated as constants.

Fubini's Theorem on a rectangle

If ff is continuous on R=[a,b]×[c,d]R=[a,b]\times[c,d], then

∬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 same conclusion holds under broader integrability hypotheses; continuity is a convenient sufficient condition for this course.

For a rectangle, reversing the order only swaps which constant interval is integrated first. For a nonrectangular region, reversing the order requires reconstructing the region's boundary equations.

Type I and Type II Regions

Type I region: vertical slices

A Type I region can be written

R={(x,y):a≤x≤b, g1(x)≤y≤g2(x)},R=\{(x,y):a\le x\le b,\ g_1(x)\le y\le g_2(x)\},

with g1(x)≤g2(x)g_1(x)\le g_2(x) throughout [a,b][a,b].

Type I Integral

Integrate with respect to y first.

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

Type II region: horizontal slices

A Type II region can be written

R={(x,y):c≤y≤d, h1(y)≤x≤h2(y)},R=\{(x,y):c\le y\le d,\ h_1(y)\le x\le h_2(y)\},

with h1(y)≤h2(y)h_1(y)\le h_2(y) throughout [c,d][c,d].

Type II Integral

Integrate with respect to x first.

∬Rf(x,y) dA=∫cd∫h1(y)h2(y)f(x,y) dx dy.\iint_R f(x,y)\,dA = \int_c^d\int_{h_1(y)}^{h_2(y)} f(x,y)\,dx\,dy.

Changing the Order of Integration

  • Sketch or analyze the region described by the original bounds.
  • Identify its projection onto the new outer-variable axis.
  • For each value of the new outer variable, determine the left/right or lower/upper boundary for the new inner integral.
  • Split the region if a single pair of boundary functions cannot describe every slice.
  • Check that the rebuilt bounds cover exactly the same planar set.

Reversed inner limits are oriented limits, not an ordinary Type I description

An expression such as

∫01∫xx2f(x,y) dy dx\int_0^1\int_x^{x^2} f(x,y)\,dy\,dx

is algebraically valid. However, for 0<x<10<x<1 we have x2<xx^2<x, so the inner limits are reversed. It equals

−∫01∫x2xf(x,y) dy dx.-\int_0^1\int_{x^2}^{x} f(x,y)\,dy\,dx.

The second form is the ordinary Type I region integral; the first is its negative because of orientation.

Polar Double Integrals

Polar coordinates are often advantageous when circles, annuli, radial functions, or angular sectors dominate the geometry. Use x=rcos⁡θx=r\cos\theta, y=rsin⁡θy=r\sin\theta, and r≥0r\ge 0.

Polar area scaling

A small polar cell has radial thickness drdr and approximate arc length r dθr\,d\theta. Its differential area is therefore

dA=r dr dθ.dA=r\,dr\,d\theta.

The factor rr is the area-scaling factor of the polar coordinate transformation. A later topic develops the full Jacobian change-of-variables theorem.

Polar Double Integral

Double integral after converting the planar region to polar coordinates.

∬Rf(x,y) dA=∫αβ∫r1(θ)r2(θ)f(rcos⁡θ,rsin⁡θ) r dr dθ.\iint_R f(x,y)\,dA = \int_{\alpha}^{\beta} \int_{r_1(\theta)}^{r_2(\theta)} f(r\cos\theta,r\sin\theta)\,r\,dr\,d\theta.

Polar Setup

  • Rewrite the region in terms of rr and θ\theta.
  • Rewrite every occurrence of xx and yy in the integrand.
  • Include the factor rr in dAdA.
  • Confirm the angular interval sweeps the intended region exactly once.
  • Check dimensions: r dr dθr\,dr\,d\theta has units of area.

Interactive Simulation

The polar simulation shows a sector and a differential annular wedge. Keep the wedge inside the active region while observing how its area grows with radius.

Polar Coordinates: Double Integral Visualizer

Concept and model scope

Explore why dA=r,dr,dθdA=r\\,dr\\,d\theta. The finite red wedge remains inside the active sector, and its area increases with radius because angular arc length is proportional to rr.

The finite cell uses fixed increments Δr=0.4\Delta r=0.4 and Δθ=0.15\Delta\theta=0.15 rad. It is a local geometric area element for visualization; the reported full integral is evaluated analytically from the selected bounds.

The active examples use f=1f=1, f=rf=r, or f=rsin⁡θf=r\sin\theta so the Jacobian factor can be separated from the integrand itself.

Integrand function
Double Integral Formulation
∫0θmax∫0Rmaxr dr dθ\int_0^{\theta_{max}}\int_0^{R_{max}} r\,dr\,d\theta
Finite wedge area:0.1200 m²
Analytical integral:14.1372 m²
Polar sweeping grid
Polar integration sectorPurple shading shows the active sector. The red annular wedge remains fully inside the radial and angular limits.dA

Wedge bounds: 1.80≤r≤2.201.80\le r\le2.20 and 0.71≤θ≤0.860.71\le\theta\le0.86.

The polar area element

In Cartesian coordinates, dA=dx,dydA=dx\\,dy is locally rectangular.

In polar coordinates, radial thickness is drdr and angular arc length is r,dθr\\,d\theta.

dA=r,dr,dθdA=r\\,dr\\,d\theta

Moving the red wedge outward increases its area while drdr and dθd\theta remain fixed. This illustrates the geometric scaling represented by the Jacobian factor rr.

The current preset reports an analytical integral for the entire purple sector, so the displayed value can be checked independently of the wedge visualization.

Applications to Area and Volume

Area of a Planar Region

Integrate one over the region.

A=∬R1 dA.A=\iint_R 1\,dA.

Signed Volume Under a Graph

Signed accumulation between z=f(x,y) and the xy-plane.

Vsigned=∬Rf(x,y) dA.V_{\mathrm{signed}}=\iint_R f(x,y)\,dA.

If f≥0f\ge 0 on RR, this signed value is also the geometric volume under the surface. If ff changes sign, geometric volume requires handling ∣f∣|f| or splitting the domain.

Mass, First Moments, and Centroids of Laminas

For a thin lamina occupying RR with surface mass density ρ(x,y)≥0\rho(x,y)\ge 0, density has units of mass per area. Integrating density over area gives mass.

Lamina Mass

Total mass of a variable-density lamina.

m=∬Rρ(x,y) dA.m=\iint_R \rho(x,y)\,dA.

First Moments

Moments about the coordinate axes.

Mx=∬Ry ρ(x,y) dA,My=∬Rx ρ(x,y) dA.M_x=\iint_R y\,\rho(x,y)\,dA, \qquad M_y=\iint_R x\,\rho(x,y)\,dA.

Center of Mass

Density-weighted average position.

xˉ=Mym,yˉ=Mxm.\bar{x}=\frac{M_y}{m}, \qquad \bar{y}=\frac{M_x}{m}.

Moment naming

MyM_y contains the lever arm xx because it is the moment about the y-axis. Likewise, MxM_x contains the lever arm yy.

Planar Moments of Inertia

Lamina Moments of Inertia

Resistance of a lamina to angular acceleration about coordinate axes.

Ix=∬Ry2ρ dA,Iy=∬Rx2ρ dA,IO=Ix+Iy=∬R(x2+y2)ρ dA.I_x=\iint_R y^2\rho\,dA, \qquad I_y=\iint_R x^2\rho\,dA, \qquad I_O=I_x+I_y=\iint_R (x^2+y^2)\rho\,dA.

Introductory Surface Area

A smooth surface can be described parametrically by a vector-valued map r(u,v)\mathbf{r}(u,v). The tangent vectors ru\mathbf{r}_u and rv\mathbf{r}_v span a small parallelogram on the surface, so the magnitude of their cross product gives the local area-scaling factor.

Parametric Surface Area

Area of a smooth parametrized surface over a parameter region D.

S=∬D∥ru×rv∥ du dv.S = \iint_D \left\lVert \mathbf{r}_u\times\mathbf{r}_v \right\rVert \,du\,dv.

For a smooth graph z=f(x,y)z=f(x,y) above a planar region RR, use the parameterization r(x,y)=⟨x,y,f(x,y)⟩\mathbf{r}(x,y)=\langle x,y,f(x,y)\rangle. The general formula then reduces to the familiar graph formula.

Surface Area of z=f(x,y)

Graph form of the parametric surface-area formula.

S=∬R1+fx2+fy2 dA.S=\iint_R \sqrt{1+f_x^2+f_y^2}\,dA.

Surface area is not signed volume

The norm ∥ru×rv∥\lVert\mathbf{r}_u\times\mathbf{r}_v\rVert is nonnegative. Surface area measures geometric area, whereas ∬Rf dA\iint_R f\,dA may be signed when ff changes sign.

Key Takeaways
  • A double integral is a signed accumulation over a two-dimensional region.
  • Fubini's Theorem converts a double integral into iterated single-variable integrals.
  • Type I and Type II descriptions encode the same region using different slice directions; changing order requires rebuilding the bounds.
  • Reversed limits change orientation and therefore change sign.
  • Polar coordinates use dA=r dr dθdA=r\,dr\,d\theta; the factor rr is essential.
  • Double integrals compute planar area, volume under nonnegative surfaces, lamina mass, centroids, planar moments of inertia, and surface area through a two-parameter description.
  • Parametric surface area uses ∥ru×rv∥\lVert\mathbf{r}_u\times\mathbf{r}_v\rVert as the local area scale; the graph formula is a special case.
  • Topic 07 extends multiple integration to three-dimensional solids and the full change-of-variables theorem.