← Back to calculus 3

Double Integrals in Polar Coordinates

Double Integrals in Polar Coordinates

Some regions and integrands have circular symmetry that makes Cartesian coordinates cumbersome. Polar coordinates (r,θ)(r, \theta) — where x=rcosθx = r\cos\theta and y=rsinθy = r\sin\theta — naturally describe disks, sectors, and annuli. This note shows how to set up and evaluate double integrals in polar coordinates.


1. The polar rectangle

Recall from Polar Coordinates that a point in the plane is described by (r,θ)(r, \theta) with r0r \ge 0 and θ[0,2π)\theta \in [0, 2\pi). Just as Cartesian rectangles [a,b]×[c,d][a,b]\times[c,d] are the building block for Cartesian integrals, polar rectangles are the natural building block for polar integrals.

A polar rectangle is the set

R={(r,θ):arb,    αθβ},R = \{\,(r,\theta) : a \le r \le b,\;\; \alpha \le \theta \le \beta\,\},

where 0a<b0 \le a < b and 0α<β2π0 \le \alpha < \beta \le 2\pi.

In Cartesian coordinates, RR is the region between two arcs of circles (radii aa and bb) and two radial lines (angles α\alpha and β\beta) — sometimes called an annular sector.

Use the sliders below to explore different polar rectangles.

Polar Rectangle: R = {(r,θ) : a ≤ r ≤ b, α ≤ θ ≤ β}

True area = ½(β−α)(b²−a²) = 4.1233
Naive ∫∫ dr dθ = (β−α)(b−a) = 2.3562 (wrong!)

2. Why we cannot use drdθdr\,d\theta naively

A natural first guess for the polar area element might be drdθdr\,d\theta (by analogy with dxdydx\,dy). Let us check this against a known area.

Take the full disk of radius RR: 0rR0 \le r \le R, 0θ2π0 \le \theta \le 2\pi. The area of a disk is πR2\pi R^2.

Naive attempt:

02π ⁣0Rdrdθ=02πRdθ=2πR.\int_0^{2\pi}\!\int_0^R dr\,d\theta = \int_0^{2\pi} R\,d\theta = 2\pi R.

This gives 2πR2\pi R, not πR2\pi R^2. The naive element drdθdr\,d\theta is wrong.

The interactive graph above displays the "true area" versus the "naive area" computed without the rr factor — you can see they differ whenever a0a \ne 0 or the sector is not degenerate.


3. The polar area element: dA=rdrdθdA = r\,dr\,d\theta

The correct area element in polar coordinates is rdrdθr\,dr\,d\theta. Here is where it comes from.

Geometry of a thin polar rectangle. Consider the polar rectangle with r[r0,r0+Δr]r \in [r_0, r_0 + \Delta r] and θ[θ0,θ0+Δθ]\theta \in [\theta_0, \theta_0 + \Delta\theta].

Its area is the area of the outer sector minus the area of the inner sector:

ΔA=12(r0+Δr)2Δθ12r02Δθ=Δθ2[(r0+Δr)2r02].\Delta A = \frac{1}{2}(r_0+\Delta r)^2\,\Delta\theta - \frac{1}{2}r_0^2\,\Delta\theta = \frac{\Delta\theta}{2}\bigl[(r_0+\Delta r)^2 - r_0^2\bigr].

Expanding:

ΔA=Δθ2[2r0Δr+(Δr)2]r0ΔrΔθas Δr0.\Delta A = \frac{\Delta\theta}{2}\bigl[2r_0\,\Delta r + (\Delta r)^2\bigr] \approx r_0\,\Delta r\,\Delta\theta \quad \text{as } \Delta r \to 0.

In the limit, ΔArdrdθ\Delta A \to r\,dr\,d\theta. The factor of rr accounts for the fact that arcs at large rr are much longer than arcs near the origin.

dA=rdrdθ.\boxed{dA = r\,dr\,d\theta.}

Verification. For the disk of radius RR:

02π ⁣0Rrdrdθ=02πR22dθ=πR2.\int_0^{2\pi}\!\int_0^R r\,dr\,d\theta = \int_0^{2\pi} \frac{R^2}{2}\,d\theta = \pi R^2. \checkmark


4. The change-of-variables formula for polar coordinates

Let f(x,y)f(x,y) be continuous on the polar rectangle R={arb,αθβ}R = \{a \le r \le b,\, \alpha \le \theta \le \beta\}. Then

Rf(x,y)dA=αβ ⁣abf(rcosθ,rsinθ)  rdrdθ.\boxed{\iint_R f(x,y)\,dA = \int_\alpha^\beta\!\int_a^b f(r\cos\theta,\, r\sin\theta)\; r\,dr\,d\theta.}

