Section 4.2 Double Integrals: Volume and Average Value
¶Consider a surface \(f(x,y)\text{.}\) In this section, we are interested in computing either the volume under \(f\) or the average function value of \(f\) over a certain area in the \(x\)\(y\)plane. You might temporarily think of this surface as representing physical topography—a hilly landscape, perhaps. What is the average height of the surface (or average altitude of the landscape) over some region?
Subsection 4.2.1 Volume and Average Value over a Rectangular Region
As with most such problems, we start by thinking about how we might approximate the answer. Suppose the region is a rectangle, \([a,b]\times[c,d]\text{.}\) We can divide the rectangle into a grid, \(m\) subdivisions in one direction and \(n\) in the other, as indicated in the 2D graph below, where \(m=6\) and \(n=4\text{.}\) We pick \(x\)values \(x_0\text{,}\) \(x_1\text{,}\)…, \(x_{m1}\) in each subdivision in the \(x\)direction, and similarly in the \(y\)direction. At each of the points \((x_i,y_j)\) in one of the smaller rectangles in the grid, we compute the height of the surface: \(f(x_i,y_j)\text{.}\) Now the average of these heights should be (depending on the fineness of the grid) close to the average height of the surface:
As both \(m\) and \(n\) go to infinity, we expect this approximation to converge to a fixed value, the actual average height of the surface. For reasonably nice functions this does indeed happen.
Using sigma notation, we can rewrite the approximation:
The two parts of this product have useful meaning: \((ba)(dc)\) is of course the area of the rectangle, and the double sum adds up \(mn\) terms of the form \(f(x_j,y_i)\Delta x\Delta y\text{,}\) which is the height of the surface at a point multiplied by the area of one of the small rectangles into which we have divided the large rectangle. In short, each term \(f(x_j,y_i)\Delta x\Delta y\) is the volume of a tall, thin, rectangular box, and is approximately the volume under the surface and above one of the small rectangles; see Figure 4.1. When we add all of these up, we get an approximation to the volume under the surface and above the rectangle \(R=[a,b]\times[c,d]\text{.}\) When we take the limit as \(m\) and \(n\) go to infinity, the double sum becomes the actual volume under the surface, which we divide by \((ba)(dc)\) to get the average height.
Double sums like this come up in many applications, so in a way it is the most important part of this example; dividing by \((ba)(dc)\) is a simple extra step that allows the computation of an average. This computation is independent of \(f\) being positive, and so as we did in the single variable case, we introduce a special notation for the limit of such a double sum, which is referred to as the double integral of \(f\) over the region \(R\text{.}\)
Definition 4.11. Double Integral over a Rectangular Region.
Let \(f(x,y)\) be a continuous function defined over the rectangle \(R=[a,b]\times[c,d]\text{,}\) then the double integral of \(f\) over \(R\) is denoted by
provided the limit exists.
Note: The notation \(dA\) indicates a small element of area, without specifying any particular order for the variables \(x\) and \(y\text{;}\) it is shorter and more generic than writing \(dx\,dy\text{.}\)
We now capture our results from the earlier calculations using the notation of the double integral.
Theorem 4.12. Average Value of a Twovariable Function.
Let \(f(x,y)\) be a continuous function defined over the rectangle \(R=[a,b]\times[c,d]\text{,}\) then the \thmfont{average value \(f_{avg}\)} of \(f\) over \(R\) is
provided the double integral exists.
Theorem 4.13. Volume Below a Surface.
Let \(f(x,y)\) be a continuous function defined over the rectangle \(R=[a,b]\times[c,d]\text{,}\) with \(f(x,y) \geq 0\text{,}\) then the \thmfont{volume \(V\) below the surface} \(f\) is
provided the double integral exists.
Subsection 4.2.2 Computing Double Integrals over Rectangular Regions
The next question, of course, is: How do we compute these double integrals? You might think that we will need some twodimensional version of the Fundamental Theorem of Calculus, but as it turns out we can get away with just the single variable version, applied twice.
Going back to the double sum, we can rewrite it to emphasize a particular order in which we want to add the terms:
In the sum in parentheses, only the value of \(x_j\) is changing; \(y_i\) is temporarily constant. As \(m\) goes to infinity, this sum has the right form to turn into an integral:
So after we take the limit as \(m\) goes to infinity, the sum is
Of course, for different values of \(y_i\) this integral has different values; in other words, it is really a function applied to \(y_i\text{:}\)
If we substitute back into the sum we get
This sum has a nice interpretation. The value \(G(y_i)\) is the area of a cross section of the region under the surface \(f(x,y)\text{,}\) namely, when \(y=y_i\text{.}\) The quantity \(G(y_i)\Delta y\) can be interpreted as the volume of a solid with face area \(G(y_i)\) and thickness \(\Delta y\text{.}\) Think of the surface \(f(x,y)\) as the top of a loaf of sliced bread. Each slice has a crosssectional area and a thickness; \(G(y_i)\Delta y\) corresponds to the volume of a single slice of bread. Adding these up approximates the total volume of the loaf. (This is very similar to the technique we used to compute volumes in Section 3.3, except that there we need the crosssections to be in some way “the same” .) Figure 4.2 shows this “sliced loaf” approximation using the same surface as shown in Figure 4.1. Nicely enough, this sum looks just like the sort of sum that turns into an integral, namely,
Let's be clear about what this means: we first will compute the inner integral, temporarily treating \(y\) as a constant. We will do this by finding an antiderivative with respect to \(x\text{,}\) then substituting \(x=a\) and \(x=b\) and subtracting, as usual. The result will be an expression with no \(x\) variable but some occurrences of \(y\text{.}\) Then the outer integral will be an ordinary onevariable problem, with \(y\) as the variable.
Example 4.14. Average Height.
Figure 4.1 shows the function \(\sin(xy)+6/5\) on \([0.5,3.5]\times[0.5,2.5]\text{.}\) Find the average height of this surface.
The volume under this surface is
The inner integral is
Unfortunately, this gives a function for which we can't find a simple antiderivative. To complete the problem we could use Sage or similar software to approximate the integral. Doing this gives a volume of approximately \(8.84\text{,}\) so the average height is approximately \(8.84/6\approx 1.47\text{.}\)
Because addition and multiplication are commutative and associative, we can rewrite the original double sum:
Now if we repeat the development above, the inner sum turns into an integral:
and then the outer sum turns into an integral:
In other words, we can compute the integrals in either order, first with respect to \(x\) then \(y\text{,}\) or vice versa. Thinking of the loaf of bread, this corresponds to slicing the loaf in a direction perpendicular to the first.
We summarize our findings by providing a general guideline for how the double integral over a rectangle, such as the one shown below, is computed.
Guideline for Computing a Double Integral over a Rectangle.
Let \(f(x,y)\) be a continuous function defined over the rectangle \(R=[a,b]\times[c,d]\text{.}\) Follow these steps to evaluate the double integral
Choose the double integral, \(\ds{\int_c^d \int_a^b f(x,y)\,dx\,dy}\) or \(\ds{\int_a^b \int_c^d f(x,y)\,dy\,dx}\text{.}\)

