9.11 Double Integrals in Polar Coordinates
INTRODUCTION
A double integral, which may be difficult or even impossible to evaluate in rectangular xy-coordinates, may become more tractable when expressed in a different coordinate system. In this section we examine double integrals in polar rθ-coordinates.
Polar Rectangles
Suppose R is a region in the plane bounded by the graphs of the polar equations r = g1(θ), r = g2(θ), and the rays θ = α, θ = β, and f is a function of r and θ that is continuous on R. In order to define the double integral of f over R, we use rays and concentric circles to partition the region into a grid of polar rectangles or subregions Rk. See FIGURE 9.11.1(a) and 9.11.1(b). The area ΔAk of typical subregion Rk, shown in Figure 9.11.1(c), is the difference of the areas of two circular sectors: ΔAk = Δθk – Δθk. Now ΔAk can be written as
where Δrk = rk+1 – rk and = (rk+1 + rk) denotes the average radius. By choosing (, ) on each Rk, the double integral of f over R is
The double integral is then evaluated by means of the iterated integral:
.(1)
On the other hand, if the region R is as given in FIGURE 9.11.2, the double integral of f over R is then
.(2)
EXAMPLE 1 Center of Mass
Find the center of mass of the lamina that corresponds to the region bounded by one leaf of the rose r = 2 sin 2θ in the first quadrant if the density at a point P in the lamina is directly proportional to the distance from the pole.
SOLUTION
By varying θ from 0 to π/2, we obtain the graph in FIGURE 9.11.3. Now, d(O, P) = |r|. Hence, ρ(r, θ) = k|r|, where k is a constant of proportionality. From (9) of Section 9.10, we have
Since x = r cos θ, we can write as
Similarly, by using x = r sin θ, we find*
Here the rectangular coordinates of the center of mass are
≡
Change of Variables: Rectangular to Polar Coordinates
In some instances a double integral f(x, y) dA that is difficult or even impossible to evaluate using rectangular coordinates may be readily evaluated when a change of variables is used. If we assume that f is continuous on the region R and if R can be described in polar coordinates as 0 ≤ g1(θ) ≤ r ≤ g2(θ), α ≤ θ ≤ β, 0 < β – α ≤ 2π, then
. (3)
Equation (3) is particularly useful when f contains the expression x2 + y2 since, in polar coordinates, we can now write
EXAMPLE 2 Changing an Integral to Polar Coordinates
Use polar coordinates to evaluate
SOLUTION
From x ≤ y ≤ , 0 ≤ x ≤ 2, we have sketched the region R of integration in FIGURE 9.11.4. Since x2 + y2 = r2, the polar description of the circle x2 + y2 = 8 is r = . Hence, in polar coordinates, the region of R is given by 0 ≤ r ≤ , π/4 ≤ θ ≤ π/2. With 1/(5 + x2 + y2) = 1/(5 + r2), the original integral becomes
EXAMPLE 3 Volume
Find the volume of the solid that is under the hemisphere z = and above the region bounded by the graph of the circle x2 + y2 – y = 0.
SOLUTION
From FIGURE 9.11.5, we see that V = dA. In polar coordinates the equations of the hemisphere and the circle become, respectively, z = and r = sin θ. Now, using symmetry, we have
Area
Note that in (1) if f(r, θ) = 1, then the area of the region R in Figure 9.11.1(a) is given by
The same observation holds for (2) and Figure 9.11.2 when f(r, θ) = 1.
REMARKS
The reader is invited to reexamine Example 3. The graph of the circle r = sin θ is obtained by varying θ from 0 to π. However, carry out the iterated integration
and see if the result is the incorrect answer π/3. What goes wrong?
9.11 Exercises Answers to selected odd-numbered problems begin on page ANS-25.
In Problems 1–4, use a double integral in polar coordinates to find the area of the region bounded by the graphs of the given polar equations.
- r = 3 + 3 sin θ
- r = 2 + cos θ
- r = 2 sin θ, r = 1, common area
- r = 8 sin 4θ, one petal
In Problems 5–10, find the volume of the solid bounded by the graphs of the given equations.
- One petal of r = 5 cos 3θ, z = 0, z = 4
- x2 + y2 = 4, z =
- Between x2 + y2 = 1 and x2 + y2 = 9,
z = - z = + y2 = 25, z = 0
- r = 1 + cos θ, z = y, z = 0, first octant
- r = cos θ, z = 2 + x2 + y2, z = 0
In Problems 11–16, find the center of mass of the lamina that has the given shape and density.
- r = 1, r = 3, x = 0, y = 0, first quadrant; ρ(r, θ) = k (constant)
- r = cos θ; density at point P directly proportional to the distance from the pole
- y = x, y = 0, x = 3; ρ(r, θ) = r2
- r = 4 cos 2θ, petal on the polar axis; ρ(r, θ) = k (constant)
- Outside r = 2 and inside r = 2 + 2 cos θ, y = 0, first quadrant; density at a point P inversely proportional to the distance from the pole
- r = 2 + 2 cos θ, y = 0, first and second quadrants; ρ(r, θ) = k (constant)
In Problems 17–20, find the indicated moment of inertia of the lamina that has the given shape and density.
- r = a; ρ(r, θ) = k (constant); Ix
- r = a; ρ(r, θ) = ; Ix
- Outside r = a and inside r = 2a cos θ; density at a point P inversely proportional to the cube of the distance from the pole; Iy
- Outside r = 1 and inside r = 2 sin 2θ, first quadrant; ρ(r, θ) = sec2θ; Iy
In Problems 21–24, find the polar moment of inertia
of the lamina that has the given shape and density.
- r = a; ρ(r, θ) = k (constant) [Hint: Use Problem 17 and the fact that Ix = Iy.]
- r = θ, 0 ≤ θ ≤ π, y = 0; density at a point P proportional to the distance from the pole
- rθ = 1, ≤ θ ≤ 1, r = 1, r = 3, y = 0; density at a point P inversely proportional to the distance from the pole [Hint: Integrate first with respect to θ.]
- r = 2a cos θ; ρ(r, θ) = k (constant)
In Problems 25–32, evaluate the given iterated integral by changing to polar coordinates.
- The liquid hydrogen tank in the space shuttle has the form of a right circular cylinder with a semi-ellipsoidal cap at each end. The radius of the cylindrical part of the tank is 4.2 m. Find the volume of the tank shown in FIGURE 9.11.6.
- Evaluate (x + y) dA over the region shown in FIGURE 9.11.7.
Discussion Problems
- The improper integral dx is important in the theory of probability, statistics, and other areas of applied mathematics. If I denotes the integral, then
Discuss how to use polar coordinates to evaluate the last integral. Find the value of I.
* We could have argued to the fact that Mx = My and hence = from the fact that the lamina and the density function are symmetric about the ray θ = π/4.