5.7 Change of Variables in Multiple Integrals

Recall from Substitution Rule the method of integration by substitution. When evaluating an integral such as 23x(x24)5dx, we substitute u=g(x)=x24. Then du=2xdx or xdx=12du and the limits change to u=g(2)=224=0 and u=g(3)=94=5. Thus the integral becomes 0512u5du and this integral is much simpler to evaluate. In other words, when solving integration problems, we make appropriate substitutions to obtain an integral that becomes much simpler than the original integral.

We also used this idea when we transformed double integrals in rectangular coordinates to polar coordinates and transformed triple integrals in rectangular coordinates to cylindrical or spherical coordinates to make the computations simpler. More generally,

abf(x)dx=cdf(g(u))g(u)du,

Where x=g(u),dx=g(u)du, and u=c and u=d satisfy c=g-1(a) and d=g-1(b).

A similar result occurs in double integrals when we substitute x=h(r,θ)=rcosθ, y=g(r,θ)=rsinθ, and dA=dxdy=rdrdθ. Then we get

Rf(x,y)dA=Sf(rcosθ,rsinθ)rdrdθ

where the domain R is replaced by the domain S in polar coordinates. Generally, the function that we use to change the variables to make the integration simpler is called a transformation or mapping.

Planar Transformations

A planar transformation T is a function that transforms a region G in one plane into a region R in another plane by a change of variables. Both G and R are subsets of 2. For example, Figure 5.71 shows a region G in the uv-plane transformed into a region R in the xy-plane by the change of variables x=g(u,v) and y=h(u,v), or sometimes we write x=x(u,v) and y=y(u,v). We shall typically assume that each of these functions has continuous first partial derivatives, which means gu,gv,hu, and hv exist and are also continuous. The need for this requirement will become clear soon.

On the left-hand side of this figure, there is a region G with point (u, v) 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 x = g(u, v) and y = h(u, v). On the right-hand side of this figure there is a region R with point (x, y) given in the Cartesian xy- plane.
Figure 5.71 The transformation of a region G in the uv-plane into a region R in the xy-plane.
Definition

A transformation T:GR, defined as T(u,v)=(x,y), is said to be a one-to-one transformation if no two points map to the same image point.

To show that T is a one-to-one transformation, we assume T(u1,v1)=T(u2,v2) and show that as a consequence we obtain (u1,v1)=(u2,v2). If the transformation T is one-to-one in the domain G, then the inverse T−1 exists with the domain R such that T−1T and TT−1 are identity functions.

Figure 5.71 shows the mapping T(u,v)=(x,y) where x and y are related to u and v by the equations x=g(u,v) and y=h(u,v). The region G is the domain of T and the region R is the range of T, also known as the image of G under the transformation T.

Determining How the Transformation Works

Suppose a transformation T is defined as T(r,θ)=(x,y) where x=rcosθ,y=rsinθ. Find the image of the polar rectangle G={(r,θ)|0<r1,0θπ/2} in the rθ-plane to a region R in the xy-plane. Show that T is a one-to-one transformation in G and find T−1(x,y).

Since r varies from 0 to 1 in the rθ-plane, we have a circular disc of radius 0 to 1 in the xy-plane. Because θ varies from 0 to π/2 in the rθ-plane, we end up getting a quarter circle of radius 1 in the first quadrant of the xy-plane (Figure 5.72). Hence R is a quarter circle bounded by x2+y2=1 in the first quadrant.

On the left-hand side of this figure, there is a rectangle G with a marked subrectangle given in the first quadrant of the Cartesian r theta-plane. Then there is an arrow from this graph to the right-hand side of the figure marked with x = r cos theta and y = r sin theta. On the right-hand side of this figure there is a quarter circle R with a marked subannulus (analogous to the rectangle in the other graph) given in the Cartesian x y-plane.
Figure 5.72 A rectangle in the rθ-plane is mapped into a quarter circle in the xy-plane.

In order to show that T is a one-to-one transformation, assume T(r1,θ1)=T(r2,θ2) and show as a consequence that (r1,θ1)=(r2,θ2). In this case, we have

T(r1,θ1)=T(r2,θ2),(x1,y1)=(x2,y2),(r1cosθ1,r1sinθ1)=(r2cosθ2,r2sinθ2),r1cosθ1=r2cosθ2r1sinθ1=r2sinθ2

Dividing, we obtain

r1cosθ1r1sinθ1=r2cosθ2r2sinθ2cosθ1sinθ1=cosθ2sinθ2cotθ1=cotθ2θ1=θ2

since the cotangent function is one-one function in the interval 0θπ/2. Also, since 0<r1, we have r1=r2,θ1=θ2. Therefore, (r1,θ1)=(r2,θ2) and T is a one-to-one transformation from G into R.

To find T−1(x,y) solve for r,θ in terms of x,y. We already know that r2=x2+y2 and tanθ=yx. Thus T−1(x,y)=(r,θ) is defined as r=x2+y2 and θ=tan−1(yx).

Finding the Image under T

Let the transformation T be defined by T(u,v)=(x,y) where x=u2v2 and y=uv. Find the image of the triangle in the uv-plane with vertices (0,0),(0,1), and (1,1).

The triangle and its image are shown in Figure 5.73. To understand how the sides of the triangle transform, call the side that joins (0,0) and (0,1) side A, the side that joins (0,0) and (1,1) side B, and the side that joins (1,1) and (0,1) side C.

On the left-hand side of this figure, there is a triangular region given in the Cartesian uv-plane with boundaries A, B, and C represented by the v axis, the line u = v, and the line v = 1, respectively. Then there is an arrow from this graph to the right-hand side of the figure marked with x = u squared minus v squared and y = u v. On the right-hand side of this figure there is a complex region given in the Cartesian x y-plane with boundaries A’, B’, and C’ given by the x axis, y axis, and a line curving from (negative 1, 0) through (0, 1), namely x = y squared minus 1, respectively.
Figure 5.73 A triangular region in the uv-plane is transformed into an image in the xy-plane.

