Green’s theorem

A closed curve in the plane bounds a region, and every vector field defined on that region secretly links two very different-looking quantities: how much the field circulates around the boundary, and how much it "curls" over the interior. Green's theorem is the precise statement of that link — it converts a line integral around a closed curve into a double integral over the enclosed region, and back again. It sits in the MA10002 vector calculus unit right after line integrals, as the flat, two-dimensional special case of the family completed later by Stokes's theorem (3D circulation) and the divergence theorem (3D flux). You are expected to state it precisely, verify it on a concrete curve and field by computing both sides independently, and use it as a computational shortcut — swapping a line integral for a double integral, or a double integral for a line integral, whichever side is easier. It also hands you a free formula for the area enclosed by any curve.

The core insight

Green's theorem says the circulation of a vector field around a positively oriented closed curve equals the double integral of its scalar curl over the region the curve encloses, ∮C(M dx+N dy)=∬R(Nx−My) dA\oint_C(M\,dx+N\,dy)=\iint_R(N_x-M_y)\,dA.

The core idea

Setup. Let CC be a positively oriented (counterclockwise), piecewise-smooth, simple closed curve in the plane, bounding a region RR. Let M(x,y)M(x,y) and N(x,y)N(x,y) have continuous partial derivatives on an open region containing RR.

The statement. Green's theorem asserts:

∮C(M dx+N dy)=∬R(∂N∂x−∂M∂y)dA.\oint_C (M\,dx + N\,dy) = \iint_R \left(\frac{\partial N}{\partial x} - \frac{\partial M}{\partial y}\right) dA.

where: CC is the boundary curve traversed counterclockwise, RR is the region it encloses, M,NM,N are the components of the vector field F⃗=(M,N)\vec F=(M,N), and ∂N/∂x−∂M/∂y\partial N/\partial x - \partial M/\partial y is the scalar (out-of-plane) curl of F⃗\vec F.

Reading it as circulation vs curl. The left side, ∮CF⃗⋅dr⃗\displaystyle\oint_C\vec F\cdot d\vec r, is the total circulation of F⃗\vec F around the boundary. The right side sums up the local rotation (curl) of F⃗\vec F at every point of RR. Green's theorem says these two totals must agree — a purely local property (curl at each point) determines a global property (net circulation), which is why the theorem is such a powerful shortcut.

Two directions of use. If the double integral is easier, use it to avoid parametrizing CC piece by piece. If the line integral is easier — or the region RR is awkward to set up — go the other way and reduce a double integral to a boundary line integral instead.

The area formula. Choosing M=−12y, N=12xM=-\tfrac12y,\ N=\tfrac12x makes Nx−My=1N_x-M_y=1, so the double integral becomes just the area of RR:

Area(R)=12∮C(x dy−y dx).\text{Area}(R) = \frac12\oint_C (x\,dy - y\,dx).

where: the right side is a line integral computed purely from the boundary curve — no need to describe the interior of RR at all. Equivalent single-term choices M=0,N=xM=0,N=x or M=−y,N=0M=-y,N=0 give Area(R)=∮Cx dy=−∮Cy dx\text{Area}(R)=\oint_C x\,dy=-\oint_C y\,dx.

Link to conservative fields. If F⃗=(M,N)\vec F=(M,N) is curl-free, Nx=MyN_x=M_y, then the right side of Green's theorem vanishes identically, so ∮CF⃗⋅dr⃗=0\displaystyle\oint_C\vec F\cdot d\vec r=0 for every closed curve in RR — exactly the zero-circulation property of conservative fields met earlier in the course, now derived as a special case.

Figure
Positively oriented boundary curve enclosing a region with the curl of the field distributed over the interior
Circulation on the boundary equals curl summed over the interior R ↻↻↻ ↻↻↻ ↻↻ local curl N_x - M_y at each point C, positively oriented (counterclockwise) — line integral of F around this boundary
Derivation AGreen's theorem as the 2D flux-of-curl identity

The goal is to show carefully, for a region RR that is both "vertically simple" (bounded above and below by graphs y=g2(x)≥y=g1(x)y=g_2(x)\ge y=g_1(x) over a≤x≤ba\le x\le b) and "horizontally simple," that the double integral ∬R(Nx−My) dA\iint_R(N_x-M_y)\,dA reduces exactly to the boundary line integral ∮C(M dx+N dy)\oint_C(M\,dx+N\,dy), by handling the MM term and the NN term separately with the fundamental theorem of calculus and matching signs to the counterclockwise orientation.

