Double Integrals in Polar Coordinates

Math 215 Examples Double Integrals in Polar Coordinates Key Concepts Polar Rectangles A polar rectangle is a region in the xy-plane defined by the inequalities a ≤ r ≤ b and α ≤ θ ≤ β in polar coordinates. For example, the unit disk can be concisely described as the polar rectangle 0 ≤ r ≤ 1 , 0 ≤ θ ≤ 2π . The area of a polar rectangle is 12 (β − α)(b2 − a2 ) , as this is the difference in area between the sectors of a radius b and a radius a circle for α ≤ θ ≤ β . In particular, if we have a polar rectangle of radial "width" Δr and angular "width" Δθ centered around (r , θ) , then β − α = Δθ and 12 (b2 − a2 ) = 12 (b + a)(b − a) = r Δr , so the area of this polar rectangle is r Δr Δθ. Integrating in Polar Coordinates Recall that, to estimate ∬ f (x, y) dA over an ordinary rectangle R , we formed the Riemann sum R ∬ f x y dA ∑ ∑ f x y m n ( , ) ≈ R ( i=1 j=1 i , j )ΔA = ∑∑f x y m n i=1 j=1 ( i , i )ΔxΔy where each (xi , yj ) is a sample point in the ij th subrectangle. Taking the limit as the number of subrectangles goes to infinity is, in fact, our definition of ∬ f (x, y) dA . R What if instead we want to integrate f (x, y) over a polar rectangle R ? We can write f (x, y) in polar coordinates as f (r cos θ, r sin θ) by using the relations x = r cos θ , y = r sin θ . Now if we subdivide R into polar subrectangles, we have ∬ f x y dA ∑ ∑ f r m n ( , ) ≈ R ( i=1 j=1 i cos θj , ri sin θj )ΔA where each (ri , θj ) is a sample point in the ij th polar subrectangle. For example, a Riemann sum over a polar rectangle with 4 radial subdivisions and 3 angular subdivisions might look graphically like If we pick the sample points in the Riemann sum to be the centers of the polar subrectangles, then as just discussed, we can write ΔA = ri ΔrΔθ , and our Riemann sum becomes ∬ f x y dA ∑ ∑ f r m n ( , ) ≈ R ( i=1 j=1 i cos θj , ri sin θj )ri ΔrΔθ (notice the ri arising from ΔA on the right hand side!). Taking the limit as m and n go to infinity turns the right hand side into an iterated integral: Double Integrals in Polar Coordinates Let R be the polar rectangle defined by 0 ≤ a ≤ r ≤ b and α ≤ θ ≤ β , with 0 ≤ β − α ≤ 2π . Then ∬ f x y dA ∫ ∫ ( , ) R = β α a b f (r cos θ, r sin θ) rdrdθ = ∫∫ b a α β f (r cos θ, r sin θ) rdθ dr. Illustrated Example Example Evaluate ∬x ( 2 R x y ) dA −3 + 2 where R is the sector of the unit disk lying in the second quadrant of the xy-plane. Worked Solution The region R is the polar rectangle 0 ≤ r ≤ 1 , π/2 ≤ θ ≤ π . Using x = r cos θ and y = r sin θ , we write the integrand in polar coordinates as x 2 x y −3 + 2 r θ r θ r θ r (cos θ + sin θ) − 3r cos θ = r − 3r cos θ. = ( cos )2 − 3 cos + ( sin )2 = 2 2 2 2 Remembering to include the extra factor of r when converting to polar coordinates, the desired integral is ∬ R x ( 2 π x y ) dA = ∫ −3 + 2 π/2 π ∫ 1 r − 3 cos ) r −3 ( 0 1 r θ rdrdθ 2 r 2 cos ) θ drdθ r 3 cos ∫ ∫ = ∫ (( 14 r = ∫ ( 14 − cos θ) dθ = π/2 π π/2 π π/2 0 ( 3 4 − 1 θ) ∣∣ ) dθ ∣0 π ( 14 θ − sin θ) ∣∣∣ π π π = ( + 0) − ( − 1) 8 π 4+ 1. = = /2 8 Visualizing the Example The following animation shows the polar Riemann sums approximating this double integral as the number of subdivisions increases. Notice that the polar rectangles closer to the origin are much narrower looking than the ones further out, so if we had two boxes in a polar Riemann sum with the same height, the one closer to the origin would contribute less to the result than the one further out. This is not true in the ordinary, non-polar Riemann sums we've looked at; in these sums, all the subrectangles have the same area ΔA = Δx Δy , so two boxes with the same height have the same volume and hence always contribute the same amount to the Riemann sum. However, in polar Riemann sums, the area of a polar subrectangle is ΔA = r Δr Δθ , which depends also on r , the distance from the origin. Thus polar subrectangles closer to the origin (with small r ) contribute less to the result than polar subrectangles further from the origin (with bigger r ). We see this graphically in the narrow rectangles near the origin, and symbolically in the extra factor of r that shows up when writing the double integral as an iterated integral in polar coordinates. Further Questions 1. Work this example again using the other order of integrals, integrating first with respect to θ then r. 2. How would our answer change if our region of integration were instead the sector of the unit disk in the first quadrant? What about the third quadrant? Fourth quadrant? Try to answer these questions without redoing the whole calculation; instead, think about which parts of the calculation would change, and how. 3. Based on your answers to Question 2, what would we get if we integrated x2 − 3x + y 2 over the entire unit disk? 4. In the box on Double Integrals in Polar Coordinates, we defined the angular range as α ≤ θ ≤ β , with 0 ≤ β − α ≤ 2π . What could go wrong with our polar rectangles if we allow β − α > 2π ? 5. In the box on Double Integrals in Polar Coordinates, we defined the radial range as 0 ≤ a ≤ r ≤ b . What could go wrong with our polar rectangles if we allow a < 0 ? Using the Mathematica Demo All graphics on this page were generated by the Mathematica notebook 15_4DoubleIntegralsInPolarCoordinates.nb. This notebook generates images and animations like those on this page of polar Riemann sums for any integrand f (r , θ) and any polar rectangle. As an exercise, use the notebook to provide clear graphical answers to Questions 2 and 3 above. Then, can you come up with an integrand f (r , θ) and bounds so that the notebook produces Riemann sums approximating the area of the unit circle, or the volume of the unit sphere? What about a cone of radius 1 and height 1?