Math graphic
📐 Concept diagram

07-06 — Green's Theorem

Phase: 7 — Calculus IV: Vector Calculus Subject: 07-06 Prerequisites: 07-05 — Line Integrals, double integrals (Phase 6) Next subject: 07-07 — Curl and Divergence


Learning Objectives

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

  1. State Green's theorem in both circulation-curl and flux-divergence forms
  2. Apply Green's theorem to convert line integrals over closed curves to double integrals, and vice versa
  3. Verify that Green's theorem holds for a given vector field and region
  4. Use Green's theorem to compute areas of planar regions as line integrals
  5. Recognize when Green's theorem applies (positively oriented, piecewise-smooth closed curve, simply-connected region, C¹ field)

Core Content

1. Statement of Green's Theorem

⚠️ CRITICAL FOUNDATION: Green's Theorem ∮_C P dx + Q dy = ∬_D (∂Q/∂x - ∂P/∂y) dA converts a line integral around a closed curve into a double integral over the enclosed region. It is the 2D precursor to Stokes' and Divergence theorems.

Let C be a positively oriented, piecewise-smooth, simple closed curve in the plane, and let D be the region bounded by C. Let F(x, y) = ⟨P(x, y), Q(x, y)⟩ be a vector field with continuous partial derivatives on an open region containing D. Then:

Circulation-Curl Form:

$∮_C P dx + Q dy = ∬_D (∂Q/∂x − ∂P/∂y) dA
$

Flux-Divergence Form (using outward normal n rather than tangent):

$∮_C F · n ds = ∬_D (∂P/∂x + ∂Q/∂y) dA = ∬_D div F dA
$

The first form is most common in calculus texts. The two forms are related by rotating the vector field 90° (replace ⟨P, Q⟩ by ⟨−Q, P⟩ or ⟨Q, −P⟩).

Orientation: "Positively oriented" means C is traversed counterclockwise, keeping the interior D on the left.

2. Circulation-Curl Form — Intuition and Examples

Green's theorem says: the circulation of F around the boundary C equals the integral of the curl (2D curl, ∂Q/∂x − ∂P/∂y) over the interior D.

Intuition: Circulation around the boundary is the sum of all the "microscopic circulation" (vorticity) within the region.

Example 1 — Verifying Green's theorem: Verify the theorem for F(x, y) = ⟨x², xy⟩ and C the boundary of the square [0, 1] × [0, 1].

Line integral (counterclockwise):
C₁ (bottom): y=0, x: 0→1. dy=0. ∫₀¹ x² dx = 1/3.
C₂ (right): x=1, y: 0→1. dx=0. ∫₀¹ 1·y dy = [y²/2]₀¹ = 1/2.
C₃ (top): y=1, x: 1→0. dy=0. ∫₁⁰ x² dx = −1/3.
C₄ (left): x=0, y: 1→0. dx=0. ∫₁⁰ 0·y dy = 0.

∮ F·dr = 1/3 + 1/2 − 1/3 + 0 = 1/2.

Green's theorem:
∂Q/∂x − ∂P/∂y = ∂(xy)/∂x − ∂(x²)/∂y = y − 0 = y.

∬_D y dA = ∫₀¹ ∫₀¹ y dy dx = ∫₀¹ [y²/2]₀¹ dx = ∫₀¹ 1/2 dx = 1/2. ✓

Example 2 — Using Green's to simplify: Evaluate ∮_C (3y + e^{√x}) dx + (5x + sin(y²)) dy where C is the circle x² + y² = 4.

Direct parametrization would be horrific (try integrating e^{√(2 cos t)}).

Use Green's theorem:
P = 3y + e^{√x}, Q = 5x + sin(y²).
∂Q/∂x = 5, ∂P/∂y = 3.
∂Q/∂x − ∂P/∂y = 2.

∮_C P dx + Q dy = ∬_D 2 dA = 2 · Area(D) = 2 · π(2)² = 8π.

The scary parts (e^{√x} and sin(y²)) disappeared! This is the power of Green's theorem.

3. Computing Area with Green's Theorem

