Jacobian Determinants — Use Jacobian determinants for coordinate transformations (intro) - Multivariable Calculus

SECTION 1

Historical Context & Motivation

Some integrals are incredibly difficult — or even impossible — to compute using standard rectangular (x, y) coordinates. Think about trying to find the area of a circle by slicing it into tiny vertical strips: the math gets messy fast. Mathematicians realized centuries ago that choosing the right coordinate system can turn a nightmare integral into something elegant. But switching coordinate systems introduces a hidden problem: the little patches of area or volume you're summing up can stretch or shrink during the transformation. The Jacobian determinant is the correction factor that accounts for this stretching.

1740s • Euler's Polar Coordinates

Leonhard Euler systematically uses polar coordinates to simplify problems involving circles and rotations, noticing that the area element changes from dx dy to r dr dθ.

1841 • Jacobi Formalizes the Determinant

Carl Gustav Jacob Jacobi publishes a comprehensive theory of determinants and introduces the Jacobian matrix. His work provides a general formula for how area and volume elements transform between any two coordinate systems.

1850s • Riemann and Integration Theory

Bernhard Riemann develops rigorous foundations for integration in multiple dimensions, and the Jacobian becomes a standard tool for changing variables in double and triple integrals.

Modern • Everyday Applications

Today the Jacobian is used everywhere — from computing gravitational fields in astrophysics to transforming images in computer graphics and mapping satellite data onto curved surfaces.

So the central question this lesson addresses is: when you switch from one set of coordinates to another, how do you correctly adjust the integral so that areas and volumes are preserved? The answer lies in the Jacobian determinant.

SECTION 2

Core Principles & Definitions

Before diving into calculations, you need a clear mental picture of what is happening when you change coordinates. A coordinate transformation is a function that takes points described in one system (like polar r, θ) and maps them to another system (like Cartesian x, y). The Jacobian determinant measures how much a tiny patch of area in the new coordinates gets stretched or compressed when mapped back to the original coordinates.

Coordinate Transformation

A set of equations that converts coordinates from one system to another, such as x = r cos θ and y = r sin θ for polar-to-Cartesian conversion.

Jacobian Matrix

A matrix of all first-order partial derivatives of the transformation. It captures how each output variable changes with respect to each input variable.

Jacobian Determinant

The determinant of the Jacobian matrix. Its absolute value gives the local scaling factor — how much an infinitesimal area (or volume) element is stretched or compressed.

Area / Volume Element

In Cartesian coordinates, the area element is dA = dx dy. After a transformation, it becomes dA = |J| du dv, where |J| is the absolute value of the Jacobian determinant.

✦

KEY TAKEAWAY

Imagine you have a rubber sheet with a grid drawn on it. If you stretch part of the sheet, the grid squares in that region get bigger. If you compress another part, those squares shrink. The Jacobian determinant tells you exactly how much each little square stretches or shrinks at every point. Without this correction factor, your integral would give the wrong answer — like measuring the area of a stretched map and claiming it's the actual land area.

SECTION 3

Visual Explanation — How Transformations Warp Area

The diagram below shows a concrete example of why the Jacobian matters. On the left you see a uniform grid in polar-like coordinates (u, v). On the right you see where those grid squares end up after a transformation into Cartesian (x, y) coordinates. Notice how the patches change size and shape — some are large, some are small. The Jacobian determinant at each point tells you the ratio of the new patch area to the original patch area.

Left: a uniform grid of equal patches in the (u, v) plane. Right: the same patches after a coordinate transformation T(u, v). Patch A₁' (violet) and A₂' (cyan) now have different sizes. The Jacobian determinant gives the local scaling factor at each point.

This is exactly what happens when you switch from Cartesian to polar coordinates. Near the origin (where r is small), the polar grid squares are tiny, so the Jacobian (which equals r) is small. Far from the origin, the patches are larger, so the Jacobian is larger. This is why the polar area element is r dr dθ instead of just dr dθ — the factor of r is the Jacobian determinant for the polar transformation.

SECTION 4

Mathematical Framework

