Skip to main content

Section 5.7 Green's theorem

Step 6: we study what happens if we integrate an exact two-form. We start by studying integrals of exact two-forms over bounded regions in R2: this gives rise to Green's theorem.

Subsection 5.7.1 Green's theorem

Let us start by recalling the corresponding statement for integrals of exact one-forms over intervals in R. Let f be a zero-form on UR. Let [a,b]U with the canonical orientation of increasing numbers. The boundary of the interval, with its induced orientation, is [a,b]={(b,+),(a,)}. Then the integral of the exact one-form ω=df along [a,b] can be evaluated as:

[a,b]df=[a,b]f=f(b)f(a).

This is nothing else but the Fundamental Theorem of Calculus, rewritten in a fancy way using one-forms and zero-forms.

We now want to do something similar for exact two-forms over closed bounded regions in R2.

We will only prove Green's theorem for recursively supported regions. (In fact, we will only write the proof for x-supported regions, but the proof for y-supported regions is analogous.) For more general closed bounded regions, one can prove Green's theorem by rewriting the region as an union of recursively supported regions.

We write η=f dx+g dy and dη=(gxfy)dxdy.

Let D be an x-supported region of the form:

D={(x,y)R2 | x[a,b],u(x)yv(x)},

for some continuous functions u,v:RR. The boundary D is a closed simple curve, with counterclockwise orientation. It can be split into four curves:

  1. C1 is the vertical line x=b between y=u(b) and y=v(b);

  2. C2 is the curve y=v(x) between x=b and x=a;

  3. C3 is the vertical line x=a between y=v(a) and y=u(a);

  4. C4 is the curve y=u(x) between x=a and x=b.

In fact, for C2 and C3, we will pick the opposite orientations (as it simplifies the parametrizations), and add a negative sign in front of the line integrals. Those four curves can be realized as parametric curves (C2 and C3 have opposite orientations):

  1. α1:[0,1]R2 with α1(t)=(b,u(b)+(v(b)u(b))t);

  2. α2:[a,b]R2 with α2(t)=(t,v(t));

  3. α3:[0,1]R2 with α3(t)=(a,u(a)+(v(a)u(a))t);

  4. α4:[a,b]R2 with α4(t)=(t,u(t));

Calculating the pullbacks, we evaluate the line integrals:

α1η=01g(b,u(b)+(v(b)u(b))t)(v(b)u(b)) dt,α2η=ab(f(t,v(t))+g(t,v(t))v(t)) dt,α3η=01g(a,u(a)+(v(a)u(a))t)(v(a)u(a)) dt,α4η=ab(f(t,u(t))+g(t,u(t))u(t)) dt.

Putting those together, remembering that we need to add a negative sign for C2 and C3, we get:

Dη=α1ηα2ηα3η+α4η=01(g(b,u(b)+(v(b)u(b))t)(v(b)u(b))g(a,u(a)+(v(a)u(a))t)(v(a)u(a)))dt+ab(f(t,u(t))f(t,v(t))+g(t,u(t))u(t)g(t,v(t))v(t)) dt.

Next, we need to evaluate the integral Ddη. Here we will do our favourite trick: the pullback. Instead of using the x-supported domain D directly, we will pullback to a rectangular domain. This will allow us to integrate the two-form. Consider the function ϕ:D2D, where

D2={(s,t)R2 | s[a,b],t[0,1]},

with

ϕ(s,t)=(s,u(s)+(v(s)u(s))t).

We see that ϕ(D2)=D, and that the function is bijective. So by Lemma 5.3.7, we know that

D2ϕ(dη)=D1dη,

so we may as well calculate the left-hand-side. In fact, we also know that ϕ(dη)=d(ϕη). So let us calculate this pullback. First,

ϕη=f(ϕ(s,t)) ds+g(ϕ(s,t))(u(s)+(v(s)u(s))t) ds+g(ϕ(s,t))(v(s)u(s)) dt.

Then

d(ϕη)=(f(ϕ(s,t))tt(g(ϕ(s,t))(u(s)+(v(s)u(s))t))+s(g(ϕ(s,t))(v(s)u(s))))dsdt.

Thus

D2ϕ(dη)=ab01(f(ϕ(s,t))tt(g(ϕ(s,t))(u(s)+(v(s)u(s))t)))dtds+01ab(s(g(ϕ(s,t))(v(s)u(s))))dsdt,