Choose P and Q such that ∂Q/∂x − ∂P/∂y = 1. Then:

$Area(D) = ∬_D 1 dA = ∮_C P dx + Q dy
$

Three common choices:

  1. P = 0, Q = x: Area = ∮_C x dy
  2. P = −y, Q = 0: Area = ∮_C −y dx = −∮_C y dx
  3. P = −y/2, Q = x/2: Area = (1/2) ∮_C (x dy − y dx) (most symmetric)

Example 3: Use Green's theorem to find the area of the ellipse x²/a² + y²/b² = 1.

$Parametrize the ellipse: x = a cos t, y = b sin t, t ∈ [0, 2π].

Use form 3: Area = (1/2) ∮_C (x dy − y dx).

x dy − y dx = (a cos t)(b cos t dt) − (b sin t)(−a sin t dt)
            = ab cos²t dt + ab sin²t dt = ab(cos²t + sin²t) dt = ab dt.

Area = (1/2) ∫₀^{2π} ab dt = (1/2) · ab · 2π = πab.
$

This is the familiar formula! Derived cleanly via Green's theorem.

Example 4 — Area of a polygon: Use form 1 (Area = ∮_C x dy) to find the area of the triangle with vertices (0,0), (2,0), (1,3).

$Break into three segments, oriented counterclockwise:
C₁: (0,0)→(2,0). y=0 → dy=0. ∫_C₁ x dy = 0.
C₂: (2,0)→(1,3). Parametrize: r(t) = ⟨2−t, 3t⟩, t∈[0,1]. y=3t, dy=3dt. x=2−t.
    ∫₀¹ (2−t)·3 dt = ∫₀¹ (6−3t) dt = [6t − 3t²/2]₀¹ = 6 − 3/2 = 9/2.
C₃: (1,3)→(0,0). r(t) = ⟨1−t, 3(1−t)⟩, t∈[0,1]. x=1−t, y=3−3t, dy=−3dt.
    ∫₀¹ (1−t)(−3) dt = −3[t − t²/2]₀¹ = −3(1 − 1/2) = −3/2.

Total Area = 0 + 9/2 − 3/2 = 6/2 = 3.

Check: base=2, height=3. Area = (1/2)(2)(3) = 3. ✓
$

4. When Green's Theorem Applies (and When It Doesn't)

Requirements: 1. C is a simple closed curve (no self-intersections) 2. C is piecewise smooth 3. C is positively oriented (counterclockwise) 4. P and Q have continuous partial derivatives on an open region containing D 5. D is a simply-connected region (no holes)

What if the curve has holes? For a region with holes (like an annulus), Green's theorem generalizes:

$∮_{C_outer} − ∮_{C_inner} = ∬_D (∂Q/∂x − ∂P/∂y) dA
$

Each boundary is oriented so D is on the left (outer CCW, inner CW).

What if ∂Q/∂x − ∂P/∂y = 0 but the line integral isn't zero? This can happen if the region has a "hole" where the field isn't defined. Example:

F(x, y) = ⟨−y/(x²+y²), x/(x²+y²)⟩

∂Q/∂x = (y²−x²)/(x²+y²)², ∂P/∂y = (y²−x²)/(x²+y²)².
∂Q/∂x − ∂P/∂y = 0 (away from origin).

But ∮_unit circle F·dr = 2π ≠ 0.
Green's theorem can't be applied because F is not defined at the origin,
which is inside the circle. The region is not simply connected with respect to F.

5. Flux-Divergence Form

If F represents a fluid velocity field, the flux form of Green's theorem measures the net flow out of D:

$Flux = ∮_C F · n ds = ∬_D div F dA
$

Where n is the outward unit normal to C, and div F = ∂P/∂x + ∂Q/∂y (the divergence).

If C is positively oriented (CCW), the outward normal is n = ⟨y'/|r'|, −x'/|r'|⟩ (right-hand rule).

The relationship between the two forms: replace ⟨P, Q⟩ by ⟨−Q, P⟩ in the circulation form, or by ⟨Q, −P⟩.

6. Proof Sketch of Green's Theorem

Green's theorem can be proved by decomposing the region into Type I and Type II subregions.

