Because we choose to cover polar coordinates in the chapter on the precalculus of multivariable functions, this section jumps right into how polar coordinates work in double integrals and provide some computational advantages. If you did not cover SectionΒ 9.9 when teaching ChapterΒ 9, you will need to cover at least the subsection on polar coordinates before teaching this section. The preview activity in this section helps students recall important elements of points, graphs, and regions in polar coordinates before we apply the same slicing and iterated integral techniques as in SectionΒ 12.4.
This section has only two activities (plus the preview activity). For instructors who introduced polar coordinates in ChapterΒ 9, it will be easy to cover this section in a single class meeting. Instead of having more activities in this section, we encourage instructors to reinforce that cylindrical coordinates are βpolar coordinates with \(z\)β so that the section on triple integrals in cylindrical coordinates can also help students build skills in this area.
Suppose that we would like to integrate \(f(x,y) = e^{x^2 + y^2}\) over the the unit disc \(D\) centered at the origin. Using the methods we have learned thus far tin the chapter, we could evaluate this double integral using the iterated integral
A two-dimensional plot of a shaded region in the \(xy\)-plane. The region is a disc of radius \(1\) centered at the origin. The upper semicircle bounding the disc is labeled \(y=\sqrt{1-x^2}\text{.}\) The lower semicircle bounding the disc is labeled \(y=-\sqrt{1-x^2}\text{.}\) There are vertical dashed lines throughout the region.
For this particular integral, we are unable to find an antiderivative of the integrand. Changing the order of integration does not help, as the region and the integrand are symmetric in \(x\) and \(y\text{.}\) Furthermore, even if we could find an antiderivative, the inner limits of integration involve relatively complicated functions. However, the region of integration has rotational symmetry about the origin. This should make you think about polar coordinates, which we introduced in SubsectionΒ 9.9.2. In this section, we will see how we can use polar coordinates to write double integrals as iterated integrals, allowing us to both streamline the writing of the limits of integration and make evaluating the iterated integral easier.
Recall that a point \(P\) in rectangular coordinates that is described by an ordered pair \((x,y)\text{,}\) where \(x\) is the displacement from \(P\) to the \(y\)-axis and \(y\) is the displacement from \(P\) to the \(x\)-axis can also be described with polar coordinates \((r,\theta)\text{,}\) where \(r\) is the distance from \(P\) to the origin and \(\theta\) is the angle formed by the line segment \(\overline{OP}\) and the positive \(x\)-axis, as shown in FigureΒ 12.6.2.
A two dimensional plot with the horizontal axis labeled \(x\text{,}\) the vertical axis labeled \(y\text{,}\) and a point in the first quadrant labeled \(P=(r,\theta)\text{.}\) There is a line segment going from the origin to \(P\text{.}\) The line segment is labeled \(r\text{.}\) An arc is shown going counterclockwise from the horizontal axis to the line segment and is labeled \(\theta\text{.}\) There is a dashed vertical line segment labeled \(y\) from the horizontal axis to \(P\text{.}\) There is a dashed horizontal line segment labeled \(x\) from the vertical axis to \(P\text{.}\)
Trigonometry and the Pythagorean Theorem allow for straightforward conversion from rectangular to polar, and vice versa. A reminder of these formulas, which were developed in SubsectionΒ 9.9.2, is included below.
The labels in the box above are useful when thinking about converting the coordinates of a point from one coordinate system to the other. As we think about iterated integrals, however, we will often need to take a function \(z=f(x,y)\) in which the inputs are described in rectangular coordinates and rewrite it using polar coordinates. This task simply requires using substitutions for \(x\) and \(y\) given in the first part of the conversions.
Recall that when considering polar coordinates, it is important to think about how the sign of \(r\) influences the position of a point. There are also multiple angle measurements in polar coordinates that can be used for a single point, so it will be helpful to think flexibly about ranges of \(\theta\) when describing regions in polar coordinates. Refer to ExampleΒ 9.9.12 and later examples for different ways to interpret and make conversions between rectangular and polar coordinates.
We can draw graphs of curves in polar coordinates in a similar way to how we do in rectangular coordinates. However, when plotting in polar coordinates, we use a grid that considers changes in angles and changes in distance from the origin. In particular, the angles \(\theta\) and distances \(r\) partition the plane into small wedges as shown in FigureΒ 12.6.3.
A two-dimensional plot of the polar coordinate grid. Concentric circles centered at the origin are shown with radius increasing in half unit increments. Rays emanante from the origin between the horizontal and vertical axes. In each quadrant, there are three rays between the axes. The middle of each of these rays is marked \(\pi/4\text{,}\)\(3\pi/4\text{,}\)\(5\pi/4\text{,}\) and \(7\pi/4\text{,}\) working counterclockwise starting in the first quadrant. The positive horizontal axis is marked \(0\text{,}\) the positive vertical axis is marked \(\pi/2\text{,}\) the negative horizontal axis is marked \(\pi\text{,}\) and the negative vertical axis is marked \(3\pi/2\text{.}\) There are tick marks for \(2\) and \(4\) along the positive portion of the horizontal axis where the circles with those radii cross it.
For each of the following polar equations, draw a plot in the plane of the curve. Write a sentence for each graph to explain why the equation is satisfied.
Sketch the region in the \(xy\)-plane with \(x\leq 0\text{,}\)\(y\leq 0\text{,}\) and inside the circle \(x^2+y^2=25\text{.}\) State inequalities for \(r\) and \(\theta\) that describe the region.
Sketch the region in the \(xy\)-plane that lies between the circles \(x^2+y^2=4\) and \(x^2+y^2=9\text{.}\) State inequalities for \(r\) and \(\theta\) that describe the region.
\begin{equation*}
\iint_D e^{x^2+y^2} \, dA
\end{equation*}
where \(D\) is the unit disk using polar coordinates, we will find a way to write this double integral as an iterated integral in the variables \(r\) and \(\theta\text{.}\) We will need to complete three steps to transform this double integral into polar coordinates:
The second and third of these steps are already things we have discussed in the context of understanding polar coordinates, and we will return to them for this specific integral shortly. However, understanding how to write the area element \(dA\) in polar coordinates is a new topic. To do this, we focus on how the area element \(dA\) can be written in terms of \(dr\) and \(d\theta\) in the context of polar coordinates.
To write \(dA\) in polar coordinates, we start by remembering what \(dA\) represents. It comes from \(\Delta A\text{,}\) which is the area of the base of a rectangular prism used to approximate a volume. In rectangular coordinates, \(\Delta A = \Delta x\, \Delta y\text{,}\) which leads to \(dA = dx\, dy\text{.}\) To think in polar coordinates, we need to consider the area of a polar rectangle with \(R_0\leq r\leq R_1 \) and \(\theta_0 \leq \theta\leq \theta_1\text{.}\)FigureΒ 12.6.4 shows two polar rectangles with the same values of \(\theta_0\) and \(\theta_1\text{.}\) However, the figure illustrates that a polar rectangle located farther from the origin has greater area than one located closer to the origin, even when \(\Delta\theta = \theta_1 - \theta_0\) and \(\Delta r = R_1 - R_0\) are the same for the two polar rectangles. This suggests that \(dA\) depends on not only \(dr\) and \(d\theta\) but also \(r\text{.}\)
A polar grid in the first quadrant of the \(xy\)-plane. Two polar regions are shown. They are bounded by the same two rays emanating from the origin. The straight sides of the regions bounded by these rays have the same lengths. The other sides of the regions are arcs of circles centered at the origin. The distance between the two circles that determine the region closer to the origin is the same as the distance between the two circles that determine the region farther from the origin. However, the region farther from the origin has area significantly larger than the one closer to the origin.
Before stating and justifying the formula for \(dA\) in polar coordinates, itβs worth also thinking about what the units on \(dA\) are. They need to be units of area, so suppose that \(x\) and \(y\) are both measured in centimeters. Then \(dA = dx\, dy\) in rectangular coordinates has units square centimeters, which is an area unit. However, in polar coordinates, \(\theta\) is measured in radians. Conventionally, we treat quantities measured in radians as lacking units, so \(d\theta\, dr\) would have units of centimeters, which is a unit of length, not area.
Remember that the Classic Calculus Approach to defining the double integral involved breaking each coordinate into smaller pieces which created a collection of smaller (rectangular) regions that we used to create our Riemann sum. We can do the same procedure in terms of polar coordinates, and we could have our smaller pieces be in terms of some small changes in \(r\) and \(\theta\text{,}\) such as the regions on the polar grid in FigureΒ 12.6.3. Notice however that, as we showed in FigureΒ 12.6.4, even though the changes in \(r\) and \(\theta\) used to define the grid are all the same, a region farther from the origin has larger area than one closer to the origin. This suggests that understanding how to measure the area of a piece described by \(R_0\leq r \leq R_1\) and \(\theta_0 \leq \theta \leq \theta_1\text{,}\) as shown in FigureΒ 12.6.7, will be more complicated than in rectangular coordinates.
We first consider the area of the annulus between \(r=R_0\) and \(r=R_{1}\text{,}\) shown in FigureΒ 12.6.6. The area of this region is the difference between the area of a disk of radius \(R_1\) and the area of a disk of radius \(R_0\text{:}\)
A plot of two circles in the plane. Both circles are centered at the origin. The circles have radii marked as \(R_0\) and \(R_1\text{,}\) with \(R_0 \lt R_1\text{.}\) The region between the circles is shaded.
We need to measure the area of the portion of the annulus satisfying \(\theta_0 \leq \theta \leq \theta_1\text{.}\) Let \(\Delta \theta=\theta_1-\theta_0\) be the change in \(\theta\) for the small piece for which we want to find the area. Because we are using radian measures for angles, the fraction of the annulus we are considering is \(\frac{\Delta \theta}{2\pi}\text{.}\) This means that the area of this small piece is
\begin{equation*}
\Delta A = \frac{\Delta \theta}{2\pi} \left(\pi\left(R_1^2-R_0^2\right)\right)\text{.}
\end{equation*}
A two-dimensional plot in the first quadrant of a region in polar coordinates. The diagram includes two dashed circular arcs with radii labeled as \(R_0\) and \(R_1\) and \(R_0\lt R_1\text{.}\) There are also two rays emanating from the origin. The measures of the angles formed by these rays and the horizontal axis are marked as \(\theta_0\) and \(\theta_1\) with \(\theta_0 \lt \theta_1\text{.}\) The angle between the rays is labled \(\Delta \theta\text{.}\) The region between the circular arcs and rays is shaded and labeled \(\Delta A\text{.}\) The side of the shaded region lying on the circular arc of radius \(R_0\) is labeled \(s_0\text{.}\) The side of the shaded region lying on the circular arc of radius \(R_1\) is labeled \(s_1\text{.}\) One side of the shaded region lying on a ray is labeled \(\Delta r\text{.}\)
where \(\Delta r = R_1-R_0\text{.}\) Note that as \(\Delta r \) gets smaller, \(R_1\) and \(R_0\) get closer to the same value, which means that
\begin{equation*}
\Delta A \approx R_1 \Delta r \, \Delta \theta\text{.}
\end{equation*}
Since we can easily convert the function \(f(x,y)\) to polar coordinates with the conversion equations to be \(f(r\cos(\theta),r\sin(\theta))\text{,}\) we can then approximate a double integral \(\iint_D f(x,y) \, dA\) with a Riemann sum of the form
Notice that here, instead of the \(\Delta x\, \Delta y\) that we had in a Riemann sum using rectangular coordinates, our integrand is multiplied by \(r_i^*\, \Delta r\, \Delta \theta\text{.}\) This leads us to the realization that we need to use \(dA = r\, dr\, d\theta\) when writing a double integral as an iterated integral in polar coordinates.
Returning briefly to our discussion of units, notice that \(dA = r\, dr\, d\theta\) has units of square centimeters if \(x\) and \(y\) (and thus \(r\)) are measured in centimeters. This is an area unit, as desired. Before getting into describing more complicated regions in polar coordinates, we return to the example of evaluating \(\iint_D e^{x^2+y^2}\, dA\) over the unit disc \(D\text{.}\)
Let \(f(x,y) = e^{x^2+y^2}\) on the disk \(D = \{(x,y) : x^2 + y^2 \leq 1\}\text{.}\) Our goal in this example is to evaluate \(\displaystyle \iint_D f(x,y) \, dA\text{.}\)
We cannot evaluate this iterated integral because \(e^{x^2 + y^2}\) does not have an elementary antiderivative with respect to either \(x\) or \(y\text{.}\) However, since \(r^2=x^2+y^2\) and the region \(D\) is circular, it is natural to wonder whether converting to polar coordinates will allow us to evaluate the new integral. To do so, we replace \(x\) with \(r \cos(\theta)\text{,}\)\(y\) with \(r \sin(\theta)\text{,}\) and \(dy \, dx\) with \(r \, dr \, d\theta\) to obtain
\begin{equation*}
\iint_D f(x,y) \, dA = \iint_D e^{r^2} \, r \, dr \, d\theta\text{.}
\end{equation*}
The disc \(D\) is described in polar coordinates by the inequalities \(0 \leq r \leq 1\) and \(0 \leq \theta \leq 2\pi\text{.}\) Therefore, it follows that
\begin{equation*}
\iint_D e^{r^2} \, r \, dr \, d\theta = \int_{\theta=0}^{\theta = 2\pi} \int_{r=0}^{r=1} e^{r^2} \, r \, dr \, d\theta\text{.}
\end{equation*}
We can evaluate the resulting iterated polar integral as follows:
We now focus on the second conversion step: convert the region of integration described by inequalities that correspond to fixing one coordinate and bounding the other in a consistent way. Similar to how we used vertically simple and horizontally simple descriptions for rectangular regions, we want to set up inequalities by holding one of our coordinates constant and setting boundary functions for the other coordinate.
A two-dimensional plot of a shaded region in the first quadrant. The region is bounded by two rays emanating from the origin as well as two wavy curves. The angles these rays make with the positive horizontal axis are labeled as \(a\) and \(b\) with \(a\lt b\text{.}\) The curve closer to the origin is labeled \(r=g_1(\theta)\text{.}\) The curver farther from the origin is labeled \(r=g_2(\theta)\text{.}\) There are dashed line segments drawn through the region in a manner that suggests they are the portion of rays emanating from the origin that lie inside the region.
A two-dimensional plot of a shaded region in the first quadrant. The region is bounded by two arcs of circles centered at the origin. The arc closer to the origin is labeled \(r=c\text{.}\) The arc farther from the origin is labeled \(r=d\text{.}\) The other boundaries of the region are wavy curves joining the circular arcs. The bounding curve with smaller \(\theta\) values is labeled \(\theta=h_1(r)\text{.}\) The other bounding curve is labeled \(\theta=h_2(r)\text{.}\) There are dashed circular arcs drawn through the region in a manner that suggests they are the portions of circles centered at the origin that lie inside the region.
Note that the unit disc \(D\) we have considered throughout this section can be considered both radially simple and angularly simple since all of the bounds of integration are constant. FigureΒ 12.6.12 shows how \(D\) can be considered radially simple.
A two-dimensional plot of a shaded region in the plane. The region is a disc of radius \(1\) centered at the origin. The circle bounding the disc is labeled \(r=1\text{.}\) The origin is labeled \(r=0\text{.}\) There are dashed rays emanating from the origin through the region.
Regions that are angularly simple but not radially simple do not arise with much frequency. Thus, while we have defined both types of simple regions for completeness, we will not spend time on an example of a region that is angularly simple but not radially simple.
Each of the types of regions above corresponds to exactly the necessary description to use an iterated integral. With a radially simple region, for any fixed value of \(\theta\) in the interval \([a,b]\text{,}\) a slice along this constant \(\theta\) value corresponds to a ray from the origin in the sector between the grey lines in FigureΒ 12.6.10. Along this ray with \(\theta=k\text{,}\) the region of interest has upper and lower bounds on the \(r\)-coordinate. In fact, \(r=g_2(\theta)\) describes the upper boundary for every value of \(\theta\) in \([a,b]\text{.}\) Similarly, \(r=g_1(\theta)\) describes the lower boundary for every value of \(\theta\) in \([a,b]\text{.}\) Geometrically, this means that along every gray dashed slice where \(\theta\) is constant, the region has the same upper and lower boundary curves, shown in red and blue, respectively.
While angularly simple regions are used less frequently than radially simple regions, a parallel structure is used. Along every slice with a constant value of \(r\) in \([c,d]\text{,}\) represented by the gray, dashed curves in FigureΒ 12.6.11, the upper bound on \(\theta\) is given by the function \(\theta=h_2(r)\) and the lower bound is given by \(\theta=h_1(r)\text{,}\) shown by the red and blue bounding curves in FigureΒ 12.6.11.
Just as in the rectangular coordinate case, especially complicated regions may need to be split into several different regions that are radially or angularly simple.
The double integral \(\iint_D f(x,y) \, dA\) in rectangular coordinates can be converted to a double integral in polar coordinates as \(\iint_D f(r\cos(\theta), r\sin(\theta)) \, r \, dr \, d\theta\text{.}\)
Note again that the order of integration corresponds to the manner in which we have sliced the region of integration. The variable being held constant along a slice is the last variable to be integrated, just as was the case with iterated integrals in rectangular coordinates.
Suppose that we have a piece of metal in the shape of a quarter annulus \(R\text{,}\) as shown in FigureΒ 12.6.14. The circular arcs bounding the region have radius \(2\) and \(3\text{,}\) and the region lies in the first quadrant. The metal is heated so that the temperature in degrees Celsius at the point \((x,y)\) is given by \(T(x,y) = \sqrt{x^2+y^2}\text{.}\) That is, the temperature is 12 times the distance between the origin and the point. In this example, we will find the average value of \(T\) over the region \(R\text{.}\)
A region in the first quadrant of the \(xy\)-plane. The region is labeled \(R\) and is bounded by the circles of radius \(2\) and \(3\text{,}\) both centered at the origin, as well as the positive \(x\)- and \(y\)-axes.
We can compute the area of \(R\) using geometry, as its area is one quarter the area of the annulus between the circles of radius \(2\) and \(3\) centered at the origin. That annulus has area \(\pi\cdot 3^2 - \pi\cdot 2^2 = 5\pi\text{,}\) so \(R\) has area \(5\pi/4\text{.}\) To compute the double integral, we will use an iterated integral. This will require finding inequalities that allow us to describe \(R\text{.}\) Notice first that \(R\) is neither horizontally simple nor vertically simple, so evaluating \(\iint_R T\, dA\) rectangular coordinates would require that we set up multiple iterated integrals. However, in polar coordinates, we see that \(R\) is both radially simple and angularly simple. That is, \(R\) is a polar rectangle, which we can describe as
Using this description of \(R\) in polar coordinates along with the facts that \(r^2 = x^2+y^2\) and \(dA = r\, dr, d\theta\text{,}\) we can evaluate the double integral:
While there is no firm rule for when polar coordinates can or should be used, they are a natural alternative anytime the domain of integration may be expressed simply in polar form and/or when the integrand involves expressions such as \(\sqrt{x^2 + y^2}\text{.}\)
In this activity, we will find the volume of the solid region \(V\) above the \(xy\)-plane, below the paraboloid \(z=f(x,y) = 25-x^2-y^2\text{,}\) and outside the cylinder \(x^2+y^2 = 1\text{.}\)
We will find the volume by computing a double integral \(\iint_R f(x,y)\, dA\) where \(R\) is an appropriate region in the \(xy\)-plane. Sketch the region \(R\text{.}\)
Write inequalities to describe \(R\) in polar coordinates. Write a couple of sentences explaining why it would be very difficult to describe \(R\) using rectangular coordinates.
The next activity asks you to find the volume between two surfaces in three-dimensional space. To do this, you will find it useful to think of the functions involved in polar coordinates from the very beginning.
Use polar coordinates to find the volume of the solid between the cone \(z=f(x,y)=\sqrt{x^2+y^2}\) and the paraboloid \(z=g(x,y)=\frac{1}{8}\left(x^2+y^2\right) + 2\text{.}\)
Perhaps the surfaces in ActivityΒ 12.6.3 brought thoughts of cylindrical coordinates to your mind. The next section will generalize the idea of a double integral to a triple integral, and after that, we will consider integration using cylindrical and spherical coordinates as well.
To convert the double integral \({\displaystyle\iint_D f(x,y) \, dA}\) to an iterated integral in polar coordinates, we substitute \(r \cos(\theta)\) for \(x\text{,}\)\(r \sin(\theta)\) for \(y\text{,}\) and \(r \, dr \, d\theta\) for \(dA\) to obtain the iterated integral (where the order of integration is determined by whether the region is described as radially or angularly simple)
\begin{equation*}
{\iint_D f(r\cos(\theta), r\sin(\theta)) \, r \, dr \, d\theta}
\end{equation*}
You can find many more problems involving the conversion of points, graphs, and region between rectangular and polar coordinates in SectionΒ 9.9. The exercises of this section will be focused on applications of these conversion to double integrals.
For each of the following, set up the integral of an arbitrary function \(f(x,y)\) over the region in whichever of rectangular or polar coordinates is most appropriate. (Use \(t\) for \(\theta\) in your expressions.)
Using polar coordinates, evaluate the integral which gives the area which lies in the first quadrant between the circles \(x^2 + y^2 = 36\) and \(x^2 - 6x + y^2 = 0\text{.}\)
(a) Graph \(r=1/(4\cos\theta)\) for \(-\pi/2\le\theta\le\pi/2\) and \(r=1\text{.}\) Then write an iterated integral in polar coordinates representing the area inside the curve \(r=1\) and to the right of \(r=1/(4\cos\theta)\text{.}\) (Use \(t\) for \(\theta\) in your work.)
Using polar coordinates, evaluate the integral \(\displaystyle \int \!\! \int_{R} \sin (x^2+y^2) dA\) where R is the region \(4 \leq x^2 + y^2 \leq 81\text{.}\)
Instructions: Please enter the integrand in the first answer box, typing theta for \(\theta\text{.}\) Depending on the order of integration you choose, enter dr and dtheta in either order into the second and third answer boxes with only one dr or dtheta in each box. Then, enter the limits of integration.
Using polar coordinates, evaluate the integral \(\displaystyle \int \!\! \int_{R} \sin (x^2+y^2) dA\) where R is the region \(16 \leq x^2 + y^2 \leq 25\text{.}\)
A cylindrical drill with radius 1 is used to bore a hole through the center of a sphere of radius 6. Find the volume of the ring shaped solid that remains.
Evaluate the double integral \(\displaystyle \iint_D x^{2}y \, dA,\) where \(D\) is the top half of the disc with center the origin and radius \(4,\) by changing to polar coordinates.
Evaluate the double integral \(\displaystyle \iint_R (2 x - y) \, dA,\) where \(R\) is the region in the first quadrant enclosed by the circle \(x^2+y^2=4\) and the lines \(x=0\) and \(y=x,\) by changing to polar coordinates.
A volcano fills the volume between the graphs \(z = 0\) and \(\displaystyle z =
\frac{1}{\left( x^{2} + y^{2} \right)^{21}}\text{,}\) and outside the cylinder \(x^{2} + y^{2} = 1\text{.}\) Find the volume of this volcano.
Find the exact center of mass of the lamina over the portion of \(D\) that lies in the first quadrant and has its mass density distribution given by \(\delta(x,y) = 1\text{.}\) (Before making any calculations, where do you expect the center of mass to lie? Why?)
Determine a polar curve in the form \(r = f(\theta)\) that traces out the circle \(x^2 + (y-1)^2 = 1\text{.}\) To do this, you should substitute \(x = r \cos(\theta)\) and \(y = r \sin(\theta)\) into the rectangular coordinate equation and solve for \(r\text{.}\)
Find the volume under the surface \(h(x,y) = x\) over the region \(D\text{,}\) where \(D\) is the region bounded above by the line \(y=x\) and below by the circle. We have provided a sketch of this as the shaded region in FigureΒ 12.6.16.
A shaded region in the \(xy\)-plane. The region is bounded by part of a circle and the portion of the line \(y=x\) inside the circle. The circle passes through the origin and the point \((0,2)\text{.}\) The circle is also symmetric about the \(y\)-axis.