Unit 4, lesson 26

Double Integrals in Polar Coordinates

Where the extra factor r in polar double integrals comes from, and how to choose radius and angle bounds that cover a disk, ring, sector or shifted circle exactly once.

  • Theory part 51
  • Practice part 52
  • Double Integrals

A round pizza is awkward to cut into a grid of little squares, yet easy to cut into slices and rings. Disks, rings and wedges are just as awkward in xx and yy, where every bound turns into a square root. In polar coordinates the same regions have constant bounds. The price is one extra factor, and getting that factor right is most of the lesson.

Key terms

Polar coordinates and their range

Definition Polar coordinates (r,θ)(r, \theta) describe the point

x=rcos⁡θ,y=rsin⁡θ,x = r\cos\theta, \qquad y = r\sin\theta,

where r≥0r \ge 0 is the distance to the origin and θ\theta is an angle in radians measured counterclockwise from the positive xx axis. Then r2=x2+y2r^2 = x^2 + y^2. To give each point other than the origin exactly one pair, restrict θ\theta to a half open interval of length 2π2\pi, such as 0≤θ<2π0 \le \theta < 2\pi or −π<θ≤π-\pi < \theta \le \pi. The origin has r=0r = 0 and any θ\theta.

In plain words The radius says how far out, the angle says which direction. Keeping rr nonnegative and the angle inside one full turn stops the same point from being named twice.

Example The point (−1,−1)(-1, -1) has r=2r = \sqrt{2}. It lies in the third quadrant, so θ=5π/4\theta = 5\pi/4 in the range [0,2π)[0, 2\pi), or θ=−3π/4\theta = -3\pi/4 in the range (−π,π](-\pi, \pi]. Note that arctan⁡(y/x)=arctan⁡1=π/4\arctan(y/x) = \arctan 1 = \pi/4 points to the wrong quadrant.

Sector area scaling

Definition A polar rectangle is a set

S={(r,θ)∣r1≤r≤r2,θ1≤θ≤θ2}\begin{aligned} S = \{ (r, \theta) \mid{} & r_1 \le r \le r_2, \\ & \theta_1 \le \theta \le \theta_2 \} \end{aligned}

with 0≤r1≤r20 \le r_1 \le r_2 and 0≤θ2−θ1≤2π0 \le \theta_2 - \theta_1 \le 2\pi. Since a sector of radius rr and angle Δθ\Delta\theta has area 12r2Δθ\tfrac{1}{2} r^2 \Delta\theta, the area of SS is

ΔA=12(r22−r12) Δθ=r∗ Δr Δθ,\begin{aligned} \Delta A &= \tfrac{1}{2}(r_2^2 - r_1^2)\,\Delta\theta \\[0.4em] &= r^* \, \Delta r \, \Delta\theta, \end{aligned}

where Δr=r2−r1\Delta r = r_2 - r_1, Δθ=θ2−θ1\Delta\theta = \theta_2 - \theta_1 and r∗=12(r1+r2)r^* = \tfrac{1}{2}(r_1 + r_2) is the average radius.

In plain words A curved patch of a ring is almost a rectangle with sides Δr\Delta r and r Δθr \, \Delta\theta. The farther out the patch sits, the longer its arc, so the area grows in proportion to rr.

Example For 1≤r≤31 \le r \le 3 and 0≤θ≤π/60 \le \theta \le \pi/6, the sector formula gives

12(9−1)⋅π6=2π3,\frac{1}{2}(9 - 1) \cdot \frac{\pi}{6} = \frac{2\pi}{3},

and the average radius formula gives 2⋅2⋅π/6=2π/32 \cdot 2 \cdot \pi/6 = 2\pi/3 as well. The same angle and radial width at 5≤r≤75 \le r \le 7 gives an area three times larger, since the average radius is 66 instead of 22.

The polar area element

Definition Let ff be continuous on the polar rectangle RR given by 0≤a≤r≤b0 \le a \le r \le b and α≤θ≤β\alpha \le \theta \le \beta, with 0≤β−α≤2π0 \le \beta - \alpha \le 2\pi. Then

∬Rf dA=∫αβ∫abf r dr dθ,\iint_R f \, dA = \int_\alpha^\beta \int_a^b f \, r \, dr \, d\theta,

where ff on the right is evaluated at (rcos⁡θ,rsin⁡θ)(r\cos\theta, r\sin\theta).

In short, dA=r dr dθdA = r \, dr \, d\theta. The factor rr is the absolute value of the Jacobian of the map (r,θ)↦(rcos⁡θ,rsin⁡θ)(r, \theta) \mapsto (r\cos\theta, r\sin\theta), which is one-to-one on the region apart from a boundary set of zero area.

