Loading [MathJax]/jax/element/mml/optable/MathOperators.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Mathematics LibreTexts

15.7: Change of Variables in Multiple Integrals

  • Gilbert Strang & Edwin “Jed” Herman
  • OpenStax

( \newcommand{\kernel}{\mathrm{null}\,}\)

Learning Objectives
  • Determine the image of a region under a given transformation of variables.
  • Compute the Jacobian of a given transformation.
  • Evaluate a double integral using a change of variables.
  • Evaluate a triple integral using a change of variables.

Recall from Substitution Rule the method of integration by substitution. When evaluating an integral such as

32x(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

5012u5du

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,

baf(x)dx=dcf(g(u))g(u)du,

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

A similar result occurs in double integrals when we substitute

  • x=f(r,θ)=rcosθ
  • y=g(r,θ)=rsinθ, and
  • dA=dxdy=rdrdθ.

Then we get

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 R^2. For example, Figure \PageIndex{1} 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 g_u, \, g_v, \, h_u, and h_v 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 \PageIndex{1}: The transformation of a region G in the uv-plane into a region R in the xy-plane.
Definition: one-to-one transformation

A transformation T: \, G \rightarrow R, 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 s a one-to-one transformation, we assume T(u_1,v_1) = T(u_2, v_2) and show that as a consequence we obtain(u_1,v_1) = (u_2, v_2). 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^{-1} \circ T and T \circ T^{-1} are identity functions.

Figure \PageIndex{2} 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.

Example \PageIndex{1A}: Determining How the Transformation Works

Suppose a transformation T is defined as T(r,\theta) = (x,y) where x = r \, \cos \, \theta, \, y = r \, \sin \, \theta. Find the image of the polar rectangle G = \{(r,\theta) | 0 \leq r \leq 1, \, 0 \leq \theta \leq \pi/2\} in the r\theta-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).

Solution

Since r varies from 0 to 1 in the r\theta-plane, we have a circular disc of radius 0 to 1 in the xy-plane. Because \theta varies from 0 to \pi/2 in the r\theta-plane, we end up getting a quarter circle of radius 1 in the first quadrant of the xy-plane (Figure \PageIndex{2}). Hence R is a quarter circle bounded by x^2 + y^2 = 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 \PageIndex{2}: A rectangle in the r\theta-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(r_1,\theta_1) = T(r_2, \theta_2) and show as a consequence that (r_1,\theta_1) = (r_2, \theta_2). In this case, we have

T(r_1,\theta_1) = T(r_2, \theta_2), \nonumber

(x_1,y_1) = (x_1,y_1), \nonumber

(r_1 \cos \, \theta_1, r_1 \sin \, \theta_1) = (r_2 \cos \, \theta_2, r_2 \sin \, \theta_2), \nonumber

r_1 \cos \, \theta_1 = r_2 \cos \, \theta_2, \, r_1 \sin \, \theta_1 = r_2 \sin \, \theta_2. \nonumber

Dividing, we obtain

\frac{r_1 \cos \, \theta_1}{r_1 \sin \, \theta_1} = \frac{ r_2 \cos \, \theta_2}{ r_2 \sin \, \theta_2} \nonumber

\frac{\cos \, \theta_1}{\sin \, \theta_1} = \frac{\cos \, \theta_2}{\sin \, \theta_2} \nonumber

\tan \, \theta_1 = \tan \, \theta_2 \nonumber

\theta_1 = \theta_2 \nonumber

since the tangent function is one-one function in the interval 0 \leq \theta \leq \pi/2. Also, since 0 \leq r \leq 1, we have r_1 = r_2, \, \theta_1 = \theta_2. Therefore, (r_1,\theta_1) = (r_2, \theta_2) and T is a one-to-one transformation from G to R.

To find T^{-1}(x,y) solve for r,\theta in terms of x,y. We already know that r^2 = x^2 + y^2 and \tan \, \theta = \frac{y}{x}. Thus T^{-1}(x,y) = (r,\theta) is defined as r = \sqrt{x^2 + y^2} and \tan^{-1} \left(\frac{y}{x}\right).

Example \PageIndex{1B}: Finding the Image under T

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

Solution

The triangle and its image are shown in Figure \PageIndex{3}. 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 \PageIndex{3}: A triangular region in the uv-plane is transformed into an image in the xy-plane.
  • For the side A: \, u = 0, \, 0 \leq v \leq 1 transforms to x = -v^2, \, y = 0 so this is the side A' that joins (-1,0) and (0,0).
  • For the side B: \, u = v, \, 0 \leq u \leq 1 transforms to x = 0, \, y = u^2 so this is the side B' that joins (0,0) and (0,1).
  • For the side C: \, 0 \leq u \leq 1, \, v = 1 transforms to x = u^2 - 1, \, y = u (hence x = y^2 - 1 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.

Exercise \PageIndex{1}

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) : \, 0 \leq u \leq 1, \, 0 \leq v \leq 2\} 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).