Compute the inner integral:
If \(\ds \int_a^b f(x,y)\,dx\text{,}\) then temporarily treat \(y\) as a constant. Find an antiderivative with respect to \(x\text{,}\) then evaluate over the integration bounds \(x = a\) and \(x = b\text{.}\)
If \(\ds \int_c^d f(x,y)\,dy\text{,}\) then temporarily treat \(x\) as a constant. Find an antiderivative with respect to \(y\text{,}\) then evaluate over the integration bounds \(y=c\) and \(y=d\text{.}\)

Compute the outer integral, which will be an ordinary onevariable integral of the form
\begin{equation*} \int_c^d G(y)\, dy \text{ or } \int_a^b G(x)\,dx\text{.} \end{equation*}
We haven't really proved that the value of a double integral is equal to the value of the corresponding two single integrals in either order of integration, but provided the function is continuous, this is true. The result is called Fubini's Theorem, which we state here without proof.
Theorem 4.15. Fubini's Theorem for Rectangular Regions of Integration.
Let \(f(x,y)\) be a continuous function defined over the rectangle \(R=[a,b]\times[c,d]\text{.}\) Then the order of integration does not matter, and so
We provide one example, where we compute the volume under a surface in two ways by switching the order of integration.
Example 4.16. Compute Volume in Two Ways.
We compute \(\ds\iint_R 1+(x1)^2+4y^2\,dA\text{,}\) where \(R=[0,3]\times[0,2]\text{,}\) in two ways.
First,
In the other order:
Subsection 4.2.3 Computing Double Integrals over any Region
Frequently we will be interested in a region that is not simply a rectangle. Let's compute the volume under the surface \(x+2y^2\) above the region described by \(0\le x\le1\) and \(0\le y\le x^2\text{,}\) shown below.
In principle there is nothing more difficult about this problem. If we imagine the threedimensional region under the surface and above the parabolic region as an oddly shaped loaf of bread, we can still slice it up, approximate the volume of each slice, and add these volumes up. For example, if we slice perpendicular to the \(x\)axis at \(x_i\text{,}\) the thickness of a slice will be \(\Delta x\) and the area of the slice will be
When we add these up and take the limit as \(\Delta x\) goes to 0, we get the double integral
We could just as well slice the solid perpendicular to the \(y\)axis, in which case we get
What is the average height of the surface over this region? As before, it is the volume divided by the area of the base, but now we need to use integration to compute the area of the base, since it is not a simple rectangle. The area is
so the average height is
Although we have not proven that the order of integration can be switched, we nonetheless capture our results and state the general version of Fubini's Theorem without proof.
Theorem 4.17. Fubini's Theorem for General Regions of Integration.
Let \(f(x,y)\) be a continuous function defined over the region \(R\text{.}\)

