Double Integrals in Polar - Calculus 3

SECTION 1

Historical Context & Motivation

The need to evaluate integrals over circular and radially symmetric regions predates even the formal development of multivariable calculus. Early astronomers and physicists confronted problems involving gravitational attraction of spherical bodies, areas of elliptical orbits, and the distribution of mass across circular plates—all of which demanded integration techniques that respected the underlying symmetry. Attempting these computations in Cartesian coordinates typically produces integrands laden with square roots and awkward limits, obscuring the geometric simplicity of the region. The polar coordinate system, which describes points by their distance from the origin and angle from a reference ray, offers a natural remedy. By aligning the coordinate system with the symmetry of the domain, double integrals in polar coordinates reduce many formidable Cartesian integrals to elegant, often closed-form expressions.

1691 • Jakob Bernoulli's Spirals

Jakob Bernoulli systematically studied curves defined in polar form, including the logarithmic spiral, demonstrating the analytical power of radial-angle descriptions for curved geometries.

1736 • Euler's Analytical Methods

Leonhard Euler developed rigorous techniques for changing variables in integrals, laying the algebraic groundwork that would eventually formalize the Jacobian determinant used in coordinate transformations.

1809 • Gauss and the Probability Integral

Carl Friedrich Gauss famously evaluated the Gaussian integral ∫

e^{−x²}

dx from −∞ to ∞ by squaring it and converting to polar coordinates—an iconic application that remains a standard example in analysis courses today.

1841 • Jacobi Formalizes the Change of Variables

Carl Gustav Jacob Jacobi published his general theory of the Jacobian determinant, providing the rigorous justification for the extra factor of r that appears in the polar area element dA = r dr dθ.

The central question that this lesson addresses is deceptively simple: when we convert a double integral from Cartesian to polar coordinates, why does the area element dx dy become r dr dθ, and how do we correctly determine the new limits of integration? Answering this question requires understanding both the geometry of polar grids and the algebra of the Jacobian, and it unlocks a powerful technique for evaluating integrals over disks, annuli, sectors, cardioids, and many other regions.

SECTION 2

Core Principles & Definitions

Before diving into computation, it is essential to internalize the foundational ideas that govern polar double integrals. The transition from rectangular to polar coordinates is more than a substitution trick; it is a genuine change of geometric perspective that reshapes both the integrand and the differential area element. The following core concepts form the backbone of every polar integration problem.

Polar Coordinate System

Every point in the plane is described by (r, θ), where r ≥ 0 is the radial distance from the origin and θ is the angle measured counterclockwise from the positive x-axis. The conversion relations are x = r cos θ and y = r sin θ.

The Polar Area Element

In Cartesian coordinates, the infinitesimal area element is dx dy. In polar coordinates, a small patch of area is not a rectangle but a curvilinear trapezoid whose area is dA = r dr dθ. The factor r accounts for arcs being longer at greater radii.

The Jacobian Determinant

The factor r arises formally as the absolute value of the Jacobian of the transformation (x, y) → (r, θ). Computing ∂(x, y)/∂(r, θ) yields J = r, confirming the geometric intuition.

Iterated Polar Integrals

A polar double integral is typically written as an iterated integral: first integrate with respect to r (radial direction), then with respect to θ (angular sweep). The limits on r may depend on θ if the boundary is not a circle centered at the origin.

When to Use Polar Coordinates

Polar coordinates are most advantageous when the region of integration is a disk, annulus, sector, or other shape naturally described by r and θ, and when the integrand involves expressions like x² + y², which simplify to r² in polar form.

✦

KEY TAKEAWAY

Think of a polar grid as a dartboard. Near the bullseye, the cells are tiny; near the outer ring, the cells are much larger even though they subtend the same angle and span the same radial increment. The factor r in r dr dθ corrects for this stretching—just as a mapmaker must scale area differently near the equator versus the poles. Without the extra r, you would systematically undercount area far from the origin and overcount area near it.

SECTION 3

Visualizing the Polar Area Element

The most important conceptual hurdle in polar double integrals is understanding why the infinitesimal area element is dA = r dr dθ rather than simply dr dθ. A Cartesian grid partitions the plane into small rectangles of uniform size dx × dy. A polar grid, by contrast, partitions the plane into curvilinear patches whose dimensions depend on the radial distance r. The diagram below illustrates a single polar patch and how its area is computed.