Suppose you have a transformation that converts coordinates (u, v) into (x, y) via two functions: x = g(u, v) and y = h(u, v). The Jacobian matrix collects all the partial derivatives of these functions into a 2×2 grid. The determinant of that matrix is the Jacobian determinant.

JACOBIAN MATRIX (2D)

J = [ ∂x/∂u ∂x/∂v ] [ ∂y/∂u ∂y/∂v ]

Each entry is a partial derivative — it measures how one output variable changes when you nudge one input variable while holding the other fixed.

JACOBIAN DETERMINANT (2D)

det(J) = (∂x/∂u)(∂y/∂v) − (∂x/∂v)(∂y/∂u)

This is the standard 2×2 determinant formula: multiply the main diagonal entries and subtract the product of the off-diagonal entries.

CHANGE-OF-VARIABLES FORMULA

∬_R f(x, y) dx dy = ∬_S f(g(u,v), h(u,v)) · |det(J)| du dv

R is the region in (x, y) coordinates, S is the corresponding region in (u, v) coordinates, and |det(J)| is the absolute value of the Jacobian determinant. You always take the absolute value because area cannot be negative.

ℹ

Why Absolute Value?

The determinant can be positive or negative depending on whether the transformation preserves or reverses orientation (think of flipping a shape in a mirror). Since we're measuring area, which is always positive, we use |det(J)|.

SECTION 5

Detailed Breakdown — The Polar Coordinate Jacobian

The most common coordinate transformation you'll encounter is the switch from Cartesian to polar coordinates. Let's derive the Jacobian for this transformation step by step, which will also show you why the area element in polar coordinates is r dr dθ.

The polar-to-Cartesian transformation is given by x = r cos θ and y = r sin θ. To build the Jacobian matrix, we compute four partial derivatives: ∂x/∂r = cos θ, ∂x/∂θ = −r sin θ, ∂y/∂r = sin θ, and ∂y/∂θ = r cos θ.

POLAR JACOBIAN MATRIX

J = [ cos θ −r sin θ ] [ sin θ r cos θ ]

Row 1 contains partial derivatives of x; row 2 contains partial derivatives of y. Columns correspond to r and θ respectively.

POLAR JACOBIAN DETERMINANT

det(J) = (cos θ)(r cos θ) − (−r sin θ)(sin θ) = r cos²θ + r sin²θ = r

Using the Pythagorean identity cos²θ + sin²θ = 1, the Jacobian simplifies beautifully to just r. This is why dA = r dr dθ in polar coordinates.

A polar grid centered at the origin. The pink patch near the center (small r) has a much smaller area than the cyan patch far from the center (large r), even though both subtend the same Δr and Δθ. The factor of r in the area element accounts for this difference.

As you can see in the diagram, when r is small the wedge-shaped patch is tiny, and when r is large the patch is much bigger. The Jacobian determinant r captures this relationship perfectly. For a region described in polar coordinates, the double integral becomes ∬ f(r cos θ, r sin θ) · r dr dθ. That extra factor of r is not something you just memorize — it's the Jacobian, and it comes directly from the determinant calculation above.

SECTION 6

Worked Example — Converting to Polar Coordinates

Let's compute the integral ∬ $e^{−x² − y²}$ dA over the disk x² + y² ≤ 4. This integral is nearly impossible in Cartesian coordinates, but polar coordinates make it straightforward.

Evaluate ∬ $e^{−x² − y²}$ dA over x² + y² ≤ 4

1
Step 1 — Identify the Transformation

Use the polar coordinate transformation: x = r cos θ, y = r sin θ. The region x² + y² ≤ 4 becomes 0 ≤ r ≤ 2, 0 ≤ θ ≤ 2π. Notice how the circular region has simple bounds in polar coordinates.

Bounds: r ∈ [0, 2], θ ∈ [0, 2π]
2
Step 2 — Substitute into the Integrand

Replace x² + y² with r² (since x² + y² = r²cos²θ + r²sin²θ = r²). The integrand $e^{−x² − y²}$ becomes $e^{−r²}$ .

Integrand: $e^{−r²}$
3
Step 3 — Apply the Jacobian