If \(R\) is bounded by \(a \leq x \leq b,\ g_1(x) \leq y \leq g_2(x)\text{,}\) with \(g_1\) and \(g_2\) continuous on \([a,b]\text{,}\) then
\begin{equation*} \iint_R f(x,y)\,dA = \int_a^b\int_{g_1(x)}^{g_2(x)} f(x,y)\, dy\,dx\text{.} \end{equation*} 
If \(R\) is bounded by \(c \leq y \leq d, \ h_1(y) \leq x \leq h_2(y)\text{,}\) with \(h_1\) and \(h_2\) continuous on \([c,d]\text{,}\) then
\begin{equation*} \iint_R f(x,y)\,dA = \int_c^d \int_{h_1(y)}^{h_2(y)} f(x,y)\, dx \,dy\text{.} \end{equation*}
We explore the order of integration with one more example of a double integral.
Example 4.18. Volume of Region.
Find the volume under the surface \(\ds z=\sqrt{1x^2}\) and above the triangle formed by \(y=x\text{,}\) \(x=1\text{,}\) and the \(x\)axis.
Let's consider the two possible ways to set this up:
Which appears easier? In the first, the first (inner) integral is easy, because we need an antiderivative with respect to \(y\text{,}\) and the entire integrand \(\ds\sqrt{1x^2}\) is constant with respect to \(y\text{.}\) Of course, the second integral may be more difficult. In the second, the first integral is mildly unpleasant—a trig substitution. So let's try the first one, since the first step is easy, and see where that leaves us.
This is quite easy, since the substitution \(u=1x^2\) works:
Then
This is a good example of how the order of integration can affect the complexity of the problem. In this case it is possible to do the other order, but it is a bit messier. In some cases one order may lead to an integral that is either difficult or impossible to evaluate; it's usually worth considering both possibilities before going very far.
Exercises for Section 4.2.
Exercise 4.2.1.
Compute the following double integrals.

