15.7: Change of Variables in Multiple Integrals
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{\!\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
Recall from Substitution Rule the method of integration by substitution. When evaluating an integral such as
\[\int_2^3 x(x^2  4)^5 dx,\]
we substitute \(u = g(x) = x^2  4\). Then \(du = 2x \, dx\) or \(x \, dx = \frac{1}{2} du\) and the limits change to \(u = g(2) = 2^2  4 = 0\) and \(u = g(3) = 9  4 = 5\). Thus the integral becomes
\[\int_0^5 \frac{1}{2}u^5 du\]
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,
\[\int_a^b f(x) dx = \int_c^d f(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,\theta) = r \, \cos \, \theta\)
 \( y = g(r, \theta) = r \, \sin \, \theta\), and
 \(dA = dx \, dy = r \, dr \, d\theta\).
Then we get
\[\iint_R f(x,y) dA = \iint_S (r \, \cos \, \theta, \, r \, \sin \, \theta)r \, dr \, d\theta\]
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.
Figure \(\PageIndex{1}\): The transformation of a region \(G\) in the \(uv\)plane into a region \(R\) in the \(xy\)plane.
Definition: onetoone transformation
A transformation \(T: \, G \rightarrow R\), defined as \(T(u,v) = (x,y)\), is said to be a onetoone transformation if no two points map to the same image point.
To show that \(T\) s a onetoone 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 onetoone 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 onetoone 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.
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 onetoone 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),\]
\[(x_1,y_1) = (x_1,y_1),\]
\[(r_1 \cos \, \theta_1, r_1 \sin \, \theta_1) = (r_2 \cos \, \theta_2, r_2 \sin \, \theta_2),\]
\[r_1 \cos \, \theta_1 = r_2 \cos \, \theta_2, \, r_1 \sin \, \theta_1 = r_2 \sin \, \theta_2.\]
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}\]
\[\frac{\cos \, \theta_1}{\sin \, \theta_1} = \frac{\cos \, \theta_2}{\sin \, \theta_2}\]
\[\tan \, \theta_1 = \tan \, \theta_2\]
\[\theta_1 = \theta_2\]
since the tangent function is oneone 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 onetoone 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\).
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 onetoone 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{3xy}{3}\) and \(v = \frac{y}{3}\)
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.\]
Note that the Jacobian is frequently denoted simply by
\[J(u,v) = \frac{\partial (x,y)}{\partial (u,v)}.\]
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} .\]
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} \nonumber \\ \dfrac{\partial y}{\partial r} & \dfrac{\partial y}{\partial \theta} \end{vmatrix} = \begin{vmatrix} \cos \theta & r\sin \theta \nonumber \\ \sin \theta & r \cos \theta \end{vmatrix} = r \, \cos^2\theta + r \, \sin^2\theta = r ( \cos^2\theta + \sin^2\theta) = r.\]
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} \nonumber \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \begin{vmatrix} 2u & v \nonumber \\ 2v & u \end{vmatrix} = 2u^2 + 2v^2.\]
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\]
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.\]
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_{ijh}) \Delta A.\]
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)\).
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.\]
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.\]
Change of Variables for Double Integrals
Let \(T(u,v) = (x,y)\) where \(x = g(u,v)\) and \(y = h(u,v)\) be a onetoone \(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.\]
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\]
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{2xx^2}} \sqrt{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.
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}\)).
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{2xx^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.\]
Exercise \(\PageIndex{3}\)
Considering the integral \(\int_0^1 \int_0^{\sqrt{1x^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\]
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,\] 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.
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{3uv}{2}\). Moreover, if we simply subtract the second equation from the first, we get \(u  v = (x  y)  (x  3y) = 2y\) and \(y = \frac{uv}{2}\).
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)\] 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} \nonumber \\ \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}\]
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.\]
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,\] 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,\] 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_{4}^4 \int_{2}^2 \frac{4}{u^2} \left(\frac{1}{2}\right) du \, dv.\]
ProblemSolving Strategy: Change of Variables
 Sketch the region given by the problem in the \(xy\)plane and then write the equations of the curves that form the boundary.
 Depending on the region or the integrand, choose the transformations \(x = g(u,v)\) and \(y = h(u,v)\).
 Determine the new limits of integration in the \(uv\)plane.
 Find the Jacobian \(J (u,v)\).
 In the integrand, replace the variables to obtain the new integrand.
 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^2y^2} dA,\] 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}\)).
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{vu}{2}\) and hence the transformation to use is \(T(u,v) = \left(\frac{u+v}{2}, \, \frac{vu}{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\}.\]
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} \nonumber \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{vmatrix} = \begin{vmatrix} 1/2 & 1/2 \nonumber \\ 1/2 & 1/2 \end{vmatrix} = \frac{1}{2}.\]
Therefore, by using the transformation \(T\), the integral changes to
\[\iint_R (x  y)e^{x^2y^2} dA = \frac{1}{2} \int_1^3 \int_{1/v}^{1/v} ue^{uv} du \, dv.\]
Doing the evaluation, we have
\[\frac{1}{2} \int_1^3 \int_{1/v}^{1/v} ue^{uv} du \, dv = \frac{4}{3e} \approx 0.490.\]
Exercise \(\PageIndex{5}\)
Using the substitutions \(x = v\) and \(y = \sqrt{u + v}\), evaluate the integral \[\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 onetoone \(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)\).
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).\]
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} \nonumber \\ \dfrac{\partial x}{\partial v} & \dfrac{\partial y}{\partial v} & \dfrac{\partial z}{\partial v} \nonumber \\ \dfrac{\partial x}{\partial w} & \dfrac{\partial y}{\partial w} & \dfrac{\partial z}{\partial w} \end{vmatrix}. \]
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} \nonumber \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} & \dfrac{\partial y}{\partial w} \nonumber \\ \dfrac{\partial z}{\partial u} & \dfrac{\partial z}{\partial v} & \dfrac{\partial z}{\partial w} \end{vmatrix}. \]
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 onetoone \(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 twodimensional case, if \(F\) is continuous on \(D\), 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 \\ &= \iiint_G H(u,v,w)  J (u,v,w)  du \, dv \, dq. \end{align}\]
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
 cylindrical and
 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} \nonumber \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \theta} & \frac{\partial y}{\partial z} \nonumber \\ \frac{\partial z}{\partial r} & \frac{\partial z}{\partial \theta} & \frac{\partial z}{\partial z} \end{vmatrix} \] \[ \begin{vmatrix} \cos \theta & r\, \sin \theta & 0 \nonumber \\ \sin \theta & r\, \cos \theta & 0 \nonumber \\ 0 & 0 & 1 \end{vmatrix} = r \, \cos^2 \theta + r \, \sin^2 \theta = r.\] 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.\]
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 \varphi\), 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} \nonumber \\ \frac{\partial y}{\partial \rho} & \frac{\partial y}{\partial \theta} & \frac{\partial y}{\partial \varphi} \nonumber \\ \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 \nonumber \\ \sin \varphi \, \sin \theta & \rho \, \sin \varphi \, \cos \theta & \rho \, \cos \varphi \, \sin \theta \nonumber \\ \cos \theta & 0 & \rho \, \sin \varphi \end{vmatrix}. \]
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 \nonumber \\ \rho \, sin \, \varphi \, \sin \theta & \rho \, \cos \varphi \, \sin \theta \end{vmatrix}  \rho \, \sin \varphi \begin{vmatrix} \sin \varphi \, \cos \theta &  \rho \, \sin \varphi \, \sin \theta \nonumber \\ \sin \varphi \, \sin \theta & \rho \, \sin \varphi \, \cos \theta \end{vmatrix} \\ &=\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 \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.\)
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.\]
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\]
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\).
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} \nonumber \\ \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} & \dfrac{\partial y}{\partial w} \nonumber \\ \dfrac{\partial z}{\partial u} & \dfrac{\partial z}{\partial v} & \dfrac{\partial z}{\partial w} \end{vmatrix} = \begin{vmatrix} 1 & 1 & 0 \nonumber \\ 0 & 2 & 0 \nonumber \\ 0 & 0 & 3 \end{vmatrix} = 6.\]
The function to be integrated becomes
\[f(x,y,z) = x + \frac{z}{3} = u + v + \frac{3w}{3} = u + v + w.\]
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 = \int_0^1 \int_0^2 \int_0^1 (u + v + w) 6du \, dv \, dw \\&= 6 \int_0^1 \int_0^2 \int_0^1 (u + v + w) \, du \, dv \, dw = 6 \int_0^1 \int_0^2 \left[ \frac{u^2}{2} + vu + wu\right]_0^1 \, dv \, dw \\ &= 6 \int_0^1 \int_0^2 \left(\frac{1}{2} + v + u\right) dv \, dw = 6 \int_0^1 \left[\frac{1}{2} v + \frac{v^2}{2} + wv \right]_0^2 dw\\ &= 6 \int_0^1 (3 + 2w)\, dw = 6[3w + w^2]_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\]
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 onetoone 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.\]
 If \(F\) is continuous on \(R\), then \[\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.\] \[= \iint_G H(u,v,w) J(u,v,w) \, du \, dv \, dw.\]