Hint

Follow the steps of Example \PageIndex{1B}.

Answer

T^{-1} (x,y) = (u,v) where u = \frac{3x-y}{3} and v = \frac{y}{3}

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 g_u, g_v, h_u and h_v exist and are also continuous. A transformation that has this property is called a C^1 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 C^1 transformation. We want to see how it transforms a small rectangular region S, \, \Delta u units by \Delta v units, in the uv-plane (Figure \PageIndex{4}).

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 \PageIndex{4}: 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 (u_0,v_0)is the coordinate of the point at the lower left corner that mapped to (x_0,y_0) = T(u_0,v_0) The line v = v_0 maps to the image curve with vector function r(u,v_0), and the tangent vector at (x_0,y_0) to the image curve is

r_u = g_u (u_0,v_0)i + h_v (u_0,v_0)j = \frac{\partial x}{\partial u}i + \frac{\partial y}{\partial u}j. \nonumber

Similarly, the line u = u_0 maps to the image curve with vector function r(u_0,v), and the tangent vector at (x_0,y_0) to the image curve is

r_v = g_v (u_0,v_0)i + h_u (u_0,v_0)j = \frac{\partial x}{\partial v}i + \frac{\partial y}{\partial v}j. \nonumber

Now, note that

r_u = \lim_{\Delta u \rightarrow 0} \frac{r (u_0 + \Delta u, v_0) - r ( u_0,v_0)}{\Delta u}\, so \, r (u_0 + \Delta u,v_0) - r(u_0,v_0) \approx \Delta u r_u. \nonumber

Similarly,

r_v = \lim_{\Delta v \rightarrow 0} \frac{r (u_0,v_0 + \Delta v) - r ( u_0,v_0)}{\Delta v}\, so \, r (u_0,v_0 + \Delta v) - r(u_0,v_0) \approx \Delta v r_v. \nonumber

This allows us to estimate the area \Delta A of the image R by finding the area of the parallelogram formed by the sides \Delta vr_v and \Delta ur_u. By using the cross product of these two vectors by adding the kth component as 0, the area \Delta A of the image R (refer to The Cross Product) is approximately |\Delta ur_u \times \Delta v r_v| = |r_u \times r_v|\Delta u \Delta v. In determinant form, the cross product is

r_u \times r_v = \begin{vmatrix} i & j & k \\ \frac{\partial x}{\partial u} & \frac{\partial y}{\partial u} & 0 \\ \frac{\partial x}{\partial v} & \frac{\partial y}{\partial v} & 0 \end{vmatrix} = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial y}{\partial u} \\ \dfrac{\partial x}{\partial v} & \dfrac{\partial y}{\partial v} \end{vmatrix} k = \left(\frac{\partial x}{\partial u} \frac{\partial y}{\partial v} - \frac{\partial x}{\partial v} \frac{\partial y}{\partial u}\right)k \nonumber

Since |k| = 1, we have

\Delta A \approx |r_u \times r_v| \Delta u \Delta v = \left( \frac{\partial x}{\partial u}\frac{\partial y}{\partial v} - \frac{\partial x}{\partial v} \frac{\partial y}{\partial u}\right) \Delta u \Delta v.

Definition: Jacobian

The Jacobian of the C^1 transformation T(u,v) = (g(u,v), \, h(u,v)) is denoted by J(u,v) and is defined by the 2 \times 2 determinant

J(u,v) = \left|\frac{\partial (x,y)}{\partial (u,v)} \right| = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial y}{\partial u} \\ \dfrac{\partial x}{\partial v} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \left( \frac{\partial x}{\partial u}\frac{\partial y}{\partial v} - \frac{\partial x}{\partial v} \frac{\partial y}{\partial u}\right). \nonumber

Using the definition, we have

\Delta A \approx J(u,v) \Delta u \Delta v = \left|\frac{\partial (x,y)}{\partial (u,v)}\right| \Delta u \Delta v. \nonumber

Note that the Jacobian is frequently denoted simply by

J(u,v) = \frac{\partial (x,y)}{\partial (u,v)}. \nonumber

Note also that

\begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial y}{\partial u} \nonumber \\ \dfrac{\partial x}{\partial v} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \left( \frac{\partial x}{\partial u}\frac{\partial y}{\partial v} - \frac{\partial x}{\partial v} \frac{\partial y}{\partial u}\right) = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} \nonumber \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix} . \nonumber

