Jacobian & Coordinate Transformations — Jacobian Concept for Coordinate Transformations (Intro) - Calculus 3

SECTION 1

Historical Context & Motivation

Many integrals in physics and engineering are naturally expressed in coordinate systems other than Cartesian. Evaluating the area of a circle, the volume of a sphere, or the gravitational potential of an ellipsoidal mass distribution all become dramatically simpler when one switches from rectangular (x, y, z) coordinates to a system—polar, cylindrical, or spherical—whose symmetry matches the geometry of the problem. Yet changing coordinates inside an integral is not merely a matter of swapping variables: one must also account for how the transformation distorts infinitesimal area or volume elements. The mathematical object that captures this distortion is the Jacobian determinant, a concept whose roots stretch back to the early nineteenth century.

1829 • Carl Gustav Jacob Jacobi

The German mathematician Carl Jacobi systematically studies functional determinants—matrices of partial derivatives—and publishes foundational work that would later bear his name. His treatise formalizes the notion of relating differential elements across coordinate systems.

1841 • "De Determinantibus Functionalibus"

Jacobi publishes his landmark paper unifying the theory of functional determinants. He proves key properties including the chain rule for Jacobians and their multiplicative behavior under composition of transformations, laying the algebraic groundwork used in modern multivariable calculus.

1869–1880s • Integration in Curvilinear Coordinates

Mathematicians including Riemann and later Lebesgue rigorously justify the change-of-variable formula for multiple integrals. The Jacobian determinant becomes the standard correction factor when transforming dA or dV, providing a bridge between abstract analysis and practical computation.

20th Century • Differential Forms & Modern Use

Élie Cartan's differential forms provide a coordinate-free framework in which the Jacobian emerges naturally as a pullback operation. Today the Jacobian appears across physics (general relativity, continuum mechanics), computer graphics (texture mapping), and probability theory (multivariate change of density).

The fundamental question the Jacobian answers is deceptively simple: when we map a tiny rectangle in one coordinate system to a tiny parallelogram in another, by what factor does its area (or volume) change? Without accounting for this scaling factor, a change of variables in a double or triple integral would yield an incorrect result. The remainder of this lesson develops the Jacobian from first principles, shows how it works visually and algebraically, and prepares you to apply it confidently in integration problems.

SECTION 2

Core Principles & Definitions

Before diving into computation, it is essential to understand the conceptual pillars underlying the Jacobian. A coordinate transformation is a smooth, invertible mapping T that takes points from one coordinate system (u, v) to another (x, y). The Jacobian matrix of T collects all first-order partial derivatives, and its determinant—the Jacobian determinant—encodes the local scaling factor for area or volume elements. Because integration accumulates infinitesimal contributions, any stretching or compression introduced by the map must be corrected by multiplying the integrand by the absolute value of this determinant.

Coordinate Transformation T

A smooth mapping T : (u, v) → (x, y) that expresses Cartesian coordinates as functions x = g(u, v) and y = h(u, v). The transformation must be one-to-one (invertible) on the region of integration for the change-of-variable formula to apply.

Jacobian Matrix

The 2 × 2 (or n × n) matrix of all first-order partial derivatives: ∂x/∂u, ∂x/∂v in the first row and ∂y/∂u, ∂y/∂v in the second. Each column describes how a unit change in one new coordinate shifts the point in the original system.

Jacobian Determinant

The determinant of the Jacobian matrix, often denoted ∂(x, y)/∂(u, v). Its absolute value gives the factor by which infinitesimal areas are stretched or compressed under T. A zero Jacobian means the map collapses an area to a line or point.

Change-of-Variable Formula

The theorem stating ∬_R f(x, y) dA = ∬_S f(g(u,v), h(u,v)) |∂(x,y)/∂(u,v)| du dv, where S is the pre-image of R under T. The absolute Jacobian replaces dA with the corrected element in (u, v) coordinates.

✦

KEY TAKEAWAY

Think of a coordinate transformation as projecting a flexible rubber sheet onto a different surface. Some regions stretch, others compress. The Jacobian determinant is the local magnification factor at every point: it tells you how much a tiny patch of area on the original sheet grows or shrinks after the mapping. Just as an engineer must account for fabric stretch when cutting a curved panel from a flat pattern, you must multiply by |J| to ensure the integral 'measures' the correct amount of area in the new coordinate system.