How to use this:

  1. Replace xrcosθx \to r\cos\theta and yrsinθy \to r\sin\theta in the integrand.
  2. Replace dArdrdθdA \to r\,dr\,d\theta.
  3. Set the limits: rr from aa to bb, θ\theta from α\alpha to β\beta.

Example. Evaluate D(x2+y2)dA\displaystyle\iint_D (x^2+y^2)\,dA where DD is the disk x2+y24x^2+y^2 \le 4.

In polar: x2+y2=r2x^2+y^2 = r^2, DD corresponds to 0r20 \le r \le 2, 0θ2π0 \le \theta \le 2\pi.

02π ⁣02r2rdrdθ=02π ⁣02r3drdθ=02πr4402dθ=02π4dθ=8π.\int_0^{2\pi}\!\int_0^2 r^2 \cdot r\,dr\,d\theta = \int_0^{2\pi}\!\int_0^2 r^3\,dr\,d\theta = \int_0^{2\pi} \frac{r^4}{4}\Bigr|_0^2\,d\theta = \int_0^{2\pi} 4\,d\theta = 8\pi.

The polar formula is a special case of the general change-of-variables theorem for double integrals.

Suppose (x,y)=T(u,v)(x,y) = \mathbf{T}(u,v) is a smooth, one-to-one transformation from a region SS in the (u,v)(u,v)-plane to a region DD in the (x,y)(x,y)-plane. The Jacobian of T\mathbf{T} is the determinant

J=(x,y)(u,v)=det(xuxvyuyv).J = \frac{\partial(x,y)}{\partial(u,v)} = \det\begin{pmatrix} x_u & x_v \\ y_u & y_v \end{pmatrix}.

The change-of-variables formula is

Df(x,y)dA=Sf(T(u,v))Jdudv.\iint_D f(x,y)\,dA = \iint_S f(\mathbf{T}(u,v))\,|J|\,du\,dv.

For polar coordinates x=rcosθx = r\cos\theta, y=rsinθy = r\sin\theta:

J=det(cosθrsinθsinθrcosθ)=rcos2θ+rsin2θ=r.J = \det\begin{pmatrix} \cos\theta & -r\sin\theta \\ \sin\theta & r\cos\theta \end{pmatrix} = r\cos^2\theta + r\sin^2\theta = r.

So J=r|J| = r, giving dA=rdrdθdA = r\,dr\,d\theta. This confirms the geometric argument above.


5. General polar domains

We are not limited to polar rectangles. A general polar domain is described by

D={(r,θ):αθβ,    h1(θ)rh2(θ)},D = \{\,(r,\theta) : \alpha \le \theta \le \beta,\;\; h_1(\theta) \le r \le h_2(\theta)\,\},

where h1h_1 and h2h_2 are continuous and 0h1(θ)h2(θ)0 \le h_1(\theta) \le h_2(\theta).

This is the polar analogue of a Type I region (sweep θ\theta from α\alpha to β\beta; for each θ\theta, rr runs from h1(θ)h_1(\theta) to h2(θ)h_2(\theta)).

Df(x,y)dA=αβ ⁣h1(θ)h2(θ)f(rcosθ,rsinθ)  rdrdθ.\boxed{\iint_D f(x,y)\,dA = \int_\alpha^\beta\!\int_{h_1(\theta)}^{h_2(\theta)} f(r\cos\theta, r\sin\theta)\; r\,dr\,d\theta.}

Example. Compute the area of one petal of the rose r=cos(2θ)r = \cos(2\theta).

The first petal lies in π/4θπ/4-\pi/4 \le \theta \le \pi/4 where cos(2θ)0\cos(2\theta) \ge 0. So h1=0h_1 = 0, h2=cos(2θ)h_2 = \cos(2\theta):

Area=π/4π/4 ⁣0cos(2θ)rdrdθ=π/4π/4cos2(2θ)2dθ=12π/4π/41+cos(4θ)2dθ=π8.\text{Area} = \int_{-\pi/4}^{\pi/4}\!\int_0^{\cos(2\theta)} r\,dr\,d\theta = \int_{-\pi/4}^{\pi/4} \frac{\cos^2(2\theta)}{2}\,d\theta = \frac{1}{2}\int_{-\pi/4}^{\pi/4}\frac{1+\cos(4\theta)}{2}\,d\theta = \frac{\pi}{8}.


6. When to use polar coordinates

Polar coordinates simplify a double integral when:

  • The domain is a disk, annulus, sector, or described by an equation like x2+y2=R2x^2+y^2 = R^2.
  • The integrand contains x2+y2x^2+y^2, x2+y2\sqrt{x^2+y^2}, arctan(y/x)\arctan(y/x), or similar combinations.
  • The domain has a circular boundary.

The classic example is the Gaussian integral: ex2dx=π\int_{-\infty}^\infty e^{-x^2}\,dx = \sqrt{\pi}.

Let I=ex2dxI = \int_{-\infty}^\infty e^{-x^2}\,dx. Then