Hence the notation J(u,v) = \frac{\partial(x,y)}{\partial(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.

Example \PageIndex{2A}: Finding the Jacobian

Find the Jacobian of the transformation given in Example \PageIndex{1A}.

Solution

The transformation in the example is T(r,\theta) = ( r \, \cos \, \theta, \, r \, \sin \, \theta) where x = r \, \cos \, \theta and y = r \, \sin \, \theta. Thus the Jacobian is

J(r, \theta) = \frac{\partial(x,y)}{\partial(r,\theta)} = \begin{vmatrix} \dfrac{\partial x}{\partial r} & \dfrac{\partial x}{\partial \theta} \\ \dfrac{\partial y}{\partial r} & \dfrac{\partial y}{\partial \theta} \end{vmatrix} = \begin{vmatrix} \cos \theta & -r\sin \theta \\ \sin \theta & r\cos\theta \end{vmatrix} = r \, \cos^2\theta + r \, \sin^2\theta = r ( \cos^2\theta + \sin^2\theta) = r. \nonumber

Example \PageIndex{2B}: Finding the Jacobian

Find the Jacobian of the transformation given in Example \PageIndex{1B}.

Solution

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

J(u,v) = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \begin{vmatrix} 2u & -2v \\ v & u \end{vmatrix} = 2u^2 + 2v^2. \nonumber

Exercise \PageIndex{2}

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

Hint

Follow the steps in the previous two examples.

Answer

J(u,v) = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} \nonumber \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \begin{vmatrix} 1 & 1 \nonumber \\ 0 & 2 \end{vmatrix} = 2 \nonumber

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 \Delta A in the xy-plane is related to the area formed by the product \Delta u \Delta v in the uv-plane by the approximation

\Delta A \approx J(u,v) \Delta u, \, \Delta v. \nonumber

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

\iint_R f(x,y)fA = \lim_{m,n \rightarrow \infty} \sum_{i=1}^m \sum_{j=1}^n f(x_{ij}, y_{ij}) \Delta A. \nonumber

Referring to Figure \PageIndex{5}, observe that we divided the region S in the uv-plane into small subrectangles S_{ij} and we let the subrectangles R_{ij} in the xy-plane be the images of S_{ij} 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 \PageIndex{5}: The subrectangles S_{ij} in the uv-plane transform into subrectangles R_{ij} in the xy-plane.

Then the double integral becomes

\iint_R = f(x,y)dA = \lim_{m,n \rightarrow \infty} \sum_{i=1}^m \sum_{j=1}^n f(x_{ij}, y_{ij}) \Delta A = \lim_{m,n \rightarrow \infty} \sum_{i=1}^m \sum_{j=1}^n f(g(u_{ij}, v_{ij}), \, h(u_{ij}, v_{ij})) | J(u_{ij}, v_{ij})| \Delta u \Delta v. \nonumber

Notice this is exactly the double Riemann sum for the integral

\iint_S f(g(u,v), \, h(u,v)) \left|\frac{\partial (x,y)}{\partial(u,v)}\right| du \, dv. \nonumber

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 C^1 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

\iint_R f(x,y) dA = \iint_S f(g(u,v), \, h(u,v)) \left|\frac{\partial (x,y)}{\partial(u,v)}\right| du \, dv. \nonumber

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 = dx \, dy = \left|\frac{\partial (x,y)}{\partial (u,v)} \right| du \, dv \nonumber

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.

Example \PageIndex{3}: Changing Variables from Rectangular to Polar Coordinates

Consider the integral

\int_0^2 \int_0^{\sqrt{2x-x^2}} \sqrt{x^2 + y^2} dy \, dx. \nonumber

Use the change of variables x = r \, \cos \, \theta and y = r \, \sin \, \theta, and find the resulting integral.

Solution

First we need to find the region of integration. This region is bounded below by y = 0 and above by y = \sqrt{2x - x^2} (Figure \PageIndex{6}).

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 \PageIndex{6}: Changing a region from rectangular to polar coordinates.

Squaring and collecting terms, we find that the region is the upper half of the circle x^2 + y^2 - 2x = 0, that is y^2 + ( x - 1)^2 = 1. In polar coordinates, the circle is r = 2 \, cos \, \theta so the region of integration in polar coordinates is bounded by 0 \leq r \leq \cos \, \theta and 0 \leq \theta \leq \frac{\pi}{2}.

The Jacobian is J(r, \theta) = r, as shown in Example \PageIndex{2A}. Since r \geq 0, we have |J(r,\theta)| = r.

The integrand \sqrt{x^2 + y^2} changes to r in polar coordinates, so the double iterated integral is

\int_0^2 \int_0^{\sqrt{2x-x^2}} \sqrt{x^2 + y^2} dy \, dx = \int_0^{\pi/2} \int_0^{2 \, cos \, \theta} r | j(r, \theta)|dr \, d\theta = \int_0^{\pi/2} \int_0^{2 \, cos \, \theta} r^2 dr \, d\theta. \nonumber

Exercise \PageIndex{3}

Considering the integral \int_0^1 \int_0^{\sqrt{1-x^2}} (x^2 + y^2) dy \, dx, use the change of variables x = r \, cos \, \theta and y = r \, sin \, \theta and find the resulting integral.

Hint

Follow the steps in the previous example.

Answer

\int_0^{\pi/2} \int_0^1 r^3 dr \, d\theta \nonumber

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.

Example \PageIndex{4}: Changing Variables

Consider the integral \iint_R (x - y) dy \, dx, \nonumber where R is the parallelogram joining the points (1,2), \, (3,4), \, (4,3), and (6,5) (Figure \PageIndex{7}). 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 \PageIndex{7}: The region of integration for the given integral.
Solution

First, we need to understand the region over which we are to integrate. The sides of the parallelogram are x - y + 1, \, x - y - 1 = 0, \, x - 3y + 5 = 0 and x - 3y + 9 = 0 (Figure \PageIndex{8}). Another way to look at them is x - y = -1, \, x - y = 1, \, x - 3y = -5, and x - 3y = 9.

Clearly the parallelogram is bounded by the lines y = x + 1, \, y = x - 1, \, y = \frac{1}{3}(x + 5), and y = \frac{1}{3}(x + 9).

Notice that if we were to make u = x - y and v = x - 3y, then the limits on the integral would be -1 \leq u \leq 1 and -9 \leq v \leq -5.

To solve for x and y, we multiply the first equation by 3 and subtract the second equation, 3u - v = (3x - 3y) - (x - 3y) = 2x. Then we have x = \frac{3u-v}{2}. Moreover, if we simply subtract the second equation from the first, we get u - v = (x - y) - (x - 3y) = 2y and y = \frac{u-v}{2}.

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 \PageIndex{8}: 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) = \left( \frac{3u - v}{2}, \, \frac{u - v}{2} \right) \nonumber and compute the Jacobian J(u,v). We have

