0% found this document useful (0 votes)
18 views11 pages

Green

DLD subject

Uploaded by

mhsinm2003
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
18 views11 pages

Green

DLD subject

Uploaded by

mhsinm2003
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 11

1

Green’s Theorem

E. L. Lady

February 14, 2000

One of the things that makes Green’s Theorem


I ZZ
∂Q ∂P
P dx + Qdy = − dx dy
C Ω ∂x ∂y

[ where C is a simple closed curve and P and Q are functions of x and y which have
continuous partial derivatives in the region enclosed by C ] look more intimidating than it
is is that it’s actually two theorems written as one:
I ZZ I ZZ
∂Q ∂P
Q dy = dx dy , P dx = − dx dy.
C Ω ∂x C Ω ∂y

Another thing that makes the basic idea not as clear as it should be is that there’s
no notation for what ought to be called a partial
Z integral, in analogy to a partial
derivative. In other words, when one writes f (x, y) dx , there are two ways of
interpreting this. In the “partial integral,” y is treated as a parameter which, for the
purposes of the integration, functions as if it were a constant.Z ZThis “partial integral” is
what one uses in the process of evaluating a double integral f (x, y) dx dy . On the
Z Ω

other hand, when doing a line integral P (x, y) dx , one does not use a partial integral:
C
here, y is not independent of x but instead represents a function of x which must be
taken into account when doing the integral. (In practice, it is often better to think of
both x and y as functions of some third variable t , but the point is that y is not held
constant during the integration.)

Now if we let the integral sign for the moment represent a “partial” integral, then this
is the opposite of the partial derivative, so that
Z
∂P (x, y)
dy = P (x, y) + C(x).
∂y
2

Here C(x) is the “constant” of integration, except that since we are only integrating y ,
it may actually be a function of x . When doing a double integral, we can ignore this
“constant” C(x) (Explain!) and we get

ZZ Z
∂P (x, y)
dx dy = P (x, y) dx,
∂y

which seems to be essentially Green’s Theorem.

Now this seems more or less plausible, but if a student is as skeptical as s/he ought to
be, this “proof” of Green’s Theorem will bother him [her] a little bit. There are in fact
several things that seem a little puzzling. It takes a while to notice all of them, but the
puzzlements are as follows:
(1) Why does the actual Green’s Theorem have a minus sign where the formula we
just derived does not? (This is the most obvious signal that the “proof” given
above is sliding over some subtleties.)
(2) Why is the single integral in Green’s Theorem a line integral?
∂P
(3) When we get P by integrating with respect to y , why does the answer still
∂y
have a y in it?

If the student is not confused at this point, s/he should stop now and take the time to
think about it until s/he becomes confused. Confusion is an essential part of the process
of learning to understand what’s really going on.

As an example, let’s see how this works out for P (x, y) = y 3 cos 5x if we let Ω be the
region bounded by the circle C around the origin with radius 2. The top boundary of Ω
p p
is the graph of the function y = 4 − x2 and the bottom is the curve y = − 4 − x2 ,
∂P
where in both cases x ranges from −2 to 2 . Then = 3y 2 cos 5x and
∂y

ZZ Z Z √
2 4−x2
∂P ∂P
dy dx = √ dy dx
Ω ∂y −2 − 4−x2 ∂y
Z 2 Z √4−x2
= √ 3y 2 cos 5x dy dx
−2 − 4−x2
Z 2 √
4−x2
= y 3 cos 5x √
y=− 4−x2
dx
−2
Z 2
=2 (4 − x2 )3/2 cos 5x dx .
−2
3

Notice how y disappears after we do the inside integral. The plausible formula which we
wrote in the abstract as ZZ Z
∂P
dy dx = P dx
Ω ∂y

turns out to be a little misleading once we fill in the details. In this example it should
actually look like

ZZ Z 2 4−x2
∂P
dy dx = P (x, y) √
dx .
Ω ∂y −2 y=− 4−x2

On the other hand, Green’s Theorem asserts that this should equal
I I
− P (x, y) dx = − y 3 cos 5x dx ,
C C

where C is the boundary of Ω , i. e the circle with radius 2 centered at the origin.

The two formulas certainly have some ressemblance, but it is far from obvious that
they are the same.

Line Integrals

In order to understand Green’s Theorem, you have to understand line integrals in the
most down-to-earth way possible. In fact, it’s going to be important to really run the
explanation into the ground.

We start with a curve C and a vector parametrization r(t) for this curve. In
general, we denote a point on the curve by (x, y) , where x = x(t) , y = y(t) . Thus
r(t) = x(t) i + y(t) j . We also have a vector field F(x, y) defined in the whole plane, or at
least enough of the plane to include the curve. (For Green’s Theorem, it is essential that
F(x, y) be defined and differentiable on the region inside the curve as well as on the curve
itself.) This is why F is written as a function of the point (x, y) and not a function of t .
It is standard to denote the components of the vector F(x, y) by P and Q :