For a region D that is both Type I and Type II:

$D: a ≤ x ≤ b, g₁(x) ≤ y ≤ g₂(x)  (Type I)
D: c ≤ y ≤ d, h₁(y) ≤ x ≤ h₂(y)  (Type II)
$

For Type I: compute ∮ Q dy along the boundary, reduce to ∬ ∂Q/∂x dA. For Type II: compute ∮ P dx along the boundary, reduce to −∬ ∂P/∂y dA.

Summing gives Green's theorem. For general regions, partition into simple subregions. Interior boundaries cancel.

Common Misconceptions

  1. "Green's theorem is just for making line integrals easier." It also works in reverse — to compute a difficult double integral by converting it to a line integral.

  2. "Green's theorem always works when ∂Q/∂x − ∂P/∂y = 0." No — if the field has a singularity inside the region, the theorem doesn't apply, and the line integral may be nonzero (as with ⟨−y, x⟩/(x²+y²)).

  3. "The orientation doesn't matter." Reversing C flips the sign of the line integral. The positive (CCW) orientation is used in the standard statement.

  4. "Green's theorem gives area only when ∂Q/∂x − ∂P/∂y = 1." Yes, that's one use, but the theorem is much more general. Any P, Q with continuous partials work.



Key Terms

Worked Examples

Example 1: Green's Theorem Simplifying a Hard Integral

Problem: Evaluate ∮_C (y + e^{cos x}) dx + (2x + ln(1 + y²)) dy where C is the circle x² + y² = 9.

Solution:

$P = y + e^{cos x}, Q = 2x + ln(1 + y²).
The non-polynomial parts are impossible to antiderivate.

Green's theorem:
∂Q/∂x = 2, ∂P/∂y = 1.
∂Q/∂x − ∂P/∂y = 2 − 1 = 1.

∮ = ∬_D 1 dA = Area(D) = π(3)² = 9π.
$

Example 2: Area of an Astroid

Problem: Find the area enclosed by the astroid x^{2/3} + y^{2/3} = a^{2/3}.

Solution:

$Parametrize: x = a cos³t, y = a sin³t, t ∈ [0, 2π].

Use Area = (1/2)∮ (x dy − y dx).

x dy = a cos³t · 3a sin²t cos t dt = 3a² cos⁴t sin²t dt
y dx = a sin³t · (−3a cos²t sin t) dt = −3a² sin⁴t cos²t dt

x dy − y dx = 3a² cos²t sin²t(cos²t + sin²t) dt = 3a² cos²t sin²t dt

∫₀^{2π} cos²t sin²t dt = (1/4)∫₀^{2π} sin²(2t) dt = (1/4)∫₀^{2π} (1−cos(4t))/2 dt
= (1/8)[t − sin(4t)/4]₀^{2π} = (1/8)(2π) = π/4.

Area = (1/2) · 3a² · (π/4) = 3πa²/8.
$

Example 3: Line Integral from Double Integral

Problem: Use Green's theorem to evaluate ∬_D e^{x²} dA over the triangle with vertices (0,0), (1,0), (0,1). (Hint: choose P and Q so that ∂Q/∂x − ∂P/∂y = e^{x²}.)

Solution:

$We need P, Q such that ∂Q/∂x − ∂P/∂y = e^{x²}.

Since e^{x²} doesn't depend on y, let P = 0 and ∂Q/∂x = e^{x²}. 
Then Q = ∫ e^{x²} dx — but this can't be expressed in elementary functions!

Alternative: Let Q = 0 and −∂P/∂y = e^{x²}. So ∂P/∂y = −e^{x²}.
Then P = −y e^{x²} (treating x constant).

So: ∬_D e^{x²} dA = ∮_C (−y e^{x²}) dx + 0 dy = −∮_C y e^{x²} dx.

Now C is the triangle boundary:
C₁: (0,0)→(1,0): y=0. −∫₀¹ 0·e^{x²} dx = 0.
C₂: (1,0)→(0,1): r(t)=⟨1−t,t⟩, t∈[0,1]. dx=−dt. y=t, x=1−t.
    −∫₀¹ t·e^{(1−t)²}·(−dt) = ∫₀¹ t e^{(1−t)²} dt.