where we used Fubini's theorem to exchange the order of two integrals in the second line. We can then evaluate the inner integrals using the Fundamental Theorem of Calculus, since they are definite integrals of derivatives. We get:

D2ϕ(dη)=ab(f(ϕ(s,0))f(ϕ(s,1))+g(ϕ(s,0))u(s)g(ϕ(s,1))v(s))ds+01(g(ϕ(b,t))(v(b)u(b))g(ϕ(a,t))(v(a)u(a))) dt.

Substituting back the expression for ϕ(s,t), we get:

D2ϕ(dη)=ab(f(s,u(s))f(s,v(s))+g(s,u(s))u(s)g(s,v(s))v(s)) ds+01(g(b,u(b)+(v(b)u(b))t)(v(b)u(b))g(a,u(a)+(v(a)u(a))t)(v(a)u(a))) dt.

Magic: this is exactly the same as the expression that we found many lines above for the line integral Dη! Woot woot! Therefore

Ddη=Dη.

Remark 5.7.2.

If D is simply connected, then D is a simple closed curve, and D is the closed curve with its interior. In this case the line integral Dη is simply the line integral of the one-form η along the simple closed curve with counterclockwise orientation.

However, in the statement of Green's theorem, D can be more general. For instance, it could be an annulus. Then D would have two components: the inner and outer circle. The line integral Dη would then be the sum of the line integrals of the one-form η along the two circles; the outer circle with counterclockwise orientation, and the inner circle with clockwise orientation, as this is the induced orientation on the boundary according to Definition 5.2.9.

Remark 5.7.3.

Following up on the previous remark, we note that we can read Green's theorem in two different ways, which is generally the case for all integral theorems of vector calculus. We could say:

  1. The integral of the exact two-form dη over the region DR2 with canonical orientation is equal to the integral of the one-form η over the boundary D with the induced orientation.

  2. The integral of the one-form η over the simple closed curve DR2 with canonical orientation is equal to the integral of its exterior derivative dη over the interior D with canonical orientation.

This is just two different readings of the same statement, depending on whether you start on the left-hand-side or the right-hand-side. Consequently, we can use Green's theorem either to evaluate integrals of exact two-forms via reading 1, or to evaluate line integrals over simple closed curves via reading 2.

In practice however, Green's theorem is generally useful mostly to evaluate line integrals by transforming them into surface integrals.

Find the line integral of the one-form ω=xy dx+(x+y) dy over the rectangle with vertices (2,1),(2,1),(2,0),(2,0), with a clockwise orientation.

We could evaluate the line integral using previous techniques, by rewriting the curve as four parametric curves for each line segment and then use the definition of line integrals. But let us instead use Green's theorem.

Let us denote the rectangular curve by C. It is the boundary of the region D consisting of the interior of the rectangle with its boundary:

D={(x,y)R2 | x[2,2],y[1,0]}.

Thus we can use Green's theorem to rewrite the line integral along C as a surface integral over D.

We have to be careful with orientation though. If D has canonical orientation, the induced orientation on the boundary rectangle C will be counterclockwise. But the question is asking us to evaluate the line integral along C with clockwise orientation. So in Green's theorem, if we relate the line integral to the surface integral over D with canonical orientation, we will need to introduce a minus sign. More precisely, let us denote C to be the rectangle with clockwise orientation, C+ the rectangle with counterclockwise orientation, and D the rectangular region with canonical (counterclockwise) orientation. Green's theorem states:

Cω=C+ω=Ddω.

So instead of calculating the line integral in the problem, we can calculate Ddω.

The two-form dω is:

dω=xdydx+dxdy=(1x)dxdy.

The integral can be evaluated:

Ddω=1022(1x) dxdy=10[xx22]22 dy=10(2(2)2+2) dy=4(0(1))=4.

Therefore, the line integral of ω along C is

Cω=Ddω=4.

It is a good exercise to check that this is the correct answer by evaluating the line integral using the standard approach with parametric curves.

While Green's theorem is mostly useful to evaluate line integrals, we can also use Green's theorem in the reverse direction, to evaluate integrals of exact two-forms. Here is an example.

Find the area enclosed by the ellipse

x2a2+y2b2=1.

