Double Integrals in Polar Coordinates

Section 11.5 Double Integrals in Polar Coordinates

Motivating Questions

What are the polar coordinates of a point in two-space?
🔗

🔗
How do we convert between polar coordinates and rectangular coordinates?
🔗

🔗
What is the area element in polar coordinates?
🔗

🔗
How do we convert a double integral in rectangular coordinates to a double integral in polar coordinates?
🔗

🔗

While we have naturally defined double integrals in the rectangular coordinate system, starting with domains that are rectangular regions, there are many of these integrals that are difficult, if not impossible, to evaluate. For example, consider the domain

D

that is the unit circle and

f (x, y) = e^{- x^{2} - y^{2}} .

To integrate

f

over

D,

we would use the iterated integral

\iint_{D} f (x, y) d A = \int_{x = - 1}^{x = 1} \int_{y = - \sqrt{1 - x^{2}}}^{y = \sqrt{1 - x^{2}}} e^{- x^{2} - y^{2}} d y d x .

🔗

For this particular integral, regardless of the order of integration, we are unable to find an antiderivative of the integrand; in addition, even if we were able to find an antiderivative, the inner limits of integration involve relatively complicated functions.

🔗

It is useful, therefore, to be able to translate to other coordinate systems where the limits of integration and evaluation of the involved integrals is simpler. In this section we provide a quick discussion of one such system — polar coordinates — and then introduce and investigate their ramifications for double integrals. The rectangular coordinate system allows us to consider domains and graphs relative to a rectangular grid. The polar coordinate system is an alternate coordinate system that allows us to consider domains less suited to rectangular coordinates, such as circles.

🔗

Preview Activity 11.5.1.

The coordinates of a point determine its location. In particular, the rectangular coordinates of a point

P

are given by an ordered pair

(x, y),

where

x

is the (signed) distance the point lies from the

y

-axis to

P

and

y

is the (signed) distance the point lies from the

x

-axis to

P .

In polar coordinates, we locate the point by considering the distance the point lies from the origin,

O = (0, 0),

and the angle the line segment from the origin to

P

forms with the positive

x

-axis.

🔗

Determine the rectangular coordinates of the following points:
🔗
1. The point $P$ that lies 1 unit from the origin on the positive $x$ -axis.
  🔗
  
  🔗
2. The point $Q$ that lies 2 units from the origin and such that $\overset{―}{O Q}$ makes an angle of $\frac{π}{2}$ with the positive $x$ -axis.
  🔗
  
  🔗
3. The point $R$ that lies 3 units from the origin such that $\overset{―}{O R}$ makes an angle of $\frac{2 π}{3}$ with the positive $x$ -axis.
  🔗
  
  🔗
🔗
Part (a) indicates that the two pieces of information completely determine the location of a point: either the traditional $(x, y)$ coordinates, or alternately, the distance $r$ from the point to the origin along with the angle $θ$ that the line through the origin and the point makes with the positive $x$ -axis. We write “ $(r, θ)$ ” to denote the point’s location in its polar coordinate representation. Find polar coordinates for the points with the given rectangular coordinates.
1. $(0, - 1)$ ii. $(- 2, 0)$ iii. $(- 1, 1)$
  🔗
  
  🔗
🔗
🔗
For each of the following points whose coordinates are given in polar form, determine the rectangular coordinates of the point.
1. $(5, \frac{π}{4})$ ii. $(2, \frac{5 π}{6})$ iii. $(\sqrt{3}, \frac{5 π}{3})$
  🔗
  
  🔗
🔗
🔗

🔗

Subsection 11.5.1 Polar Coordinates

The rectangular coordinate system is best suited for graphs and regions that are naturally considered over a rectangular grid. The polar coordinate system is an alternative that offers good options for functions and domains that have more circular characteristics. A point

P

in rectangular coordinates that is described by an ordered pair

(x, y),

where

x

is the displacement from

P

to the

y

-axis and

y

is the displacement from

P

to the

x

-axis, as seen in Preview Activity 11.5.1, can also be described with polar coordinates

(r, θ),

where

r

is the distance from

P

to the origin and

θ

is the angle formed by the line segment

