Boo-AI — Master Artificial Intelligence by Building from Scratch

Learning Objectives

By the end of this section, you will be able to:

Convert double integrals from Cartesian to polar coordinates using the transformation $x = r\cos\theta$ , $y = r\sin\theta$
Explain why the Jacobian factor $r$ appears in polar integrals geometrically and algebraically
Identify regions that are naturally described in polar coordinates (disks, sectors, annuli, cardioids, roses)
Set up appropriate limits of integration for polar regions
Evaluate double integrals over polar regions using the formula $\iint_R f(x,y)\,dA = \iint_R f(r\cos\theta, r\sin\theta)\,r\,dr\,d\theta$
Apply polar integration to compute areas, volumes, and physical quantities in circularly symmetric problems

The Big Picture: Why Polar Coordinates?

"Polar coordinates are the natural language for anything involving circles, spirals, or radial symmetry — they turn complicated Cartesian boundaries into simple constants."

In the previous sections, we evaluated double integrals over rectangular regions and general regions described by Cartesian coordinates. But consider trying to integrate over a circular disk: in Cartesian coordinates, the boundary $x^2 + y^2 = R^2$ leads to complicated limits like $-\sqrt{R^2 - x^2} \le y \le \sqrt{R^2 - x^2}$ .

Polar coordinates offer a solution: circles become simply $r = R$ (constant). The same disk that requires square roots in Cartesian becomes the elegant region $0 \le r \le R, \, 0 \le \theta \le 2\pi$ .

When to Use Polar Coordinates

Ideal for Polar ✓

Circular disks and annuli (rings)
Sectors (pie slices) of circles
Integrands containing $x^2 + y^2$
Cardioids, limaçons, rose curves, spirals
Any region with circular symmetry
Gaussian functions $e^{-(x^2+y^2)}$

Keep Cartesian ✗

Rectangles with sides parallel to axes
Regions bounded by straight lines
Integrands that are simpler in x, y form
Triangles and parallelograms
Regions without circular symmetry

The Key Transformation

When converting from Cartesian to polar coordinates, remember three things:

Replace $x = r\cos\theta$ and $y = r\sin\theta$
Replace $dA = dx\,dy$ with $dA = r\,dr\,d\theta$
Convert the region boundaries to polar form

Historical Context

The polar coordinate system has ancient roots, but its modern form emerged in the 17th century. Isaac Newton used polar-like coordinates in his Principia Mathematica (1687) to describe planetary orbits, though he did not use our modern notation.

Jacob Bernoulli (1654–1705) was among the first to systematically use polar coordinates, studying the logarithmic spiral and other curves that are naturally expressed in polar form. He famously asked to have the logarithmic spiral $r = ae^{b\theta}$ engraved on his tombstone with the inscription "Eadem mutata resurgo" (Though changed, I rise again the same) — referring to the spiral's self-similar property.

The theory of double integrals in polar coordinates was developed by Euler and Lagrange in the 18th century, who recognized that the "Jacobian factor" $r$ was essential for correct area calculations.

The Gaussian Integral

One of the most celebrated applications of polar coordinates is evaluating the Gaussian integral $\int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi}$ . This integral is impossible to evaluate using elementary Cartesian methods, but becomes tractable when we square it and convert to polar coordinates — a technique discovered by Poisson.

Review of Polar Coordinates

In polar coordinates, every point in the plane is described by two numbers:

r — the radial distance from the origin (always $r \ge 0$ for integration)
θ (theta) — the angle measured counterclockwise from the positive x-axis

Polar ↔ Cartesian Conversion

Polar → Cartesian

x = r \cos\theta

y = r \sin\theta

Cartesian → Polar

r = \sqrt{x^2 + y^2}

\theta = \arctan\left(\frac{y}{x}\right)

Key identity: $x^2 + y^2 = r^2$

Interactive Polar Coordinate Explorer

Use this interactive tool to explore how polar coordinates $(r, \theta)$ map to Cartesian coordinates $(x, y)$ . Drag the sliders to change the radius and angle, and observe how the point moves and how the conversion formulas work.

🎯Polar Coordinate System Explorer

Polar Coordinates

Radius r: 2.00

Angle θ: 45° (0.25π rad)

Coordinate Conversion

x = r cos θ =1.4142

y = r sin θ =1.4142

r = √(x² + y²) =2.0000

θ = atan2(y, x) =0.7854 rad