\(\ds \int_{0}^{2}\int_{0}^{4} 1+x \,dy\,dx\)
AnswerSolution\(16\)\begin{equation*} \begin{split} \int_0^2 \int_0^4 (1+x)\,dy\,dx \amp = \int_0^2 (1+x) y\big\vert_0^4\,dx \\ \amp= \int_0^2 4(1+x)\,dx = 4\left[x+\frac{x^2}{2}\right]_0^2 = 16.\end{split} \end{equation*} 
\(\ds \int_{1}^{1}\int_{0}^{2} x+y\,dy\,dx\)
AnswerSolution\(4\)\begin{equation*} \begin{split} \int_{1}^1\int_0^2 (x+y)\,dy\,dx \amp= \int_{1}^1 \left[xy + \frac{y^2}{2}\right]_0^2 \\ \amp= \int_{1}^2 \left(2x+2\right)\,dx = \left[x^2+2x\right]_{1}^1 = 4. \end{split} \end{equation*} 
\(\ds \int_{1}^{2}\int_{0}^{y} xy \,dx\,dy\)
AnswerSolution\(15/8\)\(\begin{aligned}\int_1^2\int_0^y xy \,dx\,dy \amp = \int_1^2 y \int_0^y x \,dx\,dy \\ \amp = \int_1^2 y\left[\frac{1}{2}x^2\right]_0^y\,dy \\ \amp = \int_1^2 \frac{1}{2}y^3\,dy = \frac{1}{2}\left[\frac{1}{4}y^4\right]_1^2 = \frac{15}{8} \end{aligned}\) 
\(\ds \int_{0}^{1}\int_{y^2/2}^{\sqrt y} \,dx\,dy\)
AnswerSolution\(1/2\)\begin{equation*} \begin{split} \int_0^1\int_{y^2/2}^{\sqrt{y}}\,dx\,dy \amp= \int_0^1 x\bigg\vert_{y^2/2}^{\sqrt{y}} \,dy \\ \amp= \int_0^1 \sqrt{y}\frac{y^2}{2}\,dy\\ \amp \left[\frac{2y^{3/2}}{3}\frac{y^3}{6}\right]_0^1 = \frac{1}{2} \end{split} \end{equation*} 
\(\ds \int_{1}^{2}\int_{1}^{x} {x^2\over y^2}\,dy\,dx\)
AnswerSolution\(5/6\)\begin{equation*} \begin{split} \int_1^2\int_1^x \frac{x^2}{y^2}\,dy\,dx \amp= \int_1^2 x^2 \left[\frac{1}{y}\right]_1^x\,dx \\ \amp= \int_1^2 x^2 \left(\frac{1}{x}+1\right)\,dx = \int_1^2 \left(x + x^2\right)\,dx\\ \amp= \left[\frac{x^2}{2} + \frac{x^3}{3}\right]_1^2 = \frac{5}{6} \end{split} \end{equation*} 
\(\ds \int_{0}^{1}\int_{0}^{x^2} y e^x\,dy\,dx\)
AnswerSolution\(\frac{9e}{2}  12\)\begin{equation*} \begin{split}\int_0^1 \int_0^{x^2} e^x y\,dy\,dx = \int_0^1 \frac{e^x y^2}{2} \big\vert_0^{x^2}\,dx = \int_0^1 \frac{e^x x^4}{2}\,dx \end{equation*}We now use integration by parts to solve the corresponding indefinite integral, with \(u=x^4\text{,}\) \(dv = e^xdx\text{,}\) and so \(v=e^x\text{,}\) and \(du=4x^3\,dx\text{:}\)\begin{equation*} \int \frac{e^x x^4}{2}\,dx = \frac{1}{2} \left[x^4e^x  4 \int e^x x^3\,dx \right] \end{equation*}We now continue by applying Integration by parts twice more with \(dv = e^x\,dx\text{:}\)\begin{equation*} \begin{split} \int \frac{e^x x^4}{2}\,dx \amp= \frac{1}{2} \left[x^4e^x  4 \int e^x x^3\,dx \right] \\ \amp= \frac{1}{2} \left[x^4 e^x  4x^3e^x + 12 \int e^x x^2\,dx \right]\\ \amp= \frac{1}{2} \left[x^4 e^x  4x^3e^x + 12x^2e^x  24 \int e^x x\,dx \right]\\ \amp= \frac{1}{2} \left[x^4 e^x  4x^3e^x + 12x^2e^x  24xe^x + 24e^x \,dx \right] \end{split} \end{equation*}Hence, we have\begin{equation*} \begin{split} \int_0^1 \int_0^{x^2} \frac{y}{e^x}\,dy\,dx \amp= \int_0^1 \frac{e^x x^4}{2}\,dx \\ \amp=\frac{1}{2} \left[x^4 e^x  4x^3e^x + 12x^2e^x  24xe^x + 24e^x \,dx \right]_0^1\\ \amp= \frac{9e}{2}  12.\end{split} \end{equation*} 
\(\ds \int_{0}^{\sqrt{\pi/2}}\int_{0}^{x^2} x\cos y\,dy\,dx\)
AnswerSolution\(1/2\)\(\begin{aligned}\int_0^1\int_{y^2/2}^{\sqrt{y}}\,dx\,dy \amp = \int_0^1 x\bigg\vert_{y^2/2}^{\sqrt{y}} \,dy \\ \amp = \int_0^1 \sqrt{y}\frac{y^2}{2}\,dy\\ \amp = \left[\frac{2y^{3/2}}{3}\frac{y^3}{6}\right]_0^1 = \frac{1}{2} \end{aligned}\) 
\(\ds \int_{0}^{\pi/2}\int_{0}^{\cos\theta}r^2(\cos\thetar) \,dr\,d\theta\)
AnswerSolution\(\pi/64\)\(\begin{aligned}\int_0^{\pi/2}\int_0^{\cos\theta} r^2 \left(\cos\theta  r\right)\,dr\,d\theta \amp =\int_0^{\pi/2} \left[\frac{r^3}{3}\cos\theta  \frac{r^4}{4}\right]_0^{\cos\theta} \,d\theta \\ \amp =\int_0^{\pi/2} \left[\frac{\cos^4\theta}{3}\frac{\cos^4\theta}{4}\right]\,d\theta = \int_0^{\pi/2} \frac{\cos^4\theta}{12} \,d\theta \end{aligned}\) Now use the identity \(\displaystyle{\cos^2\theta = \frac{1}{2}(1+\cos (2\theta))}\text{:}\)
\begin{equation*} \begin{split} \cos^4\theta \amp = (\cos^2\theta)^2 = \frac{1}{2}(1+\cos 2\theta)^2 \\ \amp = \frac{1}{4} \left(1 + \cos^22\theta + 2 \cos 2\theta\right) \\ \amp = \frac{1}{4}\left(1 + \frac{1}{2}(1+\cos 4\theta) + 2 \cos 2\theta\right) \\ \amp = \frac{3}{8} + \frac{1}{8} \cos 4\theta + \frac{1}{2}\cos 2\theta. \end{split} \end{equation*}Therefore,
\begin{equation*} \begin{split} \int_0^{\pi/2} \frac{\cos^4\theta}{12} \,d\theta \amp = \frac{1}{12}\int_0^{\pi/2} \left(\frac{3}{8} + \frac{1}{8} \cos 4\theta + \frac{1}{2}\cos 2\theta\right)\\ \amp = \frac{1}{12}\left[\frac{3\theta}{8} + \frac{\sin 4\theta}{32} + \frac{\sin 2\theta}{4}\right]_0^{\pi/2}\\ \amp = \frac{\pi}{64}. \end{split} \end{equation*}We find that
\begin{equation*} \int_0^{\pi/2}\int_0^{\cos\theta} r^2 \left(\cos\theta  r\right)\,dr\,d\theta = \frac{\pi}{64} \text{.} \end{equation*} 
\(\ds \int_0^1\int_{\sqrt{y}}^1 \sqrt{x^3+1}\,dx\,dy\)
AnswerSolution\(\frac{2}{9}\left(2\sqrt{2}1\right)\)First, we notice that we cannot evaluate the integral\begin{equation*} \int \sqrt{x^3+1} \,dx, \end{equation*}and so we try switching the order of integration. The region of integration is as follows: Therefore,\begin{equation*} \int_0^1 \int_{\sqrt{y}}^1 \sqrt{x^3+1} \,dx \,dy = \int_0^1 \int_0^{x^2} \sqrt{x^3+1} \,dy \,dx = \int_0^1 x^2 \sqrt{x^3+1} \,dx. \end{equation*}We now carry out the integration by making the following substitution: let \(u=x^3+1\text{,}\) then \(du = 3x^2\,dx\text{.}\) Then as \(x\) goes from \(0 \to 1\text{,}\) \(u\) goes from \(1 \to 2\text{:}\)\begin{equation*} \int_0^1 x^2 \sqrt{x^3+1} \,dx = \int_1^2 \frac{1}{3} \sqrt{u}\,du = \frac{2}{9} u^{3/2}\big\vert_1^2 = \frac{2}{9} \left(2\sqrt{2}1\right). \end{equation*}Therefore, we have found that\begin{equation*} \int_0^1 \int_{\sqrt{y}}^1 \sqrt{x^3+1} \,dx \,dy = \frac{2}{9}\left(2\sqrt{2}1\right). \end{equation*} 
\(\ds \int_0^1 \int_{x^2}^1 x\sqrt{1+y^2}\,dy\,dx\)
AnswerSolution\((2\sqrt21)/6\)First, notice that we cannot integrate\begin{equation*} \int x\sqrt{1+y^2}\,dy. \end{equation*}Therefore, we switch the order of integation: Therefore,\begin{equation*} \int_0^1 \int_{x^2}^1 x\sqrt{1+y^2} \,dy \,dx = \int_0^1 \int_0^{\sqrt{y}} x\sqrt{1+y^2} \,dx \,dy = \int_0^1 \frac{y}{2} \sqrt{1+y^2}\,dy. \end{equation*}Now use the substitution: \(u=1+y^2\) with \(du = 2y\,dy\text{.}\) Then as \(y\) goes from \(0 \to 1\text{,}\) we have \(u\) going from \(1 \to 2\text{.}\) Hence,\begin{equation*} \int_0^1 \frac{y}{2} \sqrt{1+y^2}\,dy = \int_1^2 \frac{\sqrt{u}}{4}\,du = \frac{u^{3/2}}{6}\big\vert_1^2 = \frac{1}{6} \left(2\sqrt{2}1\right). \end{equation*}Therefore, we have shown that\begin{equation*} \int_0^1 \int_{x^2}^1 x\sqrt{1+y^2} \,dy \,dx = \frac{1}{6} \left(2\sqrt{2}1\right). \end{equation*} 
\(\ds \int_0^1 \int_0^y {2\over\sqrt{1x^2}}\,dx\,dy\)
AnswerSolution\(\pi2\)We will find that this problem is easier if we first switch the order of integration: Therefore,\begin{equation*} \int_0^1 \int_0^y \frac{2}{\sqrt{1x^2}}\,dx \,dy = \int_0^1 \int_x^1 \frac{2}{\sqrt{1x^2}}\,dy \,dx = \int_0^1 \frac{2(1x)}{\sqrt{1x^2}}\,dx. \end{equation*}Hence, we split this integral into two:\begin{equation*} \int_0^1 \frac{2}{\sqrt{1x^2}}\,dx = 2\sin^{1}(x)\big\vert_0^1 = \pi. \end{equation*}For the second integral, we use the following substitution: \(u=1x^2\) with \(du=2x\,dx\text{.}\) Therefore,\begin{equation*} \int_0^1 \frac{2x}{\sqrt{1x^2}}\,dx = \int_0^1 \frac{1}{\sqrt{u}}\,du = 2. \end{equation*}Altogether, we find\begin{equation*} \int_0^1 \int_0^y \frac{2}{\sqrt{1x^2}}\,dx \,dy = \pi 2. \end{equation*} 
\(\ds \int_0^1 \int_{3y}^3 e^{x^2}\,dx\,dy\)
AnswerSolution\((e^91)/6\)We first switch the order of integration: Therefore,\begin{equation*} \int_0^1 \int_{3y}^3 e^{x^2}\,dx\,dy = \int_0^3\int_0^{x/3} e^{x^2}\,dy\,dy = \int_0^3 \frac{x}{3} e^{x^2}\,dx. \end{equation*}This is an integral which we can solve using substitution: let \(u=x^2\) with \(du=2x\,dx\text{.}\) Then,\begin{equation*} \int_0^3 \frac{x}{3} e^{x^2}\,dx = \int_0^9 \frac{1}{6} e^u\,du = \frac{1}{6} \left(e^91\right). \end{equation*}We have found that\begin{equation*} \int_0^1 \int_{3y}^3 e^{x^2}\,dx\,dy = \frac{1}{6} \left(e^91\right). \end{equation*} 
\(\ds \int_{1}^1\int_0^{1x^2} x^2\sqrt{y}\,dy\,dx\)
Answer\(\ds {4\over15}{\pi\over4}\) 
\(\ds \int_{0}^{\sqrt2/2}\int_{\sqrt{12x^2}}^{\sqrt{12x^2}} x\,dy\,dx\)
Answer\(1/3\)
Exercise 4.2.2.
Determine the volume of the region that is bounded as follows:

below \(z=1y\) and above the region \(1\le x\le 1\text{,}\) \(0\le y\le 1x^2\)
AnswerSolution\(4/5\)The area which bounds the surface in the \(x\)\(y\)plane is shown below. Therefore, the volume is computed as follows:\begin{equation*} \begin{split} V= \int_{1}^1 \int_{0}^{1x^2} (1y)\,dy\,dx \amp= \int_{1}^1 y\frac{y^2}{2} \big\vert_0^{1x^2} \,dx\\ \amp= \frac{1}{2}\int_{1}^1 (1x^4) \,dx = \frac{1}{2}\left[x\frac{x^5}{5}\right]_{1}^1 = \frac{4}{5}.\end{split} \end{equation*} 
enclosed by \(y=x^2\text{,}\) \(y=4\text{,}\) \(z=x^2\text{,}\) \(z=0\)
AnswerSolution\(128/15\)In the \(x\)\(y\)plane, the surface is bounded by: And so we set up the following integral to compute the volume:\begin{equation*} \begin{split} V=\int_{2}^2 \int_{x^2}^4 \int_0^{x^2} \,dz\,dy\,dx \amp= \int_{2}^2 \int_{x^2}^4 x^2\,dy\,dy \\ \amp= \int_{2}^2 x^4+4x^2\,dx = \frac{128}{15}.\end{split} \end{equation*} 
in the first octant bounded by \(y^2=4x\) and \(y=2z\)
AnswerSolution\(2\)The region is bounded by the following crosssections: Therefore, we compute the following as follows:\begin{equation*} V = \int_0^4 \int_0^{\sqrt{4x}} \frac{y}{2} \,dy\,dx=2. \end{equation*} 
in the first octant bounded by \(y^2=4x\text{,}\) \(2x+y=4\text{,}\) \(z=y\text{,}\) and \(y=0\)
AnswerSolution\(\frac{11}{3}\)In the first octant of the \(x\)\(y\)plane, the surface is bounded as follows: Therefore, we set up the integral:\begin{equation*} V = \int_0^1 \int_{2\sqrt{x}}^{42x} \int_0^y dz\,dy\,dx = \frac{11}{3}. \end{equation*} 
in the first octant bounded by \(x+y+z=9\text{,}\) \(2x+3y=18\text{,}\) and \(x+3y=9\text{.}\)
AnswerSolution\(\frac{81}{2}\)In the \(x\)\(y\)plane, the region is bounded by: Therefore, we compute the volume of the region:\begin{equation*} V = \int_0^9 \int_{(9x)/3}^{2(9x)/3} (9xy) \,dy\,dx = \frac{81}{2}. \end{equation*} 
in the first octant bounded by \(x^2+y^2=a^2\) and \(z=x+y\)
AnswerSolution\(\frac{2a^3}{3}\)In the \(x\)\(y\)plane, the region is bounded as follows: Therefore, we set up the following integral to compute the volume:\begin{equation*} V = \int_{a}^a \int_{\sqrt{a^2x^2}}^{\sqrt{a^2x^2}} (x+y)\,dy\,dx = \frac{2}{3} a^3. \end{equation*} 
bounded by \(4x^2+y^2=4z\) and \(z=2\)
AnswerSolution\(4\pi\)This volume is an elliptic paraboloid. For fixed \(z\text{,}\) the crosssections of the volume in the \(x\)\(y\) plane are ellipses. Therefore, we consider this to be the volume of the surface \(z= \frac{1}{4}\left(4x^2+y^2\right)\) below the ellipse \(\frac{x^2}{2} + \frac{y^2}{8}=1\text{:}\)\begin{equation*} V = \int_{\sqrt{2}}^{\sqrt{2}} \int_{\sqrt{84x^2}}^{\sqrt{84x^2}} \,dy\, dx = 4\pi. \end{equation*} 
bounded by \(z=x^2+y^2\) and \(z=4\)
AnswerSolution\(8\pi\)The volume bounded by the paraboloid \(z=x^2+y^2\) and the plane \(z=4\) is the volume above the surface \(z=x^2+y^2\) below the circle \(4=x^2+y^2\text{.}\) This volume can be computed by:\begin{equation*} V = \int_{2}^2 \int_{\sqrt{4y^2}}^{\sqrt{4y^2}} (x^2+y^2) \,dx \,dy = 8\pi. \end{equation*}  under the surface \(z=xy\) and above the triangle with vertices \((1,1,0)\text{,}\) \((4,1,0)\text{,}\) \((1,2,0)\) AnswerSolution\(31/8\)
The area which bounds the region in the \(x\)\(y\)plane is shown below.
Therefore, if we let \(x\) vary from \(1\) to \(4\text{,}\) we see that \(y\) is bounded between \(1\) and the line \(y=\frac{7}{3}\frac{x}{3}\text{.}\) And so the volume under the surface \(z=xy\) and above the area shown in green is calculated by\begin{equation*} \begin{split} \int_1^4\int_1^{7/3x/3} xy\,dy\,dx \amp = \int_1^4 x\left[\frac{1}{2}y^2\right]_1^{7/3x/3} \,dx\\ \amp =\int_1^4 \frac{x}{18}\left(x^214x+40\right)\,dx \\ \amp = \left[\frac{1}{18}\left(\frac{x^4}{4}\frac{14x^3}{3}+20x^2\right)\right]_1^4 =\frac{31}{8}. \end{split} \end{equation*}
Exercise 4.2.3.
Evaluate \(\ds\iint x^2\,dA\) over the region in the first quadrant bounded by the hyperbola \(xy=16\) and the lines \(y=x\text{,}\) \(y=0\text{,}\) and \(x=8\text{.}\)
AnswerExercise 4.2.4.
A swimming pool is circular with a 40 meter diameter. The depth is constant along eastwest lines and increases linearly from 2 meters at the south end to 7 meters at the north end. Find the volume of the pool.
AnswerExercise 4.2.5.
Find the average value of \(f(x,y)=e^y\sqrt{x+e^y}\) on the rectangle with vertices \((0,0)\text{,}\) \((4,0\)), \((4,1)\) and \((0,1)\text{.}\)
AnswerExercise 4.2.6.
Reverse the order of integration on each of the following integrals