\overset{―}{O P}

and the positive

x

-axis, as shown at left in Figure 11.5.1.

🔗

Figure 11.5.1. The polar coordinates of a point and the polar coordinate grid.
🔗

Trigonometry and the Pythagorean Theorem allow for straightforward conversion from rectangular to polar, and vice versa.

🔗

Converting between rectangular and polar coordinates.

Converting from rectangular to polar..

If we are given the rectangular coordinates $(x, y)$ of a point $P,$ then the polar coordinates $(r, θ)$ of $P$ satisfy

$r = \sqrt{x^{2} + y^{2}} and \tan (θ) = \frac{y}{x}, assuming x \neq 0.$

🔗

🔗
Converting from polar to rectangular..

If we are given the polar coordinates $(r, θ)$ of a point $P,$ then the rectangular coordinates $(x, y)$ of $P$ satisfy

$x = r \cos (θ) and y = r \sin (θ) .$

🔗

🔗

🔗

Note: The angle

θ

in the polar coordinates of a point is not unique. We could replace

θ

with

θ + 2 π

and still be at the same terminal point. In addition, the sign of

\tan (θ)

does not uniquely determine the quadrant in which

θ

lies, so we have to determine the value of

θ

from the location of the point. In other words, more care has to be paid when using polar coordinates than rectangular 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

θ

and distances

r

partition the plane into small wedges as shown at right in Figure 11.5.1.

🔗

Activity 11.5.2.

Most polar graphing devices can plot curves in polar coordinates of the form

r = f (θ) .

Use such a device to complete this activity.

Before plotting the polar curve $r = 1$ (where $θ$ can have any value), think about what shape it should have, in light of how $r$ is connected to $x$ and $y .$ Then use appropriate technology to draw the graph and test your intuition.
🔗

🔗
The equation $θ = 1$ does not define $r$ as a function of $θ,$ so we can’t graph this equation on many polar plotters. What do you think the graph of the polar curve $θ = 1$ looks like? Why?
🔗

🔗
Before plotting the polar curve $r = θ,$ what do you think the graph looks like? Why? Use technology to plot the curve and compare your intuition.
🔗

🔗
What does the region defined by $1 \leq r \leq 3$ (where $θ$ can have any value) look like? (Hint: Compare to your response from part (a).)
🔗

🔗
What does the region defined by $1 \leq r \leq 3$ and $π / 4 \leq θ \leq π / 2$ look like?
🔗

🔗
Consider the curve $r = \sin (θ) .$ For some values of $θ$ we will have $r < 0 .$ In these situations, we plot the point $(r, θ)$ as $(| r |, θ + π)$ (in other words, when $r < 0,$ we reflect the point through the origin). With that in mind, what do you think the graph of $r = \sin (θ)$ looks like? Plot this curve using technology and compare to your intuition.
🔗

🔗

🔗

Subsection 11.5.2 Integrating in Polar Coordinates

Consider the double integral

\iint_{D} e^{x^{2} + y^{2}} d A,

where

D

is the unit disk. While we cannot directly evaluate this integral in rectangular coordinates, a change to polar coordinates will convert it to one we can easily evaluate.

🔗

We have seen how to evaluate a double integral

\iint_{D} f (x, y) d A

as an iterated integral of the form

\int_{a}^{b} \int_{g_{1} (x)}^{g_{2} (x)} f (x, y) d y d x

in rectangular coordinates, because we know that

d A = d y d x

in rectangular coordinates. To make the change to polar coordinates, we not only need to represent the variables

x

and

y

in polar coordinates, but we also must understand how to write the area element,

d A,

in polar coordinates. That is, we must determine how the area element

d A

can be written in terms of

d r

and

d θ

in the context of polar coordinates. We address this question in the following activity.

🔗

Figure 11.5.2. Left: A polar rectangle. Right: An annulus.
🔗

Activity 11.5.3.

Consider a polar rectangle

R,

with

r

between

r_{i}

and

r_{i + 1}

and

θ

between

θ_{j}

and

θ_{j + 1}

as shown at left in Figure 11.5.2. Let

Δ r = r_{i + 1} - r_{i}

and

