Loading [MathJax]/jax/output/HTML-CSS/jax.js

We now come to the first of three important theorems that extend the Fundamental Theorem of Calculus to higher dimensions. (The Fundamental Theorem of Line Integrals has already done this in one way, but in that case we were still dealing with an essentially one-dimensional integral.) They all share with the Fundamental Theorem the following rather vague description: To compute a certain sort of integral over a region, we may do a computation on the boundary of the region that involves one fewer integrations.

Note that this does indeed describe the Fundamental Theorem of Calculus and the Fundamental Theorem of Line Integrals: to compute a single integral over an interval, we do a computation on the boundary (the endpoints) that involves one fewer integrations, namely, no integrations at all.

Theorem 18.4.1 (Green's Theorem) If the vector field F=P,Q and the region D are sufficiently nice, and if C is the boundary of D (C is a closed curve), then DQxPydA=CPdx+Qdy, provided the integration on the right is done counter-clockwise around C.

To indicate that an integral C is being done over a closed curve in the counter-clockwise direction, we usually write C. We also use the notation D to mean the boundary of D oriented in the counterclockwise direction. With this notation, C=D.

We already know one case, not particularly interesting, in which this theorem is true: If F is conservative, we know that the integral CFdr=0, because any integral of a conservative vector field around a closed curve is zero. We also know in this case that P/y=Q/x, so the double integral in the theorem is simply the integral of the zero function, namely, 0. So in the case that F is conservative, the theorem says simply that 0=0.

Example 18.4.2 We illustrate the theorem by computing both sides of Dx4dx+xydy=Dy0dA, where D is the triangular region with corners (0,0), (1,0), (0,1).

Starting with the double integral: Dy0dA=101x0ydydx=10(1x)22dx=(1x)36|10=16.

There is no single formula to describe the boundary of D, so to compute the left side directly we need to compute three separate integrals corresponding to the three sides of the triangle, and each of these integrals we break into two integrals, the "dx'' part and the "dy'' part. The three sides are described by y=0, y=1x, and x=0. The integrals are then Dx4dx+xydy=10x4dx+000dy+01x4dx+10(1y)ydy+000dx+010dy=15+015+16+0+0=16.

Alternately, we could describe the three sides in vector form as t,0, 1t,t, and 0,1t. Note that in each case, as t ranges from 0 to 1, we follow the corresponding side in the correct direction. Now Dx4dx+xydy=10t4+t0dt+10(1t)4+(1t)tdt+100+0dt=10t4dt+10(1t)4+(1t)tdt=16.

In this case, none of the integrations are difficult, but the second approach is somewhat tedious because of the necessity to set up three different integrals. In different circumstances, either of the integrals, the single or the double, might be easier to compute. Sometimes it is worthwhile to turn a single integral into the corresponding double integral, sometimes exactly the opposite approach is best.

Here is a clever use of Green's Theorem: We know that areas can be computed using double integrals, namely, D1dA computes the area of region D. If we can find P and Q so that Q/xP/y=1, then the area is also DPdx+Qdy. It is quite easy to do this: P=0,Q=x works, as do P=y,Q=0 and P=y/2,Q=x/2.

Example 18.4.3 An ellipse centered at the origin, with its two principal axes aligned with the x and y axes, is given by x2a2+y2b2=1. We find the area of the interior of the ellipse via Green's theorem. To do this we need a vector equation for the boundary; one such equation is acost,bsint, as t ranges from 0 to 2π. We can easily verify this by substitution: x2a2+y2b2=a2cos2ta2+b2sin2tb2=cos2t+sin2t=1. Let's consider the three possibilities for P and Q above: Using 0 and x gives C0dx+xdy=2π0acos(t)bcos(t)dt=2π0abcos2(t)dt. Using y and 0 gives Cydx+0dy=2π0bsin(t)(asin(t))dt=2π0absin2(t)dt. Finally, using y/2 and x/2 gives Cy2dx+x2dy=2π0bsin(t)2(asin(t))dt+acos(t)2(bcos(t))dt=2π0absin2t2+abcos2t2dt=2π0ab2dt=πab. The first two integrals are not particularly difficult, but the third is very easy, though the choice of P and Q seems more complicated.

(0,b)
(a,0)
Figure 18.4.1. A "standard'' ellipse, x2a2+y2b2=1.

Proof of Green's Theorem.
We cannot here prove Green's Theorem in general, but we can do a special case. We seek to prove that CPdx+Qdy=DQxPydA. It is sufficient to show that CPdx=DPydAandCQdy=DQxdA, which we can do if we can compute the double integral in both possible ways, that is, using dA=dydx and dA=dxdy.

For the first equation, we start with DPydA=bag2(x)g1(x)Pydydx=baP(x,g2(x))P(x,g1(x))dx. Here we have simply used the ordinary Fundamental Theorem of Calculus, since for the inner integral we are integrating a derivative with respect to y: an antiderivative of P/y with respect to y is simply P(x,y), and then we substitute g1 and g2 for y and subtract.

Now we need to manipulate CPdx. The boundary of region D consists of 4 parts, given by the equations y=g1(x), x=b, y=g2(x), and x=a. On the portions x=b and x=a, dx=0dt, so the corresponding integrals are zero. For the other two portions, we use the parametric forms x=t, y=g1(t), atb, and x=t, y=g2(t), letting t range from b to a, since we are integrating counter-clockwise around the boundary. The resulting integrals give us CPdx=baP(t,g1(t))dt+abP(t,g2(t))dt=baP(t,g1(t))dtbaP(t,g2(t))dt=baP(t,g1(t))P(t,g2(t))dt which is the result of the double integral times 1, as desired.

The equation involving Q is essentially the same, and left as an exercise.

Exercises 18.4

Here is a Sage cell if you'd like to use it.

Ex 18.4.1 Compute D2ydx+3xdy, where D is described by 0x1, 0y1. (answer)

Ex 18.4.2 Compute Dxydx+xydy, where D is described by 0x1, 0y1. (answer)

Ex 18.4.3 Compute De2x+3ydx+exydy, where D is described by 2x2, 1y1. (answer)

Ex 18.4.4 Compute Dycosxdx+ysinxdy, where D is described by 0xπ/2, 1y2. (answer)

Ex 18.4.5 Compute Dx2ydx+xy2dy, where D is described by 0x1, 0yx. (answer)

Ex 18.4.6 Compute Dxydx+x+ydy, where D is described by 1x2, 2xy4. (answer)

Ex 18.4.7 Compute D(x/y)dx+(2+3x)dy, where D is described by 1x2, 1yx2. (answer)

Ex 18.4.8 Compute Dsinydx+sinxdy, where D is described by 0xπ/2, xyπ/2. (answer)

Ex 18.4.9 Compute Dxlnydx, where D is described by 1x2, exyex2. (answer)

Ex 18.4.10 Compute D1+x2dy, where D is described by 1x1, x2y1. (answer)

Ex 18.4.11 Compute Dx2ydxxy2dy, where D is described by x2+y21. (answer)

Ex 18.4.12 Compute Dy3dx+2x3dy, where D is described by x2+y24. (answer)

Ex 18.4.13 Evaluate C(ysin(x))dx+cos(x)dy, where C is the boundary of the triangle with vertices (0,0), (1,0), and (1,2) oriented counter-clockwise. (answer)

Ex 18.4.14 Finish our proof of Green's Theorem by showing that CQdy=DQxdA.