\(\ds\int_0^9 \int_0^{\sqrt{9y}} f(x,y)\,dx \,dy\)
AnswerSolution\(\ds\int_0^3\int_0^{9x^2} f(x,y)\,dx\,dy\)The integral
\begin{equation*} \int_0^9 \int_0^{\sqrt{9y}} f(x,y)\,dx\,dy \end{equation*}gives the volume under the surface \(z=f(x,y)\) and above the region in the \(x\)\(y\)plane shown below:
So if instead we let \(x\) vary from 0 to 3, we see that \(y\) is bounded between \(y=0\) and the curve \(x=\sqrt{9y} \implies y = 9x^2\) (for \(x \in [0,3]\)). Therefore,\begin{equation*} \int_0^9 \int_0^{\sqrt{9y}} f(x,y)\,dx\,dy = \int_0^3\int_0^{9x^2} f(x,y)\,dx\,dy\text{.} \end{equation*} 
\(\ds\int_1^2 \int_0^{\ln x} f(x,y) \; dy \; dx\)
AnswerSolution\(\int_0^{\ln(2)} \int_{e^y}^{2} f(x,y)\,dx\,dy\)The area of integration is shown below: Therefore,\begin{equation*} \int_1^2 \int_0^{\ln x} f(x,y)\,dy\,dx = \int_0^{\ln(2)} \int_{e^y}^{2} f(x,y)\,dx\,dy. \end{equation*} 
\(\ds\int_0^1 \int_{\arcsin y}^{\pi/2} f(x,y) \; dx \; dy\)
AnswerSolution\(\int_0^{\pi/2} \int_0^{\sin x} f(x,y) \,dy\,dx\)The area of integration is shown below: Therefore, we see that\begin{equation*} \int_0^1 \int_{\arcsin(y)}^{\pi/2} f(x,y)\,dx\,dy = \int_0^{\pi/2} \int_0^{\sin x} f(x,y) \,dy\,dx. \end{equation*} 
\(\ds\int_0^1 \int_{4x}^{4} f(x,y) \; dy \; dx\)
AnswerSolution\(\int_0^4 \int_0^{y/4} f(x,y)\,dx \,dy\)The area of integration is shown below: Therefore, we see that\begin{equation*} \int_0^1 \int_{4x}^4 f(x,y)\,dy\,dx = \int_0^4 \int_0^{y/4} f(x,y)\,dx \,dy. \end{equation*} 
\(\ds\int_0^3 \int_{0}^{\sqrt{9y^2}} f(x,y) \; dx \; dy\)
AnswerSolution\(\int_0^3 \int_0^{\sqrt{9x^2}} f(x,y)\,dy\,dx\)Note that the circle centred at the origin with radius 3 has equation \(x^2+y^2=9\text{,}\) and in the first quadrant can be expressed as \(x=\sqrt{9y^2}\) or \(y=\sqrt{9x^2}\text{.}\) Therefore, the area of integration is: Therefore,\begin{equation*} \int_0^3\int_0^{\sqrt{9y^2}} f(x,y)\,dx\,dy = \int_0^3 \int_0^{\sqrt{9x^2}} f(x,y)\,dy\,dx. \end{equation*}