Δ θ = θ_{j + 1} - θ_{j} .

Let

Δ A

be the area of this region.

🔗

Explain why the area $Δ A$ in polar coordinates is not $Δ r Δ θ .$
🔗

🔗
Now find $Δ A$ by the following steps:
🔗
1. Find the area of the annulus (the washer-like region) between $r_{i}$ and $r_{i + 1},$ as shown at right in Figure 11.5.2. This area will be in terms of $r_{i}$ and $r_{i + 1} .$
  🔗
  
  🔗
2. Observe that the region $R$ is only a portion of the annulus, so the area $Δ A$ of $R$ is only a fraction of the area of the annulus. For instance, if $θ_{i + 1} - θ_{i}$ were $\frac{π}{4},$ then the resulting wedge would be
  
  $\frac{\frac{π}{4}}{2 π} = \frac{1}{8}$
  
  of the entire annulus. In this more general context, using the wedge between the two noted angles, what fraction of the area of the annulus is the area $Δ A ?$
  
  🔗
  
  🔗
3. Write an expression for $Δ A$ in terms of $r_{i},$ $r_{i + 1},$ $θ_{j},$ and $θ_{j + 1} .$
  🔗
  
  🔗
4. Finally, write the area $Δ A$ in terms of $r_{i},$ $r_{i + 1},$ $Δ r,$ and $Δ θ,$ where each quantity appears only once in the expression. (Hint: Think about how to factor a difference of squares.)
  🔗
  
  🔗
🔗
As we take the limit as $Δ r$ and $Δ θ$ go to 0, $Δ r$ becomes $d r,$ $Δ θ$ becomes $d θ,$ and $Δ A$ becomes $d A,$ the area element. Using your work in (iv), write $d A$ in terms of $r,$ $d r,$ and $d θ .$
🔗

🔗

🔗

From the result of Activity 11.5.3, we see when we convert an integral from rectangular coordinates to polar coordinates, we must not only convert

x

and

y

to being in terms of

r

and

θ,

but we also have to change the area element to

d A = r d r d θ

in polar coordinates. As we saw in Activity 11.5.3, the reason the additional factor of

r

in the polar area element is due to the fact that in polar coordinates, the cross sectional area element increases as

r

increases, while the cross sectional area element in rectangular coordinates is constant. So, given a double integral

\iint_{D} f (x, y) d A

in rectangular coordinates, to write a corresponding iterated integral in polar coordinates, we replace

x

with

r \cos (θ),

y

with

r \sin (θ)

and

d A

with

r d r d θ .

Of course, we need to describe the region

D

in polar coordinates as well. To summarize:

🔗

Double integrals in polar coordinates.

The double integral

\iint_{D} f (x, y) d A

in rectangular coordinates can be converted to a double integral in polar coordinates as

\iint_{D} f (r \cos (θ), r \sin (θ)) r d r d θ .

🔗

Example 11.5.3.

Let

f (x, y) = e^{x^{2} + y^{2}}

on the disk

D = {(x, y) : x^{2} + y^{2} \leq 1} .

We will evaluate

\iint_{D} f (x, y) d A .

🔗

In rectangular coordinates the double integral

\iint_{D} f (x, y) d A

can be written as the iterated integral

\iint_{D} f (x, y) d A = \int_{x = - 1}^{x = 1} \int_{y = - \sqrt{1 - x^{2}}}^{y = \sqrt{1 - x^{2}}} e^{x^{2} + y^{2}} d y d x .

🔗

We cannot evaluate this iterated integral, because

e^{x^{2} + y^{2}}

does not have an elementary antiderivative with respect to either

x

y .

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 (θ),

y

with

r \sin (θ),

and

d y d x

with

r d r d θ

to obtain

\iint_{D} f (x, y) d A = \iint_{D} e^{r^{2}} r d r d θ .

🔗

The disc

D

is described in polar coordinates by the constraints

0 \leq r \leq 1

and

0 \leq θ \leq 2 π .

Therefore, it follows that

\iint_{D} e^{r^{2}} r d r d θ = \int_{θ = 0}^{θ = 2 π} \int_{r = 0}^{r = 1} e^{r^{2}} r d r d θ .