The area is given by integrating the basic two-form ω=dxdy over the region D bounded by the ellipse with canonical orientation. We could calculate this integral directly, as a double integral. Or we can use Green's theorem, if ω is exact. But it certainly is: for instance, we can write ω=dη with η=12(x dyy dx) (we could choose other one-forms such that dη=ω as well, but this one gives a nice and easy integral). Then Green's theorem states:

Dω=Cη,

where on the right-hand-side we are evaluating the line integral of η=12(x dyy dx) along the ellipse with counterclockwise orientation.

We can parametrize the ellipse as α:[0,2π]R2 with α(t)=(acos(t),bsin(t)). The orientation is counterclockwise as required. The pullback αη is:

αη=(ab2cos2(t)ab2sin(t)(sin(t))) dt=ab2dt.

The line integral is then

Cη=02πab2 dt=πab,

which is the area of the ellipse.

Subsection 5.7.2 Vector form of Green's theorem

As always, we can translate our results to vector calculus concepts using our dictionary. It is a bit artificial though here, since Green's theorem really lives in two dimensions, while our vector calculus concepts live in three dimensions. But the idea is to think of the region D as being the trivial parametric surface in R3 where we embed it trivially in the xy-plane. That is, α:DR3, with α(x,y)=(x,y,0). The tangent vectors are Tx=(1,0,0) and Ty=(0,1,0), and the normal vector is of course n=(0,0,1)=e3.

To the one-form η=f dx+g dy, we associate the vector field F=(f,g,0). To the two-form ω=dη is then associated the vector field ×F. Using Corollary 5.6.5, we then see that we can rewrite the surface integral as

Ddη=D((×F)e3) dA.

In other words, the integrand is the z-component of the curl ×F.

As for the other side of Green's theorem, it is a line integral of the one-form η over the curve D. Assuming that we parametrize the curve with a position function r, we can rewrite the line integral as

Dη=DFdr.

The result is the vector form of Green's theorem:

D((×F)e3) dA=DFdr.

Exercises 5.7.3 Exercises

1.

Evaluate the line integral of the one-form

ω=ey dx+ex dy

along the rectangle with vertices (0,0),(3,0),(3,2),(0,2), counterclockwise, using two different methods: (a) directly, (b) using Green's theorem.

Solution.

(a) Let us call D the path along the rectangle counterclockwise. It is shown in the figure below.

The rectangle.
Figure 5.7.6. The path D along the rectangle.
To evaluate the line integral directly, we need to parametrize separately the four line segments. We write the parametrizations as αi:[0,1]R, i=1,,4, with

α1(t)=(3t,0),α2(t)=(3,2t),α3(t)=(33t,2),α4(t)=(0,22t).

We calculate the pullbacks:

α1ω=3 dt,α2ω=2e3 dt,α3ω=3e2 dt,α4ω=2 dt.

Putting all this together, the line integral becomes

Dω=α1ω+α2ω+α3ω+α4ω=01(3+2e33e22) dt=1+2e33e2.

(b) We now use Green's theorem to calculate the line integral. By Green's theorem, we know that

Dω=Ddω,

where D is the interior of the rectangle with its boundary, that is,

D={(x,y)R2 | x[0,3],y[0,2]},

with canonical orientation. We calculate the exterior derivative:

dω=ey dydx+ex dxdy=(exey) dxdy.

The integral of the two-form is:

Ddω=D(exey) dxdy=0203(exey)dxdy=02(e313ey) dy=2e323e2+3=1+2e33e2.

Therefore,

Dω=1+2e33e2,

which is the same answer as in part (a), as it should.

2.

Use Green's theorem to evaluate the line integral of the vector field F(x,y)=(x1/3+y,x4/3y1/3) along the triangle with vertices (1,3),(2,2),(1,1), clockwise.

Solution.

Let us call D the path around the triangle clockwise. It is shown in the figure below.

The triangle.
Figure 5.7.7. The path D along the triangle.
We associate to the vector field F the one-form ω=(x1/3+y) dx+(x4/3y1/3) dy. Since the path D is oriented clockwise (negative orientation), by Green's theorem we know that

Dω=Ddω,

where D is the interior of the triangle with its boundary, with canonical orientation. We can write equations for the three sides of the triangle. The blue line has equation x=1; the green line, y=x; the orange line, y=x+4. Using these equations, we can describe D as an x-supported region:

D={(x,y)R2 | x[1,2],xyx+4}.