SECTION 3

Visual Explanation — How the Jacobian Distorts Area Elements

The most intuitive way to grasp the Jacobian is to visualize what happens to a small rectangle in the (u, v) plane when mapped to the (x, y) plane. Consider a tiny rectangle with sides Δu and Δv centered at a point (u₀, v₀). Under the transformation T, each side of that rectangle becomes a vector determined by partial derivatives, and the rectangle maps approximately to a parallelogram. The area of that parallelogram equals |J| Δu Δv, where J is the Jacobian determinant evaluated at (u₀, v₀). The diagram below illustrates this mapping for a generic nonlinear transformation.

Left: a small rectangle of width Δu and height Δv in the (u, v) parameter plane. Right: the transformation T maps the rectangle to a parallelogram whose sides are the tangent vectors (∂T/∂u)Δu and (∂T/∂v)Δv. The area of this parallelogram equals |J| Δu Δv, where J is the Jacobian determinant.

In the diagram the two tangent vectors ∂T/∂u · Δu (violet) and ∂T/∂v · Δv (cyan) span the parallelogram. The magnitude of their cross product (in 2D, the absolute value of the 2 × 2 determinant formed by placing these vectors as rows or columns) yields the parallelogram's area. When you sum—that is, integrate—all such infinitesimal parallelograms over the parameter domain S, the factor |J| ensures each contribution is weighted by the correct local area scale.

SECTION 4

Mathematical Framework

We now formalize the ideas presented visually. Let T be a transformation defined by x = g(u, v) and y = h(u, v), where g and h have continuous first-order partial derivatives on an open region S in the uv-plane, and T is one-to-one on S. The Jacobian matrix of T is the 2 × 2 matrix whose entries are the partial derivatives of x and y with respect to u and v.

JACOBIAN MATRIX (2D)

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

Each column of J represents the image of a unit basis vector in the (u, v) system: the first column is (∂x/∂u, ∂y/∂u)^T and the second is (∂x/∂v, ∂y/∂v)^T.

JACOBIAN DETERMINANT (2D)

∂(x, y)/∂(u, v) = (∂x/∂u)(∂y/∂v) − (∂x/∂v)(∂y/∂u)

This scalar value is the signed area magnification factor. Its absolute value appears in the change-of-variable formula for double integrals.

CHANGE-OF-VARIABLE FORMULA (DOUBLE INTEGRAL)

∬_R f(x, y) dA = ∬_S f(g(u,v), h(u,v)) |∂(x,y)/∂(u,v)| du dv

R is the region in the xy-plane, S is its pre-image in the uv-plane under T, and |∂(x,y)/∂(u,v)| corrects the differential area element. The formula generalizes the single-variable substitution rule ∫ f(g(u)) g′(u) du.

For triple integrals the Jacobian matrix becomes 3 × 3, and the determinant is computed using cofactor expansion or the Sarrus rule. The change-of-variable formula extends naturally: dV = |∂(x, y, z)/∂(u, v, w)| du dv dw. The sign of the Jacobian determinant indicates whether the transformation preserves or reverses orientation; in integration we take the absolute value because we are measuring (positive) area or volume.

JACOBIAN DETERMINANT (3D)

∂(x,y,z)/∂(u,v,w) = det [ ∂x/∂u ∂x/∂v ∂x/∂w ] [ ∂y/∂u ∂y/∂v ∂y/∂w ] [ ∂z/∂u ∂z/∂v ∂z/∂w ]

For cylindrical coordinates (r, θ, z) → (r cos θ, r sin θ, z) the determinant evaluates to r, and for spherical coordinates (ρ, φ, θ) → (ρ sin φ cos θ, ρ sin φ sin θ, ρ cos φ) it evaluates to ρ² sin φ.

SECTION 5

Detailed Breakdown — The Polar Coordinate Jacobian