Cartesian Grid

Polar Grid

What to Explore

Notice how $x = r\cos\theta$ and $y = r\sin\theta$ change as you adjust the sliders
Set $\theta = 0$ — the point lies on the positive x-axis
Set $\theta = \pi/2$ — the point lies on the positive y-axis
Set $\theta = \pi$ — the point lies on the negative x-axis

The Polar Area Element: Why dA = r dr dθ

The most important conceptual step in polar integration is understanding why the area element changes from $dA = dx\,dy$ to $dA = r\,dr\,d\theta$ . The factor $r$ is not optional — it's essential!

Geometric Intuition

Consider a small polar "rectangle" bounded by:

Two circles of radii $r$ and $r + dr$
Two radial lines at angles $\theta$ and $\theta + d\theta$

This region is actually a curved sector (like a slice of an annulus). Its area is approximately:

\text{Area} \approx (\text{arc length}) \times (\text{radial width}) = (r\,d\theta) \times dr = r\,dr\,d\theta

The arc length at radius $r$ is $r \cdot d\theta$ (not just $d\theta$ !). This is the geometric origin of the $r$ factor.

The Jacobian: The Algebraic Explanation

The factor $r$ can also be derived rigorously using the Jacobian determinant of the coordinate transformation.

For the transformation $x = r\cos\theta, \, y = r\sin\theta$ , the Jacobian matrix is:

J = \begin{pmatrix} \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \theta} \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \theta} \end{pmatrix} = \begin{pmatrix} \cos\theta & -r\sin\theta \\ \sin\theta & r\cos\theta \end{pmatrix}

The Jacobian determinant is:

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

The Jacobian Result

dA = dx\,dy = |J|\,dr\,d\theta = r\,dr\,d\theta

The Jacobian tells us how area elements stretch or shrink under the coordinate transformation. For polar coordinates, areas stretch by a factor of $r$ .

Interactive Jacobian Demonstration

This interactive demonstration shows how polar area elements depend on the radial distance. Adjust the parameters to see how the area $dA = r\,dr\,d\theta$ changes.

📐The Jacobian: Why dA = r dr dθ

Polar Area Element

Adjust Parameters

r = 2.00 (radial distance)

θ = 45°

dr = 0.50

dθ = 0.30 rad (17°)

Area Calculations

Polar area dA:r · dr · dθ = 0.3375

Arc length at r:r · dθ = 0.6000

Radial width:dr = 0.5000

The Key Insight

The area of a polar "rectangle" is approximately (arc length) × (radial width) = r·dθ × dr. This is why we need the r factor when converting from Cartesian to polar coordinates. The Jacobian determinant captures this stretching:

J = |∂(x,y)/∂(r,θ)| = r

Compare with Cartesian

In Cartesian coordinates, a small rectangle has area dA = dx · dy. The sides are always the same size regardless of position.

In polar coordinates, the arc length r·dθ depends on the distance from the origin — farther out means longer arcs!

The Complete Conversion Formula

Double Integral in Polar Coordinates

If $R$ is a polar region described by $\alpha \le \theta \le \beta$ and $h_1(\theta) \le r \le h_2(\theta)$ , then:

\iint_R f(x, y)\,dA = \int_{\alpha}^{\beta} \int_{h_1(\theta)}^{h_2(\theta)} f(r\cos\theta, r\sin\theta) \cdot r\,dr\,d\theta

Step-by-Step Conversion Process

Sketch the region and identify if it has circular symmetry
Convert the boundaries to polar form:
- $x^2 + y^2 = R^2$ becomes $r = R$
- $y = x$ becomes $\theta = \pi/4$
- $x = 0$ (positive y-axis) becomes $\theta = \pi/2$
Replace the integrand: $f(x, y) \to f(r\cos\theta, r\sin\theta)$
Add the Jacobian: Multiply by $r$
Set up limits: Usually $\theta$ is the outer integral, $r$ is the inner
Evaluate the iterated integral

Types of Polar Regions

Polar coordinates make certain regions particularly simple to describe:

Region	Description	Polar Form	Area Formula
Full Disk	Circle of radius R	0 ≤ r ≤ R, 0 ≤ θ ≤ 2π	πR²
Sector	Pie slice from 0 to angle α	0 ≤ r ≤ R, 0 ≤ θ ≤ α	αR²/2
Annulus	Ring between r₁ and r₂	r₁ ≤ r ≤ r₂, 0 ≤ θ ≤ 2π	π(r₂² - r₁²)
Cardioid	Heart-shaped curve	0 ≤ r ≤ 1 + cos θ, 0 ≤ θ ≤ 2π	3π/2
Rose (3 petals)	r = a sin(3θ)	Depends on petal	πa²/4 per petal