J(u,v) = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \begin{vmatrix} 3/2 & -1/2 \nonumber \\ 1/2 & -1/2 \end{vmatrix} = -\frac{3}{4} + \frac{1}{4} = - \frac{1}{2} \nonumber

Therefore, |J(u,v)| = \frac{1}{2}. Also, the original integrand becomes

x - y = \frac{1}{2} [3u - v - u + v] = \frac{1}{2} [3u - u] = \frac{1}{2}[2u] = u. \nonumber

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

\iint_R (x - y) dy \, dx = \int_{-9}^{-5} \int_{-1}^1 J (u,v) u \, du \, dv = \int_{-9}^{-5} \int_{-1}^1\left(\frac{1}{2}\right) u \, du \, dv, \nonumber which is much simpler to compute.

Exercise \PageIndex{4}

Make appropriate changes of variables in the integral \iint_R \frac{4}{(x - y)^2} dy \, dx, \nonumber where R is the trapezoid bounded by the lines x - y = 2, \, x - y = 4, \, x = 0, and y = 0. Write the resulting integral.

Hint

Follow the steps in the previous example.

Answer

x = \frac{1}{2}(v + u) and y = \frac{1}{2} (v - u)

and

\int_{2}^4 \int_{-u}^u \left(\frac{1}{2}\right)\cdot\frac{4}{u^2} \,dv \, du. \nonumber

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

Problem-Solving Strategy: 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 dy \, dx or dx \, dy, whichever occurs, by J(u,v) du \, dv.

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

Example \PageIndex{5}: Evaluating an Integral

Using the change of variables u = x - y and v = x + y, evaluate the integral \iint_R (x - y)e^{x^2-y^2} dA, \nonumber where R is the region bounded by the lines x + y = 1 and x + y = 3 and the curves x^2 - y^2 = -1 and x^2 - y^2 = 1 (see the first region in Figure \PageIndex{9}).

Solution

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 \PageIndex{9}).

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 \PageIndex{9}: Transforming the region R into the region S to simplify the computation of an integral.

Given u = x - y and v = x + y, we have x = \frac{u+v}{2} and y = \frac{v-u}{2} and hence the transformation to use is T(u,v) = \left(\frac{u+v}{2}, \, \frac{v-u}{2}\right). The lines x + y = 1 and x + y = 3 become v = 1 and v = 3, respectively. The curves x^2 - y^2 = 1 and x^2 - y^2 = -1 become uv = 1 and uv = -1, respectively.