[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.
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\]
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\] 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\}\]
\[\iint_D \frac{y}{3x^2 + 1} dA, \, D = \{(x,y) 0 \leq x \leq 1, \, x \leq y \leq x\}\]
[Hide Solution]
\(0\)
\[\iint_D \sin (x^2 + y^2) dA\] where \(D\) is a disk of radius \(2\) centered at the origin \[\int_0^1 \int_0^1 xye^{x^2} dx \, dy\]
[Hide Solution]
\(\frac{1}{4}\)
\[\int_{1}^1 \int_0^z \int_0^{xz} 6dy \, dx \, dz\]
\[\iiint_R 3y \, dV,\] 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\]
\[\int_0^{2\pi} \int_0^{\pi/2} \int_1^3 \rho^2 \, \sin(\varphi) d\rho \, d\varphi, \, d\theta\]
[Hide Solution]
\(\frac{52}{3} \pi\)
\[\int_0^1 \int_{\sqrt{1x^2}}^{\sqrt{1x^2}} \int_{\sqrt{1x^2y^2}}^{\sqrt{1x^2y^2}} dz \, dy \, sx\]
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} \) ftlb
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 bestfit 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.phyastr.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};\]
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}\]
 onetoone transformation
 a transformation \(T : G \rightarrow R\) defined as \(T(u,v) = (x,y)\) is said to be onetoone 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
Contributors
Gilbert Strang (MIT) and Edwin “Jed” Herman (Harvey Mudd) with many contributing authors. This content by OpenStax is licensed with a CCBY 3/0 license. Download for free at http://cnx.org/contents/fd53eae1fa2...49835c3c@5.191.
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 onetoone \(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}\)).
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.\]
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.\]
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.\]
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.\]
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 \nonumber \\ \frac{\partial x}{\partial u} & \frac{\partial y}{\partial u} & 0 \nonumber \\ \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} \nonumber \\ \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) \]
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} \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).\]