I2=(ex2dx)(ey2dy)=R2e(x2+y2)dA.I^2 = \Bigl(\int_{-\infty}^\infty e^{-x^2}\,dx\Bigr)\Bigl(\int_{-\infty}^\infty e^{-y^2}\,dy\Bigr) = \iint_{\mathbb{R}^2} e^{-(x^2+y^2)}\,dA.

In polar coordinates x2+y2=r2x^2+y^2 = r^2 and the entire plane is 0r<0 \le r < \infty, 0θ<2π0 \le \theta < 2\pi:

I2=02π ⁣0er2rdrdθ=2π0rer2dr=2π12=π.I^2 = \int_0^{2\pi}\!\int_0^\infty e^{-r^2}\,r\,dr\,d\theta = 2\pi \int_0^\infty r\,e^{-r^2}\,dr = 2\pi \cdot \frac{1}{2} = \pi.

Therefore I=πI = \sqrt{\pi}.


7. Exercises

Use a polar integral to find the area of the disk x2+y2R2x^2+y^2 \le R^2.

Solution.

Area=02π ⁣0Rrdrdθ=02πR22dθ=πR2.\text{Area} = \int_0^{2\pi}\!\int_0^R r\,dr\,d\theta = \int_0^{2\pi}\frac{R^2}{2}\,d\theta = \pi R^2. \checkmark

Evaluate Rx2+y2dA\displaystyle\iint_R \sqrt{x^2+y^2}\,dA where RR is the annular sector 1r31 \le r \le 3, 0θπ/20 \le \theta \le \pi/2.

Solution.

x2+y2=r\sqrt{x^2+y^2} = r in polar:

0π/2 ⁣13rrdrdθ=π213r2dr=π2263=13π3.\int_0^{\pi/2}\!\int_1^3 r \cdot r\,dr\,d\theta = \frac{\pi}{2} \int_1^3 r^2\,dr = \frac{\pi}{2}\cdot\frac{26}{3} = \frac{13\pi}{3}.

Evaluate 22 ⁣4x24x2ex2+y2dydx\displaystyle\int_{-2}^{2}\!\int_{-\sqrt{4-x^2}}^{\sqrt{4-x^2}} e^{x^2+y^2}\,dy\,dx.

Solution.

The domain is the disk x2+y24x^2+y^2 \le 4 (radius 2). In polar:

02π ⁣02er2rdrdθ=2π02rer2dr=2π[er22]02=π(e41).\int_0^{2\pi}\!\int_0^2 e^{r^2}\,r\,dr\,d\theta = 2\pi \int_0^2 r\,e^{r^2}\,dr = 2\pi \cdot \Bigl[\tfrac{e^{r^2}}{2}\Bigr]_0^2 = \pi(e^4 - 1).

Find the area enclosed by the cardioid r=1+cosθr = 1 + \cos\theta.

Solution.

The cardioid traces a full loop for 0θ2π0 \le \theta \le 2\pi, with rr going from 00 to 1+cosθ1+\cos\theta:

Area=02π ⁣01+cosθrdrdθ=02π(1+cosθ)22dθ.\text{Area} = \int_0^{2\pi}\!\int_0^{1+\cos\theta} r\,dr\,d\theta = \int_0^{2\pi} \frac{(1+\cos\theta)^2}{2}\,d\theta.

Expand: (1+cosθ)2=1+2cosθ+cos2θ(1+\cos\theta)^2 = 1 + 2\cos\theta + \cos^2\theta. Integrating over [0,2π][0,2\pi]: cosθ=0\int\cos\theta = 0 and cos2θ=π\int\cos^2\theta = \pi.

=12(2π+0+π)=3π2.= \frac{1}{2}(2\pi + 0 + \pi) = \frac{3\pi}{2}.

Find the volume of the solid below z=R2x2y2z = \sqrt{R^2 - x^2 - y^2} and above the disk x2+y2R2x^2+y^2 \le R^2.

Solution.

V=x2+y2R2R2x2y2dA=02π ⁣0RR2r2rdrdθ.V = \iint_{x^2+y^2\le R^2} \sqrt{R^2-x^2-y^2}\,dA = \int_0^{2\pi}\!\int_0^R \sqrt{R^2-r^2}\,r\,dr\,d\theta.

=2π0RrR2r2dr.= 2\pi \int_0^R r\sqrt{R^2-r^2}\,dr.

Let u=R2r2u = R^2-r^2, du=2rdrdu = -2r\,dr:

=2πR20udu2=π0R2u1/2du=π23R3=2πR33.= 2\pi \int_{R^2}^0 \sqrt{u}\,\frac{-du}{2} = \pi\int_0^{R^2} u^{1/2}\,du = \pi\cdot\frac{2}{3}R^3 = \frac{2\pi R^3}{3}.

This is half the volume of a sphere of radius RR, as expected.