Thus we can describe the region S (see the second region Figure \PageIndex{9}) as

S = \left\{ (u,v) | 1 \leq v \leq 3, \, \frac{-1}{v} \leq u \leq \frac{1}{v}\right\}. \nonumber

The Jacobian for this transformation is

J(u,v) = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \begin{vmatrix} 1/2 & 1/2 \\ -1/2 & 1/2 \end{vmatrix} = \frac{1}{2}. \nonumber

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

\iint_R (x - y)e^{x^2-y^2} dA = \frac{1}{2} \int_1^3 \int_{-1/v}^{1/v} ue^{uv} du \, dv. \nonumber

Doing the evaluation, we have

\frac{1}{2} \int_1^3 \int_{-1/v}^{1/v} ue^{uv} du \, dv = \frac{2}{3e} \approx 0.245. \nonumber

Exercise \PageIndex{5}

Using the substitutions x = v and y = \sqrt{u + v}, evaluate the integral \displaystyle\iint_R y \, \sin (y^2 - x) \,dA, where R is the region bounded by the lines y = \sqrt{x}, \, x = 2 and y = 0.

Hint

Sketch a picture and find the limits of integration.

Answer

\frac{1}{2} (\sin 2 - 2)

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 \PageIndex{10}) by a one-to-one C^1 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 \PageIndex{10}: 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). \nonumber

Now we need to define the Jacobian for three variables.

Definition: Jacobian determinant

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

J(u,v,w) = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial y}{\partial u} & \dfrac{\partial z}{\partial u} \\ \dfrac{\partial x}{\partial v} & \dfrac{\partial y}{\partial v} & \dfrac{\partial z}{\partial v} \\ \dfrac{\partial x}{\partial w} & \dfrac{\partial y}{\partial w} & \dfrac{\partial z}{\partial w} \end{vmatrix}. \nonumber

This is also the same as

J(u,v,w) = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} & \dfrac{\partial x}{\partial w} \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} & \dfrac{\partial y}{\partial w} \\ \dfrac{\partial z}{\partial u} & \dfrac{\partial z}{\partial v} & \dfrac{\partial z}{\partial w} \end{vmatrix}. \nonumber

The Jacobian can also be simply denoted as \frac{\partial(x,y,z)}{\partial (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 C^1 transformation, with a nonzero Jacobian, that maps the region G in the uvw-space into the region D in the xyz-space. As in the two-dimensional case, if F is continuous on D, then

\begin{align} \iiint_D F(x,y,z) dV = \iiint_G f(g(u,v,w) \, h(u,v,w), \, k(u,v,w)) \left|\frac{\partial (x,y,z)}{\partial (u,v,w)}\right| du \, dv \, dw \\ = \iiint_G H(u,v,w) | J (u,v,w) | du \, dv \, dw. \end{align} \nonumber

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.

Example \PageIndex{6A}: Obtaining Formulas in Triple Integrals for Cylindrical and Spherical Coordinates

Derive the formula in triple integrals for

  1. cylindrical and
  2. spherical coordinates.
Solution

A.

For cylindrical coordinates, the transformation is T (r, \theta, z) = (x,y,z) from the Cartesian r\theta z-space to the Cartesian xyz-space (Figure \PageIndex{11}). Here x = r \, \cos \, \theta, \, y = r \, \sin \theta and z = z. The Jacobian for the transformation is

J(r,\theta,z) = \frac{\partial (x,y,z)}{\partial (r,\theta,z)} = \begin{vmatrix} \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \theta} & \frac{\partial x}{\partial z} \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \theta} & \frac{\partial y}{\partial z} \\ \frac{\partial z}{\partial r} & \frac{\partial z}{\partial \theta} & \frac{\partial z}{\partial z} \end{vmatrix} \nonumber

\begin{vmatrix} \cos \theta & -r\sin \theta & 0 \\ \sin \theta & r \cos \theta & 0 \\ 0 & 0 & 1 \end{vmatrix} = r \, \cos^2 \theta + r \, \sin^2 \theta = r. \nonumber

We know that r \geq 0, so |J(r,\theta,z)| = r. Then the triple integral is \iiint_D f(x,y,z)dV = \iiint_G f(r \, \cos \theta, \, r \, \sin \theta, \, z) r \, dr \, d\theta \, dz. \nonumber

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 \PageIndex{11}: The transformation from rectangular coordinates to cylindrical coordinates can be treated as a change of variables from region G in r\theta z-space to region D in xyz-space.

B.

