14 Multiple Integration 14.6 Volume Between Surfaces and Triple Integration 15 Vector Analysis

14.7 Triple Integration with Cylindrical and Spherical Coordinates

Just as polar coordinates gave us a new way of describing curves in the plane, in this section we will see how cylindrical and spherical coordinates give us new ways of describing surfaces and regions in space.

^†^†margin:
(fullscreen) Figure 14.7.1: Illustrating the principles behind cylindrical coordinates.

Cylindrical Coordinates

In short, cylindrical coordinates can be thought of as a combination of the polar and rectangular coordinate systems. One can identify a point $(x_{0},y_{0},z_{0})$ , given in rectangular coordinates, with the point $(r_{0},\theta_{0},z_{0})$ , given in cylindrical coordinates, where the $z$ -value in both systems is the same, and the point $(x_{0},y_{0})$ in the $x$ - $y$ plane is identified with the polar point $P(r_{0},\theta_{0})$ ; see Figure 14.7.1. So that each point in space that does not lie on the $z$ -axis is defined uniquely, we will restrict $r\geq 0$ and $0\leq\theta<2\pi$ .

We use the identity $z=z$ along with the identities found in Key Idea 10.4.1 to convert between the rectangular coordinate $(x,y,z)$ and the cylindrical coordinate $(r,\theta,z)$ , namely:

	$\displaystyle\text{From rectangular to cylindrical:}\qquad r=\sqrt{x^{2}+y^{2}% },\qquad\tan\theta=y/x\qquad\text{and}\qquad z=z;$
	$\displaystyle\text{From cylindrical to rectangular:}\qquad x=r\cos\theta\qquad y% =r\sin\theta\qquad\text{and}\qquad z=z.$

These identities, along with conversions related to spherical coordinates, are given later in Key Idea 14.7.1.

^†^†margin: Note: Our rectangular to polar conversion formulas used

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

, allowing for negative

r

values. Since we now restrict

r\geq 0

, we can use

r=\sqrt{x^{2}+y^{2}}

Watch the video:
Conversion From Rectangular Coordinates from https://youtu.be/w9wCdLuklw8

Example 14.7.1 Converting between rectangular and cylindrical coordinates

Convert the rectangular point $(2,-2,1)$ to cylindrical coordinates, and convert the cylindrical point $(4,3\pi/4,5)$ to rectangular.

SolutionFollowing the identities given above (also in Key Idea 14.7.1), we have that $r=\sqrt{2^{2}+(-2)^{2}}=2\sqrt{2}$ . Using $\tan\theta=y/x$ , we find that $\theta=\tan^{-1}(-2/2)=-\pi/4$ . As we restrict $\theta$ to being between $0$ and $2\pi$ , we set $\theta=7\pi/4$ . Finally, $z=1$ , giving the cylindrical point $(2\sqrt{2},7\pi/4,1)$ .

In converting the cylindrical point $(4,3\pi/4,5)$ to rectangular, we have $x=4\cos\bigl(3\pi/4\bigr)=-2\sqrt{2}$ , $y=4\sin\bigl(3\pi/4\bigr)=2\sqrt{2}$ and $z=5$ , giving the rectangular point $(-2\sqrt{2},2\sqrt{2},5)$ .

Setting each of $r$ , $\theta$ and $z$ equal to a constant defines a surface in space, as illustrated in the following example.

^†^†margin:
(fullscreen) Figure 14.7.2: Graphing the canonical surfaces in cylindrical coordinates from Example 14.7.2.

Example 14.7.2 Canonical surfaces in cylindrical coordinates

Describe the surfaces $r=1$ , $\theta=\pi/3$ and $z=2$ , given in cylindrical coordinates.

SolutionThe equation $r=1$ describes all points in space that are 1 unit away from the $z$ -axis. This surface is a “tube” or “cylinder” of radius 1, centered on the $z$ -axis, as graphed in Figure 11.1.8 (which describes the cylinder $x^{2}+y^{2}=1$ in space).

The equation $\theta=\pi/3$ describes the plane formed by extending the line $\theta=\pi/3$ , as given by polar coordinates in the $x$ - $y$ plane, parallel to the $z$ -axis.

The equation $z=2$ describes the plane of all points in space that are 2 units above the $x$ - $y$ plane. This plane is the same as the plane described by $z=2$ in rectangular coordinates.

All three surfaces are graphed in Figure 14.7.2. Note how their intersection uniquely defines the point $P=(1,\pi/3,2)$ .

Cylindrical coordinates are useful when describing certain domains in space, allowing us to evaluate triple integrals over these domains more easily than if we used rectangular coordinates.

Theorem 14.6.4 shows how to evaluate $\iiint_{D}h(x,y,z)\operatorname{d}\!V$ using rectangular coordinates. In that evaluation, we use $\operatorname{d}\!V=\operatorname{d}\!z\operatorname{d}\!y\operatorname{d}\!x$ (or one of the other five orders of integration). Recall how, in this order of integration, the bounds on $y$ are “curve to curve” and the bounds on $x$ are “point to point”: these bounds describe a region $R$ in the $x$ - $y$ plane. We could describe $R$ using polar coordinates as done in Section 14.3. In that section, we saw how we used $\operatorname{d}\!A=r\operatorname{d}\!r\operatorname{d}\!\theta$ instead of $\operatorname{d}\!A=\operatorname{d}\!y\operatorname{d}\!x$ .

Considering the above thoughts, we have $\operatorname{d}\!V=\operatorname{d}\!z\bigl(r\operatorname{d}\!r\operatorname% {d}\!\theta\bigr)=r\operatorname{d}\!z\operatorname{d}\!r\operatorname{d}\!\theta$ . We set bounds on $z$ as “surface to surface” as done in the previous section, and then use “curve to curve” and “point to point” bounds on $r$ and $\theta$ , respectively. Finally, using the identities given above, we change the integrand $h(x,y,z)$ to $h(r,\theta,z)$ .

This process should sound plausible; the following theorem states it is truly a way of evaluating a triple integral.

Theorem 14.7.1 Triple Integration in Cylindrical Coordinates

Let $w=h(r,\theta,z)$ be a continuous function on a closed, bounded region $D$ in space, bounded in cylindrical coordinates by $\alpha\leq\theta\leq\beta$ , $g_{1}(\theta)\leq r\leq g_{2}(\theta)$ and $f_{1}(r,\theta)\leq z\leq f_{2}(r,\theta)$ . Then

\iiint_{D}h(r,\theta,z)\operatorname{d}\!V=\int_{\alpha}^{\beta}\int_{g_{1}(% \theta)}^{g_{2}(\theta)}\int_{f_{1}(r,\theta)}^{f_{2}(r,\theta)}h(r,\theta,z)r% \operatorname{d}\!z\operatorname{d}\!r\operatorname{d}\!\theta.

^†^†margin:
(fullscreen) Figure 14.7.3: Visualizing the solid used in Example 14.7.3.

Example 14.7.3 Evaluating a triple integral with cylindrical coordinates