The Jacobian for polar coordinates is |det(J)| = r. Replace dA = dx dy with r dr dθ. The integral becomes ∫₀²π ∫₀² $e^{−r²}$ · r dr dθ.

∫₀²π ∫₀² r $e^{−r²}$ dr dθ
4
Step 4 — Evaluate the Inner Integral (with respect to r)

Use substitution: let u = −r², so du = −2r dr, which gives r dr = −du/2. When r = 0, u = 0; when r = 2, u = −4. The inner integral becomes ∫₀⁻⁴ eᵘ (−du/2) = (1/2)∫₋₄⁰ eᵘ du = (1/2)[eᵘ]₋₄⁰ = (1/2)(1 − e⁻⁴).

Inner integral = (1/2)(1 − e⁻⁴)
5
Step 5 — Evaluate the Outer Integral (with respect to θ)

The inner result has no θ dependence, so ∫₀²π (1/2)(1 − e⁻⁴) dθ = (1/2)(1 − e⁻⁴) × 2π = π(1 − e⁻⁴).

∬ $e^{−x² − y²}$ dA = π(1 − e⁻⁴) ≈ 3.084

💡

Notice the Power of the Jacobian

Without the Jacobian factor r, we would have integrated

e^{−r²}

dr, which has no elementary antiderivative. The presence of r dr creates the perfect substitution u = −r², making the integral solvable. The Jacobian didn't just correct the area — it also simplified the math!

SECTION 7

When to Use (and Not Use) the Jacobian

Choosing the right coordinate system is an art as much as a science. The table below compares situations where a coordinate transformation with a Jacobian helps versus situations where Cartesian coordinates may be simpler.

Scenario	Best Coordinate System	Jacobian Factor
Circular or ring-shaped region	Polar (r, θ)	\|J\| = r
Spherical region or radial symmetry	Spherical (ρ, θ, φ)	\|J\| = ρ² sin φ
Cylindrical pipe or tube	Cylindrical (r, θ, z)	\|J\| = r
Rectangular region with straight sides	Cartesian (x, y)	\|J\| = 1 (no change needed)
Parallelogram or skewed grid	Custom linear (u, v)	\|J\| = \|ad − bc\| for x = au + bv, y = cu + dv

Common coordinate systems and their Jacobians

✦

KEY TAKEAWAY

Think of choosing coordinates like choosing a camera angle for a photograph. A straight-on shot (Cartesian) works great for a rectangular building, but if you're photographing a spiral staircase (circular region), you'd want to orbit around it (polar). The Jacobian is like the lens correction that ensures the photograph (your integral) accurately represents the real scene despite the change in perspective.

SECTION 8

Connection to Advanced Theory — 3D and Beyond

Everything we've covered for double integrals extends naturally to triple integrals and even higher dimensions. In three dimensions the Jacobian matrix becomes 3×3, and its determinant gives the local volume scaling factor. The table below compares the 2D and 3D cases.

Feature	2D (Double Integrals)	3D (Triple Integrals)
Jacobian matrix size	2 × 2	3 × 3
Determinant measures	Local area scaling	Local volume scaling
Polar / Spherical Jacobian	\|J\| = r	\|J\| = ρ² sin φ
Area / Volume element	dA = \|J\| du dv	dV = \|J\| du dv dw
Number of partial derivatives	4	9

Scaling the Jacobian from 2D to 3D

In more advanced courses, you'll see the Jacobian appear in differential geometry, where it describes how manifolds (curved surfaces and spaces) are mapped and measured. In physics, it is essential for general relativity and fluid dynamics. The core idea, however, remains the same: the Jacobian determinant is the correction factor that preserves measurement when you change your description of space.

💡

Looking Ahead

Once you are comfortable computing 2×2 Jacobians, the jump to 3×3 (and the spherical/cylindrical Jacobians) is mostly bookkeeping — the same idea of partial derivatives in a matrix, just with more entries. Mastering the 2D case now gives you a strong foundation.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

In your own words, explain why switching from Cartesian to polar coordinates requires multiplying the integrand by r. What would happen to the value of the integral if you forgot this factor?

Show answer