For the side A:u=0,0v1 transforms to x=v2,y=0 so this is the side A that joins (−1,0) and (0,0).

For the side B:u=v,0u1 transforms to x=0,y=u2 so this is the side B that joins (0,0) and (0,1).

For the side C:0u1,v=1 transforms to x=u21,y=u (hence x=y21) so this is the side C that makes the upper half of the parabolic arc joining (−1,0) and (0,1).

All the points in the entire region of the triangle in the uv-plane are mapped inside the parabolic region in the xy-plane.

Let a transformation T be defined as T(u,v)=(x,y) where x=u+v,y=3v. Find the image of the rectangle G={(u,v):0u1,0v2} from the uv-plane after the transformation into a region R in the xy-plane. Show that T is a one-to-one transformation and find T−1(x,y).

T−1(x,y)=(u,v) where u=3xy3 and v=y3

Hint

Follow the steps of Example 5.65.

Jacobians

Recall that we mentioned near the beginning of this section that each of the component functions must have continuous first partial derivatives, which means that gu,gv,hu, and hv exist and are also continuous. A transformation that has this property is called a C1 transformation (here C denotes continuous). Let T(u,v)=(g(u,v),h(u,v)), where x=g(u,v) and y=h(u,v), be a one-to-one C1 transformation. We want to see how it transforms a small rectangular region S, Δu units by Δv units, in the uv-plane (see the following figure).

On the left-hand side of this figure, there is a region S with lower right corner point (u sub 0, v sub 0), 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 point (x sub 0, y sub 0) given in the Cartesian x y-plane with sides r(u, v sub 0) along the bottom and r(u sub 0, v) along the left.
Figure 5.74 A small rectangle S in the uv-plane is transformed into a region R in the xy-plane.

Since x=g(u,v) and y=h(u,v), we have the position vector r(u,v)=g(u,v)i+h(u,v)j of the image of the point (u,v). Suppose that (u0,v0) is the coordinate of the point at the lower left corner that mapped to (x0,y0)=T(u0,v0). The line v=v0 maps to the image curve with vector function r(u,v0), and the tangent vector at (x0,y0) to the image curve is

ru=gu(u0,v0)i+hu(u0,v0)j=xui+yuj.

Similarly, the line u=u0 maps to the image curve with vector function r(u0,v), and the tangent vector at (x0,y0) to the image curve is

rv=gv(u0,v0)i+hv(u0,v0)j=xvi+yvj.

Now, note that

ru=limΔu0r(u0+Δu,v0)r(u0,v0)Δusor(u0+Δu,v0)r(u0,v0)Δuru.

Similarly,

rv=limΔv0r(u0,v0+Δv)r(u0,v0)Δvsor(u0,v0+Δv)r(u0,v0)Δvrv.

This allows us to estimate the area ΔA of the image R by finding the area of the parallelogram formed by the sides Δvrv and Δuru. By using the cross product of these two vectors by adding the k component as 0, the area ΔA of the image R (refer to The Cross Product) is approximately ||Δuru×Δvrv||=||ru×rv||ΔuΔv. In determinant form, the cross product is

ru×rv=|ijkxuyu0xvyv0|=|xuyuxvyv|k=(xuyvxvyu)k.

Since ||k||=1, we have ΔA||ru×rv||ΔuΔv=(xuyvxvyu)ΔuΔv.

Definition

The Jacobian of the C1 transformation T(u,v)=(g(u,v),h(u,v)) is denoted by J(u,v) and is defined by the 2×2 determinant

J(u,v)=|(x,y)(u,v)|=|xuyuxvyv|=(xuyvxvyu).

Using the definition, we have

ΔAJ(u,v)ΔuΔv=|(x,y)(u,v)|ΔuΔv.

Note that the Jacobian is frequently denoted simply by

J(u,v)=(x,y)(u,v).

Note also that

|xuyuxvyv|=(xuyvxvyu)=|xuxvyuyv|.

Hence the notation J(u,v)=(x,y)(u,v) suggests that we can write the Jacobian determinant with partials of x in the first row and partials of y in the second row.

Finding the Jacobian

Find the Jacobian of the transformation given in Example 5.65.

The transformation in the example is T(r,θ)=(rcosθ,rsinθ) where x=rcosθ and y=rsinθ. Thus the Jacobian is

J(r,θ)=(x,y)(r,θ)=|xrxθyryθ|=|cosθrsinθsinθrcosθ|=rcos2θ+rsin2θ=r(cos2θ+sin2θ)=r.
Finding the Jacobian

Find the Jacobian of the transformation given in Example 5.66.

The transformation in the example is T(u,v)=(u2v2,uv) where x=u2v2 and y=uv. Thus the Jacobian is

J(u,v)=(x,y)(u,v)=|xuxvyuyv|=|2uv2vu|=2u2+2v2.

Find the Jacobian of the transformation: T(u,v)=(u+v,2v).

J(u,v)=(x,y)(u,v)=|xuxvyuyv|=|1102|=2

Hint

Follow the steps in the previous two examples.

Change of Variables for Double Integrals

We have already seen that, under the change of variables T(u,v)=(x,y) where x=g(u,v) and y=h(u,v), a small region ΔA in the xy-plane is related to the area formed by the product ΔuΔv in the uv-plane by the approximation

ΔAJ(u,v)Δu,Δv.

Now let’s go back to the definition of double integral for a minute:

Rf(x,y)dA=limm,ni=1mj=1nf(xij,yij)ΔA.