For spherical coordinates, the transformation is T(\rho,\theta,\varphi) from the Cartesian \rho\theta\varphi-space to the Cartesian xyz-space (Figure \PageIndex{12}). Here x = \rho \, \sin \varphi \, \cos \theta, \, y = \rho \, \sin \varphi \, \sin \theta, and z = \rho \, \cos \varphi. The Jacobian for the transformation is

J(\rho,\theta,\varphi) = \frac{\partial (x,y,z)}{\partial (\rho,\theta,\varphi)} = \begin{vmatrix} \frac{\partial x}{\partial \rho} & \frac{\partial x}{\partial \theta} & \frac{\partial x}{\partial \varphi} \\ \frac{\partial y}{\partial \rho} & \frac{\partial y}{\partial \theta} & \frac{\partial y}{\partial \varphi} \\ \frac{\partial z}{\partial \rho} & \frac{\partial z}{\partial \theta} & \frac{\partial z}{\partial \varphi} \end{vmatrix} = \begin{vmatrix} \sin \varphi \cos \theta & -\rho \sin \varphi \sin \theta & \rho \cos \varphi \cos \theta \\ \sin \varphi \sin \theta & \rho \sin \varphi \cos \theta & \rho \cos \varphi \sin \theta \\ \cos \varphi & 0 & -\rho \sin \varphi \end{vmatrix}. \nonumber

Expanding the determinant with respect to the third row:

\begin{align*} &= \cos \varphi \begin{vmatrix} -\rho \sin \varphi \sin \theta & \rho \cos \varphi \cos \theta \\ \rho \sin \varphi \cos \theta & \rho \cos \varphi \sin \theta \end{vmatrix} - \rho \sin \varphi \begin{vmatrix} \sin \varphi \cos \theta & - \rho \sin \varphi \sin \theta \\ \sin \varphi \sin \theta & \rho \sin \varphi \cos \theta \end{vmatrix} \\[4pt] &=\cos \varphi (-\rho^2 \sin \varphi \, \cos \varphi \, \sin^2 \theta - \rho^2 \, \sin \varphi \, \cos \varphi \, \cos^2\theta) \\ &\quad -\rho \sin \varphi (\rho \sin^2\varphi \cos^2\theta + \rho \sin^2\varphi \sin^2\theta) \\[4pt] &=-\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) \\[4pt] &= - \rho^2 \sin\varphi \cos^2\varphi - \rho^2 \sin \varphi \sin^2 \varphi \\[4pt] &= -\rho \sin \varphi (\cos^2 \varphi + \sin^2 \varphi) = - \rho^2 \sin \varphi. \end{align*}

Since 0 \leq \varphi \leq \pi, we must have \sin \varphi \geq 0. Thus |J(\rho,\theta, \varphi)| = |-\rho^2 \sin \varphi| = \rho^2 \sin \varphi.

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 \PageIndex{12}: The transformation from rectangular coordinates to spherical coordinates can be treated as a change of variables from region G in \rho\theta\varphi-space to region D in xyz-space.

Then the triple integral becomes

\iiint_D f(x,y,z) dV = \iiint_G f(\rho \, \sin \varphi \, \cos \theta, \, \rho \, \sin \varphi \, \sin \theta, \rho \, \cos \varphi) \rho^2 \sin \varphi \, d\rho \, d\varphi \, d\theta. \nonumber

Let’s try another example with a different substitution.

Example \PageIndex{6B}: Evaluating a Triple Integral with a Change of Variables

Evaluate the triple integral

\int_0^3 \int_0^4 \int_{y/2}^{(y/2)+1} \left(x + \frac{z}{3}\right) dx \, dy \, dz \nonumber

In xyz-space by using the transformation

u = (2x - y) /2, \, v = y/2, and w = z/3.

Then integrate over an appropriate region in uvw-space.

Solution

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 \PageIndex{13}). Clearly G in xyz-space is bounded by the planes x = y/2, \, x = (y/2) + 1, \, y = 0, \, y = 4, \, z = 0, and z = 4. We also know that we have to use u = (2x - y) /2, \, v = y/2, and w = z/3 for the transformations. We need to solve for x,y and z. Here we find that x = u + v, \, y = 2v, and z = 3w.

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 \PageIndex{13}: The region G in uvw-space is transformed to region D in xyz-space.

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 D Corresponding equations in uvw for the region G Limits for the integration in uvw
x = y/2 u + v = 2v/2 = v u = 0
x = y/2 u + v = (2v/2) + 1 = v + 1 u = 1
y = 0 2v = 0 v = 0
y = 4 2v = 4 v = 2
z = 0 3w = 0 w = 0
z = 3 3w = 3 w = 1

Now we can calculate the Jacobian for the transformation:

J(u,v,w) = \begin{vmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} & \dfrac{\partial x}{\partial w} \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} & \dfrac{\partial y}{\partial w} \\ \dfrac{\partial z}{\partial u} & \dfrac{\partial z}{\partial v} & \dfrac{\partial z}{\partial w} \end{vmatrix} = \begin{vmatrix} 1 & 1 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{vmatrix} = 6. \nonumber

The function to be integrated becomes

f(x,y,z) = x + \frac{z}{3} = u + v + \frac{3w}{3} = u + v + w. \nonumber

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

\begin{align*} \int_0^3 \int_0^4 \int_{y/2}^{(y/2)+1} \left(x + \frac{z}{3}\right) dx \, dy \, dz &= \int_0^1 \int_0^2 \int_0^1 (u + v + w) |J (u,v,w)|du \, dv \, dw \\[4pt] &= \int_0^1 \int_0^2 \int_0^1 (u + v + w) |6|du \, dv \, dw \\[4pt] &= 6 \int_0^1 \int_0^2 \int_0^1 (u + v + w) \, du \, dv \, dw \\[4pt] &= 6 \int_0^1 \int_0^2 \left[ \frac{u^2}{2} + vu + wu\right]_0^1 \, dv \, dw \\[4pt] &= 6 \int_0^1 \int_0^2 \left(\frac{1}{2} + v + u\right) dv \, dw \\[4pt] &= 6 \int_0^1 \left[\frac{1}{2} v + \frac{v^2}{2} + wv \right]_0^2 dw\\[4pt] &= 6 \int_0^1 (3 + 2w)\, dw = 6\Big[3w + w^2\Big]_0^1 = 24. \end{align*}

Exercise \PageIndex{6}

Let D be the region in xyz-space defined by 1 \leq x \leq 2, \, 0 \leq xy \leq 2, and 0 \leq z \leq 1.

Evaluate \iiint_D (x^2 y + 3xyz) \, dx \, dy \, dz by using the transformation u = x, \, v = xy, and w = 3z.

Hint

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

Answer

\int_0^3 \int_0^2 \int_1^2 \left(\frac{v}{3} + \frac{vw}{3u}\right) du \, dv \, dw = 2 + \ln 8 \nonumber

Key Concepts

  • A transformation T is a function that transforms a region G in one plane (space) into a region R. in another plane (space) by a change of variables.
  • A transformation T: G \rightarrow R defined as T(u,v) = (x,y) (or T(u,v,w) = (x,y,z))is said to be a one-to-one transformation if no two points map to the same image point.
  • If f is continuous on R, then \iint_R f(x,y) dA = \iint_S f(g(u,v), \, h(u,v)) \left|\frac{\partial(x,y)}{\partial (u,v)}\right| du \, dv. \nonumber
  • If F is continuous on R, then \begin{align*}\iiint_R F(x,y,z) \, dV &= \iiint_G F(g(u,v,w), \, h(u,v,w), \, k(u,v,w) \left|\frac{\partial(x,y,z)}{\partial (u,v,w)}\right| \,du \, dv \, dw \\[4pt] &= \iiint_G H(u,v,w) |J(u,v,w)| \, du \, dv \, dw. \end{align*}

[T] Lamé ovals (or superellipses) are plane curves of equations \left(\frac{x}{a}\right)^n + \left( \frac{y}{b}\right)^n = 1, where a, b, and n are positive real numbers.

a. Use a CAS to graph the regions R bounded by Lamé ovals for a = 1, \, b = 2, \, n = 4 and n = 6 respectively.

b. Find the transformations that map the region R bounded by the Lamé oval x^4 + y^4 = 1 also called a squircle and graphed in the following figure, into the unit disk.

A square of side length 2 with rounded corners.

c. Use a CAS to find an approximation of the area A (R) of the region R bounded by x^4 + y^4 = 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 \left(\frac{x}{a}\right)^n + \left( \frac{y}{b}\right)^n = 1 with \frac{a}{b} = \frac{9}{7} 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 yards.

[Hide Solution]

A(R) \simeq 83,999.2

Chapter Review Exercises

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

\int_a^b \int_c^d f(x,y) \, dy \, dx = \int_c^d \int_a^b f(x,y) \, dy \, dx \nonumber

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

[Hide solution]

True.

The integral \int_0^{2\pi} \int_0^1 \int_0^1 dz \, dr \, d\theta \nonumber represents the volume of a right cone.

The Jacobian of the transformation for x = u^2 - 2v, \, y = 3v - 2uv is given by -4u^2 + 6u + 4v.

[Hide Solution]

False.

Evaluate the following integrals.

\iint_R (5x^3y^2 - y^2) \, dA, \, R = \{(x,y)|0 \leq x \leq 2, \, 1 \leq y \leq 4\} \nonumber

\iint_D \frac{y}{3x^2 + 1} dA, \, D = \{(x,y) |0 \leq x \leq 1, \, -x \leq y \leq x\} \nonumber

[Hide Solution]

0

\iint_D \sin (x^2 + y^2) dA \nonumber where D is a disk of radius 2 centered at the origin \int_0^1 \int_0^1 xye^{x^2} dx \, dy \nonumber

[Hide Solution]

\frac{1}{4}

\int_{-1}^1 \int_0^z \int_0^{x-z} 6dy \, dx \, dz \nonumber

\iiint_R 3y \, dV, \nonumber where R = \{(x,y,z) |0 \leq x \leq 1, \, 0 \leq y \leq x, \, 0 \leq z \leq \sqrt{9 - y^2}\}

[Hide Solution]

1.475

\int_0^2 \int_0^{2\pi} \int_r^1 r \, dz \, d\theta \, dr \nonumber

\int_0^{2\pi} \int_0^{\pi/2} \int_1^3 \rho^2 \, \sin(\varphi) d\rho \, d\varphi, \, d\theta \nonumber

[Hide Solution]

\frac{52}{3} \pi

\int_0^1 \int_{-\sqrt{1-x^2}}^{\sqrt{1-x^2}} \int_{-\sqrt{1-x^2-y^2}}^{\sqrt{1-x^2-y^2}} dz \, dy \, sx \nonumber

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

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

[Hide Solution]

\frac{\pi}{16}

The volume of the solid that lies between the paraboloid z = 2x^2 + 2y^2 and the plane z = 8.

The volume of the solid bounded by the cylinder x^2 + y^2 = 16 and from z = 1 to z + x = 2.

[Hide Solution]

93.291

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.

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

[Hide Solution]

\left(\frac{8}{15}, \frac{8}{15}\right)

\rho(x,y) = (y + 1) \sqrt{x} in the region bounded by y = e^x, \, y = 0, and x = 1.

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

\left(0,0,\frac{8}{5}\right)

The volume an ice cream cone that is given by the solid above z = \sqrt{(x^2 + y^2)} and below z^2 + x^2 + y^2 = 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 400 \, lb/ft^3, find the amount of work required to build the mountain.

[Hide Solution]

1.452 \pi \times 10^{15} 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 400 \, lb/ft^3 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)

