Change of Variables and Jacobians

Welcome, students! Today’s lesson is all about a powerful tool in multivariable calculus: the change of variables formula and Jacobians. By the end of this lesson, you’ll understand how to transform complicated integrals into simpler ones by using new coordinate systems. We’ll explore how Jacobians help us account for area and volume distortions in these transformations. Get ready to unlock a key technique used in everything from physics to engineering to computer graphics! 🚀

What is a Change of Variables?

Let’s start with the big idea. Imagine you’re trying to solve a double integral over a complicated region—maybe it’s a wavy shape or an oddly rotated rectangle. Integrating directly in $x$ and $y$ coordinates might be tough. What if we could stretch, rotate, or skew our coordinate system to make the region simpler? That’s exactly what the change of variables technique does!

A change of variables means we introduce new variables, often named $u$ and $v$ (or $r$ and $\theta$ in polar coordinates). These new variables are related to the old ones by some transformation. Once we transform both the region and the function we’re integrating, the integral often becomes much simpler.

Mathematically, the general idea looks like this for two variables:

x = x(u, v), \quad y = y(u, v)

When we switch from $(x, y)$ to $(u, v)$, we need to adjust for how areas (or volumes in 3D) change. This is where the Jacobian comes in.

Introducing the Jacobian

The Jacobian is a special determinant that measures how a transformation stretches or compresses areas (or volumes). Think of it as a scaling factor.

Let’s say we have a transformation from $(x, y)$ to $(u, v)$. The Jacobian $J$ is defined as the determinant of the matrix of partial derivatives:

$J = \begin{vmatrix}$

\frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\

\frac{\partial y}{\partial u} & \frac{\partial y}{\partial v}

$\end{vmatrix}$

This determinant tells us how much an infinitesimal area in the $u$-$v$ plane gets scaled when we map it into the $x$-$y$ plane. If $|J| > 1$, the area gets stretched. If $0 < |J| < 1$, the area gets shrunk.

In 3D, the Jacobian is a $3 \times 3$ determinant:

$J = \begin{vmatrix}$

\frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} & \frac{\partial x}{\partial w} \\

\frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} & \frac{\partial y}{\partial w} \\

\frac{\partial z}{\partial u} & \frac{\partial z}{\partial v} & \frac{\partial z}{\partial w}

$\end{vmatrix}$

We use $|J|$ to adjust the integral when we switch variables. This ensures that the total “amount” being integrated—like mass, probability, or energy—stays the same.

The Change of Variables Formula

Here’s the grand finale: the change of variables formula for double integrals. If $f(x, y)$ is a function we want to integrate over a region $R$ in the $x$-$y$ plane, and we use a transformation to new variables $(u, v)$, then the integral transforms as follows:

$\iint_{R}$ f(x, y) \, dx\, dy = $\iint_{S}$ f(x(u, v), y(u, v)) \, |J| \, du\, dv

Here, $S$ is the corresponding region in the $u$-$v$ plane. Notice that we multiply by $|J|$ to account for how the area changes under the transformation.

In 3D, the formula is similar:

$\iiint_{R}$ f(x, y, z) \, dx\, dy\, dz = $\iiint_{S}$ f(x(u, v, w), y(u, v, w), z(u, v, w)) \, |J| \, du\, dv\, dw

This formula is a game-changer. It allows us to solve integrals that would otherwise be extremely difficult. Let’s dive into some examples to see this in action!

Example: Polar Coordinates Transformation

One of the most common examples of a change of variables is switching from Cartesian coordinates $(x, y)$ to polar coordinates $(r, \theta)$. You’ve probably seen this one before, but let’s break it down step by step.

In polar coordinates:

x = r $\cos($$\theta)$, \quad y = r $\sin($$\theta)$

Let’s compute the Jacobian for this transformation. We’ll find the partial derivatives:

\frac{\partial x}{\partial r} = $\cos($$\theta)$, \quad \frac{\partial x}{\partial \theta} = -r $\sin($$\theta)$

\frac{\partial y}{\partial r} = $\sin($$\theta)$, \quad \frac{\partial y}{\partial \theta} = r $\cos($$\theta)$

Now let’s form the Jacobian matrix:

$J = \begin{vmatrix}$

$\cos(\theta) & -r \sin(\theta) \\$

$\sin(\theta) & r \cos(\theta)$

$\end{vmatrix}$

We compute the determinant:

J = $\cos($$\theta)$ $\cdot$ (r $\cos($$\theta)$) - (-r $\sin($$\theta)$) $\cdot$ $\sin($$\theta)$ = r ($\cos^2$($\theta)$ + $\sin^2$($\theta)$) = r

This is a super important result: the Jacobian for the polar transformation is $r$. That means when we switch from $(x, y)$ to $(r, \theta)$, we need to multiply by $r$ to account for the area scaling.

So the integral in polar coordinates looks like:

$\iint_{R}$ f(x, y) \, dx\, dy = $\int_{\theta_1}$^{$\theta_2$} $\int_{r_1}$^{r_2} f(r $\cos($$\theta)$, r $\sin($$\theta)$) \, r \, dr\, d$\theta$

This explains why in polar integrals, you’ll always see that extra $r$ factor. It’s the Jacobian!

Real-World Application: Calculating Area of a Circle

Let’s use polar coordinates to find the area of a circle of radius $a$. We know the area should be $\pi a^2$, but let’s prove it using integration.

We want to integrate the constant function $f(x, y) = 1$ over the circle $x^2 + y^2 \leq a^2$. In polar coordinates, the region is $0 \leq r \leq a$, $0 \leq \theta \leq 2\pi$. The integral becomes:

$\int_{0}$^{$2\pi$} $\int_{0}$^{a} $1 \cdot$ r \, dr\, d$\theta$

First, integrate with respect to $r$:

$\int_{0}$^{a} r \, dr = $\frac{r^2}{2}$ $\Big|_{0}$^{a} = $\frac{a^2}{2}$

Now integrate with respect to $\theta$:

$\int_{0}$^{$2\pi$} $\frac{a^2}{2}$ \, d$\theta$ = $\frac{a^2}{2}$ $\cdot 2$$\pi$ = $\pi$ a^2

Voilà! We’ve found the area of the circle using the change of variables technique.

Example: Transformation to Elliptical Coordinates

Let’s tackle a more complex example: transforming an integral over an ellipse. Suppose we want to integrate over the region defined by:

$\frac{x^2}{a^2} + \frac{y^2}{b^2} \leq 1$

This is an ellipse with semi-major axis $a$ and semi-minor axis $b$. Integrating directly in $x$ and $y$ might be tricky. Let’s use a change of variables to transform this ellipse into a circle.

We define a new set of variables $(u, v)$ by:

x = a u, \quad y = b v

In these new coordinates, the ellipse becomes the unit circle:

$\frac{(a u)^2}{a^2}$ + $\frac{(b v)^2}{b^2}$ = u^2 + v^$2 \leq 1$

Now let’s find the Jacobian. We compute the partial derivatives:

\frac{\partial x}{\partial u} = a, \quad \frac{\partial x}{\partial v} = 0

\frac{\partial y}{\partial u} = 0, \quad \frac{\partial y}{\partial v} = b

So the Jacobian matrix is:

$J = \begin{vmatrix}$

a & 0 \\

0 & b

$\end{vmatrix}$

The determinant is:

J = a $\cdot$ b = ab

This means that areas in the $u$-$v$ plane are scaled by a factor of $ab$ when we go back to the $x$-$y$ plane. So the integral transforms like this:

$\iint_{R}$ f(x, y) \, dx\, dy = $\iint_{S}$ f(a u, b v) \, ab \, du\, dv