We calculate the exterior derivative dω:

dω=dydx+43x1/3dxdy=(43x1/31)dxdy.

Its integral along D is:

Ddω=12xx+4(43x1/31) dydx=12(43x1/31)[y]xx+4 dx=12(43x1/31)(x+4x) dx=12(83x4/3+2x+163x1/34) dx=(87x7/3+x2+4x4/34x)x=1x=2=32721/3+4+821/38+8714+4=37(821/39).

Therefore,

Dω=Ddω=37(821/39).

3.

Use Green's theorem to find the work done by the force

F(x,y)=(x2+y2,xy)

while moving an object first along the parabola y=x2 from the origin to the point (2,4), then along a horizontal line back to the y-axis, and then back to the origin along a vertical line.

Solution.

We denote the closed path by D. It is shown in the figure below.

The path.
Figure 5.7.8. The closed path D.
From the description of the path we know that the object is moving counterclockwise along the path. Therefore, if we associate the one-form ω=(x2+y2) dx+xy dy to the vector field F, by Green's theorem we know that

Dω=Ddω,

with D being the region enclosed by the path (with its boundary), with canonical orientation. The region D can be described as an x-supported region:

D={(x,y)R2 | x[0,2],x2y4}.

The exterior derivative dω is:

dω=2y dydx+y dxdy=y dxdy

Its integral along D is:

Ddω=Dy dxdy=02x24y dydx=02[y22]y=x2y=4 dx=02(8x42) dx=(16165)=645.

Therefore, by Green's theorem, the work done by the force F while moving an object along the path D is:

Dω=645.

4.

Suppose that a polygon has vertices (x1,y1),(x2,y2),,(xn,yn), in counterclockwise order. The area A of the polygon is given by integrating the basic two-form ω=dxdy over the polygon. Use Green's theorem to show that

A=12((x1y2x2y1)+(x2y3x3y2)++(xn1ynxnyn1)+(xny1x1yn)).

This is a well known formula for the area of an arbitrary polygon, see for instance https://en.wikipedia.org/wiki/Polygon#Area.

Solution.

We denote by D the closed path that goes around the polygon counterclockwise. It consists of n line segments Li, i=1,,n, with Li being the line segment between the points (xi,yi) and (xi+1,yi+1) (for simplicity of notation, we define (xn+1,yn+1):=(x1,y1)).

Further, we notice that the two ω=dxdy is exact. We write it as ω=dη for the one-form

η=x dy.

Note that this is certainly not the only choice of η, but any η such that dη=ω will work.

Then, by Green's theorem, we know that the area A of the polygon is given by

A=Ddxdy=Dη=i=1nαiη,

where in the last line αi is a parametrization of the line segment Li. More explicitly, since the line segment Li joins the points (xi,yi) and (xi+1,yi+1), we can parametrize it as αi:[0,1]R2 with

αi(t)=(xi+t(xi+1xi),yi+t(yi+1yi)).

Then

αiη=(xi+t(xi+1xi))(yi+1yi) dt.

The line integral along αi is then:

αiη=01(xi+t(xi+1xi))(yi+1yi) dt=(yi+1yi)(xi+12(xi+1xi))=12(yi+1yi)(xi+xi+1)=12(xiyi+1xi+1yi+xi+1yi+1xiyi).

Finally, we add up these line integrals to get the area of the polygon. We get:

A=i=1nαiη=12i=1n(xiyi+1xi+1yi+xi+1yi+1xiyi)=12((x1y2x2y1+x2y2x1y1)+(x2y3x3y2+x3y3x2y2)++(xny1x1yn+x1y1xnyn)),

where in the last line we used the fact that xn+1=x1 and yn+1=y1. We see that the terms x1y1,x2y2, all cancel out, and we are left with

A=12((x1y2x2y1)+(x2y3x3y2)++(xny1x1yn)).

Bingo!

5.

We already studied the one-form

η=yx2+y2 dx+xx2+y2 dy,