The most widely encountered coordinate transformation in two dimensions is from polar coordinates (r, θ) to Cartesian (x, y). Because this example is both concrete and practically important, it serves as the ideal test case for the Jacobian machinery. The transformation is x = r cos θ, y = r sin θ, and you have likely already used the area element dA = r dr dθ in earlier courses—possibly without a rigorous justification. The Jacobian provides exactly that justification.

Left: a small rectangle [r₀, r₀+Δr] × [θ₀, θ₀+Δθ] in the polar parameter plane. Right: the same patch mapped to the Cartesian plane becomes a curved annular wedge whose area is approximately r₀ Δr Δθ. The factor r₀ is exactly the absolute value of the Jacobian determinant evaluated at (r₀, θ₀).

Let us compute the Jacobian explicitly. With x = r cos θ and y = r sin θ, the partial derivatives are ∂x/∂r = cos θ, ∂x/∂θ = −r sin θ, ∂y/∂r = sin θ, and ∂y/∂θ = r cos θ. The Jacobian determinant is then:

POLAR JACOBIAN

∂(x, y)/∂(r, θ) = (cos θ)(r cos θ) − (−r sin θ)(sin θ) = r cos²θ + r sin²θ = r

Since r ≥ 0 in polar coordinates, |J| = r. This confirms the familiar area element dA = r dr dθ.

💡

Why r, not r²?

A common source of confusion: the polar Jacobian is r (not r²) because we are in two dimensions. The factor r² sin φ that appears in spherical coordinates arises in three dimensions, where the 3 × 3 Jacobian determinant accounts for an additional angular variable φ. Always match the dimension of the Jacobian matrix to the number of variables.

SECTION 6

Worked Example — Evaluating a Double Integral via Jacobian

Consider the integral ∬_R (x² + y²) dA, where R is the disk x² + y² ≤ 4. We will convert to polar coordinates, apply the Jacobian, and evaluate.

Integrating x² + y² Over a Disk of Radius 2

1
Step 1 — Define the Transformation

Use x = r cos θ and y = r sin θ. The disk x² + y² ≤ 4 maps to the rectangle S: 0 ≤ r ≤ 2, 0 ≤ θ ≤ 2π in the (r, θ) plane.

S = { (r, θ) : 0 ≤ r ≤ 2, 0 ≤ θ ≤ 2π }
2
Step 2 — Rewrite the Integrand

In polar coordinates, x² + y² = r². Therefore the integrand becomes f(r, θ) = r².

f = r²
3
Step 3 — Attach the Jacobian

We computed |∂(x,y)/∂(r,θ)| = r. The area element transforms as dA = r dr dθ. Multiplying the integrand by the Jacobian gives the full integrand in polar: r² · r = r³.

Integrand with Jacobian: r³ dr dθ
4
Step 4 — Set Up and Evaluate the Iterated Integral

∫₀²π ∫₀² r³ dr dθ. Evaluate the inner integral first: ∫₀² r³ dr = [r⁴/4]₀² = 16/4 = 4. Then the outer integral: ∫₀²π 4 dθ = 4 · 2π = 8π.

∬_R (x² + y²) dA = 8π
5
Step 5 — Verify Reasonableness

The integrand x² + y² ranges from 0 at the origin to 4 on the boundary. The area of the disk is 4π. A rough average value of the integrand is about 2 (the midpoint), giving an estimate of 2 × 4π = 8π, which matches exactly. This confirms our result.

Estimate check: average value × area ≈ 2 × 4π = 8π ✓

SECTION 7

Strengths, Limitations & When to Use the Jacobian

The Jacobian-based change of variables is one of the most powerful tools in multivariable calculus, but it is not without subtleties and limitations. Knowing when a coordinate transformation simplifies an integral—and when it may not—is a skill that develops with practice. The table below contrasts the benefits and pitfalls.