This makes it much easier to integrate over the ellipse, because we’ve turned it into a simple unit circle in the $u$-$v$ plane.

Real-World Application: Elliptical Orbits in Physics

In physics, elliptical orbits are common—planets orbit the sun in ellipses. Calculating areas or integrals over elliptical regions is often necessary. The change of variables technique we just used is the same idea behind Kepler’s Second Law, which says that a planet sweeps out equal areas in equal times. The Jacobian helps us understand how these areas change as we switch between coordinate systems.

Example: 3D Change of Variables—Spherical Coordinates

Now let’s move up to three dimensions. One of the most powerful transformations in 3D is the switch to spherical coordinates. This is super useful for problems with spherical symmetry, like electric fields or gravitational fields around a planet.

Spherical coordinates $(\rho, \phi, \theta)$ are defined by:

x = $\rho$ $\sin($$\phi)$ $\cos($$\theta)$, \quad y = $\rho$ $\sin($$\phi)$ $\sin($$\theta)$, \quad z = $\rho$ $\cos($$\phi)$

Here, $\rho$ is the distance from the origin, $\phi$ is the angle from the positive $z$-axis (often called the polar angle), and $\theta$ is the angle in the $xy$-plane (the azimuthal angle).

Let’s find the Jacobian. We compute the partial derivatives (this is a bit more involved, but hang in there!):

\frac{\partial x}{\partial \rho} = $\sin($$\phi)$ $\cos($$\theta)$, \quad \frac{\partial x}{\partial \phi} = $\rho$ $\cos($$\phi)$ $\cos($$\theta)$, \quad \frac{\partial x}{\partial \theta} = -$\rho$ $\sin($$\phi)$ $\sin($$\theta)$

\frac{\partial y}{\partial \rho} = $\sin($$\phi)$ $\sin($$\theta)$, \quad \frac{\partial y}{\partial \phi} = $\rho$ $\cos($$\phi)$ $\sin($$\theta)$, \quad \frac{\partial y}{\partial \theta} = $\rho$ $\sin($$\phi)$ $\cos($$\theta)$

\frac{\partial z}{\partial \rho} = $\cos($$\phi)$, \quad \frac{\partial z}{\partial \phi} = -$\rho$ $\sin($$\phi)$, \quad \frac{\partial z}{\partial \theta} = 0

Now we form the $3 \times 3$ Jacobian matrix:

$J = \begin{vmatrix}$

$\sin($$\phi)$$\cos($$\theta)$ & $\rho$ $\cos($$\phi)$$\cos($$\theta)$ & -$\rho$ $\sin($$\phi)$$\sin($$\theta)$ \\

$\sin($$\phi)$$\sin($$\theta)$ & $\rho$ $\cos($$\phi)$$\sin($$\theta)$ & $\rho$ $\sin($$\phi)$$\cos($$\theta)$ \\

$\cos(\phi) & -\rho \sin(\phi) & 0$

$\end{vmatrix}$

After a bit of algebra (and maybe some patience!), the determinant simplifies to:

$J = \rho^2 \sin(\phi)$

So the Jacobian for spherical coordinates is $\rho^2 \sin(\phi)$. This means that when we change variables from $(x, y, z)$ to $(\rho, \phi, \theta)$, we multiply by $\rho^2 \sin(\phi)$ to account for the volume scaling.

The integral in spherical coordinates looks like:

$\iiint_{R}$ f(x, y, z) \, dx\, dy\, dz = $\int_{\theta_1}$^{$\theta_2$} $\int_{\phi_1}$^{$\phi_2$} $\int_{\rho_1}$^{$\rho_2$} f($\rho$, $\phi$, $\theta)$ \, $\rho^2$ $\sin($$\phi)$ \, d$\rho$\, d$\phi$\, d$\theta$

Real-World Application: Electric Field of a Sphere