Referring to Figure 5.75, observe that we divided the region S in the uv-plane into small subrectangles Sij and we let the subrectangles Rij in the xy-plane be the images of Sij under the transformation T(u,v)=(x,y).

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.
Figure 5.75 The subrectangles Sij in the uv-plane transform into subrectangles Rij in the xy-plane.

Then the double integral becomes

Rf(x,y)dA=limm,ni=1mj=1nf(xij,yij)ΔA=limm,ni=1mj=1nf(g(uij,vij),h(uij,vij))|J(uij,vij)|ΔuΔv.

Notice this is exactly the double Riemann sum for the integral

Sf(g(u,v),h(u,v))|(x,y)(u,v)|dudv.
Change of Variables for Double Integrals

Let T(u,v)=(x,y) where x=g(u,v) and y=h(u,v) be a one-to-one C1 transformation, with a nonzero Jacobian on the interior of the region S in the uv-plane; it maps S into the region R in the xy-plane. If f is continuous on R, then

Rf(x,y)dA=Sf(g(u,v),h(u,v))|(x,y)(u,v)|dudv.

With this theorem for double integrals, we can change the variables from (x,y) to (u,v) in a double integral simply by replacing

dA=dxdy=|(x,y)(u,v)|dudv

when we use the substitutions x=g(u,v) and y=h(u,v) and then change the limits of integration accordingly. This change of variables often makes any computations much simpler.

Changing Variables from Rectangular to Polar Coordinates

Consider the integral

0202xx2x2+y2dydx.

Use the change of variables x=rcosθ and y=rsinθ, and find the resulting integral.

First we need to find the region of integration. This region is bounded below by y=0 and above by y=2xx2 (see the following figure).

A semicircle in 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)
Figure 5.76 Changing a region from rectangular to polar coordinates.

Squaring and collecting terms, we find that the region is the upper half of the circle x2+y22x=0, that is, y2+(x1)2=1. In polar coordinates, the circle is r=2cosθ so the region of integration in polar coordinates is bounded by 0r2 cosθ and 0θπ2.

The Jacobian is J(r,θ)=r, as shown in Example 5.67. Since r0, we have |J(r,θ)|=r.

The integrand x2+y2 changes to r in polar coordinates, so the double iterated integral is

0202xx2x2+y2dydx=0π/202cosθr|J(r,θ)|drdθ=0π/202cosθr2drdθ.

Considering the integral 0101x2(x2+y2)dydx, use the change of variables x=rcosθ and y=rsinθ, and find the resulting integral.

0π/201r3drdθ

Hint

Follow the steps in the previous example.

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.

Changing Variables

Consider the integral R(xy)dydx, where R is the parallelogram joining the points (1,2), (3,4),(4,3), and (6,5) (Figure 5.77). Make appropriate changes of variables, and write the resulting integral.

A parallelogram R with corners (1, 2), (3, 4), (6, 5), and (4, 3).
Figure 5.77 The region of integration for the given integral.

First, we need to understand the region over which we are to integrate. The sides of the parallelogram are xy+1=0,xy1=0, x3y+5=0,andx3y+9=0 (Figure 5.78). Another way to look at them is xy=−1,xy=1, x3y=−5, and x3y=9.

Clearly the parallelogram is bounded by the lines y=x+1,y=x1,y=13(x+5), and y=13(x+9).

Notice that if we were to make u=xy and v=x3y, then the limits on the integral would be −1u1 and −9v5.

To solve for x and y, we multiply the first equation by 3 and subtract the second equation, 3uv=(3x3y)(x3y)=2x. Then we have x=3uv2. Moreover, if we simply subtract the second equation from the first, we get uv=(xy)(x3y)=2y and y=uv2.

A parallelogram R with corners (1, 2), (3, 4), (6, 5), and (4, 3) formed by the lines y = x + 1, y = x minus 1, y = (x + 9)/3, and y = (x + 5)/3.
Figure 5.78 A parallelogram in the xy-plane that we want to transform by a change in variables.

Thus, we can choose the transformation

T(u,v)=(3uv2,uv2)

and compute the Jacobian J(u,v). We have

J(u,v)=(x,y)(u,v)=|xuxvyuyv|=|3/21/21/21/2|=34+14=12.

Therefore, |J(u,v)|=12. Also, the original integrand becomes

xy=12[3uvu+v]=12[3uu]=12[2u]=u.

Therefore, by the use of the transformation T, the integral changes to

R(xy)dydx=−9−5−11J(u,v)ududv=−9−5−11(12)ududv,

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 R4(xy)2dydx, where R is the trapezoid bounded by the lines xy=2,xy=4,x=0,andy=0. Write the resulting integral.

x=12(v+u) and y=12(vu) and 24−uu4u2(12)dvdu.

Hint

Follow the steps in the previous example.

We are ready to give a problem-solving strategy for change of variables.

Change of Variables
  1. Sketch the region given by the problem in the xy-plane and then write the equations of the curves that form the boundary.
  2. Depending on the region or the integrand, choose the transformations x=g(u,v) and y=h(u,v).
  3. Determine the new limits of integration in the uv-plane.
  4. Find the Jacobian J(u,v).
  5. In the integrand, replace the variables to obtain the new integrand.
  6. Replace dydx or dxdy, whichever occurs, by J(u,v)dudv.

In the next example, we find a substitution that makes the integrand much simpler to compute.

Evaluating an Integral

Using the change of variables u=xy and v=x+y, evaluate the integral

R(xy)ex2y2dA,

where R is the region bounded by the lines x+y=1 and x+y=3 and the curves x2y2=−1 and x2y2=1 (see the first region in Figure 5.79).

As before, first find the region R and picture the transformation so it becomes easier to obtain the limits of integration after the transformations are made (Figure 5.79).