C₃: (0,1)→(0,0): x=0. −∫₁⁰ y·e⁰·0 = 0.

So ∬_D e^{x²} dA = ∫₀¹ t e^{(1−t)²} dt.
Let u = 1−t, du = −dt. = ∫₁⁰ (1−u) e^{u²} (−du) = ∫₀¹ (1−u) e^{u²} du.
= ∫₀¹ e^{u²} du − ∫₀¹ u e^{u²} du.
= ∫₀¹ e^{u²} du − [e^{u²}/2]₀¹
= ∫₀¹ e^{u²} du − (e−1)/2.

Unfortunately, ∫₀¹ e^{u²} du is the imaginary error function — not elementary.
But either representation (the double integral or the line integral) is valid.
$

Quiz

Q1: Green's Theorem relates a line integral around a closed curve C to:

A) A triple integral over the solid bounded by C B) A double integral over the region D enclosed by C C) A surface integral over the surface bounded by C D) A line integral along any other curve

Correct: B)


Q2: To apply Green's Theorem, the curve C must be:

A) Any curve, open or closed B) A simple closed curve, oriented counterclockwise (positively) C) A straight line D) Parametrized by arc length

Correct: B)


Q3: Green's Theorem can be used to compute area by setting:

A) P = 0, Q = 0 B) P = −y/2, Q = x/2 C) P = 1, Q = 1 D) P = x, Q = y

Correct: B)


Q4: For F = ⟨−y, x⟩ around the unit circle (counterclockwise), ∮_C F · dr equals:

A) 0 B) π C) 2π D) 1

Correct: C)


Q5: If ∮_C F · dr = 0 for every closed curve C in a region, then in that region:

A) F = 0 everywhere B) ∂Q/∂x − ∂P/∂y = 0 everywhere (F is irrotational) C) F is constant D) The region must be simply connected

Correct: B)


Q6: Green's Theorem fails for a region with a hole because:

A) The integral is always zero B) The curve is not simple closed C) The region is not simply connected — you need to add cuts or use the extended version D) The theorem only works for circles

Correct: C)


Practice Problems

(Answers are below. Try each problem before checking.)

Problem 1: Verify Green's theorem for F(x, y) = ⟨x²y, xy²⟩ over the square [0, 1] × [0, 1].

Problem 2: Use Green's theorem to evaluate ∮_C (xy dx + x² dy) where C is the triangle with vertices (0,0), (1,0), (0,2).

Problem 3: Use Green's theorem to evaluate ∮_C (y² dx − x² dy) where C is the circle x² + y² = 1.

Problem 4: Find the area of one arch of the cycloid x = t − sin t, y = 1 − cos t for t ∈ [0, 2π] and the x-axis.

Problem 5: Use Green's theorem to evaluate ∬_D (x² + y²) dA where D is the region bounded by the cardioid r = 1 + cos θ. (Hint: find a line integral representation.)

Problem 6: Evaluate ∮_C ⟨−y³, x³⟩ · dr where C is the circle x² + y² = 4.

Problem 7: Use the area formula (1/2)∮(x dy − y dx) to find the area inside the ellipse x²/16 + y²/9 = 1.