In physics, spherical coordinates are perfect for problems involving spheres. For example, suppose we want to find the total electric charge inside a uniformly charged sphere. The charge density might be constant, say $\rho_0$. We want to integrate this density over the volume of the sphere.

In spherical coordinates, the sphere of radius $R$ is described by $0 \leq \rho \leq R$, $0 \leq \phi \leq \pi$, $0 \leq \theta \leq 2\pi$. So the total charge $Q$ is:

Q = $\iiint_{R}$ $\rho_0$ \, dx\, dy\, dz = $\int_{0}$^{$2\pi$} $\int_{0}$^{$\pi$} $\int_{0}$^{R} $\rho_0$ \, $\rho^2$ $\sin($$\phi)$ \, d$\rho$\, d$\phi$\, d$\theta$

First, integrate with respect to $\rho$:

$\int_{0}^{R} \rho^2 \, d\rho = \frac{R^3}{3}$

Next, integrate with respect to $\phi$:

$\int_{0}^{\pi} \sin(\phi) \, d\phi = 2$

Finally, integrate with respect to $\theta$:

$\int_{0}^{2\pi} \, d\theta = 2\pi$

Putting it all together:

Q = $\rho_0$ $\cdot$ $\frac{R^3}{3}$ $\cdot 2$ $\cdot 2$$\pi$ = $\frac{4\pi \rho_0 R^3}{3}$

This is the total charge inside the sphere. Without spherical coordinates, this integral would have been much more complicated.

A Quick Note on Inverse Transformations

It’s important to remember that when we do a change of variables, we need to make sure the transformation is invertible. That means we can go back from $(u, v)$ to $(x, y)$ uniquely. If the transformation isn’t one-to-one, we might run into trouble. For example, polar coordinates aren’t one-to-one at the origin, so we have to be careful when dealing with integrals there.

Conclusion

Congratulations, students! You’ve just unlocked a powerful tool in multivariable calculus: the change of variables formula and the Jacobian. We explored how transformations can simplify integrals by switching to new coordinate systems, and we saw how the Jacobian helps us adjust for area and volume changes. From polar coordinates to elliptical transformations to spherical coordinates, these techniques have wide-ranging applications in math, physics, and engineering. Keep practicing, and you’ll soon be a master of multivariable integration! 🎯

Study Notes

The change of variables formula for double integrals:

$\iint_{R}$ f(x, y) \, dx\, dy = $\iint_{S}$ f(x(u, v), y(u, v)) \, |J| \, du\, dv

The Jacobian for a transformation $(x, y) = (x(u, v), y(u, v))$ is:

$ J = \begin{vmatrix}$

\frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\

\frac{\partial y}{\partial u} & \frac{\partial y}{\partial v}

$ \end{vmatrix}$

For polar coordinates $(r, \theta)$:

x = r $\cos($$\theta)$, \quad y = r $\sin($$\theta)$

Jacobian: $J = r$
Integral formula:

$\iint_{R}$ f(x, y) \, dx\, dy = $\int_{\theta_1}$^{$\theta_2$} $\int_{r_1}$^{r_2} f(r $\cos($$\theta)$, r $\sin($$\theta)$) \, r \, dr\, d$\theta$

For elliptical transformations $(u, v)$:

x = a u, \quad y = b v

Jacobian: $J = ab$
Integral formula:

$\iint_{R}$ f(x, y) \, dx\, dy = $\iint_{S}$ f(a u, b v) \, ab \, du\, dv

For spherical coordinates $(\rho, \phi, \theta)$:

x = $\rho$ $\sin($$\phi)$ $\cos($$\theta)$, \quad y = $\rho$ $\sin($$\phi)$ $\sin($$\theta)$, \quad z = $\rho$ $\cos($$\phi)$

Jacobian: $J = \rho^2 \sin(\phi)$
Integral formula:

$\iiint_{R}$ f(x, y, z) \, dx\, dy\, dz = \