Now let’s go back to the definition of double integral for a minute:
Referring to Figure 1, observe that we divided the region [latex]S[/latex] in the [latex]uv[/latex]-plane into small subrectangles [latex]S_[/latex] and we let the subrectangles [latex]R_[/latex] in the [latex]xy[/latex]-plane be the images of [latex]S_[/latex] under the transformation [latex]T(u, v)=(x, y)[/latex].
On the left-hand side of this figure, there is a rectangle S with an inscribed red oval and a subrectangle with lower right corner point (u sub ij, v sub ij), height Delta v, and length Delta u given in the Cartesian u v-plane. Then there is an arrow from this graph to the right-hand side of the figure marked with T. On the right-hand side of this figure there is a region R with inscribed (deformed) red oval and a subrectangle R sub ij with corner point (x sub ij, y sub ij) given in the Cartesian x y-plane. The subrectangle is blown up and shown with vectors pointing along the edge from the corner point." width="765" height="321" />
Figure 1. The subrectangles [latex]S_[/latex] in the [latex]uv[/latex]-plane transform into subrectangles [latex]R_[/latex] in the [latex]xy[/latex]-plane.
Then the double integral becomes
Notice this is exactly the double Riemann sum for the integral
Let [latex]T(u, v)=(x, y)[/latex] where [latex]x=g(u, v)[/latex] and [latex]y=h(u, v)[/latex] be a one-to-one [latex]C1[/latex] transformation, with a nonzero Jacobian on the interior of the region [latex]S[/latex] in the [latex]uv[/latex]-plane; it maps [latex]S[/latex] into the region [latex]R[/latex] in the [latex]xy[/latex]-plane. If [latex]f[/latex] is continuous on [latex]R[/latex], then
With this theorem for double integrals, we can change the variables from [latex](x, y)[/latex] to [latex](u, v)[/latex] in a double integral simply by replacing
when we use the substitutions [latex]x=g(u, v)[/latex] and [latex]y=h(u, v)[/latex] and then change the limits of integration accordingly. This change of variables often makes any computations much simpler.
Consider the integral
Use the change of variables [latex]x=r\cos\theta[/latex] and [latex]y=r\sin\theta[/latex], and find the resulting integral.
Show SolutionFirst we need to find the region of integration. This region is bounded below by [latex]y=0[/latex] and above by [latex]y=\sqrt[/latex] (see the following figure).
the first quadrant of the xy plane with radius 1 and center (1, 0). The equation for this curve is given as y = the square root of (2x minus x squared)" width="290" height="234" />
Figure 2. Changing a region from rectangular to polar coordinates.
Squaring and collecting terms, we find that the region is the upper half of the circle [latex]x^+y^-2x=0[/latex] that is, [latex]y^+(x-1)^[/latex]. In polar coordinates, the circle is [latex]r=2\cos\theta[/latex] so the region of integration in polar coordinates is bounded by [latex]0\leq\leq\cos\theta[/latex] and [latex]0\leq\theta\leq\frac<\pi>2[/latex].
The Jacobian is [latex]J(r,\theta)=r[/latex], as shown in Example “Finding the Jacobian”. Since [latex]r\geq0[/latex], we have [latex]|J(r,\theta)|=r[/latex].
The integrand [latex]\sqrt[/latex] changes to [latex]r[/latex] in polar coordinates, so the double iterated integral is
[latex]\displaystyle\int_0^2\displaystyle\int_0^>\sqrtdy \ dx=\displaystyle\int_0^<\pi/2>\displaystyle\int_0^r|J(r,\theta)|dr \ d\theta=\displaystyle\int_0^<\pi/2>\displaystyle\int_0^r^2 \ dr \ d\theta[/latex].
Considering the integral [latex]\displaystyle\int_0^1\displaystyle\int_0^>(x^2+y^2)dy \ dx[/latex], use the change of variables [latex]x=r\cos\theta[/latex] and [latex]y=\sin\theta[/latex], and find the resulting integral.
Show Solution[latex]\displaystyle\int_0^<\pi/2>\displaystyle\int_0^1r^3 \ dr \ d\theta[/latex]
Notice in the next example that the region over which we are to integrate may suggest a suitable transformation for the integration. This is a common and important situation.
Consider the integral [latex]\underset<\displaystyle\iint>(x-y)dy \ dx[/latex], where [latex]R[/latex] is the parallelogram joining the points [latex](1, 2)[/latex], [latex](3, 4)[/latex], [latex](4, 3)[/latex], and [latex](6, 5)[/latex] (Figure 3). Make appropriate changes of variables, and write the resulting integral.
the first quadrant of the xy plane with radius 1 and center (1, 0). The equation for this curve is given as y = the square root of (2x minus x squared)" width="328" height="272" />
Figure 3. The region of integration for the given integral.
Show SolutionFirst, we need to understand the region over which we are to integrate. The sides of the parallelogram are [latex]x-y+1=0[/latex], [latex]x-y-1=0[/latex], [latex]x-3y+5=0[/latex], and [latex]x-3y+9=0[/latex] (Figure 4). Another way to look at them is [latex]x-y=-1[/latex], [latex]x-y=1[/latex], [latex]x-3y=-5[/latex], and [latex]x-3y=-9[/latex].
Clearly the parallelogram is bounded by the lines [latex]y=x+1[/latex], [latex]y=x-1[/latex], [latex]y=\frac13(x+5)[/latex], and [latex]y=\frac13(x+9)[/latex].
Notice that if we were to make [latex]u=x-y[/latex] and [latex]v=x-3y[/latex], then the limits on the integral would be [latex]-1\leq\leq1[/latex] and [latex]-9\leq\leq-5[/latex].
To solve for [latex]x[/latex] and [latex]y[/latex], we multiply the first equation by [latex]3[/latex] and subtract the second equation, [latex]3u-v=(3x-3y)-(x-3y)=2x[/latex]. Then we have [latex]x=\frac2[/latex]. Moreover, if we simply subtract the second equation from the first, we get [latex]u-v=(x-y)-(x-3y)=2y[/latex] and [latex]y=\frac2[/latex].
Figure 4. A parallelogram in the [latex]xy[/latex]-plane that we want to transform by a change in variables.
Thus, we can choose the transformation
and compute the Jacobian [latex]J(u, v)[/latex]. We have
Therefore, [latex]|J(u,v)|=\frac12[/latex]. Also, the original integrand becomes
Therefore, by the use of the transformation [latex]T[/latex], the integral changes to
[latex]\large<\underset
which is much simpler to compute. In fact, it is easily found to be zero. And this is just one example of why we transform integrals like this.
Make appropriate changes of variables in the integral [latex]\underset<\displaystyle\iint>\frac4dy \ dx[/latex], where [latex]R[/latex] is the trapezoid bounded by the lines [latex]x-y=2[/latex], [latex]x-y=4[/latex], [latex]x=0[/latex], and [latex]y=0[/latex]. Write the resulting integral.
Show Solution[latex]x=\frac12(v+u)[/latex] and [latex]y=\frac12(u-v)[/latex] and [latex]\displaystyle\int_2^4\displaystyle\int_^u\frac4\left(\frac12\right)dv \ du[/latex].
We are ready to give a problem-solving strategy for change of variables.
In the next example, we find a substitution that makes the integrand much simpler to compute.
Using the change of variables [latex]u=x-y[/latex] and [latex]v=x+y[/latex], evaluate the integral
where [latex]R[/latex] is the region bounded by the lines [latex]x+y=1[/latex] and [latex]x+y=3[/latex] and the curves [latex]x^-y^=-1[/latex] and [latex]x^-y^=1[/latex] (see the first region in Figure 5).
Show SolutionAs before, first find the region [latex]R[/latex] and picture the transformation so it becomes easier to obtain the limits of integration after the transformations are made (Figure 5).
Figure 5. Transforming the region [latex]R[/latex] into the region [latex]S[/latex] to simplify the computation of an integral.
Given [latex]u=x-y[/latex] and [latex]v=x+y[/latex], we have [latex]x=\frac2[/latex] and [latex]y=\frac2[/latex] and hence the transformation to use is [latex]T(u,v)=\left(\frac2,\frac2\right)[/latex]. The lines [latex]x+y=1[/latex] and [latex]x+y=3[/latex] become [latex]v=1[/latex] and [latex]v=3[/latex], respectively. The curves [latex]x^-y^=1[/latex] and [latex]x^-y^=-1[/latex] become [latex]uv=1[/latex] and [latex]uv=-1[/latex], respectively.
Thus we can describe the region [latex]S[/latex] (see the second region Figure 5) as
The Jacobian for this transformation is
Therefore, by using the transformation [latex]T[/latex], the integral changes to
Doing the evaluation, we have
Using the substitutions [latex]x=v[/latex] and [latex]y=\sqrt[/latex], evaluate the integral [latex]\underset<\displaystyle\iint>y\sin(y^2-x)dA[/latex] where [latex]R[/latex] is the region bounded by the lines [latex]y=\sqrt,\ x=2[/latex] and [latex]y=0[/latex].
Show SolutionWatch the following video to see the worked solution to the above Try It
Changing variables in triple integrals works in exactly the same way. Cylindrical and spherical coordinate substitutions are special cases of this method, which we demonstrate here.
Suppose that [latex]G[/latex] is a region in [latex]uvw[/latex]-space and is mapped to [latex]D[/latex] in [latex]xyz[/latex]-space (Figure 6) by a one-to-one [latex]C1[/latex] transformation [latex]T(u, v, w)=(x, y, z)[/latex] where [latex]x=g(u, v, w)[/latex], [latex]y=h(u, v, w)[/latex], and [latex]z=k(u, v, w)[/latex].
On the left-hand side of this figure, there is a region G in u v w space. Then there is an arrow from this graph to the right-hand side of the figure marked with x = g(u, v, w), y = h(u, v, w), and z = k(u, v, w). On the right-hand side of this figure there is a region D in xyz space." width="930" height="350" />
Figure 6. A region [latex]G[/latex] in [latex]uvw[/latex]-space mapped to a region [latex]D[/latex] in [latex]xyz[/latex]-space.
Then any function [latex]F(x, y, z)[/latex] defined on [latex]D[/latex] can be thought of as another function [latex]H(u, v, w)[/latex] that is defined on [latex]G[/latex]:
Now we need to define the Jacobian for three variables.
The Jacobian determinant [latex]J(u, v, w)[/latex] in three variables is defined as follows:
This is also the same as
The Jacobian can also be simply denoted as [latex]\frac<\partial(x,y,z)><\partial(u,v,w)>[/latex].
With the transformations and the Jacobian for three variables, we are ready to establish the theorem that describes change of variables for triple integrals.
Let [latex]T(u, v, w)=(x, y, z)[/latex] where [latex]x=g(u, v, w)[/latex], [latex]y=h(u, v, w)[/latex], and [latex]z=k(u, v, w)[/latex], be a one-to-one [latex]C1[/latex] transformation, with a nonzero Jacobian, that maps the region [latex]G[/latex] in the [latex]uvw[/latex]-plane into the region [latex]D[/latex] in the [latex]xyz[/latex]-plane. As in the two-dimensional case, if [latex]F[/latex] is continuous on [latex]D[/latex], then
Let us now see how changes in triple integrals for cylindrical and spherical coordinates are affected by this theorem. We expect to obtain the same formulas as in Triple Integrals in Cylindrical and Spherical Coordinates.
Derive the formula in triple integrals for
On the left-hand side of this figure, there is a cube G with sides parallel to the coordinate axes in cylindrical coordinate space. Then there is an arrow from this graph to the right-hand side of the figure marked with x = r cos theta, y = r sin theta, and z = z. On the right-hand side of this figure there is a region D in x y z space that is a thick annulus. The top is labeled z = constant, the flat vertical side is labeled theta = constant, and the outermost side is labeled r = constant." width="898" height="391" />
Figure 7. The transformation from rectangular coordinates to cylindrical coordinates can be treated as a change of variables from region [latex]G[/latex] in [latex]r\theta[/latex]-space to region [latex]D[/latex] in [latex]xyz[/latex]-space.
Expanding the determinant with respect to the third row:
[latex]\hspace\begin &=\cos\varphi\begin-\rho\sin\varphi\sin\theta&\rho\cos\varphi\cos\theta \\ \rho\sin\varphi\sin\theta&\rho\cos\varphi\sin\theta\end-\rho\sin\varphi\begin\sin\varphi\cos\theta&-\rho\sin\varphi\sin\theta \\ \sin\varphi\sin\theta&\rho\sin\varphi\cos\theta\end \\ &=\cos\varphi(-\rho^2\sin\varphi\cos\varphi\sin^2\theta-\rho^2\sin\varphi\cos\varphi\cos^2\theta)-\rho\sin\varphi(\rho\sin^2\varphi\cos^2\theta+\rho\sin^2\varphi\sin^2\theta) \\ &=-\rho^2\sin\varphi\cos^2\varphi(\sin^2\theta+\cos^2\theta)-\rho^2\sin\varphi\sin^2\varphi(\sin^2\theta+\cos^2\theta) \\ &=-\rho^2\sin\varphi\cos^2\varphi-\rho^2\sin\varphi\sin^2\varphi \\ &=-\rho^2\sin\varphi(\cos^2\varphi+\sin^2\varphi)=-\rho^2\sin\varphi\end[/latex]
Since [latex]0\leq\varphi\leq\pi[/latex], we must have [latex]\sin\varphi\geq0[/latex]. Thus [latex]|J(\rho,\theta,\varphi)|=|-\rho^2\sin\varphi|=\rho^2\sin\varphi[/latex].
On the left-hand side of this figure, there is a cube G with sides parallel to the coordinate axes in rho phi theta space. Then there is an arrow from this graph to the right-hand side of the figure marked with x = rho sin phi cos theta, y = rho sin phi sin theta, and z = rho cos phi. On the right-hand side of this figure there is a region D in xyz space that is a thick annulus and has the point (x, y, z) shown as being equal to (rho, phi, theta). The top is labeled phi = constant, the flat vertical side is labeled theta = constant, and the outermost side is labeled rho = constant." width="895" height="391" />
Figure 8. The transformation from rectangular coordinates to cylindrical coordinates can be treated as a change of variables from region [latex]G[/latex] in [latex]\rho\theta\varphi[/latex]-space to region [latex]D[/latex] in [latex]xyz[/latex]-space.
Let’s try another example with a different substitution.
Evaluate the triple integral
in [latex]xyz[/latex]-space by using the transformation
Then integrate over an appropriate region in [latex]uvw[/latex]-space
Show SolutionAs before, some kind of sketch of the region [latex]G[/latex] in [latex]xyz[/latex]-space over which we have to perform the integration can help identify the region [latex]D[/latex] in [latex]uvw[/latex]-space (Figure 9). Clearly [latex]G[/latex] in [latex]xyz[/latex]-space is bounded by the planes [latex]x=y/2[/latex], [latex]x=(y/2)+1[/latex], [latex]y=0[/latex],[latex]y=4[/latex], [latex]z=0[/latex], and [latex]z=4[/latex]. We also know that we have to use [latex]u=(2x-y)/2[/latex], [latex]v=y/2[/latex] and [latex]w=z/3[/latex] for the transformations. We need to solve for [latex]x[/latex], [latex]y[/latex], and [latex]z[/latex]. Here we find that [latex]x=u+v[/latex], [latex]y=2v[/latex], and [latex]z=3w[/latex].
Using elementary algebra, we can find the corresponding surfaces for the region [latex]G[/latex] and the limits of integration in [latex]uvw[/latex]-space. It is convenient to list these equations in a table.
Equations in [latex]xyz[/latex] for the region [latex]D[/latex] | Corresponding equations in [latex]uvw[/latex] for the region [latex]G[/latex] | Limits for the integration in [latex]uvw[/latex] |
---|---|---|
[latex]x=y/2[/latex] | [latex]u+v=2v/2=v[/latex] | [latex]u=0[/latex] |
[latex]x=y/2[/latex] | [latex]u+v=(2v/2)+1=v+1[/latex] | [latex]u=1[/latex] |
[latex]y=0[/latex] | [latex]2v=0[/latex] | [latex]v=0[/latex] |
[latex]y=4[/latex] | [latex]2v=4[/latex] | [latex]v=2[/latex] |
[latex]z=0[/latex] | [latex]3w=0[/latex] | [latex]w=0[/latex] |
[latex]z=3[/latex] | [latex]3w=3[/latex] | [latex]w=1[/latex] |
On the left-hand side of this figure, there is a box G with sides 1, 2, and 1 along the u, v, and w axes, respectively. Then there is an arrow from this graph to the right-hand side of the figure marked with x = u + v, y = 2v, and z = 3w. On the right-hand side of this figure there is a region D in xyz space that is a rotated box with sides 1, 4, and 3 along the x, y, and z axes. The rear plane is marked x = y/2 or y = 2x. The front plane is marked x = y/2 + 1 or y = 2x minus 2." width="894" height="435" />
Figure 9. The region [latex]G[/latex] in [latex]uvw[/latex]-space is transformed to region [latex]D[/latex] in [latex]xyz[/latex]-space
Now we can calculate the Jacobian for the transformation:
The function to be integrated becomes
We are now ready to put everything together and complete the problem.
[latex]\hspace\large &\hspace\displaystyle\int_0^3\displaystyle\int_0^4\displaystyle\int_^<(y/2)+1>\left(x+\frac3\right)dx \ dy \ dz \\ &=\displaystyle\int_0^1\displaystyle\int_0^2\displaystyle\int_0^1(u+v+w)|J(u,v,w)|du \ dv \ dw=\displaystyle\int_0^1\displaystyle\int_0^2\displaystyle\int_0^1(u+v+w)|6|du \ dv \ dw \\ &=6\displaystyle\int_0^1\displaystyle\int_0^2\displaystyle\int_0^1(u+v+w)du \ dv \ dw=6\displaystyle\int_0^1\displaystyle\int_0^2\left[\frac2+vu+wu\right]_0^1dv \ dw \\ &=6\displaystyle\int_0^1\displaystyle\int_0^2\left(\frac12+v+w\right)dv \ dw=6\displaystyle\int_0^1\left[\frac12v+\frac2+wv\right]_0^2dw \\ &=6\displaystyle\int_0^1(3+2w)dw=6\left[3w+w^2\right]_0^1=24 \end>[/latex].
Let [latex]D[/latex] be the region in [latex]xyz[/latex]-space defined by [latex]1\leq\leq2, \ 0\leq\leq2, \text \ 0\leq\leq1[/latex].
Evaluate [latex]\underset<\displaystyle\iiint>(x^2y+3xyz)dx \ dy \ dz[/latex] by using the transformation [latex]u=x[/latex], [latex]v=xy[/latex], and [latex]w=3z[/latex].
Show Solution[latex]\displaystyle\int_0^3\displaystyle\int_0^2\displaystyle\int_1^2\left(\frac3+\frac\right)dy \ dv \ dw=2+\ln8 = 2+ 3\ln2[/latex].
Watch the following video to see the worked solution to the above Try It
Licenses and Attributions CC licensed content, Original