Answers (click to expand) **Problem 1:** Line integral (CCW): C₁ (bottom, y=0): ∫₀¹ 0 dx = 0. C₂ (right, x=1): ∫₀¹ 1·y² dy = [y³/3]₀¹ = 1/3. C₃ (top, y=1): ∫₁⁰ x²·1 dx = −[x³/3]₁⁰ = −1/3. C₄ (left, x=0): ∫₁⁰ 0 dy = 0. ∮ = 1/3 − 1/3 = 0. Green's: ∂Q/∂x − ∂P/∂y = ∂(xy²)/∂x − ∂(x²y)/∂y = y² − x². ∬_D (y² − x²) dA = ∫₀¹∫₀¹ (y²−x²) dy dx = ∫₀¹ [y³/3 − x²y]₀¹ dx = ∫₀¹ (1/3 − x²) dx = [x/3 − x³/3]₀¹ = 1/3 − 1/3 = 0. ✓ **Problem 2:** P = xy, Q = x². ∂Q/∂x − ∂P/∂y = 2x − x = x. ∬_D x dA over triangle: x ∈ [0,1], y ∈ [0, 2(1−x)]. ∬ x dA = ∫₀¹ ∫₀^{2(1−x)} x dy dx = ∫₀¹ 2x(1−x) dx = 2[x²/2 − x³/3]₀¹ = 2(1/2 − 1/3) = 2(1/6) = 1/3. **Problem 3:** P = y², Q = −x². ∂Q/∂x − ∂P/∂y = −2x − 2y. ∮ = ∬_D (−2x − 2y) dA over unit disk. By symmetry: ∬ x dA = 0, ∬ y dA = 0. ∮ = 0. **Problem 4:** The cycloid arch: closed curve C = arch + line segment x-axis from 2π back to 0. Area = −∮_C y dx (use the −y dx formula). Along the x-axis: y=0, so ∫ 0 dx = 0. Along the arch: x = t−sin t, y = 1−cos t, t: 0→2π. dx = (1−cos t) dt. ∫_{arch} y dx = ∫₀^{2π} (1−cos t)(1−cos t) dt = ∫₀^{2π} (1−cos t)² dt. = ∫₀^{2π} (1 − 2cos t + cos²t) dt = [t − 2sin t]₀^{2π} + ∫₀^{2π} (1+cos 2t)/2 dt = 2π + [t/2 + sin(2t)/4]₀^{2π} = 2π + π = 3π. But the boundary orientation needs care: the standard CCW orientation for the area formula traverses the arch from (2π,0) to (0,0), then the x-axis back. Let's use the (1/2)∮(x dy − y dx): Param for arch: t: 0→2π, but orientation is from (2π,0) to (0,0) along the arch. Better: area = ∮_C x dy (traverse CCW). Arch: x = t−sin t, dy = sin t dt. ∫₀^{2π} (t−sin t) sin t dt. = ∫₀^{2π} (t sin t − sin²t) dt = [−t cos t + sin t]₀^{2π} − π = (−2π·1+0) − (−0+0) − π = −2π − π = −3π. Then x-axis back: dy=0, so contributes 0. Total ∮ x dy = −3π. But since we traversed CW? Let me use (1/2)∮(x dy − y dx) more carefully. The cycloid from (0,0) to (2π,0) makes a CCW closed curve with the x-axis from (2π,0) back to (0,0). Along the x-axis segment (reversed): x goes 2π→0, y=0, dy=0. x dy = 0, y dx = 0. Along arch (t: 0→2π): x = t−sin t, y = 1−cos t. dx = (1−cos t)dt, dy = sin t dt. x dy − y dx = (t−sin t) sin t − (1−cos t)(1−cos t) = t sin t − sin²t − (1 − 2cos t + cos²t) = t sin t − sin²t − 1 + 2cos t − cos²t = t sin t − (sin²t+cos²t) − 1 + 2cos t = t sin t − 1 − 1 + 2cos t = t sin t + 2cos t − 2. ∫₀^{2π} (t sin t + 2cos t − 2) dt = [−t cos t + sin t + 2sin t − 2t]₀^{2π} = (−2π·1 + 0 + 0 − 4π) − (0) = −6π. (1/2)(−6π) = −3π. Negative because orientation is CW relative to the usual formula. The absolute area is 3π. ✓ (Classic result: area under one arch of cycloid = 3πr² = 3π.) **Problem 5:** Choose P, Q such that ∂Q/∂x − ∂P/∂y = x² + y². Example: P = −y³/3, Q = x³/3. Then ∂Q/∂x − ∂P/∂y = x² − (−y²) = x² + y². Perfect. ∬_D (x²+y²) dA = ∮_C (−y³/3 dx + x³/3 dy). Cardioid: r = 1+cos θ. x = (1+cos θ)cos θ, y = (1+cos θ)sin θ, θ: 0→2π. dx = (−sin θ cos θ − (1+cos θ)sin θ) dθ... this is tedious in Cartesian. Better to compute the double integral directly in polar: ∬_D r² · r dr dθ = ∫₀^{2π} ∫₀^{1+cos θ} r³ dr dθ = ∫₀^{2π} [r⁴/4]₀^{1+cos θ} dθ = (1/4)∫₀^{2π} (1+cos θ)⁴ dθ. Expand: (1+cos θ)⁴ = 1 + 4cos θ + 6cos²θ + 4cos³θ + cos⁴θ. ∫₀^{2π} cos θ dθ = 0. ∫₀^{2π} cos³θ dθ = 0. ∫₀^{2π} cos²θ dθ = π. ∫₀^{2π} cos⁴θ dθ = 3π/4. (1/4)[2π + 6π + 0 + 0 + 3π/4] = (1/4)[8π + 3π/4] = (1/4)(35π/4) = 35π/16. **Problem 6:** F = ⟨−y³, x³⟩. P = −y³, Q = x³. ∂Q/∂x − ∂P/∂y = 3x² − (−3y²) = 3(x² + y²). Disk D: x²+y² ≤ 4. ∮ = ∬_D 3(x²+y²) dA = 3∫₀^{2π} ∫₀² r² · r dr dθ = 3·2π·[r⁴/4]₀² = 6π·4 = 24π. **Problem 7:** x = 4 cos t, y = 3 sin t, t ∈ [0, 2π]. x dy − y dx = (4 cos t)(3 cos t) − (3 sin t)(−4 sin t) = 12 cos²t + 12 sin²t = 12. Area = (1/2)∫₀^{2π} 12 dt = 6·2π = 12π. Check: πab = π·4·3 = 12π. ✓