Interactive Polar Region Explorer

Explore different types of polar regions and see how they are partitioned for integration. Notice how the partition elements (polar "rectangles") are actually curved sectors whose areas depend on their distance from the origin.

🌀Polar Region Types for Double Integrals

Select Region Type:

Circular Sector

A pie-slice region from 0 to 60°

∫₀^(π/3) ∫₀^2 f(r,θ) r dr dθ

Partition Segments: 20

More segments = finer approximation of the area

Why the r Factor?

Notice how polar "rectangles" are actually wedge-shaped sectors. The area element dA = r dr dθ accounts for the fact that outer sectors are larger than inner ones — they sweep a larger arc for the same angle dθ.

Region Bounds

θ: [0.00π, 0.33π]

r: [r₁(θ), r₂(θ)]

Worked Examples

Example 1: Area of a Disk

Find the area of a disk of radius $R$ centered at the origin.

Solution: We integrate $f = 1$ over the disk.

A = \iint_D 1\,dA = \int_0^{2\pi} \int_0^R r\,dr\,d\theta

= \int_0^{2\pi} \left[\frac{r^2}{2}\right]_0^R d\theta = \int_0^{2\pi} \frac{R^2}{2}\,d\theta

= \frac{R^2}{2} \cdot 2\pi = \pi R^2

This confirms the well-known formula for the area of a circle!

Example 2: Integral of $x^2 + y^2$

Evaluate $\iint_D (x^2 + y^2)\,dA$ where $D$ is the disk $x^2 + y^2 \le 4$ .

Solution: In polar coordinates, $x^2 + y^2 = r^2$ .

\iint_D (x^2 + y^2)\,dA = \int_0^{2\pi} \int_0^2 r^2 \cdot r\,dr\,d\theta

= \int_0^{2\pi} \int_0^2 r^3\,dr\,d\theta = \int_0^{2\pi} \left[\frac{r^4}{4}\right]_0^2 d\theta

= \int_0^{2\pi} 4\,d\theta = 8\pi

Example 3: The Gaussian Integral (Famous Result)

Evaluate $I = \int_{-\infty}^{\infty} e^{-x^2}\,dx$ .

Solution: This is one of the most beautiful applications of polar coordinates. We cannot evaluate $I$ directly, but we can evaluate $I^2$ .

I^2 = \int_{-\infty}^{\infty} e^{-x^2}\,dx \cdot \int_{-\infty}^{\infty} e^{-y^2}\,dy = \iint_{\mathbb{R}^2} e^{-(x^2+y^2)}\,dx\,dy

Now convert to polar coordinates. Since $x^2 + y^2 = r^2$ :

I^2 = \int_0^{2\pi} \int_0^{\infty} e^{-r^2} \cdot r\,dr\,d\theta

The inner integral is perfect for substitution. Let $u = r^2$ , so $du = 2r\,dr$ :

\int_0^{\infty} e^{-r^2} r\,dr = \frac{1}{2}\int_0^{\infty} e^{-u}\,du = \frac{1}{2}[-e^{-u}]_0^{\infty} = \frac{1}{2}

Therefore:

I^2 = \int_0^{2\pi} \frac{1}{2}\,d\theta = \pi

So $I = \sqrt{\pi}$ !

Why This Matters

The Gaussian integral appears throughout probability theory (the normal distribution), quantum mechanics (wave functions), and machine learning (Gaussian processes, softmax normalization). This polar coordinate trick is the standard way to derive its value.

Applications in Science and Engineering

Physics: Moments of Inertia

For a flat disk of radius $R$ and uniform density $\rho$ , the moment of inertia about the center is:

I = \iint_D \rho r^2\,dA = \rho \int_0^{2\pi} \int_0^R r^2 \cdot r\,dr\,d\theta = \rho \cdot 2\pi \cdot \frac{R^4}{4} = \frac{1}{2}\pi\rho R^4

If the total mass is $M = \pi\rho R^2$ , then $I = \frac{1}{2}MR^2$ .

Probability: Circular Normal Distributions

The standard bivariate normal distribution has density:

f(x, y) = \frac{1}{2\pi} e^{-\frac{x^2 + y^2}{2}}

Verifying it integrates to 1 is natural in polar coordinates:

\int_0^{2\pi} \int_0^{\infty} \frac{1}{2\pi} e^{-r^2/2} \cdot r\,dr\,d\theta = \frac{1}{2\pi} \cdot 2\pi \cdot 1 = 1

Engineering: Antenna Patterns

Many antenna radiation patterns are described in polar coordinates. The power radiated in different directions often follows formulas like $P(\theta) = P_0 \cos^2\theta$ , naturally integrated in polar form.

Applications in Machine Learning

Gaussian Processes

Gaussian process priors often use radial basis function (RBF) kernels $k(x, x') = \exp(-\|x-x'\|^2/2\sigma^2)$ . Understanding polar integrals helps analyze kernel properties.

Normalization Constants

Many ML models require computing normalization constants for distributions. The softmax temperature scaling and Gaussian normalizations all trace back to polar integral techniques.

Image Processing

Circular and radial filters for edge detection (e.g., Laplacian of Gaussian) are naturally analyzed in polar coordinates. FFT-based methods often leverage radial symmetry.

Variational Autoencoders

VAEs use Gaussian latent spaces where the KL divergence involves Gaussian integrals. Understanding these distributions deeply requires polar coordinate analysis.

Python Implementation

Numerical Polar Integration

Numerical Double Integrals in Polar Coordinates

🐍polar_integration.py

Explanation(5)

Code(72)

3Polar Double Integration

This function numerically computes double integrals in polar coordinates. The key is including the Jacobian factor r in the integrand.

17The Jacobian Factor

When converting to polar coordinates, we must multiply by r. This is the Jacobian determinant that accounts for the stretching of area elements.

29Area of a Disk

To find area, we integrate f = 1. The polar integral ∫₀^(2π) ∫₀^R r dr dθ = ∫₀^(2π) [r²/2]₀^R dθ = ∫₀^(2π) R²/2 dθ = πR².

39Cardioid Area

The cardioid r = 1 + cos(θ) is a classic polar curve. Its area is ∫₀^(2π) ∫₀^(1+cosθ) r dr dθ = 3π/2, a beautiful result.

48Converting Cartesian Functions

When the integrand involves x² + y², we use the polar identity x² + y² = r². This often simplifies integrals dramatically.

67 lines without explanation

1import numpy as np
2from scipy import integrate
3import matplotlib.pyplot as plt
4
5def polar_double_integral(f, r_inner, r_outer, theta_min, theta_max,
6                          nr=50, ntheta=50):
7    """
8    Compute a double integral in polar coordinates using numerical methods.
9
10    The integral is: ∫∫ f(r,θ) r dr dθ
11
12    Parameters:
13    - f: function of (r, theta) to integrate
14    - r_inner, r_outer: functions of theta giving inner/outer r bounds
15    - theta_min, theta_max: angle bounds
16    - nr, ntheta: number of integration points
17
18    Returns: approximate value of the integral
19    """
20    def integrand(theta, r):
21        # The integrand includes the Jacobian factor r
22        return f(r, theta) * r
23
24    def inner_integral(theta):
25        r_low = r_inner(theta)
26        r_high = r_outer(theta)
27        result, _ = integrate.quad(lambda r: integrand(theta, r),
28                                    r_low, r_high)
29        return result
30
31    result, error = integrate.quad(inner_integral, theta_min, theta_max)
32    return result
33
34# Example 1: Area of a disk of radius R
35R = 2.0
36f_constant = lambda r, theta: 1  # f = 1 gives area
37
38area = polar_double_integral(
39    f=f_constant,
40    r_inner=lambda theta: 0,
41    r_outer=lambda theta: R,
42    theta_min=0,
43    theta_max=2*np.pi
44)
45print(f"Area of disk (radius {R}): {area:.6f}")
46print(f"Expected (πR²): {np.pi * R**2:.6f}")
47
48# Example 2: Area inside cardioid r = 1 + cos(θ)
49cardioid_area = polar_double_integral(
50    f=f_constant,
51    r_inner=lambda theta: 0,
52    r_outer=lambda theta: 1 + np.cos(theta),
53    theta_min=0,
54    theta_max=2*np.pi
55)
56print(f"\nArea of cardioid: {cardioid_area:.6f}")
57print(f"Expected (3π/2): {3*np.pi/2:.6f}")
58
59# Example 3: ∫∫ (x² + y²) dA over disk of radius 2
60# In polar: x² + y² = r²
61def f_r_squared(r, theta):
62    return r**2  # This is x² + y² in polar
63
64integral_r2 = polar_double_integral(
65    f=f_r_squared,
66    r_inner=lambda theta: 0,
67    r_outer=lambda theta: 2,
68    theta_min=0,
69    theta_max=2*np.pi
70)
71print(f"\n∫∫(x² + y²) dA over disk radius 2: {integral_r2:.6f}")
72print(f"Expected (8π): {8*np.pi:.6f}")