Aspect	Strengths	Limitations / Pitfalls
Geometric matching	Choosing coordinates that match the region's symmetry can convert complicated boundaries into simple rectangles, drastically simplifying limits of integration.	If no natural symmetry exists, the transformed region may have boundaries just as complex as the original, negating the advantage.
Integrand simplification	Expressions like x² + y² become r² in polar, eliminating square roots and reducing algebraic complexity.	The Jacobian factor (e.g., r, or ρ² sin φ) adds a multiplicative term that may complicate the integrand, especially if the original function was already simple.
One-to-one requirement	When satisfied, the formula provides an exact, rigorous equivalence between the two integrals.	At isolated points where the Jacobian vanishes (e.g., r = 0 in polar), the transformation is not locally invertible. These sets of measure zero do not affect the integral but can cause conceptual confusion.
Higher dimensions	The framework generalizes seamlessly: 2D, 3D, and even n-dimensional Jacobians follow the same determinant recipe.	Computing an n × n determinant grows in complexity; for dimensions ≥ 4, alternative techniques (wedge products, differential forms) are often preferred.

When coordinate transformations help — and when they don't

🧭

WHEN TO TRANSFORM

A good rule of thumb: use a coordinate transformation when the boundary of the region or the integrand contains an expression (such as x² + y², or x/y) that simplifies dramatically in another coordinate system. If neither the region nor the integrand simplifies, the Jacobian overhead is not worth it.

SECTION 8

Connection to Advanced Theory — Differential Forms & Higher Dimensions

The Jacobian determinant you have learned in this lesson is the introductory face of a much deeper structure. In modern mathematics the change-of-variable formula is understood through the language of differential forms. A differential 2-form dx ∧ dy transforms under a map T by 'pulling back' through the Jacobian, producing (∂(x,y)/∂(u,v)) du ∧ dv automatically. This perspective unifies the substitution rule across all dimensions and extends naturally to curved manifolds encountered in differential geometry and general relativity.

Feature	Introductory Jacobian (This Lesson)	Differential Forms (Advanced)
Objects	Jacobian matrix and its determinant computed explicitly from partial derivatives	Exterior algebra: k-forms, wedge products, pullback maps
Setting	Euclidean ℝ² or ℝ³ with explicit coordinate formulas	Smooth manifolds of arbitrary dimension, coordinate-free
Orientation	Handled by taking absolute value of determinant	Built in: signed forms track orientation naturally
Generalization	Extends to n × n determinants with increasing computational cost	Stokes' theorem subsumes divergence, Green's, and Gauss's theorems in one statement

Introductory Jacobian vs. Differential Forms

For now the explicit Jacobian determinant is all you need for evaluating double and triple integrals via substitution. As you advance into courses on real analysis, differential geometry, or mathematical physics, you will encounter the pullback of forms as the natural generalization. Mastering the matrix-and-determinant approach here builds the concrete intuition that makes the abstract theory accessible later.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain in your own words why the Jacobian determinant is needed when changing variables in a double integral. What would go wrong if you simply replaced x with g(u, v) and y with h(u, v) without including |∂(x,y)/∂(u,v)|?

Show answer

Without the Jacobian, the integral would still sum contributions over the parameter domain S, but each infinitesimal rectangle du dv in S does not in general correspond to the same amount of area in the xy-plane. The transformation T stretches and compresses different parts of the domain by different amounts. The Jacobian determinant provides the local area magnification factor: it tells us that the true area of the image of a small du × dv rectangle is approximately |J| du dv, not just du dv. Omitting |J| means every patch is treated as if it contributes equally, which produces an incorrect total—analogous to computing a weighted sum with all weights set to one instead of their actual values.

Question 2: BASIC CALCULATION

Compute the Jacobian determinant ∂(x, y)/∂(u, v) for the transformation x = u² − v², y = 2uv.

Show answer

Compute the four partial derivatives: ∂x/∂u = 2u, ∂x/∂v = −2v, ∂y/∂u = 2v, ∂y/∂v = 2u. The Jacobian determinant is (2u)(2u) − (−2v)(2v) = 4u² + 4v² = 4(u² + v²). Observe that this Jacobian is always non-negative and equals zero only at the origin, confirming the transformation is locally invertible everywhere except (0, 0). Note also that x + iy = (u + iv)², so this is the squaring map in the complex plane—the Jacobian 4|w|² reflects how the squaring function doubles angles and squares moduli.

Question 3: INTERMEDIATE

Let T be the transformation x = u + v, y = u − v. Compute the Jacobian determinant and use it to evaluate ∬_R 1 dA, where R is the square with vertices (0,0), (1,1), (2,0), (1,−1) in the xy-plane. (Hint: determine what R looks like in the uv-plane.)