🔗

We can evaluate the resulting iterated polar integral as follows:

\begin{aligned} \int_{θ = 0}^{θ = 2 π} \int_{r = 0}^{r = 1} e^{r^{2}} r d r d θ & = \int_{θ = 0}^{2 π} (\frac{1}{2} e^{r^{2}} |_{r = 0}^{r = 1}) d θ \\ = \frac{1}{2} \int_{θ = 0}^{θ = 2 π} (e - 1) d θ \\ = \frac{1}{2} (e - 1) \int_{θ = 0}^{θ = 2 π} d θ \\ = \frac{1}{2} (e - 1) [θ] |_{θ = 0}^{θ = 2 π} \\ = π (e - 1) . \end{aligned}

🔗

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}} .

🔗

Activity 11.5.4.

Let

f (x, y) = x + y

and

D = {(x, y) : x^{2} + y^{2} \leq 4} .

Sketch the region $D$ and then write the double integral of $f$ over $D$ as an iterated integral in rectangular coordinates.
🔗

🔗
Write the double integral of $f$ over $D$ as an iterated integral in polar coordinates.
🔗

🔗
Evaluate one of the iterated integrals. Why is the final value you found not surprising?
🔗

🔗

🔗

Activity 11.5.5.

Consider the circle given by

x^{2} + (y - 1)^{2} = 1

as shown in Figure 11.5.4.

Figure 11.5.4. The graphs of $y = x$ and $x^{2} + (y - 1)^{2} = 1,$ for use in Activity 11.5.5.
🔗

Determine a polar curve in the form $r = f (θ)$ that traces out the circle $x^{2} + (y - 1)^{2} = 1 .$ (Hint: Recall that a circle centered at the origin of radius $r$ can be described by the equations $x = r \cos (θ)$ and $y = r \sin (θ) .$ )
🔗

🔗
Find the exact average value of $g (x, y) = \sqrt{x^{2} + y^{2}}$ over the interior of the circle $x^{2} + (y - 1)^{2} = 1 .$
🔗

🔗
Find the volume under the surface $h (x, y) = x$ over the region $D,$ where $D$ is the region bounded above by the line $y = x$ and below by the circle (this is the shaded region in Figure 11.5.4).
🔗

🔗
Explain why in both (b) and (c) it is advantageous to use polar coordinates.
🔗

🔗

🔗

Subsection 11.5.3 Summary

The polar representation of a point $P$ is the ordered pair $(r, θ)$ where $r$ is the distance from the origin to $P$ and $θ$ is the angle the ray through the origin and $P$ makes with the positive $x$ -axis.
🔗

🔗
The polar coordinates $r$ and $θ$ of a point $(x, y)$ in rectangular coordinates satisfy

$r = \sqrt{x^{2} + y^{2}} and \tan (θ) = \frac{y}{x};$

the rectangular coordinates $x$ and $y$ of a point $(r, θ)$ in polar coordinates satisfy

$x = r \cos (θ) and y = r \sin (θ) .$

🔗

🔗
The area element $d A$ in polar coordinates is determined by the area of a slice of an annulus and is given by

$d A = r d r d θ .$

🔗

🔗
To convert the double integral $\iint_{D} f (x, y) d A$ to an iterated integral in polar coordinates, we substitute $r \cos (θ)$ for $x,$ $r \sin (θ)$ for $y,$ and $r d r d θ$ for $d A$ to obtain the iterated integral

$\iint_{D} f (r \cos (θ), r \sin (θ)) r d r d θ .$

🔗

🔗

🔗

Exercises 11.5.4 Exercises

1.

For each set of Polar coordinates, match the equivalent Cartesian coordinates.

🔗

2.

(a) The Cartesian coordinates of a point are

(- 1, - \sqrt{3}) .

🔗

(i) Find polar coordinates

(r, θ)

of the point, where

r > 0

and

0 \leq θ < 2 π .

🔗

r =

🔗

θ =

🔗

(ii) Find polar coordinates

(r, θ)

of the point, where

r < 0

and

0 \leq θ < 2 π .

🔗

r =

🔗

θ =

🔗