When you switch to polar coordinates, the tiny area elements near the origin are much smaller than those far from the origin, even if they span the same Δr and Δθ. The factor r is the Jacobian determinant, and it corrects for this varying patch size. If you forgot the factor r, you would underweight contributions far from the origin and overweight those near the origin, giving an incorrect (too small) result. Essentially, you'd be treating all polar patches as if they had the same area, which they don't.

Question 2: BASIC CALCULATION

Find the Jacobian determinant for the transformation x = 3u + v, y = u − 2v.

Show answer

Compute the four partial derivatives: ∂x/∂u = 3, ∂x/∂v = 1, ∂y/∂u = 1, ∂y/∂v = −2. The Jacobian determinant is det(J) = (3)(−2) − (1)(1) = −6 − 1 = −7. Taking the absolute value, |J| = 7. This means that every 1×1 square in the (u, v) plane maps to a parallelogram with area 7 in the (x, y) plane.

Question 3: INTERMEDIATE

Consider the transformation x = u², y = u + v. Compute the Jacobian determinant and determine for which values of u the transformation is non-degenerate (i.e., |J| ≠ 0).

Show answer

The partial derivatives are ∂x/∂u = 2u, ∂x/∂v = 0, ∂y/∂u = 1, ∂y/∂v = 1. The Jacobian determinant is det(J) = (2u)(1) − (0)(1) = 2u. The transformation is non-degenerate whenever 2u ≠ 0, that is, whenever u ≠ 0. At u = 0, the Jacobian is zero, meaning the transformation collapses area to zero locally — it squishes a patch of area in (u, v) into a line segment in (x, y). You should avoid using this transformation at u = 0.

Question 4: APPLIED

Use polar coordinates and the Jacobian to evaluate ∬ (x² + y²) dA over the annular region (ring) 1 ≤ x² + y² ≤ 9.

Show answer

In polar coordinates, x² + y² = r², and the region becomes 1 ≤ r ≤ 3, 0 ≤ θ ≤ 2π. The Jacobian is r, so the integral is ∫₀²π ∫₁³ r² · r dr dθ = ∫₀²π ∫₁³ r³ dr dθ. The inner integral gives [r⁴/4]₁³ = 81/4 − 1/4 = 80/4 = 20. The outer integral gives ∫₀²π 20 dθ = 40π. So the answer is 40π ≈ 125.66. Without the Jacobian factor r, you would have integrated r² dr and gotten a different (incorrect) result.

Question 5: CRITICAL THINKING

Suppose you create a custom transformation x = u + v, y = uv. Compute the Jacobian determinant and explain geometrically what it means that the Jacobian depends on both u and v (not a constant). What does this imply about the way this transformation distorts shapes?

Show answer

The partial derivatives are ∂x/∂u = 1, ∂x/∂v = 1, ∂y/∂u = v, ∂y/∂v = u. So det(J) = (1)(u) − (1)(v) = u − v. The absolute value |u − v| is not constant — it varies across the (u, v) plane. This means the transformation distorts different regions by different amounts. Near the line u = v, the Jacobian approaches zero, meaning the transformation compresses area dramatically (parallelogram patches become nearly degenerate). Far from that line, |u − v| is large, meaning patches get stretched significantly. Unlike a linear transformation (which has a constant Jacobian), this nonlinear map creates uneven distortion, and the Jacobian tells you exactly how the distortion varies from point to point.

SUMMARY

Lesson Summary

Summary

A coordinate transformation converts an integral from one coordinate system to another — for example, from Cartesian (x, y) to polar (r, θ). When you change coordinates, the tiny patches of area or volume you're summing can stretch or compress. The Jacobian determinant is the correction factor that accounts for this distortion. You compute it by building a matrix of partial derivatives (the Jacobian matrix) and taking its determinant. The absolute value |det(J)| multiplies the new area element to ensure the integral gives the correct result.

For polar coordinates, the Jacobian simplifies to r, explaining why dA = r dr dθ. For spherical coordinates, it becomes ρ² sin φ. The key idea is universal: whenever you change variables in a multiple integral, compute the Jacobian determinant and include its absolute value as a scaling factor to preserve the true area or volume being measured.