A polar area element (shaded patch) bounded by radii r and r + dr and angles θ and θ + dθ. The inner arc has length r dθ, the outer arc has length (r + dr) dθ, and the radial sides have length dr. The area of this patch is approximately r dr dθ.

In the diagram above, the shaded region represents a single infinitesimal polar patch. Its inner boundary is an arc of radius r and its outer boundary is an arc of radius r + dr. The two radial sides are separated by the angle dθ. Because the arc length at radius r is r dθ, the patch closely resembles a rectangle with dimensions dr (radially) and r dθ (tangentially), yielding an area of dA ≈ r dr dθ. As dr and dθ shrink to zero, this approximation becomes exact. This geometric argument provides the intuitive underpinning for the Jacobian factor r that appears in every polar double integral.

SECTION 4

Mathematical Framework

We now formalize the conversion of a double integral from Cartesian to polar coordinates. Let R be a region in the xy-plane, and suppose we wish to evaluate ∬ᵣ f(x, y) dA. The change-of-variables theorem for multiple integrals tells us how to rewrite this integral in terms of new coordinates (r, θ), provided we include the absolute value of the Jacobian determinant of the transformation.

The Coordinate Transformation

POLAR-TO-CARTESIAN

x = r cos θ, y = r sin θ

where r ≥ 0 is the radial distance and θ ∈ [0, 2π) is the polar angle. The inverse relations are r = √(x² + y²) and θ = arctan(y/x) (with appropriate quadrant adjustment).

Computing the Jacobian

The Jacobian of the transformation T(r, θ) = (r cos θ, r sin θ) is the determinant of the matrix of first-order partial derivatives:

JACOBIAN MATRIX

∂(x, y)/∂(r, θ) = | ∂x/∂r ∂x/∂θ | = | cos θ −r sin θ | = r cos²θ + r sin²θ = r | ∂y/∂r ∂y/∂θ | | sin θ r cos θ |

The absolute value |J| = |r| = r (since r ≥ 0). This confirms the geometric observation that the area element scales linearly with r.

The General Conversion Formula

DOUBLE INTEGRAL IN POLAR

∬ᵣ f(x, y) dA = ∫ᵅᵝ ∫ₕ₁(θ)ʰ²⁽θ⁾ f(r cos θ, r sin θ) · r dr dθ

Here α ≤ θ ≤ β defines the angular sweep, and h₁(θ) ≤ r ≤ h₂(θ) defines the radial extent, which may vary with θ. The crucial factor r multiplying dr dθ must never be omitted.

⚠

Common Mistake

The single most frequent error students make when converting to polar coordinates is forgetting the factor of r in the area element. Remember: dA = r dr dθ, not dr dθ. Without this factor, the integral will be dimensionally incorrect and will yield a wrong numerical answer.

USEFUL IDENTITY

x² + y² = r²

This identity dramatically simplifies integrands containing x² + y². For example, $e^{−(x²+y²)}$ becomes $e^{−r²}$ , and 1/√(x² + y²) becomes 1/r.

SECTION 5

Setting Up Limits of Integration

Correctly determining the limits of integration is often the most challenging step. Unlike Cartesian double integrals where the region is sliced by horizontal or vertical strips, polar integrals slice the region with radial rays emanating from the origin. For each fixed angle θ, the variable r ranges from the inner boundary to the outer boundary of the region. The angle θ then sweeps through all directions that intersect the region. The table below catalogs the most common region types and their corresponding polar limits.

Region Type	Description	Polar Limits
Full Disk	x² + y² ≤ a²	0 ≤ r ≤ a, 0 ≤ θ ≤ 2π
Annulus	a² ≤ x² + y² ≤ b²	a ≤ r ≤ b, 0 ≤ θ ≤ 2π
Sector	Wedge-shaped region of a disk	0 ≤ r ≤ a, α ≤ θ ≤ β
Cardioid	Interior of r = 1 + cos θ	0 ≤ r ≤ 1 + cos θ, 0 ≤ θ ≤ 2π
Semicircle (upper)	x² + y² ≤ a², y ≥ 0	0 ≤ r ≤ a, 0 ≤ θ ≤ π
Off-center circle	x² + y² = 2ax (circle through origin)	0 ≤ r ≤ 2a cos θ, −π/2 ≤ θ ≤ π/2