defined on U=R2{(0,0)}, since it is the typical example of a one-form that is closed but not exact. For instance, we proved in Exercise 3.4.3.2 that it is not exact by showing that its line integral along a circle of radius one around the origin is non-vanishing. In this problem we show that the line integral of η is in fact non-vanishing for every simple closed curve that encloses the origin, and vanishing for every simple closed curve that does not pass through or enclose the origin.

  1. Consider an arbitrary simple closed curve C0 with canonical orientation that does not pass through or enclose the origin. Use Green's theorem to show that

    C0η=0.

  2. Let C1 be an arbitrary simple closed curve with canonical orientation that encloses the origin. Explain why the argument of (a) does not apply here. So we need to do something else to study the line integral of ω along C1.

  3. As in (b), let C1 be an arbitrary simple closed curve with canonical orientation that encloses the origin. Suppose that K is a circle centered at the origin, with a radius small enough that K lies completely inside C1. Give K a counterclockwise orientation. Use Green's theorem to show that

    C1η=Kη.

  4. Using part (c), show that it implies that

    C1η=2π.

You have showed that the line integral of η along an arbitrary simple closed curve that encloses the origin is non-vanishing (and is in fact equal to 2π), while the line integral along an arbitrary simple closed curve that does not pass through or enclose the origin is vanishing. Isn't it neat?

Solution.

The key in this problem is to be very careful with the domain of definition of the one-form η. The one-form is define on the open set U=R2{(0,0)}.

As a starting point, recall that the one-form η is closed, that is, dω=0. Indeed,

dη=(x2+y22y2(x2+y2)2)dydx+(x2+y22x2(x2+y2)2)dxdy=(x2y2)+(y2x2)(x2+y2)2 dxdy=0.

This will be useful in this problem.

(a) We assume that C0 is a simple closed curve that does not pass through or enclose the origin, with canonical orientation. Therefore, C0U. Moreover, if we denote by D0 the region consisting of the closed curve C0 and its interior, then D0U as well, since the origin is not in D0. Therefore, by Green's theorem,

C0η=D0dη=0,

since dη=0.

(b) If C1 is a simple closed curve that encloses the origin, we cannot apply the argument in (a). The reason is that, if we denote by D1 the region consistin of the closed curve C1 and its interior, then D1 is not a subset of U, since D1 contains the origin. In other words, the one-form η is not defined on D1, which is one of the assumptions in the statement of Green's theorem. Therefore, Green's theorem does not apply. (Indeed, as we will see, if Green's theorem applied, it would imply that the line integral of η along C1 vanishes, which is not true, as we will see.)

(c) We suppose that C1 is an arbitrary simple closed curve with canonical orientation that encloses the origin, and that K is a circle centered at the origin that lies completely inside C1, also oriented counterclockwise. For instance, C1 could be the orange curve in the figure below, and K the blue circle.

An example.
Figure 5.7.9. An example of curves C1 and K.
We let D1 be the region bounded by the two closed curves C1 and K, and we give D1 a canonical orientation. By definition, the origin is not in D1, so D1U and the one-form is well defined on D1. We can then apply Green's theorem. The boundary of D1 has two components, D1=C1K, and if D1 is oriented counterclockwise, then the induced orientation on the boundary is counterclockwise for the outer boundary C1 and clockwise for the inner boundary K. To clarify the notation, we write (C1)+ for C1 with counterclockwise orientation, K+ for K with counterclockwise orientation, and K for K with clockwise orientation. So we can write the oriented boundary of D1 as D1=(C1)+K. By Green's theorem, we have:

D1dη=D1η=(C1)+η+Kη=(C1)+ηK+η,

where we used the fact that line integrals pick a sign if we change the orientation. But the left-hand-side of the above equation is necessarily zero, since dη=0. Therefore,

(C1)+η=K+η,

which is the statement that we were trying to prove, with both curves canonically oriented.

(d) The power of what we did in (c) is to replace the evaluation of the line integral of η along an arbitrary simple closed curve C1 that encloses the origin by the evaluation of the line integral of η along a circle K centered at the origin of a certain non-zero radius. But we can evaluate the latter explicitly! Indeed, if K has radius R, for some RR>0, then we can parametrize K by α:[0,2π]R2 with

α(θ)=(Rcos(θ),Rsin(θ)).

The pullback of η is

αη=1R2(Rsin(θ)(Rsin(θ))+Rcos(θ)(Rcos(θ)))dθ=dθ.

Therefore,

αη=[0,2π]αη=02πdθ=2π.

So the line integral of η along any circle K centered at the origin is always non-zero and equal to 2π. By (c), we thus conclude that

C1η=2π

for arbitrary simple closed curves C1 that enclose the origin. Very powerful result!