(b) The Cartesian coordinates of a point are

(- 2, 3) .

🔗

(i) Find polar coordinates

(r, θ)

of the point, where

r > 0

and

0 \leq θ < 2 π .

🔗

r =

🔗

θ =

🔗

(ii) Find polar coordinates

(r, θ)

of the point, where

r < 0

and

0 \leq θ < 2 π .

🔗

r =

🔗

θ =

🔗

3.

(a) You are given the point

(1, π / 2)

in polar coordinates.

🔗

(i) Find another pair of polar coordinates for this point such that

r > 0

and

2 π \leq θ < 4 π .

🔗

r =

🔗

θ =

🔗

(ii) Find another pair of polar coordinates for this point such that

r < 0

and

0 \leq θ < 2 π .

🔗

r =

🔗

θ =

🔗

(b) You are given the point

(- 2, π / 4)

in polar coordinates.

🔗

(i) Find another pair of polar coordinates for this point such that

r > 0

and

2 π \leq θ < 4 π .

🔗

r =

🔗

θ =

🔗

(ii) Find another pair of polar coordinates for this point such that

r < 0

and

- 2 π \leq θ < 0.

🔗

r =

🔗

θ =

🔗

(3, 2)

in polar coordinates.

🔗

(i) Find another pair of polar coordinates for this point such that

r > 0

and

2 π \leq θ < 4 π .

🔗

r =

🔗

θ =

🔗

(ii) Find another pair of polar coordinates for this point such that

r < 0

and

0 \leq θ < 2 π .

🔗

r =

🔗

θ =

🔗

4.

Decide if the points given in polar coordinates are the same. If they are the same, enter T . If they are different, enter F .

🔗

a.)

(7, \frac{π}{3}), (- 7, \frac{- π}{3})

🔗

b.)

(2, \frac{51 π}{4}), (2, - \frac{51 π}{4})

🔗

c.)

(0, 7 π), (0, \frac{6 π}{4})

🔗

d.)

(1, \frac{117 π}{4}), (- 1, \frac{π}{4})

🔗

e.)

(17, \frac{44 π}{3}), (- 17, \frac{- π}{3})

🔗

f.)

(7, 7 π), (- 7, 7 π)

🔗

5.

A curve with polar equation

🔗

r = \frac{9}{8 \sin θ + 25 \cos θ}

🔗

represents a line. Write this line in the given Cartesian form.

🔗

y =

🔗

Note: Your answer should be a function of

x

🔗

6.

Find a polar equation of the form

r = f (θ)

for the curve represented by the Cartesian equation

x = - y^{2} .

🔗

Note: Since

θ

is not a symbol on your keyboard, use

t

in place of

θ

in your answer.

🔗

r =

🔗

7.

By changing to polar coordinates, evaluate the integral

\iint_{D} (x^{2} + y^{2})^{5 / 2} d x d y

where

D

is the disk

x^{2} + y^{2} \leq 4 .

🔗

Answer =

🔗

8.

Convert the integral

\int_{0}^{\sqrt{2}} \int_{- x}^{x} d y d x

to polar coordinates and evaluate it (use

t

for

θ

🔗

With

a =

b =

c =

and

d =

🔗

\int_{0}^{\sqrt{2}} \int_{- x}^{x} d y d x = \int_{a}^{b} \int_{c}^{d}

d r d t

🔗

= \int_{a}^{b}

d t

🔗

=

|_{a}^{b}

🔗

=

🔗

9.

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

θ

in your expressions.)

🔗

(a) The region

🔗

With

a =

b =

🔗

c =

, and

d =

🔗

integral =

\int_{a}^{b} \int_{c}^{d}

d

d

🔗

(b) The region

🔗

With

a =

b =

🔗

c =

, and

d =

🔗

integral =

\int_{a}^{b} \int_{c}^{d}

d

d

🔗

10.

A Cartesian equation for the polar equation

r = 4

can be written as:

🔗

x^{2} + y^{2} =

🔗

11.

Using polar coordinates, evaluate the integral which gives the area which lies in the first quadrant between the circles

x^{2} + y^{2} = 100

and

x^{2} - 10 x + y^{2} = 0 .