F(x, y) = P (x, y) i + Q(x, y) j.

Now we want to define the concept of the line integral of F along C . The most
conceptual notation for this is Z
F · dr.
C

(“Conceptual” in this case means that the formula does not refer to the coordinates
x and y . Thus the formula does not change when we make a change of coordinates. In
4

general, the more “conceptual” a formula is, the less practical it is for doing an actual
calculation.)

Now I want to put aside all considerations of what this formula really means (in other
words, all “conceptual” considerations) and simply try to understand this line integral in
terms of a simple integral by substitution.

We have the following formulas:


F(x, y) = P (x, y) i + Q(x, y) j
r(t) = x(t) i + y(t) j
dr(t) = r0 (t) dt
= ( x0 (t) i + y 0 (t) j ) dt
F(x, y) · dr = [P (x, y)x0(t) + Q(x, y)y 0 (t)] dt.

We can put this all together, and also write

dx = x0 (t) dt , dy = y 0 (t) dt

and thus get


Z Z t1
F · dr = [P (t) x0 (t) + Q(t) y 0 (t)] dt
C t0
Z
= P dx + Q dy,
C

where in the first line we have used a notation common in physics but generally
avoided by mathematicians by writing P (t) as an abbreviation for P (x(t), y(t)) and
Q(t) = Q(x(t), y(t)) .

The P dx + Q dy notation is very common and often convenient. However it looks


rather strange because one is not used to seeing a differential
Z Zthat is the sum of two parts.
One could make it less strange by instead writing P dx + Q dy . Although this is
C C
quite correct, and corresponds more closely to the actual process of doing a calculation, it
is conceptually less desirable, since it makes the line integral seem like two things rather
than a single unified thing.
Z
In practice, in any case, in order to actually calculate P (x, y) dx one usually needs
C
to use a parameter t so that
Z Z t1
P (x, y) dx = P (x(t), y(t)) x0 (t) dt.
C t0
5
Z
In any case, it is important to remember that P (x, y) dx is not a partial integral.
C
The variable y in this integral is not an independent parameter, it is a function that
depends on the point on the curve, just as x does.
Z
To understand how the line integral P dx is different from an ordinary integral in
C
Calculus I, suppose that P (x, y) = y and C is the ellipse

x2 y2
+ = 1.
4 9

The easiest way to do this problem is to parametrize the ellipse as x(t) = 2 cos t ,
y(t) = 3 sin t . But it’s more enlightening to use a method that exposes us to some possible
pitfalls.

First of all, suppose we want to do the integral as written:


I Z x1
P dx = y dx.
C x0

Remembering that this is not a partial integral, it will be necessary to substitute for y
the appropriate function of x . But first there’s a more intriguing problem: what should
x0 and x1 be?

The values x0 and x1 should be the x-coordinates of the beginning and end of
the curve. But since the ellipse is a closed curve, it starts and ends at the same point.
It might make sense, for instance, to think of it starting and ending at (2, 0) . Thus
x0 = x1 = 2 , so should we write
I Z 2
y dx = y(x) dx ?
C 2

But this integral would be 0. And maybe, in fact, 0 is the right answer. But there are a
couple of reasons why we should be skeptical of this.

(1) It’s always wise to be skeptical of an answer that comes too easily. Sometimes,
of course, there really is an extremely easy way to do a problem. But when you get an
answer this easily, it’s always a good idea to stop and think for a while, to make sure that
it really makes sense.

(2) In getting this answer, we never used the formula for P or the formula for C . So
the logic of what we’ve done would say that every line integral around any closed curve
would be zero. But that doesn’t seem like it ought to be true. If it were true, it would
certainly make Green’s Theorem completely unnecessary.
6

Since it doesn’t work to take x0 = x1 , we need to rethink the process. One idea would
be that since (−2, 0) is the point on the ellipse furthest to the left, and (2, 0) is the point
furthest to the right, we should take
I Z 2
y dx = y(x) dx (?)
C −2

But this also turns out to be incorrect.

In learning mathematics, it’s very important to make mistakes. Often only the process
of trying to figure out what’s wrong when you’ve made a mistake can cause you to think
about what’s happening deeply enough to get real insight.

In this case, what’s happening is that the line integral has to go all the way around
the ellipse. When it does this, x will start at a certain x0 and end at the same x0 , but in
between the point must travel over both the top and the bottom of the ellipse.

Now x alone is inadequate to determine whether one is one the top or bottom of the
ellipse. To know this, we have to know y .

We can write I Z Z
−2 2
y dx = y dx + y dx.
C 2 −2