In plain words Replace xx and yy by their polar forms, and never forget the extra rr. It is the stretching factor from the sector area above.

Example For the disk x2+y2≤4x^2 + y^2 \le 4,

∬e−(x2+y2) dA=∫02π∫02e−r2 r dr dθ=2π⋅1−e−42=π(1−e−4).\begin{aligned} &\iint e^{-(x^2 + y^2)} \, dA \\[0.4em] &\quad = \int_0^{2\pi} \int_0^2 e^{-r^2} \, r \, dr \, d\theta \\[0.4em] &\quad = 2\pi \cdot \frac{1 - e^{-4}}{2} = \pi \left( 1 - e^{-4} \right). \end{aligned}

Without the rr, the inner integral ∫02e−r2 dr\int_0^2 e^{-r^2} \, dr has no elementary antiderivative.

Polar regions

Definition A polar region of the general kind is

D={(r,θ)∣α≤θ≤β,h1(θ)≤r≤h2(θ)},\begin{aligned} D = \{ (r, \theta) \mid{} & \alpha \le \theta \le \beta, \\ & h_1(\theta) \le r \le h_2(\theta) \}, \end{aligned}

where 0≤β−α≤2π0 \le \beta - \alpha \le 2\pi and h1h_1, h2h_2 are continuous on [α,β][\alpha, \beta] with 0≤h1≤h20 \le h_1 \le h_2. If ff is continuous on DD, then

∬Df dA=∫αβ∫h1(θ)h2(θ)f r dr dθ,\iint_D f \, dA = \int_\alpha^\beta \int_{h_1(\theta)}^{h_2(\theta)} f \, r \, dr \, d\theta,

with ff evaluated at (rcos⁡θ,rsin⁡θ)(r\cos\theta, r\sin\theta).

In plain words For each direction θ\theta, a ray from the origin enters the region at r=h1(θ)r = h_1(\theta) and leaves at r=h2(θ)r = h_2(\theta). That ray plays the part of the vertical strip from rectangular coordinates.

Example Take the upper half of the ring 1≤x2+y2≤41 \le x^2 + y^2 \le 4, so 0≤θ≤π0 \le \theta \le \pi and 1≤r≤21 \le r \le 2. Then

∬Dy dA=∫0π∫12(rsin⁡θ) r dr dθ=∫0πsin⁡θ dθ⋅∫12r2 dr=2⋅73=143.\begin{aligned} \iint_D y \, dA &= \int_0^\pi \int_1^2 (r\sin\theta) \, r \, dr \, d\theta \\[0.4em] &= \int_0^\pi \sin\theta \, d\theta \cdot \int_1^2 r^2 \, dr \\[0.4em] &= 2 \cdot \frac{7}{3} = \frac{14}{3}. \end{aligned}

Bounds from the geometry

Definition To set up polar bounds, convert each boundary curve to polar form using x=rcos⁡θx = r\cos\theta, y=rsin⁡θy = r\sin\theta and x2+y2=r2x^2 + y^2 = r^2, and keep only the branch with r≥0r \ge 0. Then find the angles of the rays that touch the region, which give α\alpha and β\beta, and for each such angle read off where the ray enters and leaves. Intersection angles come from setting the two polar equations equal.

In plain words Stand at the origin, sweep a ray through every direction that hits the region, and note where the ray goes in and comes out.

Example The circle x2+y2=2xx^2 + y^2 = 2x has center (1,0)(1, 0) and radius 11. In polar form r2=2rcos⁡θr^2 = 2r\cos\theta, so r=2cos⁡θr = 2\cos\theta, which is nonnegative for −π/2≤θ≤π/2-\pi/2 \le \theta \le \pi/2. The volume under the cone z=x2+y2z = \sqrt{x^2 + y^2} above this disk is

∫−π/2π/2∫02cos⁡θr⋅r dr dθ=∫−π/2π/283cos⁡3θ dθ=83⋅43=329.\begin{aligned} &\int_{-\pi/2}^{\pi/2} \int_0^{2\cos\theta} r \cdot r \, dr \, d\theta \\[0.4em] &\quad = \int_{-\pi/2}^{\pi/2} \frac{8}{3}\cos^3\theta \, d\theta \\[0.4em] &\quad = \frac{8}{3} \cdot \frac{4}{3} = \frac{32}{9}. \end{aligned}

For a region between two curves, take the part inside r=2cos⁡θr = 2\cos\theta and outside r=1r = 1. The curves meet where 2cos⁡θ=12\cos\theta = 1, at θ=±π/3\theta = \pm\pi/3, so its area is