🔗

12.

(a) Graph

r = 1 / (6 \cos θ)

for

- π / 2 \leq θ \leq π / 2

and

r = 1 .

Then write an iterated integral in polar coordinates representing the area inside the curve

r = 1

and to the right of

r = 1 / (6 \cos θ) .

(Use

t

for

θ

in your work.)

🔗

With

a =

b =

🔗

c =

, and

d =

🔗

area =

\int_{a}^{b} \int_{c}^{d}

d

d

🔗

(b) Evaluate your integral to find the area.

🔗

area =

🔗

13.

Using polar coordinates, evaluate the integral

\int \int_{R} \sin (x^{2} + y^{2}) d A

where R is the region

9 \leq x^{2} + y^{2} \leq 49 .

🔗

14.

Sketch the region of integration for the following integral.

🔗

\int_{0}^{π / 4} \int_{0}^{6 / \cos (θ)} f (r, θ) r d r d θ

🔗

The region of integration is bounded by

🔗

$y = 0, x = \sqrt{36 - y^{2}}, and y = 6$
🔗

🔗
$y = 0, y = x, and x = 6$
🔗

🔗
$y = 0, y = x, and y = 6$
🔗

🔗
$y = 0, y = \sqrt{36 - x^{2}}, and x = 6$
🔗

🔗
None of the above
🔗

🔗

🔗

15.

Use the polar coordinates to find the volume of a sphere of radius 10.

🔗

16.

Consider the solid under the graph of

z = e^{- x^{2} - y^{2}}

above the disk

x^{2} + y^{2} \leq a^{2},

where

a > 0 .

🔗

(a) Set up the integral to find the volume of the solid.

🔗

Instructions: Please enter the integrand in the first answer box, typing theta for

θ .

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.

🔗

\int_{A}^{B} \int_{C}^{D}

A =

B =

C =

D =

(b) Evaluate the integral and find the volume. Your answer will be in terms of

a .

🔗

Volume V =

🔗

a \to \infty ?

🔗

lim_{a \to \infty} V =

🔗

17.

Consider the iterated integral

I = \int_{- 3}^{0} \int_{- \sqrt{9 - y^{2}}}^{0} \frac{y}{x^{2} + y^{2} + 1} d x d y .

Sketch (and label) the region of integration.
🔗

🔗
Convert the given iterated integral to one in polar coordinates.
🔗

🔗
Evaluate the iterated integral in (b).
🔗

🔗
State one possible interpretation of the value you found in (c).
🔗

🔗

🔗

18.

Let

D

be the region that lies inside the unit circle in the plane.

Set up and evaluate an iterated integral in polar coordinates whose value is the area of $D .$
🔗

🔗
Determine the exact average value of $f (x, y) = y$ over the upper half of $D .$
🔗

🔗
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 $δ (x, y) = 1 .$ (Before making any calculations, where do you expect the center of mass to lie? Why?)
🔗

🔗
Find the exact volume of the solid that lies under the surface $z = 8 - x^{2} - y^{2}$ and over the unit disk, $D .$
🔗

🔗

🔗

19.

For each of the following iterated integrals,

sketch and label the region of integration,
🔗

🔗
convert the integral to the other coordinate system (if given in polar, to rectangular; if given in rectangular, to polar), and
🔗

🔗
choose one of the two iterated integrals to evaluate exactly.
🔗

🔗

$\int_{π}^{3 π / 2} \int_{0}^{3} r^{3} d r d θ$
🔗

🔗
$\int_{0}^{2} \int_{- \sqrt{1 - (x - 1)^{2}}}^{\sqrt{1 - (x - 1)^{2}}} \sqrt{x^{2} + y^{2}} d y d x$
🔗

🔗
$\int_{0}^{π / 2} \int_{0}^{\sin (θ)} r \sqrt{1 - r^{2}} d r d θ .$
🔗

🔗
$\int_{0}^{\sqrt{2} / 2} \int_{y}^{\sqrt{1 - y^{2}}} \cos (x^{2} + y^{2}) d x d y .$
🔗

🔗

🔗

You have attempted 1 of 17 activities on this page.

🔗

Prev Top Next