This is correct, but misleading because it is overly abbreviated. In this formula, it seems
that the two formulas are negatives of each other and hence cancel. But this doesn’t
happen, because it’s not the same y in the two formulas. The first integral is on the top
2 2 p
of the ellipse x4 + y9 = 1 , so that y = 13 1 − (x2 /4) and the second integral is on the
p
bottom, so y = − 13 1 − (x2 /4) .

Thus we finally get the correct formula for the integral around the ellipse:
I Z r Z r
−2 2
1 x2 1 x2
y dx = 1− dx + − 1− dx
C 2 3 4 −2 3 4

Z r
2
2 x2
=− 1− dx.
3 −2 4
7

Proof of Green’s Theorem

Suppose now we have a situation similar to the ellipse just discussed, namely a simple
closed curve C consisting of a bottom half, given by the graph of a function y = ϕ(x) for
x between x0 and x1 and a top half given by y = ψ(x) , for x between x1 and x0 . We
assume that this curve is oriented counterclockwise, as will be the case if x0 < x1 . Also
continue to suppose that P and Q are functions of x and y having continuous partial
derivatives throughout the region enclosed by C .

Then using the logic above, we get


I Z x1 Z x0
P dx = P (x, ϕ(x)) dx + P (x, ψ(x)) dx
C x0 x1

Z x1 Z x1
= P (x, ϕ(x)) dx − P (x, ψ(x)) dx .
x0 x0

Now let Ω be the area enclosed by C . Using the usual rules for double integrals,
ZZ Z x1 Z ψ(x)
∂P ∂P (x, y)
dx dy = dy dx .
Ω ∂y x0 ϕ(x) ∂y
Notice that, as previously discussed, these integrals, unlike the line integral, are partial
∂P
integrals, so when we integral with respect to dy we get P :
∂y
Z ψ(x) y=ψ(x)
∂P
dy = P (x, y)
ϕ(x) ∂y y=ϕ(x)

= P (x, ψ(x)) − P (x, ϕ(x)) .


Thus
ZZ Z x1 Z ψ(x)
∂P ∂P (x, y)
dx dy = dy dx
Ω ∂y x0 ϕ(x) ∂y
Z x1
= P (x, ψ(x)) − P (x, ϕ(x)) dx .
x0

On the other hand, as we have seen,


I Z x1 Z x1
P (x, y) dx = P (x, ϕ(x)) dx − P (x, ψ(x)) dx .
C x0 x0

It is now obvious that I ZZ


∂P
P (x, y) dx = − dx dy .
C Ω ∂y
8

Now we have the proof of Green’s Theorem. It is worthwhile at this point to go back
to the very beginning of the explanation and answer the questions that arose.

It seemed at first that the formula


ZZ I
∂P
dy dx = − P (x, y) dx
Ω ∂y C

∂P
was almost obvious from the fact that when you integrate with respect to dy you
∂y
should get P . However we noticed a few puzzling questions.
(1) Why does this formula have a minus sign in it? (This is the most obvious signal
that the “proof” originally given was sliding over some subtleties.)
(2) Why is right-hand side of this formula a line integral?
∂P
(3) When we get P by integrating with respect to y , why does the answer still
∂y
have a y in it?

We can now answer these questions.

(2) The right-hand side of the formula in (2) is a line integral over C because in
evaluating the double integral, P has to be integrated between the bottom boundary and
the top boundary of Ω and these boundaries are what make up C .

(3) When we integrate P (x, y) with respect to y , the answer only looks like it still
has a y in it. When we do the line integral, y is not an independent variable but gets
replaced by a function of x .

(1) The minus sign on the right-hand side of the formula occurs because the line
integral is done counter-clockwise, and this means that we go backwards over the top of
the boundary of Ω and forwards over the bottom, which is exactly the opposite of what is
required for the double integral.

Finally, we should note that this proof only applies to a special case. There are
many regions Ω where the top and bottom boundaries will not be graphs for functions
y = ψ(x) and y = ϕ(x) . But this special case basically shows the logic that makes
Green’s Theorem true. Most calculus books show how to generalize the proof for regions
with more complicated boundaries.

We trust that the student will easily figure out for him/herself to prove the other half
needed for Green’s Theorem, namely
I ZZ
∂Q
Q dy = dx dy.
C Ω ∂x
9

Example. Let F(x, y) = −y i + x j and let C be the circle of radius a centered at the
origin. We can parametrize C in the usual way:

r(θ) = a cos θ i + a sin θ j.

We then see that r0 (θ) = −a sin θ i + a cos θ j = −y i + x j = F(r(θ)), so that


I Z 2π Z 2π Z 2π
0 2
F · dr = F · r dθ = F · F dθ = a dθ = 2πa2 .
C 0 0 0

On the other hand,


ZZ ZZ
∂Q ∂P
− dx dy = 1 − (−1) dx dy = 2πa2 ,
Ω ∂x ∂y Ω