Find the mass of the solid represented by the region in space bounded by $z=0$ , $z=\sqrt{4-x^{2}-y^{2}}+3$ and the cylinder $x^{2}+y^{2}=4$ (as shown in Figure 14.7.3), with density function $\delta(x,y,z)=x^{2}+y^{2}+z+1$ , using a triple integral in cylindrical coordinates. Distances are measured in centimeters and density is measured in grams/cm³.

SolutionWe begin by describing this region of space with cylindrical coordinates. The plane $z=0$ is left unchanged; with the identity $r=\sqrt{x^{2}+y^{2}}$ , we convert the hemisphere of radius 2 to the equation $z=\sqrt{4-r^{2}}$ ; the cylinder $x^{2}+y^{2}=4$ is converted to $r^{2}=4$ , or, more simply, $r=2$ . We also convert the density function: $\delta(r,\theta,z)=r^{2}+z+1$ .

To describe this solid with the bounds of a triple integral, we bound $z$ with $0\leq z\leq\sqrt{4-r^{2}}+3$ ; we bound $r$ with $0\leq r\leq 2$ ; we bound $\theta$ with $0\leq\theta\leq 2\pi$ .

Using Definition 14.6.4 and Theorem 14.7.1, we have the mass of the solid is

	$\displaystyle M=\iiint_{D}\delta(x,y,z)\operatorname{d}\!V$	$\displaystyle=\int_{0}^{2\pi}\int_{0}^{2}\int_{0}^{\sqrt{4-r^{2}}+3}\bigl(r^{2% }+z+1\bigr)r\operatorname{d}\!z\operatorname{d}\!r\operatorname{d}\!\theta$
		$\displaystyle=\int_{0}^{2\pi}\int_{0}^{2}\bigl((r^{3}+4r)\sqrt{4-r^{2}}+\frac{% 5}{2}r^{3}+\frac{19}{2}r\bigr)\operatorname{d}\!r\operatorname{d}\!\theta$
		$\displaystyle=\frac{1318\pi}{15}\approx 276.04\text{ gm},$

where we leave the details of the remaining double integral to the reader.

Example 14.7.4 Finding the center of mass using cylindrical coordinates

Find ^†^†margin:
(fullscreen) Figure 14.7.4: Visualizing the solid used in Example 14.7.4. the center of mass of the solid with constant density whose base can be described by the polar curve $r=\cos(3\theta)$ and whose top is defined by the plane $z=1-x+0.1y$ , where distances are measured in feet, as seen in Figure 14.7.4. (The volume of this solid was found in Example 14.3.5.)