Layer Depth from center (km) Temperature ^oC
Rocky Crust 0 to 40 0
Upper Mantle 40 to 150 870
Mantle 400 to 650 870
Inner Mantel 650 to 2700 870
Molten Outer Core 2890 to 5150 4300
Inner Core 5150 to 6378 7200

Source: http://www.enchantedlearning.com/sub...h/Inside.shtml

[Hide Solution]

y = -1.238 \times 10^{-7} x^3 + 0.001196 x^2 - 3.666x + 7208; average temperature approximately 2800 ^oC

[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.

Layer Depth from center (km) Density (g/cm^3)
Inner Core 0 12.95
Outer Core 1228 11.05
Mantle 3488 5.00
Upper Mantle 6338 3.90
Crust 6378 2.55

Source: http://hyperphysics.phy-astr.gsu.edu...rthstruct.html

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 V_x = 2\pi A \bar{y}, and when you revolve around the y-axis the volume is given by V_y = 2\pi A \bar{x}, where A is the area of R. Consider the region bounded by x^2 + y^2 = 1 and above y = x + 1.

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

[Hide Solution]

\frac{\pi}{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 \times 2 determinant:

J(u,v) = \begin{vmatrix} \frac{\partial x}{\partial u} \frac{\partial y}{\partial u} \nonumber \\ \frac{\partial x}{\partial v} \frac{\partial y}{\partial v} \end{vmatrix}; \nonumber

the Jacobian J (u,v,w) in three variables is a 3 \times 3 determinant:

J(u,v,w) = \begin{vmatrix} \frac{\partial x}{\partial u} \frac{\partial y}{\partial u} \frac{\partial z}{\partial u} \nonumber \\ \frac{\partial x}{\partial v} \frac{\partial y}{\partial v} \frac{\partial z}{\partial v} \nonumber \\ \frac{\partial x}{\partial w} \frac{\partial y}{\partial w} \frac{\partial z}{\partial w}\end{vmatrix} \nonumber

one-to-one transformation
a transformation T : G \rightarrow R 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 GG in one plane into a region RR in another plane by a change of variables

This page titled 15.7: Change of Variables in Multiple Integrals is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Gilbert Strang & Edwin “Jed” Herman (OpenStax) via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?