∫−π/3π/3∫12cos⁡θr dr dθ=∫−π/3π/34cos⁡2θ−12 dθ=π3+32.\begin{aligned} &\int_{-\pi/3}^{\pi/3} \int_1^{2\cos\theta} r \, dr \, d\theta \\[0.4em] &\quad = \int_{-\pi/3}^{\pi/3} \frac{4\cos^2\theta - 1}{2} \, d\theta \\[0.4em] &\quad = \frac{\pi}{3} + \frac{\sqrt{3}}{2}. \end{aligned}

Avoiding overlap

Definition The polar formula counts every point of DD once only if the map (r,θ)↦(rcos⁡θ,rsin⁡θ)(r, \theta) \mapsto (r\cos\theta, r\sin\theta) is one-to-one on the parameter region, apart from a set of zero area there. Typical exceptions are the edge r=0r = 0, which all goes to the origin, and the two edges θ=α\theta = \alpha and θ=β\theta = \beta when β−α=2π\beta - \alpha = 2\pi, which go to the same ray. So the angle interval must not exceed length 2π2\pi, every bound must satisfy 0≤h1(θ)≤h2(θ)0 \le h_1(\theta) \le h_2(\theta), and when a curve r=h(θ)r = h(\theta) traces itself more than once, the angle interval must cover only one tracing.

In plain words Running the angle too far paints part of the region twice, and the integral counts it twice. Check how much angle the curve needs to be drawn once.

Example The circle r=2sin⁡θr = 2\sin\theta has center (0,1)(0, 1) and radius 11, and it is traced once as θ\theta runs from 00 to π\pi. Its area is

∫0π∫02sin⁡θr dr dθ=∫0π2sin⁡2θ dθ=π,\begin{aligned} \int_0^\pi \int_0^{2\sin\theta} r \, dr \, d\theta &= \int_0^\pi 2\sin^2\theta \, d\theta \\[0.4em] &= \pi, \end{aligned}

as it should be. Now let θ\theta run from 00 to 2π2\pi by mistake. For π<θ<2π\pi < \theta < 2\pi the upper bound 2sin⁡θ2\sin\theta is negative, so there is no rr with 0≤r≤2sin⁡θ0 \le r \le 2\sin\theta and those directions contain no point of the disk. Yet the inner integral computed blindly is still 2sin⁡2θ>02\sin^2\theta > 0, and the total comes out as 2π2\pi. If negative rr were allowed, those angles would trace the same circle a second time and sweep the disk again, and that second sweep is exactly what the extra π\pi counts.

Check yourself

  1. Use polar coordinates to evaluate the integral of x2x^2 over the unit disk x2+y2≤1x^2 + y^2 \le 1.
Show the answer

With x2=r2cos⁡2θx^2 = r^2\cos^2\theta and dA=r dr dθdA = r \, dr \, d\theta,

∫02π∫01r3cos⁡2θ dr dθ=14∫02πcos⁡2θ dθ=π4.\begin{aligned} &\int_0^{2\pi} \int_0^1 r^3 \cos^2\theta \, dr \, d\theta \\[0.4em] &\quad = \frac{1}{4} \int_0^{2\pi} \cos^2\theta \, d\theta = \frac{\pi}{4}. \end{aligned}
  1. Convert to polar coordinates and evaluate
∫02∫04−x2xy dy dx.\int_0^2 \int_0^{\sqrt{4 - x^2}} xy \, dy \, dx.
Show the answer

The region is the quarter disk of radius 22 in the first quadrant, so 0≤θ≤π/20 \le \theta \le \pi/2 and 0≤r≤20 \le r \le 2. The integrand times rr is r3cos⁡θsin⁡θr^3\cos\theta\sin\theta, and

∫0π/2cos⁡θsin⁡θ dθ⋅∫02r3 dr=12⋅4=2.\begin{aligned} &\int_0^{\pi/2} \cos\theta\sin\theta \, d\theta \cdot \int_0^2 r^3 \, dr \\[0.4em] &\quad = \frac{1}{2} \cdot 4 = 2. \end{aligned}
  1. Find the area of one petal of the rose r=cos⁡2θr = \cos 2\theta, and say which angle range to use.
Show the answer

The petal on the positive xx axis is traced once for −π/4≤θ≤π/4-\pi/4 \le \theta \le \pi/4, where cos⁡2θ≥0\cos 2\theta \ge 0. Its area is

∫−π/4π/4∫0cos⁡2θr dr dθ=∫−π/4π/4cos⁡22θ2 dθ=π8.\begin{aligned} &\int_{-\pi/4}^{\pi/4} \int_0^{\cos 2\theta} r \, dr \, d\theta \\[0.4em] &\quad = \int_{-\pi/4}^{\pi/4} \frac{\cos^2 2\theta}{2} \, d\theta = \frac{\pi}{8}. \end{aligned}

How this lesson connects

It also uses Polar coordinate conversion with nonnegative radius and Double integrals and elementary trigonometric integration.