On the left-hand side of this figure, there is a complex region R in the Cartesian x y-plane bounded by x squared minus y squared = negative 1, x squared minus y squared = 1, x + y = 3, and x + y = 1. Then there is an arrow from this graph to the right-hand side of the figure marked with x = (u + v)/2 and y = (v minus u)/2. On the right-hand side of this figure there is a simpler region S in the Cartesian u v-plane bounded by u v = negative 1, u v = 1, v = 1, and v = 3.
Figure 5.79 Transforming the region R into the region S to simplify the computation of an integral.

Given u=xy and v=x+y, we have x=u+v2 and y=vu2 and hence the transformation to use is T(u,v)=(u+v2,vu2). The lines x+y=1 and x+y=3 become v=1 and v=3, respectively. The curves x2y2=1 and x2y2=−1 become uv=1 and uv=−1, respectively.

Thus we can describe the region S (see the second region Figure 5.79) as

S={(u,v)|1v3,−1vu1v}.

The Jacobian for this transformation is

J(u,v)=(x,y)(u,v)=|xuxvyuyv|=|1/21/21/21/2|=12.

Therefore, by using the transformation T, the integral changes to

R(xy)ex2y2dA=1213−1/v1/vueuvdudv.

Doing the evaluation, we have

1213−1/v1/vueuvdudv=23e0.245.

Using the substitutions x=v and y=u+v, evaluate the integral Rysin(y2x)dA where R is the region bounded by the lines y=x,x=2,andy=0.

12(sin22)

Hint

Sketch a picture and find the limits of integration.

Change of Variables for Triple Integrals

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 G is a region in uvw-space and is mapped to D in xyz-space (Figure 5.80) by a one-to-one C1 transformation T(u,v,w)=(x,y,z) where x=g(u,v,w), y=h(u,v,w), and z=k(u,v,w).

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.
Figure 5.80 A region G in uvw-space mapped to a region D in xyz-space.

Then any function F(x,y,z) defined on D can be thought of as another function H(u,v,w) that is defined on G:

F(x,y,z)=F(g(u,v,w),h(u,v,w),k(u,v,w))=H(u,v,w).

Now we need to define the Jacobian for three variables.

Definition

The Jacobian determinant J(u,v,w) in three variables is defined as follows:

J(u,v,w)=|xuyuzuxvyvzvxwywzw|.

This is also the same as

J(u,v,w)=|xuxvxwyuyvywzuzvzw|.

The Jacobian can also be simply denoted as (x,y,z)(u,v,w).

With the transformations and the Jacobian for three variables, we are ready to establish the theorem that describes change of variables for triple integrals.

Change of Variables for Triple Integrals

Let T(u,v,w)=(x,y,z) where x=g(u,v,w),y=h(u,v,w), and z=k(u,v,w), be a one-to-one C1 transformation, with a nonzero Jacobian, that maps the region G in the uvw-space into the region R in the xyz-space. As in the two-dimensional case, if F is continuous on R, then

RF(x,y,z)dV=GF(g(u,v,w),h(u,v,w),k(u,v,w))|(x,y,z)(u,v,w)|dudvdw=GH(u,v,w)|J(u,v,w)|dudvdw.

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.

Obtaining Formulas in Triple Integrals for Cylindrical and Spherical Coordinates

Derive the formula in triple integrals for

  1. cylindrical and
  2. spherical coordinates.
  1. For cylindrical coordinates, the transformation is T(r,θ,z)=(x,y,z) from the Cartesian rθz-plane to the Cartesian xyz-plane (Figure 5.81). Here x=rcosθ, y=rsinθ, and z=z. The Jacobian for the transformation is
    J(r,θ,z)=(x,y,z)(r,θ,z)=|xrxθxzyryθyzzrzθzz|=|cosθrsinθ0sinθrcosθ0001|=rcos2θ+rsin2θ=r(cos2θ+sin2θ)=r.

    We know that r0, so |J(r,θ,z)|=r. Then the triple integral is
    Df(x,y,z)dV=Gf(rcosθ,rsinθ,z)rdrdθdz.

    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.
    Figure 5.81 The transformation from rectangular coordinates to cylindrical coordinates can be treated as a change of variables from region G in rθz-space to region D in xyz-space.
  2. For spherical coordinates, the transformation is T(ρ,θ,φ)=(x,y,z) from the Cartesian pθφ-plane to the Cartesian xyz-plane (Figure 5.82). Here x=ρsinφcosθ, y=ρsinφsinθ, and z=ρcosφ. The Jacobian for the transformation is
    J(ρ,θ,φ)=(x,y,z)(ρ,θ,φ)=|xρxθxφyρyθyφzρzθzφ|=|sinφcosθρsinφsinθρcosφcosθsinφsinθρsinφcosθρcosφsinθcosϕ0ρsinφ|.

    Expanding the determinant with respect to the third row:
    =cosφ|ρsinφsinθρcosφcosθρsinφsinθρcosφsinθ|ρsinφ|sinφcosθρsinφsinθsinφsinθρsinφcosθ|=cosφ(ρ2sinφcosφsin2θρ2sinφcosφcos2θ)ρsinφ(ρsin2φcos2θ+ρsin2φsin2θ)=ρ2sinφcos2φ(sin2θ+cos2θ)ρ2sinφsin2φ(sin2θ+cos2θ)=ρ2sinφcos2φρ2sinφsin2φ=ρ2sinφ(cos2φ+sin2φ)=ρ2sinφ.

    Since 0φπ, we must have sinφ0. Thus |J(ρ,θ,φ)|=|ρ2sinφ|=ρ2sinφ.
    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.
    Figure 5.82 The transformation from rectangular coordinates to spherical coordinates can be treated as a change of variables from region G in ρθφ-space to region D in xyz-space.

    Then the triple integral becomes
    Df(x,y,z)dV=Gf(ρsinφcosθ,ρsinφsinθ,ρcosφ)ρ2sinφdρdφdθ.