Summary

  1. Green's theorem ∮_C P dx + Q dy = ∬_D (∂Q/∂x − ∂P/∂y) dA connects the circulation around a closed curve to the integral of the curl over the enclosed region — it's the 2D analog of the Fundamental Theorem of Calculus
  2. The theorem transforms difficult line integrals into simple double integrals (especially when ∂Q/∂x − ∂P/∂y is constant) and can also work in reverse to compute double integrals via line integrals
  3. Three area formulas emerge as special cases: Area = ∮ x dy = −∮ y dx = (1/2)∮ (x dy − y dx); the last form is particularly useful for parametric curves
  4. Green's theorem requires P and Q to have continuous partial derivatives on an open region containing D; if the field has a singularity inside D, the theorem does not apply
  5. The flux-divergence form ∮ F · n ds = ∬ div F dA measures net outward flow and is the 2D precursor to the Divergence Theorem

Pitfalls

  1. Applying Green's theorem when the field has a singularity inside the region. If P or Q is undefined at any point within D, the theorem's hypotheses are violated. The classic case is ⟨−y, x⟩/(x² + y²) around the unit circle — curl is zero everywhere except at the origin, yet the line integral is 2π, not 0.

  2. Using the wrong orientation. Green's theorem in its standard form requires counterclockwise (positive) orientation. If the curve is traversed clockwise, the result is the negative. Always check orientation before applying the formula.

  3. Forgetting to verify that ∂Q/∂x − ∂P/∂y simplifies before integrating. The power of Green's theorem comes when the curl simplifies to a constant or simple function. Students sometimes set up the double integral without first computing and simplifying the integrand.

  4. Confusing the circulation-curl and flux-divergence forms. The circulation form uses ∂Q/∂x − ∂P/∂y (the 2D curl); the flux form uses ∂P/∂x + ∂Q/∂y (the 2D divergence). Using one when the problem calls for the other will produce a wrong answer.

  5. Misapplying the area formulas. The formulas Area = (1/2)∮(x dy − y dx), Area = ∮ x dy, and Area = −∮ y dy work for any simple closed curve, but only when the integrand is correctly paired with the differential. Using, say, ∮ x dx — which always gives zero around a closed curve — is a common mistake.



Next Steps

Move on to 07-07 — Curl and Divergence to explore the physical meaning and computation of these two fundamental differential operators in 3D, and to understand how they characterize vector fields.