Visualization of Polar Partitions

Visualizing Polar Area Elements

🐍polar_visualization.py

Explanation(4)

Code(106)

6Polar Partition Visualization

This function creates a visual representation of how a polar region is divided into small sectors. Notice how sectors farther from the origin are larger!

27Wedge Sectors

Each polar 'rectangle' is actually a wedge-shaped sector. The area depends on the radial distance — this is why we need the r factor.

35Color by Area

The coloring shows that outer sectors have larger areas. This visualization makes it intuitive why dA = r dr dθ, not just dr dθ.

56Linear Growth

The bar chart shows that area elements grow linearly with r. At twice the distance, the arc length (and thus area) doubles.

102 lines without explanation

1import numpy as np
2import matplotlib.pyplot as plt
3from matplotlib.patches import Wedge
4from matplotlib.collections import PatchCollection
5
6def visualize_polar_partition(r_max=3, n_r=5, n_theta=12):
7    """
8    Visualize how a polar region is partitioned for integration.
9    Each polar 'rectangle' is actually a curved sector.
10    """
11    fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 6))
12
13    # Left: Show the partition
14    ax1.set_aspect('equal')
15    ax1.set_xlim(-r_max - 0.5, r_max + 0.5)
16    ax1.set_ylim(-r_max - 0.5, r_max + 0.5)
17
18    dr = r_max / n_r
19    dtheta = 2 * np.pi / n_theta
20
21    patches = []
22    colors = []
23
24    for i in range(n_r):
25        for j in range(n_theta):
26            r_inner = i * dr
27            r_outer = (i + 1) * dr
28            theta_start = j * dtheta * 180 / np.pi
29            theta_end = (j + 1) * dtheta * 180 / np.pi
30
31            # Create a wedge (sector)
32            wedge = Wedge((0, 0), r_outer, theta_start, theta_end,
33                          width=dr)
34            patches.append(wedge)
35
36            # Color based on area (proportional to r)
37            r_mid = (r_inner + r_outer) / 2
38            colors.append(r_mid / r_max)
39
40    collection = PatchCollection(patches, cmap='viridis', alpha=0.7)
41    collection.set_array(np.array(colors))
42    ax1.add_collection(collection)
43
44    ax1.axhline(y=0, color='k', linewidth=0.5)
45    ax1.axvline(x=0, color='k', linewidth=0.5)
46    ax1.set_title('Polar Partition: dA = r·dr·dθ', fontsize=14)
47    ax1.set_xlabel('x')
48    ax1.set_ylabel('y')
49
50    # Add colorbar
51    cbar = plt.colorbar(collection, ax=ax1)
52    cbar.set_label('Relative area (proportional to r)')
53
54    # Right: Show how area elements grow with r
55    r_values = np.linspace(0.5, r_max, 20)
56    areas = r_values * dr * dtheta  # Area = r * dr * dθ
57
58    ax2.bar(r_values, areas, width=dr*0.8, color='steelblue', alpha=0.7)
59    ax2.set_xlabel('Radial distance r')
60    ax2.set_ylabel('Area element dA = r·dr·dθ')
61    ax2.set_title('Area Elements Grow Linearly with r', fontsize=14)
62    ax2.grid(True, alpha=0.3)
63
64    # Add annotation
65    ax2.annotate('Outer sectors have\nlarger areas!',
66                 xy=(r_max*0.8, areas[-1]*0.9),
67                 fontsize=10, ha='center')
68
69    plt.tight_layout()
70    plt.savefig('polar_partition.png', dpi=150, bbox_inches='tight')
71    plt.show()
72
73visualize_polar_partition()
74
75# Demonstrate area calculation for specific regions
76def area_by_polar_integration():
77    """
78    Calculate areas of various polar regions.
79    """
80    results = []
81
82    # 1. Full disk of radius R
83    R = 2
84    area_disk = np.pi * R**2
85    results.append(("Full disk (r ≤ 2)", area_disk))
86
87    # 2. Annulus (ring) 1 ≤ r ≤ 3
88    r1, r2 = 1, 3
89    area_annulus = np.pi * (r2**2 - r1**2)
90    results.append((f"Annulus ({r1} ≤ r ≤ {r2})", area_annulus))
91
92    # 3. Sector (pie slice) r ≤ 2, 0 ≤ θ ≤ π/4
93    R, theta = 2, np.pi/4
94    area_sector = 0.5 * R**2 * theta
95    results.append(("Sector (r ≤ 2, θ ≤ π/4)", area_sector))
96
97    # 4. Cardioid r = 1 + cos(θ)
98    area_cardioid = 3 * np.pi / 2
99    results.append(("Cardioid r = 1 + cos(θ)", area_cardioid))
100
101    print("Area Calculations:")
102    print("=" * 40)
103    for name, area in results:
104        print(f"{name}: {area:.4f}")
105
106area_by_polar_integration()