Step 1. Handle the MM part first: show ∬R(−∂M∂y)dA=∮CM dx\displaystyle\iint_R\left(-\frac{\partial M}{\partial y}\right)dA=\oint_C M\,dx. Treat RR as vertically simple, a≤x≤ba\le x\le b, g1(x)≤y≤g2(x)g_1(x)\le y\le g_2(x). Then ∬R(−∂M∂y)dA=−∫ab∫g1(x)g2(x)∂M∂y dy dx=−∫ab[M(x,g2(x))−M(x,g1(x))] dx\displaystyle\iint_R\left(-\frac{\partial M}{\partial y}\right)dA=-\int_a^b\int_{g_1(x)}^{g_2(x)}\frac{\partial M}{\partial y}\,dy\,dx=-\int_a^b\big[M(x,g_2(x))-M(x,g_1(x))\big]\,dx, using the fundamental theorem of calculus on the inner integral.

Step 2. Recognize the right side as a line integral. Split CC into the lower piece C1C_1 (y=g1(x)y=g_1(x), xx increasing from aa to bb) and the upper piece C2C_2 (y=g2(x)y=g_2(x), xx decreasing from bb to aa, since counterclockwise traversal goes along the top backwards). Then ∫abM(x,g1(x)) dx=∫C1M dx\displaystyle\int_a^b M(x,g_1(x))\,dx=\int_{C_1}M\,dx and −∫abM(x,g2(x)) dx=∫C2M dx\displaystyle-\int_a^b M(x,g_2(x))\,dx=\int_{C_2}M\,dx (the minus sign appears because C2C_2 runs from bb to aa). Adding these two pieces gives exactly ∮CM dx=−∬R∂M∂y dA\displaystyle\oint_C M\,dx=-\iint_R\frac{\partial M}{\partial y}\,dA.

Step 3. Repeat symmetrically for the NN part, treating RR as horizontally simple (c≤y≤dc\le y\le d, h1(y)≤x≤h2(y)h_1(y)\le x\le h_2(y)), to obtain ∮CN dy=∬R∂N∂x dA\displaystyle\oint_C N\,dy=\iint_R\frac{\partial N}{\partial x}\,dA by the same fundamental-theorem argument along xx instead of yy.

Step 4. Add the two results from Step 2 and Step 3. ∮CM dx+∮CN dy=−∬R∂M∂y dA+∬R∂N∂x dA\displaystyle\oint_C M\,dx+\oint_C N\,dy=-\iint_R\frac{\partial M}{\partial y}\,dA+\iint_R\frac{\partial N}{\partial x}\,dA, which rearranges to ∮C(M dx+N dy)=∬R(∂N∂x−∂M∂y)dA\displaystyle\oint_C(M\,dx+N\,dy)=\iint_R\left(\frac{\partial N}{\partial x}-\frac{\partial M}{\partial y}\right)dA — precisely Green's theorem. General regions are handled by cutting them into finitely many pieces that are each both vertically and horizontally simple; the interior cuts cancel in pairs because each interior boundary is traversed once in each direction.

Worked example

Problem. Verify Green's theorem for F⃗=(M,N)=(xy, x2)\vec F=(M,N)=(xy,\ x^2) on the region RR bounded by y=x2y=x^2 and y=xy=x between x=0x=0 and x=1x=1, with CC the positively oriented boundary.

Step 1. Double-integral side. Compute Nx−My=2x−x=xN_x-M_y=2x-x=x. On RR, for 0≤x≤10\le x\le1, yy runs from y=x2y=x^2 (lower) to y=xy=x (upper). So ∬Rx dA=∫01∫x2xx dy dx=∫01x(x−x2) dx=∫01(x2−x3) dx=[x33−x44]01=13−14=112\displaystyle\iint_R x\,dA=\int_0^1\int_{x^2}^{x}x\,dy\,dx=\int_0^1 x(x-x^2)\,dx=\int_0^1(x^2-x^3)\,dx=\left[\frac{x^3}3-\frac{x^4}4\right]_0^1=\frac13-\frac14=\frac1{12}.