SolutionWe convert the equation of the plane to use cylindrical coordinates: $z=1-r\cos\theta+0.1r\sin\theta$ . Thus the region is space is bounded by $0\leq z\leq 1-r\cos\theta+0.1r\sin\theta$ , $0\leq r\leq\cos(3\theta)$ , $0\leq\theta\leq\pi$ (recall that the rose curve $r=\cos(3\theta)$ is traced out once on $[0,\pi]$ .

Since density is constant, we set $\delta=1$ and finding the mass is equivalent to finding the volume of the solid. We set up the triple integral to compute this but do not evaluate it; we leave it to the reader to confirm it evaluates to the same result found in Example 14.3.5.

M=\iiint_{D}\delta\operatorname{d}\!V=\int_{0}^{\pi}\int_{0}^{\cos(3\theta)}% \int_{0}^{1-r\cos\theta+0.1r\sin\theta}r\operatorname{d}\!z\operatorname{d}\!r% \operatorname{d}\!\theta=\frac{\pi}{4}.

From Definition 14.6.4 we set up the triple integrals to compute the moments about the three coordinate planes. The computation of each is left to the reader (using technology is recommended):

	$\displaystyle M_{yz}=\iiint_{D}x\operatorname{d}\!V$	$\displaystyle=\int_{0}^{\pi}\int_{0}^{\cos(3\theta)}\int_{0}^{1-r\cos\theta+0.% 1r\sin\theta}(r\cos\theta)r\operatorname{d}\!z\operatorname{d}\!r\operatorname% {d}\!\theta$
		$\displaystyle=-\frac{3\pi}{64}.$
	$\displaystyle M_{xz}=\iiint_{D}y\operatorname{d}\!V$	$\displaystyle=\int_{0}^{\pi}\int_{0}^{\cos(3\theta)}\int_{0}^{1-r\cos\theta+0.% 1r\sin\theta}(r\sin\theta)r\operatorname{d}\!z\operatorname{d}\!r\operatorname% {d}\!\theta$
		$\displaystyle=\frac{3\pi}{640}.$
	$\displaystyle M_{xy}=\iiint_{D}z\operatorname{d}\!V$	$\displaystyle=\int_{0}^{\pi}\int_{0}^{\cos(3\theta)}\int_{0}^{1-r\cos\theta+0.% 1r\sin\theta}(z)r\operatorname{d}\!z\operatorname{d}\!r\operatorname{d}\!\theta$
		$\displaystyle=\frac{1903\pi}{12800}.$

The center of mass is located at

(-3/16,3/160,1903/3200)=(-0.1875,0.01875,0.5946875),

which lies outside the bounds of the solid.

Spherical Coordinates

^†^†margin:
(fullscreen) Figure 14.7.5: Illustrating the principles behind spherical coordinates.

In short, spherical coordinates can be thought of as a “double application” of the polar coordinate system. In spherical coordinates, a point $P$ is identified with $(\rho,\theta,\varphi)$ , where $\rho$ is the distance from the origin to $P$ , $\theta$ is the same angle as would be used to describe $P$ in the cylindrical coordinate system, and $\varphi$ is the angle between the positive $z$ -axis and the ray from the origin to $P$ ; see Figure 14.7.5. So that each point in space that does not lie on the $z$ -axis is defined uniquely, we will restrict $\rho\geq 0$ , $0\leq\theta<2\pi$ and $0\leq\varphi\leq\pi$ .

The following Key Idea gives conversions to/from our three spatial coordinate systems.

^†^†margin: Note: The symbol

\rho

is the Greek letter “rho.” Traditionally it is used in the spherical coordinate system, while

r

is used in the polar and cylindrical coordinate systems.

Key Idea 14.7.1 Converting Between Rectangular, Cylindrical and Spherical Coordinates

Rectangular and Cylindrical $\displaystyle r^{2}$ $\displaystyle=x^{2}+y^{2},$ $\displaystyle\tan\theta$ $\displaystyle=y/x,$ $\displaystyle z$ $\displaystyle=z$ $\displaystyle x$ $\displaystyle=r\cos\theta,$ $\displaystyle y$ $\displaystyle=r\sin\theta,$ $\displaystyle z$ $\displaystyle=z$ Rectangular and Spherical $\displaystyle\rho$ $\displaystyle=\sqrt{x^{2}+y^{2}+z^{2}},$ $\displaystyle\tan\theta$ $\displaystyle=\frac{y}{x},$ $\displaystyle\cos\varphi$ $\displaystyle=\frac{z}{\sqrt{x^{2}+y^{2}+z^{2}}}$ $\displaystyle x$ $\displaystyle=\rho\sin\varphi\cos\theta,$ $\displaystyle y$ $\displaystyle=\rho\sin\varphi\sin\theta,$ $\displaystyle z$ $\displaystyle=\rho\cos\varphi$ Cylindrical and Spherical $\displaystyle\rho$ $\displaystyle=\sqrt{r^{2}+z^{2}},$ $\displaystyle\theta$ $\displaystyle=\theta,$ $\displaystyle\tan\varphi$ $\displaystyle=r/z$ $\displaystyle r$ $\displaystyle=\rho\sin\varphi,$ $\displaystyle\theta$ $\displaystyle=\theta,$ $\displaystyle z$ $\displaystyle=\rho\cos\varphi$

^†^†margin: Note: The role of

\theta

and

\varphi

in spherical coordinates differs between mathematicians and physicists. When reading about physics in spherical coordinates, be careful to note how that particular author uses these variables and recognize that these identities may be rearranged.

Example 14.7.5 Converting between rectangular and spherical coordinates

Convert the rectangular point $(2,-2,1)$ to spherical coordinates, and convert the spherical point $(6,\pi/3,\pi/2)$ to rectangular and cylindrical coordinates.

SolutionThis rectangular point is the same as used in Example 14.7.1. Using Key Idea 14.7.1, we find $\rho=\sqrt{2^{2}+(-1)^{2}+1^{2}}=3$ . Using the same logic as in Example 14.7.1, we find $\theta=7\pi/4$ . Finally, $\cos\varphi=1/3$ , giving $\varphi=\cos^{-1}(1/3)\approx 1.23$ , or about $70.53^{\circ}$ . Thus the spherical coordinates are approximately $(3,7\pi/4,1.23)$ .

Converting the spherical point $(6,\pi/3,\pi/2)$ to rectangular, we have

	$\displaystyle x$	$\displaystyle=6\sin(\pi/2)\cos(\pi/3)=3,$
	$\displaystyle y$	$\displaystyle=6\sin(\pi/2)\sin(\pi/3)=3\sqrt{3},\text{ and}$
	$\displaystyle z$	$\displaystyle=6\cos(\pi/2)=0.$

Thus the rectangular coordinates are $(3,3\sqrt{3},0)$ .

To convert this spherical point to cylindrical, we have $r=6\sin(\pi/2)=6$ , $\theta=\pi/3$ and $z=6\cos(\pi/2)=0$ , giving the cylindrical point $(6,\pi/3,0)$ .

Example 14.7.6 Canonical surfaces in spherical coordinates

Describe the surfaces $\rho=1$ , $\theta=\pi/3$ and $\varphi=\pi/6$ , given in spherical coordinates. ^†^†margin:
(fullscreen) Figure 14.7.6: Graphing the canonical surfaces in spherical coordinates from Example 14.7.6.

SolutionThe equation $\rho=1$ describes all points in space that are 1 unit away from the origin: this is the sphere of radius 1, centered at the origin.

The equation $\theta=\pi/3$ describes the same surface in spherical coordinates as it does in cylindrical coordinates: beginning with the line $\theta=\pi/3$ in the $x$ - $y$ plane as given by polar coordinates, extend the line parallel to the $z$ -axis, forming a plane.

The equation $\varphi=\pi/6$ describes all points $P$ in space where the ray from the origin to $P$ makes an angle of $\pi/6$ with the positive $z$ -axis. This describes a cone, with the positive $z$ -axis its axis of symmetry, with point at the origin.

All three surfaces are graphed in Figure 14.7.6. Note how their intersection uniquely defines the point $P=(1,\pi/3,\pi/6)$ .

Spherical coordinates are useful when describing certain domains in space, allowing us to evaluate triple integrals over these domains more easily than if we used rectangular coordinates or cylindrical coordinates. The crux of setting up a triple integral in spherical coordinates is appropriately describing the “small amount of volume,” $d V$ , used in the integral.

^†^†margin:
(fullscreen) Figure 14.7.7: Approximating the volume of a standard region in space using spherical coordinates.

Considering Figure 14.7.7, we can make a small “spherical wedge” by varying $\rho$ , $\theta$ and $\varphi$ each a small amount, $\Delta\rho$ , $\Delta\theta$ and $\Delta\varphi$ , respectively. This wedge is approximately a rectangular solid when the change in each coordinate is small, giving a volume of about

\Delta V\approx\Delta\rho\ \times\ \rho\Delta\varphi\ \times\ \rho\sin(\varphi% )\Delta\theta.

Given a region $D$ in space, we can approximate the volume of $D$ with many such wedges. As the size of each of $\Delta\rho$ , $\Delta\theta$ and $\Delta\varphi$ goes to zero, the number of wedges increases to infinity and the volume of $D$ is more accurately approximated, giving

\operatorname{d}\!V=\operatorname{d}\!\rho\ \times\ \rho\operatorname{d}\!% \varphi\ \times\ \rho\sin(\varphi)\operatorname{d}\!\theta=\rho^{2}\sin(% \varphi)\operatorname{d}\!\rho\operatorname{d}\!\theta\operatorname{d}\!\varphi.

Again, this development of $\operatorname{d}\!V$ should sound reasonable, and the following theorem states it is the appropriate manner by which triple integrals are to be evaluated in spherical coordinates.

^†^†margin: Note: It is often most intuitive to evaluate the triple integral in Theorem 14.7.2 by integrating with respect to

\rho

first; it often does not matter whether we next integrate with respect to

\theta

\varphi

. Different texts present different standard orders, some preferring

\operatorname{d}\!\varphi\operatorname{d}\!\theta

instead of

\operatorname{d}\!\theta\operatorname{d}\!\varphi

. As the bounds for these variables are usually constants in practice, it generally is a matter of preference.

Theorem 14.7.2 Triple Integration in Spherical Coordinates

Let $w=h(\rho,\theta,\varphi)$ be a continuous function on a closed, bounded region $D$ in space, bounded in spherical coordinates by $\alpha_{1}\leq\varphi\leq\alpha_{2}$ , $\beta_{1}\leq\theta\leq\beta_{2}$ and $f_{1}(\theta,\varphi)\leq\rho\leq f_{2}(\theta,\varphi)$ . Then

\iiint_{D}h(\rho,\theta,\varphi)\operatorname{d}\!V=\int_{\alpha_{1}}^{\alpha_% {2}}\int_{\beta_{1}}^{\beta_{2}}\int_{f_{1}(\theta,\varphi)}^{f_{2}(\theta,% \varphi)}h(\rho,\theta,\varphi)\rho^{2}\sin(\varphi)\operatorname{d}\!\rho% \operatorname{d}\!\theta\operatorname{d}\!\varphi.

Example 14.7.7 Establishing the volume of a sphere

Let $D$ be the region in space bounded by the sphere, centered at the origin, of radius $r$ . Use a triple integral in spherical coordinates to find the volume $V$ of $D$ .

SolutionThe sphere of radius $r$ , centered at the origin, has equation $\rho=r$ . To obtain the full sphere, the bounds on $\theta$ and $\varphi$ are $0\leq\theta\leq 2\pi$ and $0\leq\varphi\leq\pi$ . This leads us to:

	$\displaystyle V$	$\displaystyle=\iiint_{D}\operatorname{d}\!V$
		$\displaystyle=\int_{0}^{\pi}\int_{0}^{2\pi}\int_{0}^{r}\bigl(\rho^{2}\sin(% \varphi)\bigr)\operatorname{d}\!\rho\operatorname{d}\!\theta\operatorname{d}\!\varphi$
		$\displaystyle=\int_{0}^{\pi}\int_{0}^{2\pi}\left(\frac{1}{3}\rho^{3}\sin(% \varphi)\Bigr\|_{\rho=0}^{\rho=r}\right)\operatorname{d}\!\theta\operatorname{d% }\!\varphi$
		$\displaystyle=\int_{0}^{\pi}\int_{0}^{2\pi}\left(\frac{1}{3}r^{3}\sin(\varphi)% \right)\operatorname{d}\!\theta\operatorname{d}\!\varphi$
		$\displaystyle=\int_{0}^{\pi}\left(\frac{2\pi}{3}r^{3}\sin(\varphi)\right)% \operatorname{d}\!\varphi$
		$\displaystyle=\left.\left(-\frac{2\pi}{3}r^{3}\cos(\varphi)\right)\right\|_{% \varphi=0}^{\varphi=\pi}$
		$\displaystyle=\frac{4\pi}{3}r^{3},$

the familiar formula for the volume of a sphere. Note how the integration steps were easy, not using square-roots nor integration steps such as Substitution.

^†^†margin:
(fullscreen) Figure 14.7.8: Graphing the solid, and its center of mass, from Example 14.7.8.

Example 14.7.8 Finding the center of mass using spherical coordinates

Find the center of mass of the solid with constant density enclosed above by $\rho=4$ and below by $\varphi=\pi/6$ , as illustrated in Figure 14.7.8.

SolutionWe will set up the four triple integrals needed to find the center of mass (i.e., to compute $M$ , $M_{yz}$ , $M_{xz}$ and $M_{xy}$ ) and leave it to the reader to evaluate each integral. Because of symmetry, we expect the $x$ - and $y$ - coordinates of the center of mass to be 0.

While the surfaces describing the solid are given in the statement of the problem, to describe the full solid $D$ , we use the following bounds: $0\leq\rho\leq 4$ , $0\leq\theta\leq 2\pi$ and $0\leq\varphi\leq\pi/6$ . Since density $\delta$ is constant, we assume $\delta=1$ .

The mass of the solid:

	$\displaystyle M$	$\displaystyle=\iiint_{D}\operatorname{d}\!m=\iiint_{D}\operatorname{d}\!V$
		$\displaystyle=\int_{0}^{\pi/6}\int_{0}^{2\pi}\int_{0}^{4}\bigl(\rho^{2}\sin(% \varphi)\bigr)\operatorname{d}\!\rho\operatorname{d}\!\theta\operatorname{d}\!\varphi$
		$\displaystyle=\frac{64}{3}\bigl(2-\sqrt{3}\bigr)\pi\approx 17.958.$

To compute $M_{yz}$ , the integrand is $x$ ; using Key Idea 14.7.1, we have $x=\rho\sin\varphi\cos\theta$ . This gives:

	$\displaystyle M_{yz}$	$\displaystyle=\iiint_{D}x\operatorname{d}\!m$
		$\displaystyle=\int_{0}^{\pi/6}\int_{0}^{2\pi}\int_{0}^{4}\bigl((\rho\sin(% \varphi)\cos(\theta))\rho^{2}\sin(\varphi)\bigr)\operatorname{d}\!\rho% \operatorname{d}\!\theta\operatorname{d}\!\varphi$
		$\displaystyle=\int_{0}^{\pi/6}\int_{0}^{2\pi}\int_{0}^{4}\bigl(\rho^{3}\sin^{2% }(\varphi)\cos(\theta)\bigr)\operatorname{d}\!\rho\operatorname{d}\!\theta% \operatorname{d}\!\varphi$
		$\displaystyle=0,$

which we expected as we expect $\overline{x}=0$ .

To compute $M_{xz}$ , the integrand is $y$ ; using Key Idea 14.7.1, we have $y=\rho\sin\varphi\sin\theta$ . This gives:

	$\displaystyle M_{xz}$	$\displaystyle=\iiint_{D}y\operatorname{d}\!m$
		$\displaystyle=\int_{0}^{\pi/6}\int_{0}^{2\pi}\int_{0}^{4}\bigl((\rho\sin(% \varphi)\sin(\theta))\rho^{2}\sin(\varphi)\bigr)\operatorname{d}\!\rho% \operatorname{d}\!\theta\operatorname{d}\!\varphi$
		$\displaystyle=\int_{0}^{\pi/6}\int_{0}^{2\pi}\int_{0}^{4}\bigl(\rho^{3}\sin^{2% }(\varphi)\sin(\theta)\bigr)\operatorname{d}\!\rho\operatorname{d}\!\theta% \operatorname{d}\!\varphi$
		$\displaystyle=0,$

which we also expected as we expect $\overline{y}=0$ .

To compute $M_{xy}$ , the integrand is $z$ ; using Key Idea 14.7.1, we have $z=\rho\cos\varphi$ . This gives:

	$\displaystyle M_{xy}$	$\displaystyle=\iiint_{D}z\operatorname{d}\!m$
		$\displaystyle=\int_{0}^{\pi/6}\int_{0}^{2\pi}\int_{0}^{4}\bigl((\rho\cos(% \varphi))\rho^{2}\sin(\varphi)\bigr)\operatorname{d}\!\rho\operatorname{d}\!% \theta\operatorname{d}\!\varphi$
		$\displaystyle=\int_{0}^{\pi/6}\int_{0}^{2\pi}\int_{0}^{4}\bigl(\rho^{3}\cos(% \varphi)\sin(\varphi)\bigr)\operatorname{d}\!\rho\operatorname{d}\!\theta% \operatorname{d}\!\varphi$
		$\displaystyle=16\pi\approx 50.266.$

Thus the center of mass is $(0,0,M_{xy}/M)=(0,0,3(2+\sqrt{3})/4)\approx(0,0,2.799)$ , as indicated in Figure 14.7.8.

Change of Variables in Multiple Integrals

Given the difficulty of evaluating multiple integrals, the reader may be wondering if it is possible to simplify those integrals using a suitable substitution for the variables. The answer is yes, though it is a bit more complicated than the substitution method which you learned in single-variable calculus.

Recall that if you are given, for example, the definite integral

\int_{1}^{2}x^{3}\sqrt{x^{2}-1}\operatorname{d}\!x,

then you would make the substitution

	$\displaystyle u$	$\displaystyle=x^{2}-1\Rightarrow x^{2}=u+1$
	$\displaystyle\operatorname{d}\!u$	$\displaystyle=2x\operatorname{d}\!x$
which changes the limits of integration
	$\displaystyle x$	$\displaystyle=1\Rightarrow u=0$
	$\displaystyle x$	$\displaystyle=2\Rightarrow u=3$

so that we get

	$\displaystyle\int_{1}^{2}x^{3}\sqrt{x^{2}-1}\operatorname{d}\!x$	$\displaystyle=\int_{1}^{2}\frac{1}{2}x^{2}\cdot 2x\sqrt{x^{2}-1}\operatorname{% d}\!x$
		$\displaystyle=\int_{0}^{3}\frac{1}{2}(u+1)\sqrt{u}\operatorname{d}\!u$
		$\displaystyle=\frac{1}{2}\int_{0}^{3}\left(u^{3/2}+u^{1/2}\right)\operatorname% {d}\!u$
		$\displaystyle=\frac{14\sqrt{3}}{5}.$

Let us take a different look at what happened when we did that substitution, which will give some motivation for how substitution works in multiple integrals. First, we let $u=x^{2}-1$ . On the interval of integration $[1,2]$ , the function $x\mapsto x^{2}-1$ is strictly increasing (and maps $[1,2]$ onto $[0,3]$ ) and hence has an inverse function (defined on the interval $[0,3]$ ). That is, on $[0,3]$ we can define $x$ as a function of $u$ , namely

x=g(u)=\sqrt{u+1}.

Then substituting that expression for $x$ into the function $f(x)=x^{3}\sqrt{x^{2}-1}$ gives

f(x)=f(g(u))=(u+1)^{3/2}\sqrt{u},

and we see that

	$\displaystyle\frac{\operatorname{d}\!x}{\operatorname{d}\!u}=g\mkern 1.35mu^{% \prime}(u)\Rightarrow\operatorname{d}\!x$	$\displaystyle=g\mkern 1.35mu^{\prime}(u)\operatorname{d}\!u$
	$\displaystyle\operatorname{d}\!x$	$\displaystyle=\frac{1}{2}(u+1)^{-1/2}\operatorname{d}\!u,$

so since

	$\displaystyle g(0)=1\Rightarrow 0=g^{-1}(1)$
	$\displaystyle g(3)=2\Rightarrow 3=g^{-1}(2)$

then performing the substitution as we did earlier gives

	$\displaystyle\int_{1}^{2}f(x)\operatorname{d}\!x$	$\displaystyle=\int_{1}^{2}x^{3}\sqrt{x^{2}-1}\operatorname{d}\!x$
		$\displaystyle=\int_{0}^{3}\frac{1}{2}(u+1)\sqrt{u}\operatorname{d}\!u,\text{ % which can be written as}$
		$\displaystyle=\int_{0}^{3}(u+1)^{3/2}\sqrt{u}\,\cdot\frac{1}{2}(u+1)^{-1/2}% \operatorname{d}\!u,\text{ which means}$
	$\displaystyle\int_{1}^{2}f(x)\operatorname{d}\!x$	$\displaystyle=\int_{g^{-1}(1)}^{g^{-1}(2)}f(g(u))\,g\mkern 1.35mu^{\prime}(u)% \operatorname{d}\!u.$

In general, if $x=g(u)$ is a one-to-one, differentiable function from an interval $[c,d]$ (which you can think of as being on the “ $u$ -axis”) onto an interval $[a,b]$ (on the $x$ -axis), which means that $g\mkern 1.35mu^{\prime}(u)\neq 0$ on the interval $(c,d)$ , so that $a=g(c)$ and $b=g(d)$ , then $c=g^{-1}(a)$ and $d=g^{-1}(b)$ , and

\int_{a}^{b}f(x)\operatorname{d}\!x=\int_{g^{-1}(a)}^{g^{-1}(b)}f(g(u))\,g% \mkern 1.35mu^{\prime}(u)\operatorname{d}\!u.

This is called the change of variable formula for integrals of single-variable functions, and it is what you were implicitly using when doing integration by substitution. This formula turns out to be a special case of a more general formula which can be used to evaluate multiple integrals. We will state the formulas for double and triple integrals involving real-valued functions of two and three variables, respectively. We will assume that all the functions involved are continuously differentiable and that the regions and solids involved all have “reasonable” boundaries. The proof of the following theorem is beyond the scope of the text.

Theorem 14.7.3 Change of Variables Formula for Multiple Integrals

Let $x=x(u,v)$ and $y=y(u,v)$ define a one-to-one mapping of a region $R^{\prime}$ in the $u v$ -plane onto a region $R$ in the $x y$ -plane such that the determinant

J(u,v)=\begin{vmatrix}\dfrac{\partial x}{\partial u}&\dfrac{\partial x}{% \partial v}\vskip 5.69054pt\\ \dfrac{\partial y}{\partial u}&\dfrac{\partial y}{\partial v}\end{vmatrix}

(14.7.1)

is never $0$ in $R^{\prime}$ . Then

\iint_{R}f(x,y)\operatorname{d}\!A(x,y)=\iint_{R^{\prime}}f(x(u,v),y(u,v))\,% \left\lvert J(u,v)\right\rvert\operatorname{d}\!A(u,v).

(14.7.2)

We use the notation $\operatorname{d}\!A(x,y)$ and $\operatorname{d}\!A(u,v)$ to denote the area element in the $(x,y)$ and $(u,v)$ coordinates, respectively.

Similarly, if $x=x(u,v,w)$ , $y=y(u,v,w)$ and $z=z(u,v,w)$ define a one-to-one mapping of a solid $S^{\prime}$ in $u v w$ -space onto a solid $S$ in $x y z$ -space such that the determinant

J(u,v,w)=\begin{vmatrix}\dfrac{\partial x}{\partial u}&\dfrac{\partial x}{% \partial v}&\dfrac{\partial x}{\partial w}\vskip 5.69054pt\\ \dfrac{\partial y}{\partial u}&\dfrac{\partial y}{\partial v}&\dfrac{\partial y% }{\partial w}\vskip 5.69054pt\\ \dfrac{\partial z}{\partial u}&\dfrac{\partial z}{\partial v}&\dfrac{\partial z% }{\partial w}\end{vmatrix}

(14.7.3)

is never $0$ in $S^{\prime}$ , then

\iiint_{S}f(x,y,z)\operatorname{d}\!V(x,y,z)=\\ \iiint_{S^{\prime}}f(x(u,v,w),y(u,v,w),z(u,v,w))\,\left\lvert J(u,v,w)\right% \rvert\operatorname{d}\!V(u,v,w).

(14.7.4)

^†^†margin: There are two different uses of

\left\lvert\cdot\right\rvert

in these equations. The first

\left\lvert\cdot\right\rvert

refers to the Jacobian, which is the determinant of the matrix. We then take the absolute value of this determinant.

The determinant $J(u,v)$ in Equation (14.7.1) is called the Jacobian of $x$ and $y$ with respect to $u$ and $v$ , and is sometimes written as

J(u,v)=\left\lvert\frac{\partial(x,y)}{\partial(u,v)}\right\rvert.

Similarly, the Jacobian $J(u,v,w)$ of three variables is sometimes written as

J(u,v,w)=\left\lvert\frac{\partial(x,y,z)}{\partial(u,v,w)}\right\rvert.

Notice that Equation (14.7.2) is saying that $\operatorname{d}\!A(x,y)=\left\lvert J(u,v)\right\rvert\operatorname{d}\!A(u,v)$ , which you can think of as a two-variable version of the relation $\operatorname{d}\!x=g\mkern 1.35mu^{\prime}(u)\operatorname{d}\!u$ in the single-variable case.

Watch the video:
Jacobian from https://youtu.be/Bw5yEqwMjQU

The following example shows how the change of variables formula is used.

Example 14.7.9 Change of Variables

Evaluate $\displaystyle\iint_{R}e^{\frac{x-y}{x+y}}\operatorname{d}\!A$ , where $R=\{(x,y):x\geq 0,y\geq 0,x+y\leq 1\}$ .

SolutionFirst, note that evaluating this double integral without using substitution is probably impossible, at least in a closed form. By looking at the numerator and denominator of the exponent of $e$ , we will try the substitution $u=x-y$ and $v=x+y$ . To use the change of variables Equation (14.7.2), we need to write both $x$ and $y$ in terms of $u$ and $v$ . So solving for $x$ and $y$ gives $x=\frac{1}{2}(u+v)$ and $y=\frac{1}{2}(v-u)$ . In Figure 14.7.9 below, we see how the mapping $x=x(u,v)=\frac{1}{2}(u+v)$ , $y=y(u,v)=\frac{1}{2}(v-u)$ maps the region $R^{\prime}$ onto $R$ in a one-to-one manner.

Figure 14.7.9: The regions

R

and

R^{\prime}

Now we see that

J(u,v)=\begin{vmatrix}\dfrac{\partial x}{\partial u}&\dfrac{\partial x}{% \partial v}\vskip 5.69054pt\\ \dfrac{\partial y}{\partial u}&\dfrac{\partial y}{\partial v}\end{vmatrix}=% \begin{vmatrix}\frac{1}{2}&\frac{1}{2}\vskip 5.69054pt\\ -\frac{1}{2}&\frac{1}{2}\end{vmatrix}=\frac{1}{2}\Rightarrow\left\lvert J(u,v)% \right\rvert=\left\lvert\frac{1}{2}\right\rvert=\frac{1}{2},

so using horizontal slices in $R^{\prime}$ , we have

	$\displaystyle\iint_{R}e^{\frac{x-y}{x+y}}\operatorname{d}\!A$	$\displaystyle=\iint_{R^{\prime}}f(x(u,v),y(u,v))\,\left\lvert J(u,v)\right% \rvert\operatorname{d}\!A$
		$\displaystyle=\int_{0}^{1}\int_{-v}^{v}e^{\frac{u}{v}}\,\frac{1}{2}% \operatorname{d}\!u\operatorname{d}\!v$
		$\displaystyle=\int_{0}^{1}\left(\left.\frac{v}{2}e^{\frac{u}{v}}\,\right\|_{u=-% v}^{u=v}\right)\operatorname{d}\!v$
		$\displaystyle=\int_{0}^{1}\frac{v}{2}(e-e^{-1})\operatorname{d}\!v$
		$\displaystyle=\left.\frac{v^{2}}{4}(e-e^{-1})\,\right\|_{0}^{1}=\frac{1}{4}% \left(e-\frac{1}{e}\right)=\frac{e^{2}-1}{4e}.$

The change of variables formula can be used to evaluate double integrals in polar coordinates. Letting

x=x(r,\theta)=r\cos\theta\quad\text{and}\quad y=y(r,\theta)=r\sin\theta,

we have

J(u,v)=\begin{vmatrix}\dfrac{\partial x}{\partial r}&\dfrac{\partial x}{% \partial\theta}\vskip 5.69054pt\\ \dfrac{\partial y}{\partial r}&\dfrac{\partial y}{\partial\theta}\end{vmatrix}% =\begin{vmatrix}\cos\theta&-r\sin\theta\vskip 5.69054pt\\ \sin\theta&r\cos\theta\end{vmatrix}=r\cos^{2}\theta+r\sin^{2}\theta=r\\ \Rightarrow\left\lvert J(u,v)\right\rvert=\left\lvert r\right\rvert=r,

which verifies Key Idea 14.3.1. In a similar fashion, Exercises 51 and 52 ask you to verify Theorems 14.7.1 and 14.7.2 by computing the appropriate Jacobian.

This section has provided a brief introduction into two new coordinate systems useful for identifying points in space. Each can be used to define a variety of surfaces in space beyond the canonical surfaces graphed as each system was introduced.

However, the usefulness of these coordinate systems does not lie in the variety of surfaces that they can describe nor the regions in space these surfaces may enclose. Rather, cylindrical coordinates are mostly used to describe cylinders and spherical coordinates are mostly used to describe spheres. These shapes are of special interest in the sciences, especially in physics, and computations on/inside these shapes is difficult using rectangular coordinates. For instance, in the study of electricity and magnetism, one often studies the effects of an electrical current passing through a wire; that wire is essentially a cylinder, described well by cylindrical coordinates.

This chapter investigated the natural follow-on to partial derivatives: iterated integration. We learned how to use the bounds of a double integral to describe a region in the plane using both rectangular and polar coordinates, then later expanded to use the bounds of a triple integral to describe a region in space. We used double integrals to find volumes under surfaces, surface area, and the center of mass of lamina; we used triple integrals as an alternate method of finding volumes of space regions and also to find the center of mass of a region in space.

Integration does not stop here. We could continue to iterate our integrals, next investigating “quadruple integrals” whose bounds describe a region in 4-dimensional space (which are very hard to visualize). We can also look back to “regular” integration where we found the area under a curve in the plane. A natural analogue to this is finding the “area under a curve,” where the curve is in space, not in a plane. These are just two of many avenues to explore under the heading of “integration.”

Exercises

Terms and Concepts

1.

Explain the difference between the roles $r$ , in cylindrical coordinates, and $\rho$ , in spherical coordinates, play in determining the location of a point.
2.

Why are points on the $z$ -axis not determined uniquely when using cylindrical and spherical coordinates?
3.

What surfaces are naturally defined using cylindrical coordinates?
4.

What surfaces are naturally defined using spherical coordinates?

Problems

In Exercises 5–8, points are given in either the rectangular, cylindrical or spherical coordinate systems. Find the coordinates of the points in the other systems.

5.
- Points in rectangular coordinates:
  $(2,2,1)$ and $(-\sqrt{3},1,0)$
  
  Points in cylindrical coordinates:
  $(2,\pi/4,2)$ and $(3,3\pi/2,-4)$
  
  Points in spherical coordinates:
  $(2,\pi/4,\pi/4)$ and $(1,0,0)$
6.
- Points in rectangular coordinates:
  $(0,1,1)$ and $(-1,0,1)$
  
  Points in cylindrical coordinates:
  $(0,\pi,1)$ and $(2,4\pi/3,0)$
  
  Points in spherical coordinates:
  $(2,\pi/6,\pi/2)$ and $(3,\pi,\pi)$
7.
- Points in rectangular coordinates:
  $(2,2\sqrt{3},-1)$ and $(-5,5,6)$
  
  Points in cylindrical coordinates:
  $(1,\frac{\pi}{3},-2)$ and $(2,\frac{5\pi}{6},3)$
  
  Points in spherical coordinates:
  $(4,\frac{\pi}{3},\frac{\pi}{4})$ and $(1,0,\frac{\pi}{2})$
8.
- Points in rectangular coordinates:
  $(\sqrt{21},-\sqrt{7},0)$ and $(0,\sqrt{2},2)$
  
  Points in cylindrical coordinates:
  $(\frac{1}{2},\frac{\pi}{2},1)$ and $(6,\frac{5\pi}{3},\sqrt{7})$
  
  Points in spherical coordinates:
  $(3,\frac{3\pi}{2},\frac{5\pi}{6})$ and $(2,\frac{7\pi}{6},\frac{3\pi}{4})$

In Exercises 9–10, describe the curve, surface or region in space determined by the given bounds.

9.
Bounds in cylindrical coordinates:
- $r=1$ , $0\leq\theta\leq 2\pi$ , $0\leq z\leq 1$
  
  $1\leq r\leq 2$ , $0\leq\theta\leq\pi$ , $0\leq z\leq 1$
Bounds in spherical coordinates:
- $\rho=3$ , $0\leq\theta\leq 2\pi$ , $0\leq\varphi\leq\pi/2$
  
  $2\leq\rho\leq 3$ , $0\leq\theta\leq 2\pi$ , $0\leq\varphi\leq\pi$
10.
Bounds in cylindrical coordinates:
- $1\leq r\leq 2$ , $\theta=\pi/2$ , $0\leq z\leq 1$
  
  $r=2$ , $0\leq\theta\leq 2\pi$ , $z=5$
Bounds in spherical coordinates:
- $0\leq\rho\leq 2$ , $0\leq\theta\leq\pi$ , $\varphi=\pi/4$
  
  $\rho=2$ , $0\leq\theta\leq 2\pi$ , $\varphi=\pi/6$

In Exercises 11–12, standard regions in space, as defined by cylindrical and spherical coordinates, are shown. Set up the triple integral that integrates the given function over the graphed region.

11.

Cylindrical coordinates, integrating $h(r,\theta,z)$ :

(fullscreen)
12.

Spherical coordinates, integrating $h(\rho,\theta,\varphi)$ :

(fullscreen)

In Exercises 13–18, a triple integral in cylindrical coordinates is given. Describe the region in space defined by the bounds of the integral.

13.

$\displaystyle\int_{0}^{\pi/2}\int_{0}^{2}\int_{0}^{2}r\operatorname{d}\!z% \operatorname{d}\!r\operatorname{d}\!\theta$
14.

$\displaystyle\int_{0}^{2\pi}\int_{3}^{4}\int_{0}^{5}r\operatorname{d}\!z% \operatorname{d}\!r\operatorname{d}\!\theta$
15.

$\displaystyle\int_{0}^{2\pi}\int_{0}^{1}\int_{0}^{1-r}r\operatorname{d}\!z% \operatorname{d}\!r\operatorname{d}\!\theta$
16.

$\displaystyle\int_{0}^{\pi}\int_{0}^{1}\int_{0}^{2-r}r\operatorname{d}\!z% \operatorname{d}\!r\operatorname{d}\!\theta$
17.

$\displaystyle\int_{0}^{\pi}\int_{0}^{3}\int_{0}^{\sqrt{9-r^{2}}}r\operatorname% {d}\!z\operatorname{d}\!r\operatorname{d}\!\theta$
18.

$\displaystyle\int_{0}^{2\pi}\int_{0}^{a}\int_{0}^{\sqrt{a^{2}-r^{2}}+b}r% \operatorname{d}\!z\operatorname{d}\!r\operatorname{d}\!\theta$

In Exercises 19–24, a triple integral in spherical coordinates is given. Describe the region in space defined by the bounds of the integral.

19.

$\displaystyle\int_{0}^{\pi/2}\int_{0}^{\pi}\int_{0}^{1}\rho^{2}\sin(\varphi)% \operatorname{d}\!\rho\operatorname{d}\!\varphi\operatorname{d}\!\theta$
20.

$\displaystyle\int_{0}^{\pi}\int_{0}^{\pi}\int_{1}^{1.1}\rho^{2}\sin(\varphi)% \operatorname{d}\!\rho\operatorname{d}\!\varphi\operatorname{d}\!\theta$
21.

$\displaystyle\int_{0}^{2\pi}\int_{0}^{\pi/4}\int_{0}^{2}\rho^{2}\sin(\varphi)% \operatorname{d}\!\rho\operatorname{d}\!\varphi\operatorname{d}\!\theta$
22.

$\displaystyle\int_{0}^{2\pi}\int_{\pi/6}^{\pi/4}\int_{0}^{2}\rho^{2}\sin(% \varphi)\operatorname{d}\!\rho\operatorname{d}\!\varphi\operatorname{d}\!\theta$
23.

$\displaystyle\int_{0}^{2\pi}\int_{0}^{\pi/6}\int_{0}^{\sec\varphi}\rho^{2}\sin% (\varphi)\operatorname{d}\!\rho\operatorname{d}\!\varphi\operatorname{d}\!\theta$
24.

$\displaystyle\int_{0}^{2\pi}\int_{0}^{\pi/6}\int_{0}^{a\sec\varphi}\rho^{2}% \sin(\varphi)\operatorname{d}\!\rho\operatorname{d}\!\varphi\operatorname{d}\!\theta$

In Exercises 25–28, a solid is described along with its density function. Find the mass of the solid using cylindrical coordinates.

25.

Bounded by the cylinder $x^{2}+y^{2}=4$ and the planes $z=0$ and $z=4$ with density function $\delta(x,y,z)=\sqrt{x^{2}+y^{2}}+1$ .
26.

Bounded by the cylinders $x^{2}+y^{2}=4$ and $x^{2}+y^{2}=9$ , between the planes $z=0$ and $z=10$ with density function $\delta(x,y,z)=z$ .
27.

Bounded by $y\geq 0$ , the cylinder $x^{2}+y^{2}=1$ , and between the planes $z=0$ and $z=4-y$ with density function $\delta(x,y,z)=1$ .
28.

The upper half of the unit ball, bounded between $z=0$ and $z=\sqrt{1-x^{2}-y^{2}}$ , with density function $\delta(x,y,z)=1$ .

In Exercises 29–32, a solid is described along with its density function. Find the center of mass of the solid using cylindrical coordinates. (Note: these are the same solids and density functions as found in Exercises 25–28.)

29.

Bounded by the cylinder $x^{2}+y^{2}=4$ and the planes $z=0$ and $z=4$ with density function $\delta(x,y,z)=\sqrt{x^{2}+y^{2}}+1$ .
30.

Bounded by the cylinders $x^{2}+y^{2}=4$ and $x^{2}+y^{2}=9$ , between the planes $z=0$ and $z=10$ with density function $\delta(x,y,z)=z$ .
31.

Bounded by $y\geq 0$ , the cylinder $x^{2}+y^{2}=1$ , and between the planes $z=0$ and $z=4-y$ with density function $\delta(x,y,z)=1$ .
32.

The upper half of the unit ball, bounded between $z=0$ and $z=\sqrt{1-x^{2}-y^{2}}$ , with density function $\delta(x,y,z)=1$ .

In Exercises 33–36, a solid is described along with its density function. Find the mass of the solid using spherical coordinates.

33.

The upper half of the unit ball, bounded between $z=0$ and $z=\sqrt{1-x^{2}-y^{2}}$ , with density function $\delta(x,y,z)=1$ .
34.

The spherical shell bounded between $x^{2}+y^{2}+z^{2}=16$ and $x^{2}+y^{2}+z^{2}=25$ with density function $\delta(x,y,z)=\sqrt{x^{2}+y^{2}+z^{2}}$ .
35.

The conical region bounded above $z=\sqrt{x^{2}+y^{2}}$ and below the sphere $x^{2}+y^{2}+z^{2}=1$ with density function $\delta(x,y,z)=z$ .
36.

The cone bounded above $z=\sqrt{x^{2}+y^{2}}$ and below the plane $z=1$ with density function $\delta(x,y,z)=z$ .

In Exercises 37–40, a solid is described along with its density function. Find the center of mass of the solid using spherical coordinates. (Note: these are the same solids and density functions as found in Exercises 33–36.)

37.

The upper half of the unit ball, bounded between $z=0$ and $z=\sqrt{1-x^{2}-y^{2}}$ , with density function $\delta(x,y,z)=1$ .
38.

The spherical shell bounded between $x^{2}+y^{2}+z^{2}=16$ and $x^{2}+y^{2}+z^{2}=25$ with density function $\delta(x,y,z)=\sqrt{x^{2}+y^{2}+z^{2}}$ .
39.

The conical region bounded above $z=\sqrt{x^{2}+y^{2}}$ and below the sphere $x^{2}+y^{2}+z^{2}=1$ with density function $\delta(x,y,z)=z$ .
40.

The cone bounded above $z=\sqrt{x^{2}+y^{2}}$ and below the plane $z=1$ with density function $\delta(x,y,z)=z$ .

In Exercises 41–44, a region is space is described. Set up the triple integrals that find the volume of this region using rectangular, cylindrical and spherical coordinates, then comment on which of the three appears easiest to evaluate.

41.

The region enclosed by the unit sphere, $x^{2}+y^{2}+z^{2}=1$ .
42.

The region enclosed by the cylinder $x^{2}+y^{2}=1$ and planes $z=0$ and $z=1$ .
43.

The region enclosed by the cone $z=\sqrt{x^{2}+y^{2}}$ and plane $z=1$ .
44.

The cube enclosed by the planes $x=0$ , $x=1$ , $y=0$ , $y=1$ , $z=0$ and $z=1$ . (Hint: in spherical, use order of integration $\operatorname{d}\!\rho\operatorname{d}\!\varphi\operatorname{d}\!\theta$ .)

45.

Find the center and radius of the sphere given by the spherical equation

$\rho=4\sin\phi\cos\theta+6\sin\phi\sin\theta-2\cos\phi.$
46.

Show that for $a\neq 0$ , the equation $\rho=2a\sin\phi\,\cos\theta$ in spherical coordinates describes a sphere centered at $(a,0,0)$ with radius $\left\lvert a\right\rvert$ .
47.

Let $P=(a,\theta,\phi)$ be a point in spherical coordinates, with $a>0$ and $0<\phi<\pi$ . Then $P$ lies on the sphere $\rho=a$ . Since $0<\phi<\pi$ , the line segment from the origin to $P$ can be extended to intersect the cylinder given by $r=a$ (in cylindrical coordinates). Find the cylindrical coordinates of that point of intersection.
48.

Let $P_{1}$ and $P_{2}$ be points whose spherical coordinates are $(\rho_{1},\theta_{1},\phi_{1})$ and $(\rho_{2},\theta_{2},\phi_{2})$ , respectively. Let $\vec{v}_{1}$ be the vector from the origin to $P_{1}$ , and let $\vec{v}_{2}$ be the vector from the origin to $P_{2}$ . For the angle $\gamma$ between $\vec{v}_{1}$ and $\vec{v}_{2}$ , show that

$\cos\gamma=\cos\phi_{1}\cos\phi_{2}+\sin\phi_{1}\sin\phi_{2}\cos(\theta_{2}-% \theta_{1}).$

This formula is used in electrodynamics to prove the addition theorem for spherical harmonics, which provides a general expression for the electrostatic potential at a point due to a unit charge.
49.

Show that the distance $d$ between the points $P_{1}$ and $P_{2}$ with cylindrical coordinates $(r_{1},\theta_{1},z_{1})$ and $(r_{2},\theta_{2},z_{2})$ , respectively, is

$d=\sqrt{r_{1}^{2}+r_{2}^{2}-2r_{1}\,r_{2}\cos(\,\theta_{2}-\theta_{1}\,)+(z_{2% }-z_{1})^{2}}\,.$
50.

Show that the distance $d$ between the points $P_{1}$ and $P_{2}$ with spherical coordinates $(\rho_{1},\theta_{1},\phi_{1})$ and $(\rho_{2},\theta_{2},\phi_{2})$ , respectively, is $d=$

$\sqrt{\rho_{1}^{2}+\rho_{2}^{2}-2\rho_{1}\,\rho_{2}[\sin\phi_{1}\,\sin\phi_{2}% \,\cos(\,\theta_{2}-\theta_{1}\,)+\cos\phi_{1}\,\cos\phi_{2}]}\,.$
51.

Prove Theorem 14.7.1 by finding the Jacobian of the cylindrical coordinate transformation.
52.

Prove Theorem 14.7.2 by finding the Jacobian of the spherical coordinate transformation.
53.

Evaluate $\displaystyle\iint_{R}\sin\left(\frac{x+y}{2}\right)\,\cos\left(\frac{x-y}{2}% \right)\operatorname{d}\!A$ , where $R$ is the triangle with vertices $(0,0)$ , $(2,0)$ and $(1,1)$ . (Hint: Use the change of variables $u=(x+y)/2$ , $v=(x-y)/2$ .)