Common polar regions and their integration limits

Three common polar regions: a full disk with constant radial bound, a sector with restricted angular range, and a cardioid where the upper limit on r depends on θ.

When setting up limits, always begin by sketching the region and drawing a typical radial ray from the origin at a generic angle θ. Identify where the ray enters the region (this gives the lower limit on r) and where it exits (the upper limit on r). Then determine the full angular range over which such rays intersect the region. This systematic approach prevents errors, especially for regions like off-center circles and rose curves where the limits on r are functions of θ.

SECTION 6

Worked Example: Area Under $e^{−(x²+y²)}$ Over a Disk

Let us evaluate the integral ∬ᵣ $e^{−(x²+y²)}$ dA, where R is the disk x² + y² ≤ 9 (a disk of radius 3 centered at the origin). This integral arises naturally in probability theory and physics, and it is practically impossible to evaluate in Cartesian coordinates because the antiderivative of $e^{−x²}$ cannot be expressed in terms of elementary functions. Polar coordinates transform the problem into a straightforward computation.

Evaluate ∬ᵣ $e^{−(x²+y²)}$ dA over x² + y² ≤ 9

1
Step 1 — Identify the Region and Convert

The region R is a disk of radius 3 centered at the origin. In polar coordinates, this is described by 0 ≤ r ≤ 3 and 0 ≤ θ ≤ 2π. The integrand becomes $e^{−(x²+y²)}$ = $e^{−r²}$ , and the area element is dA = r dr dθ.

∬ᵣ $e^{−(x²+y²)}$ dA = ∫₀²π ∫₀³ $e^{−r²}$ · r dr dθ
2
Step 2 — Evaluate the Inner Integral (with respect to r)

The inner integral is ∫₀³ r $e^{−r²}$ dr. Use the substitution u = −r², so du = −2r dr, which gives r dr = −du/2. When r = 0, u = 0; when r = 3, u = −9. The integral becomes −(1/2) ∫₀⁻⁹ eᵘ du = (1/2)[−eᵘ]₀⁻⁹ = (1/2)(1 − e⁻⁹).

∫₀³ r $e^{−r²}$ dr = (1/2)(1 − e⁻⁹)
3
Step 3 — Evaluate the Outer Integral (with respect to θ)

Since the result of the inner integral does not depend on θ, the outer integral is simply ∫₀²π (1/2)(1 − e⁻⁹) dθ = (1/2)(1 − e⁻⁹) · [θ]₀²π = π(1 − e⁻⁹).

∬ᵣ $e^{−(x²+y²)}$ dA = π(1 − e⁻⁹)
4
Step 4 — Interpret the Result

Since e⁻⁹ ≈ 0.000123, the value is approximately π(1 − 0.000123) ≈ π(0.999877) ≈ 3.1412. Note that as the radius of the disk grows to infinity, the integral approaches π, recovering the famous result ∫₋∞^∞ ∫₋∞^∞ $e^{−(x²+y²)}$ dx dy = π, which is the square of the Gaussian integral √π.

π(1 − e⁻⁹) ≈ 3.1412

💡

KEY INSIGHT

The magic of this example lies in the substitution u = −r². In Cartesian coordinates, the integrand

e^{−x²}

has no elementary antiderivative. But when we convert to polar, the extra factor of r from the Jacobian pairs perfectly with the

e^{−r²}

to create r

e^{−r²}

, which is trivially integrable by u-substitution. This is a prime example of polar coordinates not merely simplifying an integral but making it possible.

SECTION 7

When to Choose Polar Over Cartesian

Not every double integral benefits from a switch to polar coordinates. The choice depends on both the geometry of the region and the structure of the integrand. The table below compares the two coordinate systems across several important criteria, helping you develop the judgment to choose wisely in any given problem.

Criterion	Cartesian (x, y)	Polar (r, θ)
Best region shapes	Rectangles, triangles, regions bounded by lines or parabolas	Disks, annuli, sectors, cardioids, rose curves, spirals
Integrand simplification	Polynomials in x, y; products of single-variable functions	Functions of x² + y² (= r²); trigonometric dependence on angle
Area element	dA = dx dy (no extra factor)	dA = r dr dθ (Jacobian factor r required)
Typical difficulty	Limits are usually linear functions; integrands may involve √(a² − x²)	Limits may involve trigonometric functions of θ; integrand is often simpler
Key pitfall	Awkward limits when region has curved boundaries	Forgetting the Jacobian factor r; wrong angular limits