Let’s try another example with a different substitution.

Evaluating a Triple Integral with a Change of Variables

Evaluate the triple integral

0304y/2(y/2)+1(x+z3)dxdydz

in xyz-space by using the transformation

u=(2xy)/2,v=y/2,andw=z/3.

Then integrate over an appropriate region in uvw-space.

As before, some kind of sketch of the region G in xyz-space over which we have to perform the integration can help identify the region D in uvw-space (Figure 5.83). Clearly G in xyz-space is bounded by the planes x=y/2,x=(y/2)+1,y=0, y=4, z=0,andz=4. We also know that we have to use u=(2xy)/2,v=y/2,andw=z/3 for the transformations. We need to solve for x,y,andz. Here we find that x=u+v, y=2v, and z=3w.

Using elementary algebra, we can find the corresponding surfaces for the region G and the limits of integration in uvw-space. It is convenient to list these equations in a table.

Equations in xyz for the region DCorresponding equations in uvw for the region GLimits for the integration in uvw
x=y/2u+v=2v/2=vu=0
x=y/2u+v=(2v/2)+1=v+1u=1
y=02v=0v=0
y=42v=4v=2
z=03w=0w=0
z=33w=3w=1

This table consists of three columns and seven rows. The first row is a header row and reads Equations in x y z for the region D, Corresponding equations in u v w for the region G, and Limits for the integration in u v w. The first column down reads x = y/2, x = (y/2) + 1, y =, y = 4, z = 0, and z = 3. The second column down reads u + v = 2v/2 = v, u + v = (2v/2) + 1 = v + 1, 2v = 0, 2v = 4, 3w = 0, and 3w = 3. The last column reads u = 0, u =1, v = 0, v = 2, w = 0, and w = 1.

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.
Figure 5.83 The region G in uvw-space is transformed to region D in xyz-space.

Now we can calculate the Jacobian for the transformation:

J(u,v,w)=|xuxvxwyuyvywzuzvzw|=|110020003|=6.

The function to be integrated becomes

f(x,y,z)=x+z3=u+v+3w3=u+v+w.

We are now ready to put everything together and complete the problem.

0304y/2(y/2)+1(x+z3)dxdydz=010201(u+v+w)|J(u,v,w)|dudvdw=010201(u+v+w)|6|dudvdw=6010201(u+v+w)dudvdw=60102[u22+vu+wu]01dvdw=60102(12+v+w)dvdw=601[12v+v22+wv]02dw=601(3+2w)dw=6[3w+w2]01=24.

Let D be the region in xyz-space defined by 1x2,0xy2,and0z1.

Evaluate D(x2y+3xyz)dxdydz by using the transformation u=x,v=xy, and w=3z.

030212(v3+vw3u)dudvdw=2+ln8

Hint

Make a table for each surface of the regions and decide on the limits, as shown in the example.

Key Concepts

In the following exercises, the function T:SR,T(u,v)=(x,y) on the region S={(u,v)|0u1,0v1} bounded by the unit square is given, where R2 is the image of S under T.

  1. Justify that the function T is a C1 transformation.
  2. Find the images of the vertices of the unit square S through the function T.
  3. Determine the image R of the unit square S and graph it.

x=2u,y=3v

x=u2,y=v3

a. T(u,v)=(g(u,v),h(u,v)),x=g(u,v)=u2 and y=h(u,v)=v3. The functions g and h are continuous and differentiable, and the partial derivatives gu(u,v)=12, gv(u,v)=0,hu(u,v)=0andhv(u,v)=13 are continuous on S; b. T(0,0)=(0,0), T(1,0)=(12,0),T(0,1)=(0,13), and T(1,1)=(12,13); c. R is the rectangle of vertices (0,0),(12,0),(12,13),and(0,13) in the xy-plane; the following figure.

A rectangle with one corner at the origin, horizontal length 0.5, and vertical height 0.34.

x=uv,y=u+v

x=2uv,y=u+2v

a. T(u,v)=(g(u,v),h(u,v)),x=g(u,v)=2uv, and y=h(u,v)=u+2v. The functions g and h are continuous and differentiable, and the partial derivatives gu(u,v)=2, gv(u,v)=−1, hu(u,v)=1, and hv(u,v)=2 are continuous on S; b. T(0,0)=(0,0), T(1,0)=(2,1), T(0,1)=(−1,2), and T(1,1)=(1,3); c. R is the sqaure of vertices (0,0),(2,1),(1,3),and(−1,2) in the xy-plane; see the following figure.

A square of side length square root of 5 with one corner at the origin and another at (2, 1).

x=u2,y=v2

x=u3,y=v3

a. T(u,v)=(g(u,v),h(u,v)),x=g(u,v)=u3, and y=h(u,v)=v3. The functions g and h are continuous and differentiable, and the partial derivatives gu(u,v)=3u2, gv(u,v)=0, hu(u,v)=0, and hv(u,v)=3v2 are continuous on S; b. T(0,0)=(0,0), T(1,0)=(1,0), T(0,1)=(0,1), and T(1,1)=(1,1); c. R is the unit square in the xy-plane; see the following figure.

A Cartesian coordinate graph with axes ranging from 0 to 1 on both x (horizontal) and y (vertical) axes, displaying an empty plot area with no data or functions presented.

In the following exercises, determine whether the transformations T:SR are one-to-one or not.

x=u2,y=v2,whereS is the rectangle of vertices (−1,0),(1,0),(1,1),and(−1,1).

x=u4,y=u2+v,whereS is the triangle of vertices (−2,0),(2,0),and(0,2).

T is not one-to-one: two points of S have the same image. Indeed, T(−2,0)=T(2,0)=(16,4).