Step 2. Line-integral side, piece C1C_1: y=x2y=x^2, x:0→1x:0\to1, dy=2x dxdy=2x\,dx. ∫C1(xy dx+x2dy)=∫01[x⋅x2+x2⋅2x]dx=∫013x3 dx=34\displaystyle\int_{C_1}(xy\,dx+x^2dy)=\int_0^1\big[x\cdot x^2+x^2\cdot2x\big]dx=\int_0^1 3x^3\,dx=\frac34.

Step 3. Piece C2C_2: y=xy=x, traversed back from (1,1)(1,1) to (0,0)(0,0), so x:1→0x:1\to0, dy=dxdy=dx. ∫C2(xy dx+x2dy)=∫10[x2+x2]dx=∫102x2 dx=2[x33]10=−23\displaystyle\int_{C_2}(xy\,dx+x^2dy)=\int_1^0\big[x^2+x^2\big]dx=\int_1^0 2x^2\,dx=2\left[\frac{x^3}3\right]_1^0=-\frac23.

Step 4. Sum the two pieces: ∮CF⃗⋅dr⃗=34−23=912−812=112\displaystyle\oint_C\vec F\cdot d\vec r=\frac34-\frac23=\frac9{12}-\frac8{12}=\frac1{12}. Both sides equal 112\dfrac1{12}, confirming Green's theorem on this curve and field.

Competitive stretch

Applying the area formula to the ellipse x2/a2+y2/b2=1x^2/a^2+y^2/b^2=1, parametrized by x=acos⁡t, y=bsin⁡tx=a\cos t,\ y=b\sin t for t∈[0,2π]t\in[0,2\pi], gives x dy−y dx=(acos⁡t)(bcos⁡t dt)−(bsin⁡t)(−asin⁡t dt)=ab dtx\,dy-y\,dx=(a\cos t)(b\cos t\,dt)-(b\sin t)(-a\sin t\,dt)=ab\,dt, so Area=12∫02πab dt=πab\text{Area}=\frac12\int_0^{2\pi}ab\,dt=\pi ab — the familiar ellipse area, obtained entirely from a boundary line integral, with no need to set up a double integral over the interior at all.

Common traps

  • Using the wrong orientation: Green's theorem requires CC positively oriented (counterclockwise, region on the left as you walk along CC); a clockwise traversal flips the sign of the answer.
  • Forgetting that M,NM,N must be continuously differentiable throughout the closed region RR, not just on CC — a field with a singularity inside RR (such as 1/(x2+y2)1/(x^2+y^2) at the origin) invalidates a direct application.
  • Sign slip in Nx−MyN_x-M_y: it is NxN_x minus MyM_y, not MyM_y minus NxN_x; reversing the order flips the sign of the whole double integral.
  • Setting up the double integral with the wrong bounds when the two curves bounding RR cross more than once — always sketch RR first to identify the correct upper and lower (or left and right) functions, as in the worked example above.

Ask the tutor: Want to see Green's theorem verified on a curve with three straight sides instead of a parabola and a line, or work through what changes when the region has a hole in it?

Why does Green's theorem need C to be positively oriented?

The proof builds the double integral by splitting CC into a lower piece traversed left-to-right and an upper piece traversed right-to-left (see the derivation above), and the signs from the fundamental theorem of calculus only line up correctly for that counterclockwise traversal. Reverse the direction and every sign in the derivation flips, so the line integral picks up an overall minus sign relative to the double integral. Practically: always check that you are walking along CC with the enclosed region on your left before applying the formula directly.

Can I use Green's theorem if the field has a singularity inside the region?

Not directly — the theorem requires MM and NN to have continuous partial derivatives on an open set containing the whole closed region RR, including the interior. If the field blows up at a point inside RR (a common example is a field like 1/(x2+y2)1/(x^2+y^2) at the origin), you have to exclude a small disk around the singularity, apply Green's theorem on the resulting punctured region, and handle the small circle around the singularity separately. Simply plugging into the formula as if nothing were wrong gives a wrong answer.

What is the fastest way to decide whether to compute the line integral or the double integral?

Look at both sides before committing. If Nx−MyN_x-M_y simplifies to something simple like a constant or a low-degree polynomial, and RR has clean bounds, the double integral is usually faster — this is exactly what happens in the ellipse area example. If instead the curve is simple to parametrize (a line, a circle, an ellipse) but the region's bounds would be awkward, or if Nx−MyN_x-M_y is messy, compute the line integral directly. There is no universal rule; a quick glance at Nx−MyN_x-M_y versus the shape of CC usually settles it in seconds.

Before this

Where this leads