Comparison of Cartesian and polar approaches for double integrals

✦

DECISION HEURISTIC

Ask two questions before choosing coordinates. First: is the region naturally described by circles or angles? If yes, polar is likely better. Second: does the integrand simplify when x² + y² is replaced by r²? If the answer to either question is yes, convert to polar. If both answers are no—say, you are integrating a polynomial over a rectangle—stay in Cartesian.

💡

Hybrid Regions

Some regions are described most easily with a mix of Cartesian and polar descriptions. For instance, the region above the line y = 1 inside the circle x² + y² = 4 requires careful conversion: the line y = 1 becomes r sin θ = 1, i.e., r = csc θ, and the circle is r = 2. You would integrate with limits r from csc θ to 2 and θ from π/6 to 5π/6.

SECTION 8

Connection to Triple Integrals & Beyond

The idea of converting to coordinates that align with the symmetry of the problem extends far beyond polar double integrals. In three dimensions, the same philosophy leads to cylindrical coordinates (r, θ, z) and spherical coordinates (ρ, φ, θ), each with their own Jacobian factors. Understanding the polar Jacobian r thoroughly prepares you for these generalizations, because the underlying principle—the Jacobian captures how volumes distort under a change of variables—is identical.

Coordinate System	Dimensions	Volume/Area Element	Jacobian Factor
Polar	2D	dA = r dr dθ	r
Cylindrical	3D	dV = r dr dθ dz	r
Spherical	3D	dV = ρ² sin φ dρ dφ dθ	ρ² sin φ

Comparison of coordinate systems and their Jacobian factors

Notice that cylindrical coordinates are essentially polar coordinates with an extra Cartesian z-dimension appended, so the Jacobian remains r. Spherical coordinates introduce a more complex Jacobian ρ² sin φ because the transformation involves two angular variables. In all cases, the Jacobian is derived from the same determinant procedure you learned in this lesson. As you progress through Calculus 3, you will also encounter the general change-of-variables theorem for arbitrary coordinate transformations, which unifies all of these special cases under a single framework.

ℹ

Looking Ahead

Surface integrals, line integrals, and the theorems of Green, Stokes, and Gauss all build upon the change-of-variables machinery introduced here. Mastering the polar double integral—and deeply understanding why the Jacobian appears—provides the conceptual foundation for virtually all of multivariable and vector calculus.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain geometrically why the polar area element is dA = r dr dθ rather than simply dA = dr dθ. What would go wrong if the factor r were omitted?

Show answer

A polar grid divides the plane into curvilinear patches bounded by arcs and radial lines. At radius r, the tangential width of such a patch is r dθ (arc length), not simply dθ. Since the radial width is dr, the area of the patch is approximately (r dθ)(dr) = r dr dθ. If the factor r were omitted, the integral would treat all patches as having the same tangential width regardless of their distance from the origin. This would systematically undercount area at large radii and overcount near the origin, yielding an incorrect result. Formally, r is the Jacobian determinant of the transformation from (r, θ) to (x, y), and omitting it violates the change-of-variables theorem.

Question 2: BASIC CALCULATION

Evaluate ∬ᵣ (x² + y²) dA, where R is the disk x² + y² ≤ 4.

Show answer

Convert to polar: x² + y² = r², dA = r dr dθ, and the region is 0 ≤ r ≤ 2, 0 ≤ θ ≤ 2π. The integral becomes ∫₀²π ∫₀² r² · r dr dθ = ∫₀²π ∫₀² r³ dr dθ. The inner integral is [r⁴/4]₀² = 16/4 = 4. The outer integral is ∫₀²π 4 dθ = 4 · 2π = 8π. Therefore the answer is 8π. This result can be interpreted as the second moment of area (with respect to the origin) of a unit-density disk of radius 2.

Question 3: INTERMEDIATE

Evaluate ∫₀² ∫₀^(√(4−x²)) √(x² + y²) dy dx by converting to polar coordinates.