x=2u,y=3v,whereS is the square of vertices (−1,1),(−1,−1),(1,−1),and(1,1).

T(u,v)=(2uv,u), where S is the triangle of vertices (−1,1),(−1,−1),and(1,−1).

T is one-to-one: We argue by contradiction. T(u1,v1)=T(u2,v2) implies 2u1v1=2u2v2 and u1=u2. Thus, u1=u2 and v1=v2.

x=u+v+w,y=u+v,z=w, where S=R=3.

x=u2+v+w,y=u2+v,z=w, where S=R=3.

T is not one-to-one: T(1,v,w)=(−1,v,w)

In the following exercises, the transformations T:SR are one-to-one. Find their related inverse transformations T−1:RS.

x=4u,y=5v, where S=R=2.

x=u+2v,y=u+v, where S=R=2.

u=x2y3,v=x+y3

x=e2u+v,y=euv, where S=2 and R={(x,y)|x>0,y>0}

x=lnu,y=ln(uv), where S={(u,v)|u>0,v>0} and R=2.

u=ex,v=ex+y

x=u+v+w,y=3v,z=2w, where S=R=3.

x=u+v,y=v+w,z=u+w, where S=R=3.

u=xy+z2,v=x+yz2,w=x+y+z2

In the following exercises, the transformation T:SR,T(u,v)=(x,y) and the region R2 are given. Find the region S2.

x=au,y=bv,R={(x,y)|x2+y2a2b2}, where a,b>0

x=au,y=bv,R={(x,y)|x2a2+y2b21}, where a,b>0

S={(u,v)|u2+v21}

x=ua,y=vb,z=wc, R={(x,y)|x2+y2+z21}, where a,b,c>0

x=au,y=bv,z=cw,R={(x,y)|x2a2y2b2z2c21,z>0}, where a,b,c>0

R={(u,v,w)|u2v2w21,w>0}

In the following exercises, find the Jacobian J of the transformation.

x=u+2v,y=u+v

x=u32,y=vu2

32

x=e2uv,y=eu+v

x=uev,y=ev

−1

x=ucos(ev),y=usin(ev)

x=vsin(u2),y=vcos(u2)

2uv

x=ucoshv,y=usinhv,z=w

x=vcosh(1u),y=vsinh(1u),z=u+w2

2vwu2or2wu+w2u2

x=u+v,y=v+w,z=u

x=uv,y=u+v,z=u+v+w

2

The triangular region R with the vertices (0,0),(1,1),and(1,2) is shown in the following figure.

A triangle with corners at the origin, (1, 1), and (1, 2).
  1. Find a transformation T:SR, T(u,v)=(x,y)=(au+bv,cu+dv), where a,b,c, and d are real numbers with adbc0 such that T−1(0,0)=(0,0),T−1(1,1)=(1,0), and T−1(1,2)=(0,1).
  2. Use the transformation T to find the area A(R) of the region R.

The triangular region R with the vertices (0,0),(2,0),and(1,3) is shown in the following figure.

A triangle with corners at the origin, (2, 0), and (1, 3).
  1. Find a transformation T:SR, T(u,v)=(x,y)=(au+bv,cu+dv), where a,b,c and d are real numbers with adbc0 such that T−1(0,0)=(0,0), T−1(2,0)=(1,0), and T−1(1,3)=(0,1).
  2. Use the transformation T to find the area A(R) of the region R.

a. T(u,v)=(2u+v,3v); b. The area of R is
A(R)=03y/3(6y)/3dxdy=0101u|(x,y)(u,v)|dvdu=0101u6dvdu=3.

In the following exercises, use the transformation u=yx,v=y, to evaluate the integrals on the parallelogram R of vertices (0,0),(1,0),(2,1),and(1,1) shown in the following figure.

A rhombus with corners at the origin, (1, 0), (1, 1), and (2, 1).

R(yx)dA

R(y2xy)dA

14

In the following exercises, use the transformation yx=u,x+y=v to evaluate the integrals on the square R determined by the lines y=x,y=x+2,y=x+2, and y=x shown in the following figure.

A square with side lengths square root of 2 rotated 45 degrees with one corner at the origin and another at (1, 1).

Rex+ydA

Rsin(xy)dA

−1+cos2

In the following exercises, use the transformation x=u,5y=v to evaluate the integrals on the region R bounded by the ellipse x2+25y2=1 shown in the following figure.

An ellipse with center at the origin, major axis 2, and minor 0.4.

Rx2+25y2dA

R(x2+25y2)2dA

π15

In the following exercises, use the transformation u=x+y,v=xy to evaluate the integrals on the trapezoidal region R determined by the points (1,0),(2,0),(0,2),and(0,1) shown in the following figure.

A trapezoid with corners at (1, 0), (0, 1), (0, 2), and (2, 0).

R(x22xy+y2)ex+ydA

R(x3+3x2y+3xy2+y3)dA

315

The circular annulus sector R bounded by the circles 4x2+4y2=1 and 9x2+9y2=64, the line x=y3, and the y-axis is shown in the following figure. Find a transformation T from a rectangular region S in the rθ-plane to the region R in the xy-plane. Graph S.

In the first quadrant, a section of an annulus described by an inner radius of 0.5, outer radius slightly more than 2.5, and center the origin. There is a line dividing this annulus that comes from approximately a 30 degree angle. The portion corresponding to 60 degrees is shaded.

The solid R bounded by the circular cylinder x2+y2=9 and the planes z=0,z=1, x=0,andy=0 is shown in the following figure. Find a transformation T from a cylindrical box S in rθz-space to the solid R in xyz-space.

A quarter of a cylinder with height 1 and radius 3. The center axis is the z axis.

T(r,θ,z)=(rcosθ,rsinθ,z);S=[0,3]×[0,π2]×[0,1] in the rθz-space