Common Mistakes to Avoid

Mistake 1: Forgetting the r Factor

The most common error is writing $dA = dr\,d\theta$ instead of $dA = r\,dr\,d\theta$ . This will give you the wrong answer every time! The $r$ factor is the Jacobian and is essential.

Mistake 2: Wrong Limits for a Circle

For a full circle, $\theta$ goes from $0$ to $2\pi$ , not from $0$ to $\pi$ (which gives only a semicircle).

Mistake 3: Not Converting the Integrand

If the Cartesian integrand is $f(x, y) = x$ , you must write $f = r\cos\theta$ , not just leave it as $x$ .

Mistake 4: Negative r Values

In integration, $r \ge 0$ always. We don't use negative radii. If a curve like $r = \sin\theta$ gives negative values, we need to carefully set up the limits.

Pro Tip: Check with Known Results

After computing a polar integral, verify with known results:

Area of a disk of radius 2 = 4π
Area of a semicircle of radius 2 = 2π
Area inside cardioid r = 1 + cos θ = 3π/2

Test Your Understanding

📝Test Your Understanding

0 / 8

Q1When converting a double integral from Cartesian to polar coordinates, what factor must be included in the integrand?

Q2What is the area of a region bounded by the cardioid r = 1 + cos(θ)?

Q3To integrate over a disk of radius 3 centered at the origin, what are the limits of integration?

Q4What is ∫∫_D (x² + y²) dA where D is the disk x² + y² ≤ 4?

Q5Which type of region is most naturally described in polar coordinates?

Q6To find the area inside r = 2sin(θ), what are the correct limits for θ?

Q7In the integral ∫₀^(π/4) ∫₁^(sec θ) f(r,θ) r dr dθ, what Cartesian region is being integrated over?

Q8What is the Jacobian determinant for the transformation x = r cos θ, y = r sin θ?

Summary

Double integrals in polar coordinates are powerful tools for evaluating integrals over regions with circular symmetry. The key insight is understanding why and how the area element transforms.

Key Formulas

Concept	Formula	Notes
Coordinate transform	x = r cos θ, y = r sin θ	Basic conversion
Area element	dA = r dr dθ	The r is the Jacobian
Jacobian	\|∂(x,y)/∂(r,θ)\| = r	Measures area stretching
Key identity	x² + y² = r²	Simplifies many integrands
Full disk	∫₀^(2π) ∫₀^R f r dr dθ	Complete circle
Gaussian	∫∫ e^(−r²) r dr dθ = π	Famous result

Key Takeaways

Polar coordinates simplify circular regions — what requires square roots in Cartesian becomes simple constant bounds
The Jacobian factor r is essential — it accounts for the stretching of area elements as we move away from the origin
Convert everything: the integrand, the limits, and the area element
The Gaussian integral $\int e^{-x^2} dx = \sqrt{\pi}$ is derived using polar coordinates — a technique appearing throughout ML and physics
Common curves like cardioids, roses, and limaçons have elegant polar descriptions

The Essence of Polar Integration:

"Circles become constants, and the Jacobian r reminds us that outer arcs are longer than inner ones."

Coming Next: In the next section, we'll explore Applications of Double Integrals, including computing volumes, mass, centroids, and moments of inertia for various 2D and 3D objects.