Show answer

The Cartesian limits describe the quarter-disk in the first quadrant: 0 ≤ x ≤ 2, 0 ≤ y ≤ √(4 − x²), which corresponds to the region x² + y² ≤ 4 with x ≥ 0 and y ≥ 0. In polar coordinates, this is 0 ≤ r ≤ 2, 0 ≤ θ ≤ π/2. The integrand √(x² + y²) = r, and dA = r dr dθ. The integral becomes ∫₀^(π/2) ∫₀² r · r dr dθ = ∫₀^(π/2) ∫₀² r² dr dθ. The inner integral is [r³/3]₀² = 8/3. The outer integral is (8/3)(π/2) = 4π/3. The answer is 4π/3 ≈ 4.189.

Question 4: APPLIED

A circular plate of radius R meters has a temperature distribution given by T(x, y) = T₀

e^{−(x²+y²)/σ²}

degrees Celsius, where T₀ = 100 °C and σ = 2 m. Find the average temperature over the plate when R = 3 m.

Show answer

The average temperature is T̄ = (1/A) ∬ᵣ T dA, where A = πR² = 9π. Converting to polar: T = 100

e^{−r²/4}

, dA = r dr dθ, with 0 ≤ r ≤ 3, 0 ≤ θ ≤ 2π. Thus T̄ = (1/9π) ∫₀²π ∫₀³ 100

e^{−r²/4}

· r dr dθ = (100/9π) · 2π ∫₀³ r

e^{−r²/4}

dr = (200/9) ∫₀³ r

e^{−r²/4}

dr. Use u = −r²/4, du = −r/2 dr: the integral becomes (200/9) · [−2

e^{−r²/4}

]₀³ = (200/9)(2)(1 −

e^{−9/4}

) = (400/9)(1 −

e^{−2.25}

) ≈ (44.44)(1 − 0.1054) ≈ 44.44 × 0.8946 ≈ 39.76 °C. The average temperature is approximately 39.8 °C.

Question 5: CRITICAL THINKING

Derive the formula for the area enclosed by one petal of the rose curve r = cos(2θ) using a double integral in polar coordinates. Explain why the limits on θ must be chosen carefully and identify the correct angular bounds.

Show answer

The rose r = cos(2θ) has four petals. One petal lies where cos(2θ) ≥ 0. For the petal along the positive x-axis, cos(2θ) ≥ 0 when −π/4 ≤ θ ≤ π/4. The area is A = ∬ dA = ∫₋

(π/4)^{π/4}

∫₀^(cos 2θ) r dr dθ. The inner integral gives [r²/2]₀^(cos 2θ) = cos²(2θ)/2. So A = (1/2) ∫₋

(π/4)^{π/4}

cos²(2θ) dθ. Using the identity cos²(2θ) = (1 + cos(4θ))/2, we get A = (1/4) ∫₋

(π/4)^{π/4}

(1 + cos 4θ) dθ = (1/4)[θ + sin(4θ)/4]₋

(π/4)^{π/4}

= (1/4)[(π/4 + 0) − (−π/4 + 0)] = (1/4)(π/2) = π/8. The careful choice of angular limits is essential: if we naively integrated from 0 to 2π, the radial function cos(2θ) takes negative values in some intervals, and the formula r = cos(2θ) does not trace a valid curve there. Integrating only where cos(2θ) ≥ 0 isolates exactly one petal. The area of a single petal is π/8.

SUMMARY

Lesson Summary

Summary

Converting a double integral to polar coordinates involves three essential steps: substituting x = r cos θ and y = r sin θ into the integrand, replacing the Cartesian area element dx dy with the polar area element r dr dθ, and re-expressing the region of integration in terms of radial and angular bounds. The factor r in the area element arises from the Jacobian determinant of the coordinate transformation and geometrically reflects the fact that arcs grow longer as the radial distance increases.

This technique is most powerful when the region is a disk, annulus, sector, or polar curve, and when the integrand involves x² + y² = r². The classic example is the Gaussian integral, where the Jacobian factor r creates a perfect u-substitution opportunity with $e^{−r²}$ . The same change-of-variables philosophy extends to cylindrical and spherical coordinates in three dimensions, making this lesson a gateway to all subsequent coordinate-transformation techniques in multivariable calculus.