Show that Rf(x23+y23)dA=2π1501f(ρ)ρdρ, where f is a continuous function on [0,1] and R is the region bounded by the ellipse 5x2+3y2=15.

Show that Rf(16x2+4y2+z2)dV=π201f(ρ)ρ2dρ, where f is a continuous function on [0,1] and R is the region bounded by the ellipsoid 16x2+4y2+z2=1.

[T] Find the area of the region bounded by the curves xy=1,xy=3,y=2x, and y=3x by using the transformation u=xy and v=yx. Use a computer algebra system (CAS) to graph the boundary curves of the region R.

[T] Find the area of the region bounded by the curves xy=2,xy=3,y=x, and y=2x by using the transformation u=xy and v=yx.

The area of the region is  ln22.  A graph of the region is:

Four lines are drawn, namely, y = 3, y = 2, y = 3/(x squared), and y = 2/(x squared). The lines y = 3 and y = 2 are parallel to each other. The lines y = 3/(x squared) and y = 2/(x squared) are curves that run somewhat parallel to each other.

Evaluate the triple integral 0112zz+1(y+1)dxdydz by using the transformation u=xz, v=3y,andw=z2.

Evaluate the triple integral 02463z3z+2(54y)dxdzdy by using the transformation u=x3z,v=4y,andw=z.

8

A transformation T:22,T(u,v)=(x,y) of the form x=au+bv,y=cu+dv, where a,b,c,andd are real numbers, is called linear. Show that a linear transformation for which adbc0 maps parallelograms to parallelograms.

The transformation Tθ:22,Tθ(u,v)=(x,y), where x=ucosθvsinθ, y=usinθ+vcosθ, is called a rotation of angle θ. Show that the inverse transformation of Tθ satisfies Tθ−1=Tθ, where Tθ is the rotation of angle θ.

[T] Use the results of the previous exercise to find the region S in the uv-plane whose image through a rotation of angle π4 is the region R enclosed by the ellipse x2+4y2=1. Use a CAS to answer the following questions.

  1. Graph the region S.
  2. Evaluate the integral Se−2uvdudv. Round your answer to two decimal places.

[T] The transformations Ti:22, i=1,…,4, defined by T1(u,v)=(u,v), T2(u,v)=(u,v),T3(u,v)=(u,v), and T4(u,v)=(v,u) are called reflections about the x-axis,y-axis, origin, and the line y=x, respectively.

  1. Find the image of the region S={(u,v)|u2+v22u4v+10} in the xy-plane through the transformation T1T2T3T4.
  2. Use a CAS to graph the image of S.
  3. Evaluate the integral Ssin(u2)dudv by using a CAS. Round your answer to two decimal places.

a. R={(x,y)|y2+x22y4x+10}; b. R is graphed in the following figure;

A circle with radius 2 and center (2, 1).


c. 3.16

[T] The transformation Tk,1,1:33,Tk,1,1(u,v,w)=(x,y,z) of the form x=ku, y=v,z=w, where k1 is a positive real number, is called a stretch if k>1 and a compression if 0<k<1 in the x-direction. Use a CAS to evaluate the integral Se(4x2+9y2+25z2)dxdydz on the solid S={(x,y,z)|4x2+9y2+25z21} by considering the compression T2,3,5(u,v,w)=(x,y,z) defined by x=u2,y=v3, and z=w5. Round your answer to four decimal places.

[T] The transformation Ta,0:22,Ta,0(u,v)=(u+av,v), where a0 is a real number, is called a shear in the x-direction. The transformation, T0,b:22,T0,b(u,v)=(u,bu+v), where b0 is a real number, is called a shear in the y-direction.

  1. Find transformations T0,2T3,0.
  2. Find the image R of the trapezoidal region S bounded by u=0,v=0,v=1, and v=2u through the transformation T0,2T3,0.
  3. Use a CAS to graph the image R in the xy-plane.
  4. Find the area of the region R by using the area of region S.

a. T0,2T3,0(u,v)=(u+3v,2u+7v); b. The image R is the quadrilateral of vertices (0,0),(3,7),(2,4),and(4,9); c. S is graphed in the following figure;

A four-sided figure with points the origin, (2, 4), (4, 9), and (3, 7).


d. 32

Use the transformation, x=au,y=av,z=cw and spherical coordinates to show that the volume of a region bounded by the spheroid x2+y2a2+z2c2=1 is 4πa2c3.

Find the volume of a football whose shape is a spheroid x2+y2a2+z2c2=1 whose length from tip to tip is 11 inches and circumference at the center is 22 inches. Round your answer to two decimal places.

26623π282.45in3

[T] Lamé ovals (or superellipses) are plane curves of equations (xa)n+(yb)n=1, where a, b, and n are positive real numbers.

  1. Use a CAS to graph the regions R bounded by Lamé ovals for a=1,b=2,n=4 and n=6, respectively.
  2. Find the transformations that map the region R bounded by the Lamé oval x4+y4=1, also called a squircle and graphed in the following figure, into the unit disk.
    A square of side length 2 with rounded corners.
  3. Use a CAS to find an approximation of the area A(R) of the region R bounded by x4+y4=1. Round your answer to two decimal places.

[T] Lamé ovals have been consistently used by designers and architects. For instance, Gerald Robinson, a Canadian architect, has designed a parking garage in a shopping center in Peterborough, Ontario, in the shape of a superellipse of the equation (xa)n+(yb)n=1 with ab=97 and n=e. Use a CAS to find an approximation of the area of the parking garage in the case a=900 yards, b=700 yards, and n=2.72.

A(R)=2,177,216 yd2

Chapter Review Exercises

True or False? Justify your answer with a proof or a counterexample.