since the area of Ω is πa2 . Thus the two calculations give the same result, as predicted
by Green’s Theorem. X

There are still a few comments on Green’s Theorem that needed to be added. To start
with, an extremely important one:

Green’s Theorem is not the formula we have proved.

Instead, Green’s Theorem consists of this formula together with the words that
surround it. Now it is understandable that one should give slight attention to these
words, regarding them as essential but not very important wrapping paper enveloping the
assertion of real interest. But to take this point of view is the make a very serious error.
The most important part of the wrapping paper here is the assertion in Green’s
Theorem that P and Q must be differentiable functions of (x, y) within the area enclosed
by C and their partial derivatives must be continuous. The assertion that C be a simple
closed curve is also important, but there is less of a temptation to get oneself into trouble
by forgetting it.

Consider the following example:


−y x
P = , Q= .
x2 + y2 x2 + y2

These functions are mildly unpleasant to differentiate, but if one takes the effort one finds
that
∂Q ∂P y 2 − x2 y 2 − x2
− = 2 − = 0.
∂x ∂y (x + y 2 )2 (x2 + y 2 )2
10

Now let C be a circle with radius a around the origin. We can parametrize C in the usual
way:
r(θ) = a cos θ i + a sin θ j , −0 ≤ θ ≤ 2π.
(We have denoted the parameter by θ rather than t in order to notice a relationship to
polar coordinates.)
I
If we now attempt to calculate P dx + Q dy using Green’s Theorem, we get
ZZ C

0 dx dy = 0 . However, it is easy to do the integral directly and thus see that the

correct answer is different.
In fact, notice that in polar coordinates

−y −r sin θ − sin θ cos θ


P = = = , Q= ,
x2+y 2 r 2 r r
p
where, as usual, r = x2 + y 2 , so that r = a on the circle C . Furthermore, we see
immediately that
r0 (θ) = −a sin θ i + a cos θ j,
so that
Z Z Z 2π
0 a(sin2 θ + cos2 θ)
P dx + Q dy = (P i + Q j) · r dθ = dθ = 2π.
C C 0 a

Green’s Theorem fails for this example, because P (x, y) and Q(x, y) are not
continuous at the origin, where P and Q have zero in the denominator.
It is at first somewhat amazing that the failure of P and Q to be continuous at a
single point could cause Green’s Theorem to go bad here. But it is indeed the case, and is
∂Q ∂P
a little less surprising if one realizes that and both blow up as (x, y) approaches
ZZ ZZ ∂x ∂y
∂Q ∂P
the origin, so that dx dy and dx dy become improper integrals and, as
Ω ∂x Ω ∂y
far as I know, don’t converge.
It is interesting that in the line integral in this example, the result, 2π , is independent
of the radius a of the circle. In fact, by doing the calculation slightly differently, we can
see that much more is true. Let now C be any simple closed curve going around the
origin counterclockwise and which can be described in polar coordinates: r = r(θ). This
gives a parametrization of the curve with θ as the parameter:

x = r(θ) cos θ, y = r(θ) sin θ.


11

From this, we get by the product rule,


dx
dx = dθ = (r 0 cos θ − r sin θ) dθ

dy
dy = dθ = (r 0 sin θ + r cos θ) dθ.

Changing P and Q to polar coordinates as before, we get
Z Z 2π
−r sin θ(r 0 cos θ − r sin θ) + r cos θ(r 0 sin θ − r sin θ)
P dx + Q dy = dθ
C 0 r2
Z 2π 2
r (sin2 θ + cos2 θ)
= dθ = 2π.
0 r2
Thus the answer is the same for any closed curve going around the origin
counterclockwise. (As the curve goes around the origin, it sweeps through all four
quadrants of the plane, so that θ goes from 0 to 2π .)
This fact can actually be seen to be a consequence of Green’s Theorem.

I also want to say that it’s a little disingenous of me to state that Green’s Theorem
is actually two separate theorems, one for P and one for Q . Although this is a useful
way to think of it for purposes of finding a proof, it really puts too much emphasis on the
coordinates x and y and obscures the fact that Green’s Theorem has a significance that
transcends a particular choice of a coordinate system.
Furthermore, the expression
∂Q ∂P

∂x ∂y
actually has an intrinsic significance for any vector field. It is essentially what in three
dimensional vector analysis is called the “curl,” and measures the extent to which there is
a rotation in the vector field. (The easiest way of understanding this is actually through
the use of Green’s Theorem.)
For the moment, let’s just point out a simple case showing the significance of the curl.

Theorem. If P (x, y) and Q(x, y) are continuously differentiable functions throughout a


certain region and if
∂Q ∂P
− =0
∂x ∂y
through out that region, then for any simple closed curve C whose interior is completely
enclosed within the region, I
P dx + Q dy = 0.
C

proof: This is an immediate consequence of Green’s Theorem. X

You might also like