Section14.7Triple 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 desribing surfaces and regions in space.
Subsection14.7.1Cylindrical 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)\text{,}\) given in rectangular coordinates, with the point \((r_0,\theta_0,z_0)\text{,}\) given in cylindrical coordinates, where the \(z\)-value in both systems is the same, and the point \((x_0,y_0)\) in the \(xy\)-plane is identified with the polar point \(P(r_0,\theta_0)\text{;}\) see Figure 14.7.2. 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\leq 2\pi\text{.}\)
We use the identity \(z=z\) along with the identities found in Key Idea 10.4.6 to convert between the rectangular coordinate \((x,y,z)\) and the cylindrical coordinate \((r,\theta,z)\text{,}\) namely:
\begin{equation*}
\begin{array}{l} \text{ From rectangular to cylindrical: } r=\sqrt{x^2+y^2}, \tan(\theta) = y/x \text{ and } z=z;\\
\text{ From cylindrical to rectangular: } x=r\cos(\theta), y=r\sin(\theta) \text{ and } z=z. \end{array}
\end{equation*}
These identities, along with conversions related to spherical coordinates, are given later in Key Idea 14.7.16.
Example14.7.4.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.
Following the identities given above (and, later in Key Idea 14.7.16), we have \(r = \sqrt{2^2+(-2)^2} = 2\sqrt{2}\text{.}\) Using \(\tan(\theta) = y/x\text{,}\) we find \(\theta = \tan^{-1}(-2/2) =-\pi/4\text{.}\) As we restrict \(\theta\) to being between \(0\) and \(2\pi\text{,}\) we set \(\theta = 7\pi/4\text{.}\) Finally, \(z = 1\text{,}\) giving the cylindrical point \((2\sqrt2,7\pi/4,1)\text{.}\)
In converting the cylindrical point \((4,3\pi/4,5)\) to rectangular, we have \(x = 4\cos\big(3\pi/4\big) = -2\sqrt{2}\text{,}\)\(y = 4\sin\big(3\pi/4\big) = 2\sqrt{2}\) and \(z=5\text{,}\) giving the rectangular point \((-2\sqrt{2},2\sqrt{2},5)\text{.}\)
Setting each of \(r\text{,}\)\(\theta\) and \(z\) equal to a constant defines a surface in space, as illustrated in the following example.
Example14.7.5.Canonical surfaces in cylindrical coordinates.
Describe the surfaces \(r=1\text{,}\)\(\theta = \pi/3\) and \(z=2\text{,}\) given in cylindrical coordinates.
The 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.17 (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\text{,}\) as given by polar coordinates in the \(xy\)-plane, parallel to the \(z\)-axis.
The equation \(z=2\) describes the plane of all points in space that are 2 units above the \(xy\)-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.6. Note how their intersection uniquely defines the point \(P=(1,\pi/3,2)\text{.}\)
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.29 shows how to evaluate \(\iiint_Dh(x,y,z)\, dV\) using rectangular coordinates. In that evaluation, we use \(dV = dz\,dy\,dx\) (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 \(xy\)-plane. We could describe \(R\) using polar coordinates as done in Section 14.3. In that section, we saw how we used \(dA = r\,dr\,d\theta\) instead of \(dA = dy\,dx\text{.}\)
Considering the above thoughts, we have \(dV = dz\big(r\,dr\,d\theta\big) = r\,dz\,dr\,d\theta\text{.}\) 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\text{,}\) respectively. Finally, using the identities given above, we change the integrand \(h(x,y,z)\) to \(h(r,\theta,z)\text{.}\)
This process should sound plausible; the following theorem states it is truly a way of evaluating a triple integral.
Theorem14.7.7.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\text{,}\)\(g_1(\theta)\leq r \leq g_2(\theta)\) and \(f_1(r,\theta) \leq z \leq f_2(r,\theta)\text{.}\) Then
\begin{equation*}
\iiint_D h(r,\theta,z)\, dV = \int_\alpha^\beta\int_{g_1(\theta)}^{g_2(\theta)}\int_{f_1(r,\theta)}^{f_2(r,\theta)}h(r,\theta,z) r\,dz\,dr\,d\theta\text{.}
\end{equation*}
Example14.7.9.Evaluating a triple integral with cylindrical coordinates.
Find the mass of the solid represented by the region in space bounded by \(z=0\text{,}\)\(z=\sqrt{4-x^2-y^2}+3\) and the cylinder \(x^2+y^2=4\) (as shown in Figure 14.7.10), with density function \(\delta(x,y,z) = x^2+y^2+z+1\text{,}\) using a triple integral in cylindrical coordinates. Distances are measured in centimeters and density is measured in grams per cubic centimeter.
We 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}\text{,}\) we convert the hemisphere of radius 2 to the equation \(z=\sqrt{4-r^2}\text{;}\) the cylinder \(x^2+y^2=4\) is converted to \(r^2=4\text{,}\) or, more simply, \(r=2\text{.}\) We also convert the density function: \(\delta(r,\theta,z) = r^2+z+1\text{.}\)
To describe this solid with the bounds of a triple integral, we bound \(z\) with \(0\leq z\leq \sqrt{4-r^2}+3\text{;}\) we bound \(r\) with \(0 \leq r \leq 2\text{;}\) we bound \(\theta\) with \(0 \leq \theta \leq 2\pi\text{.}\)
where we leave the details of the remaining double integral to the reader.
Example14.7.11.Finding the center of mass using cylindrical coordinates.
Find 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\text{,}\) where distances are measured in feet, as seen in Figure 14.7.12. (The volume of this solid was found in Example 14.3.13.)
We convert the equation of the plane to use cylindrical coordinates: \(z= 1-r\cos(\theta)+0.1r\sin(\theta)\text{.}\) Thus the region is space is bounded by \(0 \leq z \leq 1-r\cos(\theta) + 0.1r\sin(\theta)\text{,}\)\(0 \leq r \leq \cos(3\theta)\text{,}\)\(0 \leq \theta \leq \pi\) (recall that the rose curve \(r=\cos(3\theta)\) is traced out once on \([0,\pi]\text{.}\)
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.13.
\begin{equation*}
M = \iiint_D\delta \, dV = \int_0^{\pi}\int_0^{\cos(3\theta)}\int_0^{1-r\cos(\theta)+0.1r\sin(\theta)} r\,dz\,dr\,d\theta = \frac{\pi}{4}\text{.}
\end{equation*}
From Definition 14.6.30 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):
The center of mass in rectangular coordinates, found by dividing the respective moments by the mass, is approximately located at \((-0.188,0.019,0.595)\text{,}\) which lies outside the bounds of the solid.
Subsection14.7.2Spherical 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)\text{,}\) where \(\rho\) is the distance from the origin to \(P\text{,}\)\(\theta\) is the same angle as would be used to describe \(P\) in the cylindrical coordinate system, and \(\varphi\) is the angle between the \(xy\)-plane and the ray from the origin to \(P\text{;}\) see Figure 14.7.14. So that each point in space that does not lie on the \(z\)-axis is defined uniquely, we will restrict \(\rho \geq 0\text{,}\)\(0 \leq \theta \leq 2\pi\) and \(-\pi/2 \leq \varphi \leq \pi/2\text{.}\)
Convention14.7.13.
Note that most mathematics textbooks define \(\varphi\) to be measured from the positive \(z\)-axis, with values in \([0,\pi]\text{,}\) rather than from the \(xy\)-plane.
We have chosen our convention with a number of considerations in mind:
The coordinates \((\rho,\theta,\varphi)\) form a right-handed coordinate system: one in which the orientation matches that of our usual \((x,y,z)\) coordinates, where the “right-hand rule” applies. If \(\varphi\) is measured from the \(z\)-axis, the order \((\rho, \varphi, \theta)\) is needed to get a right-handed system.
Points of the form \((a,\alpha,0)\) are the same in both cylindrical and spherical coordinates.
Some integration problems become slightly easier: we will see soon that the volume element in spherical coordinates involves \(\cos(\varphi)\text{,}\) which integrates to \(\sin(\varphi)\text{.}\) In the usual convention, the volume element involves \(\sin(\varphi)\text{,}\) which integrates to \(-\cos(\varphi)\) – a source of many common sign errors.
Students of Physics will encounter yet another convention. In Physics, the variable \(r\) is preferred as the radial coordinate, and spherical coordinates are given as \((r,\theta,\varphi)\text{;}\) however, in Physics, \(\varphi\) becomes the angle in the \(xy\)-plane, while \(\theta\) is the angle measured from the positive \(z\)-axis.
Note that the angle in the \(xy\)-plane (\(\theta\text{,}\) in our case) is known as the azimuthal angle. Our angle \(\varphi\) is known as the elevation angle. The angle used in other conventions that is measured from the positive \(z\)-axis (often identified with the north pole) is known as the polar angle. For further discussion, the Wikipedia article 1 is quite useful.
The following Key Idea gives conversions to/from our three spatial coordinate systems.
Key Idea14.7.16.Converting Between Rectangular, Cylindrical and Spherical Coordinates.
Example14.7.17.Converting between rectangular and spherical coordinates.
Convert the rectangular point \((2,-2,1)\) to spherical coordinates, and convert the spherical point \((6,\pi/3,0)\) to rectangular and cylindrical coordinates.
This rectangular point is the same as used in Example 14.7.4. Using Key Idea 14.7.16, we find \(\rho = \sqrt{2^2+(-1)^2+1^2} = 3\text{.}\) Using the same logic as in Example 14.7.4, we find \(\theta = 7\pi/4\text{.}\) Finally, \(\sin(\varphi) = 1/3\text{,}\) giving \(\varphi = \sin^{-1}(1/3) \approx 0.34\text{,}\) or about \(19.47^\circ\text{.}\) Thus the spherical coordinates are approximately \((3,7\pi/4,0.34)\text{.}\)
Converting the spherical point \((6,\pi/3,0)\) to rectangular, we have \(x = 6\cos(0)\cos(\pi/3) = 3\text{,}\)\(y = 6\cos(0)\sin(\pi/3) = 3\sqrt{3}\) and \(z = 6\sin(0) = 0\text{.}\) Thus the rectangular coordinates are \((3,3\sqrt{3},0)\text{.}\)
To convert this spherical point to cylindrical, we have \(r = 6\cos(0) = 6\text{,}\)\(\theta = \pi/3\) and \(z = 6\sin(0) =0\text{,}\) giving the cylindrical point \((6,\pi/3,0)\text{.}\)
Example14.7.18.Canonical surfaces in spherical coordinates.
Describe the surfaces \(\rho=1\text{,}\)\(\theta = \pi/3\) and \(\varphi = \pi/3\text{,}\) given in spherical coordinates.
The 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 \(xy\)-plane as given by polar coordinates, extend the line parallel to the \(z\)-axis, forming a plane.
The equation \(\varphi=\pi/3\) describes all points \(P\) in space where the ray from the origin to \(P\) makes an angle of \(\pi/3\) with the \(xy\)-plane. 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.19. Note how their intersection uniquely defines the point \(P=(1,\pi/3,\pi/6)\text{.}\)
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,” \(dV\text{,}\) used in the integral.
Considering Figure 14.7.20, we can make a small “spherical wedge” by varying \(\rho\text{,}\)\(\theta\) and \(\varphi\) each a small amount, \(\Delta\rho\text{,}\)\(\Delta\theta\) and \(\Delta\varphi\text{,}\) respectively. This wedge is approximately a rectangular solid when the change in each coordinate is small, giving a volume of about
\begin{equation*}
\Delta V \approx \Delta\rho\ \times\ \rho\Delta\varphi\ \times\ \rho\cos(\varphi)\Delta\theta\text{.}
\end{equation*}
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\text{,}\)\(\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
Again, this development of \(dV\) should sound reasonable, and the following theorem states it is the appropriate manner by which triple integrals are to be evaluated in spherical coordinates.
Theorem14.7.21.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\text{,}\)\(\beta_1 \leq \theta \leq \beta_2\) and \(f_1(\theta,\varphi) \leq \rho \leq f_2(\theta,\varphi)\text{.}\) Then
Example14.7.22.Establishing the volume of a sphere.
Let \(D\) be the region in space bounded by the sphere, centered at the origin, of radius \(r\text{.}\) Use a triple integral in spherical coordinates to find the volume \(V\) of \(D\text{.}\)
The sphere of radius \(r\text{,}\) centered at the origin, has equation \(\rho = r\text{.}\) To obtain the full sphere, the bounds on \(\theta\) and \(\varphi\) are \(0\leq \theta \leq 2\pi\) and \(-\pi/2 \leq \varphi \leq \pi/2\text{.}\) This leads us to:
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.
Example14.7.24.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/3\text{,}\) as illustrated in Figure 14.7.25.
We will set up the four triple integrals needed to find the center of mass (i.e., to compute \(M\text{,}\)\(M_{yz}\text{,}\)\(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\text{,}\) we use the following bounds: \(0 \leq \rho \leq 4\text{,}\)\(0 \leq \theta \leq 2\pi\) and \(0 \leq \varphi \leq \pi/3\text{.}\) Since density \(\delta\) is constant, we assume \(\delta =1\text{.}\)
Thus the center of mass is \((0,0,M_{xy}/M) \approx (0,0,2.799)\text{,}\) as indicated in Figure 14.7.25.
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.”
Exercises14.7.3Exercises
Terms and Concepts
1.
Explain the difference between the roles \(r\text{,}\) in cylindrical coordinates, and \(\rho\text{,}\) 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
Exercise Group.
In the following exercises, points are given in either the rectangular, cylindrical or spherical coordinate systems. Find the coordinates of the points in the other systems.
5.
(a)
Points in rectangular coordinates: \((2,2,1)\) and \((-\sqrt{3},1,0)\)
(b)
Points in cylindrical coordinates: \((2,\pi/4,2)\) and \((3,3\pi/2,-4)\)
(c)
Points in spherical coordinates: \((2,\pi/4,\pi/4)\) and \((1,0,0)\)
6.
(a)
Points in rectangular coordinates: \((0,1,1)\) and \((-1,0,1)\)
(b)
Points in cylindrical coordinates: \((0,\pi,1)\) and \((2,4\pi/3,0)\)
(c)
Points in spherical coordinates: \((2,\pi/6,0)\) and \((3,\pi,-\pi/2)\)
Exercise Group.
In the following exercises, describe the curve, surface or region in space determined by the given bounds in cylindrical coordinates.
This is a curve, a circle of radius 1 centered at \((0,0,\sqrt 3)\text{,}\) lying parallel to the \(xy\)-plane.
Exercise Group.
In the following exercises, 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.
In the following exercises, 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\text{.}\)
26.
Bounded by the cylinders \(x^2+y^2=4\) and \(x^2+y^2=9\text{,}\) between the planes \(z=0\) and \(z=10\) with density function \(\delta(x,y,z) =z\text{.}\)
27.
Bounded by \(y\geq 0\text{,}\) the cylinder \(x^2+y^2=1\text{,}\) and between the planes \(z=0\) and \(z=4-y\) with density function \(\delta(x,y,z) =1\text{.}\)
28.
The upper half of the unit ball, bounded between \(z= 0\) and \(z=\sqrt{1-x^2-y^2}\text{,}\) with density function \(\delta(x,y,z) =1\text{.}\)
Exercise Group.
In the following exercises, 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 14.7.3.25 through Exercise 14.7.3.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\text{.}\)
30.
Bounded by the cylinders \(x^2+y^2=4\) and \(x^2+y^2=9\text{,}\) between the planes \(z=0\) and \(z=10\) with density function \(\delta(x,y,z) =z\text{.}\)
31.
Bounded by \(y\geq 0\text{,}\) the cylinder \(x^2+y^2=1\text{,}\) and between the planes \(z=0\) and \(z=4-y\) with density function \(\delta(x,y,z) =1\text{.}\)
32.
The upper half of the unit ball, bounded between \(z= 0\) and \(z=\sqrt{1-x^2-y^2}\text{,}\) with density function \(\delta(x,y,z) =1\text{.}\)
Exercise Group.
In the following exercises, 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}\text{,}\) with density function \(\delta(x,y,z) =1\text{.}\)
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}\text{.}\)
35.
The conical region bounded below by \(z=\sqrt{x^2+y^2}\) and above by the sphere \(x^2+y^2+z^2=1\) with density function \(\delta(x,y,z) = z\text{.}\)
36.
The cone that lies above the cone \(z=\sqrt{x^2+y^2}\) and below the plane \(z=1\) with density function \(\delta(x,y,z) = z\text{.}\)
Exercise Group.
In the following exercises, 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 14.7.3.33 through Exercise 14.7.3.36.)
37.
The upper half of the unit ball, bounded between \(z= 0\) and \(z=\sqrt{1-x^2-y^2}\text{,}\) with density function \(\delta(x,y,z) =1\text{.}\)
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}\text{.}\)
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\text{.}\)
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\text{.}\)
Exercise Group.
In the following exercises, 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\text{.}\)
42.
The region enclosed by the cylinder \(x^2+y^2=1\) and planes \(z=0\) and \(z=1\text{.}\)
43.
The region enclosed by the cone \(z=\sqrt{x^2+y^2}\) and plane \(z=1\text{.}\)
44.
The cube enclosed by the planes \(x=0\text{,}\)\(x=1\text{,}\)\(y=0\text{,}\)\(y=1\text{,}\)\(z=0\) and \(z=1\text{.}\) (Hint: in spherical, use order of integration \(d\rho\, d\varphi\, d\theta\text{.}\))