abcdf(x,y)dydx=cdabf(x,y)dydx

Fubini’s theorem can be extended to three dimensions, as long as f is continuous in all variables.

True.

The integral 02π01r1dzdrdθ represents the volume of a right cone.

The Jacobian of the transformation x=u22v,y=3v2uv is given by −4u2+6u+4v.

False.

Evaluate the following integrals.

R(5x3y2y2)dA,R={(x,y)|0x2,1y4}

Dy3x2+1dA,D={(x,y)|0x1,xyx}

0

Dsin(x2+y2)dA where D is a disk of radius 2 centered at the origin

01y1xyex2dxdy

14

−110z0xz6dydxdz

R3ydV, where R={(x,y,z)|0x1,0yx,0z9y2}

1.475

0202πr1rdzdθdr

02π0π/213ρ2sin(φ)dρdφdθ

523π

011x21x21x2y21x2y2dzdydx

For the following problems, find the specified area or volume.

The area of region enclosed by one petal of r=cos(4θ).

π16

The volume of the solid that lies between the paraboloid z=2x2+2y2 and the plane z=8.

The volume of the solid bounded by the cylinder x2+y2=16 and from z=1 and z+x=2.

16π50.265

The volume of the intersection between two spheres of radius 1, the top whose center is (0,0,0.25) and the bottom, which is centered at (0,0,0).

For the following problems, find the center of mass of the region.

ρ(x,y)=xy on the circle with radius 1 in the first quadrant only.

(815,815)

ρ(x,y)=(y+1)x in the region bounded by y=ex, y=0, and x=1.

ρ(x,y,z)=z on the inverted cone with radius 2 and height 2.

(0,0,45)

The volume an ice cream cone that is given by the solid above z=(x2+y2) and below z2+x2+y2=z.

The following problems examine Mount Holly in the state of Michigan. Mount Holly is a landfill that was converted into a ski resort. The shape of Mount Holly can be approximated by a right circular cone of height 1100 ft and radius 6000 ft.

If the compacted trash used to build Mount Holly on average has a density 400lb/ft3, find the amount of work required to build the mountain.

1.452π×1015 ft-lb

In reality, it is very likely that the trash at the bottom of Mount Holly has become more compacted with all the weight of the above trash. Consider a density function with respect to height: the density at the top of the mountain is still density 400lb/ft3 and the density increases. Every 100 feet deeper, the density doubles. What is the total weight of Mount Holly?

The following problems consider the temperature and density of Earth’s layers.

[T] The temperature of Earth’s layers is exhibited in the table below. Use your calculator to fit a polynomial of degree 3 to the temperature along the radius of the Earth. Then find the average temperature of Earth. (Hint: begin at 0 in the inner core and increase outward toward the surface)

Source: http://www.enchantedlearning.com/subjects/astronomy/planets/earth/Inside.shtml
LayerDepth from center (km)Temperature °C
Rocky Crust0 to 400
Upper Mantle40 to 150870
Mantle400 to 650870
Inner Mantel650 to 2700870
Molten Outer Core2890 to 51504300
Inner Core5150 to 63787200

The layers from crust to core, starting at the top: Rock crust, 0 to 40 km, 0 degrees C; Upper mantle, 40 to 150 km, 870 degrees C. These two layers are grouped together as the Lithosphere, and the line between them is marked as the Mohorovicic Discontinuity. Then the next layer is the Mantle, 400 to 650 km, 870 degrees C; Inner mantle (semi-rigid), 650 to 2700 km, 870 degrees C. Then there is a short region between 2700 and 2890 km marked the Gutenberg Discontinuity, which has temperature 3700 degrees C. The next layer is the Molten Outer Core (Iron/nickel) from 2890 to 5150, 4300 degrees C; and finally there is the inner core (Iron/nickel) from 5150 to 6378, 7200 degrees C.

y=−1.238×10−7x3+0.001196x23.666x+7208; average temperature approximately 2800°C

[T] The density of Earth’s layers is displayed in the table below. Using your calculator or a computer program, find the best-fit quadratic equation to the density. Using this equation, find the total mass of Earth.

Source: http://hyperphysics.phy-astr.gsu.edu/hbase/geophys/earthstruct.html
LayerDepth from center (km)Density (g/cm3)
Inner Core012.95
Outer Core122811.05
Mantle34885.00
Upper Mantle63383.90
Crust63782.55

This table consists of three columns and six rows. The first row is a header row, and from left to right, it reads Layer, Depth from center (km), and Density (g/cm3). The first informational row reads Inner Core, 0, 12.95. The second row reads Outer Core, 1228, 11.05. The third row reads Mantle, 2488, 5.00. The fourth row reads Upper Mantle, 6338, 3.90. And the final row reads Crust, 6378, 2.55.

The following problems concern the Theorem of Pappus (see Moments and Centers of Mass for a refresher), a method for calculating volume using centroids. Assuming a region R, when you revolve around the x-axis the volume is given by Vx=2πAy, and when you revolve around the y-axis the volume is given by Vy=2πAx, where A is the area of R. Consider the region bounded by x2+y2=1 and above y=x+1.

Find the volume when you revolve the region around the x-axis.

π3

Find the volume when you revolve the region around the y-axis.

Glossary

Jacobian
the Jacobian J(u,v) in two variables is a 2×2 determinant:
J(u,v)=|xuyuxvyv|;

the Jacobian J(u,v,w) in three variables is a 3×3 determinant:
J(u,v,w)=|xuyuzuxvyvzvxwywzw|
one-to-one transformation
a transformation T:GR defined as T(u,v)=(x,y) is said to be one-to-one if no two points map to the same image point
planar transformation
a function T that transforms a region G in one plane into a region R in another plane by a change of variables
transformation
a function that transforms a region G in one plane into a region R in another plane by a change of variables