Show answer

First, the Jacobian: ∂x/∂u = 1, ∂x/∂v = 1, ∂y/∂u = 1, ∂y/∂v = −1, so J = (1)(−1) − (1)(1) = −2, and |J| = 2. The inverse transformation gives u = (x + y)/2 and v = (x − y)/2. Mapping the four vertices: (0,0) → (0,0), (1,1) → (1,0), (2,0) → (1,1), (1,−1) → (0,1). So R maps to the unit square S = [0,1] × [0,1] in the uv-plane. The integral becomes ∬_S 1 · |J| du dv = ∬_S 2 du dv = 2 × 1 = 2. We can verify: the area of the diamond R in the xy-plane is ½ × diagonal₁ × diagonal₂ = ½ × 2 × 2 = 2. ✓

Question 4: APPLIED

The mass of a thin plate occupying the region bounded by x² + y² = 9 with density ρ(x, y) =

e^{−(x²+y²)}

is given by M = ∬_R

e^{−(x²+y²)}

dA. Convert to polar coordinates using the Jacobian and evaluate the integral.

Show answer

Under polar coordinates, x² + y² = r² and the Jacobian gives dA = r dr dθ. The region R maps to S: 0 ≤ r ≤ 3, 0 ≤ θ ≤ 2π. Thus M = ∫₀²π ∫₀³

e^{−r²}

· r dr dθ. The inner integral: let w = r², dw = 2r dr, so ∫₀³ r

e^{−r²}

dr = ½ ∫₀⁹

e^{−w}

dw = ½[−

e^{−w}

]₀⁹ = ½(1 −

e^{−9}

). The outer integral contributes a factor of 2π, yielding M = 2π · ½(1 −

e^{−9}

) = π(1 −

e^{−9}

). Since

e^{−9}

≈ 0.000123, the mass is approximately π(0.999877) ≈ 3.1412. This integral would be extremely difficult—arguably impossible in closed form—without the polar substitution, demonstrating the practical power of the Jacobian.

Question 5: CRITICAL THINKING

Suppose T₁ : (u, v) → (s, t) has Jacobian determinant J₁, and T₂ : (s, t) → (x, y) has Jacobian determinant J₂. Prove that the Jacobian determinant of the composite transformation T₂ ∘ T₁ : (u, v) → (x, y) equals J₁ · J₂. (Hint: use the chain rule for partial derivatives and the multiplicative property of determinants.)

Show answer

The Jacobian matrix of the composite T₂ ∘ T₁ is found via the multivariable chain rule. If we denote the Jacobian matrices as matrices A = [∂(s,t)/∂(u,v)] for T₁ and B = [∂(x,y)/∂(s,t)] for T₂, then by the chain rule ∂(x,y)/∂(u,v) = B · A (matrix product). Taking determinants of both sides and applying the property det(BA) = det(B) det(A), we obtain ∂(x,y)/∂(u,v) = J₂ · J₁. This is sometimes called the chain rule for Jacobians. It mirrors the single-variable chain rule (f ∘ g)′ = f′ · g′ and extends to arbitrary dimension. A practical consequence: if you compose two changes of variables, you can compute each Jacobian separately and multiply them, rather than deriving the composite transformation's Jacobian from scratch.

SUMMARY

Lesson Summary

Summary

A coordinate transformation T : (u, v) → (x, y) maps regions between coordinate systems, and the Jacobian matrix collects all first-order partial derivatives of that map into a single matrix. Its determinant—the Jacobian determinant ∂(x, y)/∂(u, v)—measures the local area (or volume) scaling factor at every point. In the change-of-variable formula for multiple integrals, the absolute value |J| replaces the differential area element dA (or dV) to account for how the transformation distorts infinitesimal patches.

For polar coordinates (r, θ), the Jacobian is r, yielding dA = r dr dθ; for spherical coordinates (ρ, φ, θ), it is ρ² sin φ. The key to successful application is choosing a transformation whose symmetry matches the region of integration or the integrand, and always remembering to multiply by |J|. This determinant-based framework extends to arbitrary dimension and, through differential forms, to integration on curved manifolds in advanced mathematics and physics.