Skip to main content
Logo image

Coordinated Multivariable Calculus

Section 4.6 Green’s Theorem

We know from Section 4.4 that a vector field is path-independent if and only if the circulation around every closed curve in its domain is \(0\text{.}\) It is probably not surprising that, given the multitude of names we have for path-independent vector fields, they are important vector fields that arise frequently. However, not every vector field is path-independent, and many times we will want to calculate the circulation around a closed curve in a vector field that is not path-independent. This section explores a connection between line integrals and double integrals that you may find surprising.

Preview Activity 4.6.1.

We will consider the vector field \(\vF = \langle 2y,3x^2 y\rangle\text{,}\) which is defined on the entire plane. Suppose that we want to calculate the circulation of \(\vF\) around the circle \(C\) of radius \(2\text{,}\) centered at \((0,0)\text{,}\) and oriented counterclockwise.

(a)

Verify that \(\vF\) is not path-independent by calculating the circulation of \(\vF\) around the circle \(C\text{.}\) The SageMath cell below is set up to assist you with this, but you will need to supply a parametrization of \(C\) on line 3.
Answer.
\(-8\pi\)
Solution.
If we parametrize \(C\) as \(\vr(t) = \langle 2\cos(t),2\sin(t)\rangle\text{,}\) then \(\oint_C \vF\cdot d\vr = -8\pi\text{.}\) To use the SageMath cell to do this, make line 3 read C = vector([2*cos(t),2*sin(t)]) and execute. Because the circulation around \(C\) is not \(0\text{,}\) we know that \(\vF\) is not path-independent.

(b)

Consider \(\vF = \langle F_1(x,y),F_2(x,y)\rangle\text{.}\) In Section 4.5, we discussed the circulation density of a two-dimensional vector field at a given point and determined that it can be algebraically determined by \(\displaystyle \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y}\text{.}\) Calculate the circulation density for this vector field.
Answer.
No; \(6xy-2\text{.}\)
Solution.
Here we have that \(\partial F_2/\partial x =6xy\) and \(\partial F_1/\partial y = 2\text{,}\) so we have that \(\partial F_2/\partial x - \partial F_1/\partial y= 6xy-2\text{.}\)

(c)

Sketch the curve \(C\) and shade the region inside it. Describe the region bounded by \(C\) in both rectangular and polar coordinates.
Answer.
Polar: \(0\leq r\leq 2\text{,}\) \(0\leq \theta\leq 2\pi\text{.}\) Rectangular: \(-2\leq x\leq 2\text{,}\) \(-\sqrt{4-x^2}\leq y\leq \sqrt{4-x^2}\text{.}\)
Solution.
Since we know the radius of \(C\) is \(2\text{,}\) the region the circle bounds is covered by all values of \(r\) with \(0\leq r\leq 2\text{.}\) For \(\theta\text{,}\) since we have the full disk, the bounds are \(0\leq \theta\leq 2\pi\text{.}\) In rectangular coordinates, the circle is \(x^2+y^2=4\text{.}\) The range for \(x\) is \(-2\leq x\leq 2\text{,}\) and solving the equation of the circle for \(y\) gives us \(-\sqrt{4-x^2}\leq y\leq \sqrt{4-x^2}\) as the bounds for \(y\text{.}\)

(d)

Calculate the double integral of \(\displaystyle \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y}\) over the region inside the circle \(C\text{.}\) This integral is not the most fun to do by hand, so a SageMath cell has been provided to assist you.
Hint.
If you choose polar coordinates, use \(\theta\) for the first variable and \(r\) for the second variable.
Answer.
\(-8\pi\)
Solution.
Denote by \(D\) the region inside the circle \(C\text{.}\) We can describe \(D\) readily in polar coordinates with the inequalities \(0\leq r \leq 2\) and \(0\leq \theta\leq 2\pi\text{.}\) Thus, we need to calculate
\begin{align*} \iint_D \left(\frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y}\right)\, dA &= \iint_D (6xy-2)\, dA\\ &= \int_0^2 \int_0^{2\pi} (6r^2\cos(\theta)\sin(\theta) - 2)r\, d\theta\, dr\text{.} \end{align*}
Using SageMath to calculate this, we obtain \(-8\pi\text{.}\)
A natural question after completing this Preview Activity is if there is something special about the vector field \(\vF\) or the curve \(C\) that led to the results you obtained, and investigating this will be our principal task in this section

Subsection 4.6.1 Circulation

In Preview Activity 4.6.1, you integrated the circulation density of a smooth vector field over a disk and found the result was equal to the circulation of the vector field along the region’s circular boundary. This relationship is the theme of this section. To see why this would make sense, consider the region \(R\) bounded by the curve \(C\) shown in Figure 4.6.1. We have placed a square grid inside \(R\) to suggest the idea of breaking \(R\) up into many smaller regions, most of which are square. This idea should make you think of the methods we have already seen of breaking up a region into smaller and smaller regions for Riemann sums.
described in detail following the image
A bean-shaped region in the plane that is bounded by a closed curve. The curve is oriented so that when walking along the curve, the region is always on the left. There is a square grid dividing the up the interior of the region.
Figure 4.6.1. A region bounded by a curve
The critical idea here is that if we integrate the circulation density over each of the small regions and add those up, this is the same as integrating the circulation density over the entirety of the region \(R\) because of the fundamental properties of integrals. Also, integrating circulation per unit area over a two-dimensional region should be related to the total circulation on that region. Now look at Figure 4.6.2 and think of these four square regions as being four of the square regions inside \(R\) in Figure 4.6.1.
described in detail following the image
A square region in the plane that is bounded by a closed curve oriented counterclockwise. The square is divided into four squares of equal size, which are labeled as \(R_1\text{,}\) \(R_2\text{,}\) \(R_3\text{,}\) and \(R_4\text{.}\) There are circular arrows inside each of the smaller regions indicating counterclockwise orientations.
Figure 4.6.2. A square region divided into four smaller squares
We orient the boundary \(C_i\) of square region \(R_i\) in the manner suggested by the circular arrows. This means that the vertical boundary in common between \(R_1\) and \(R_2\) is oriented up when we calculate the line integral \(\oint_{C_1}\vF\cdot d\vr\) and oriented down when we evaluate \(\oint_{C_2}\vF\cdot d\vr\text{.}\) Thus, this line segment does not contribute to the sum \(\oint_{C_1}\vF\cdot d\vr + \oint_{C_2}\vF\cdot d\vr\text{.}\) The same idea applies to the other three interior line segments so that
\begin{equation*} \oint_{C_1}\vF\cdot d\vr + \oint_{C_2}\vF\cdot d\vr + \oint_{C_3}\vF\cdot d\vr + \oint_{C_4}\vF\cdot d\vr \end{equation*}
is equal to the line integral of \(\vF\) along the boundary of the large square.
Returning to Figure 4.6.1, if we find the circulation along the boundary of each of the smaller regions, the line integrals along the boundaries that lie inside the region \(R\) will all offset. Thus the value should equal \(\oint_C\vF\cdot d\vr\text{.}\) If we make our grid fine enough, all of the smaller regions into which \(R\) is divided will be very close to rectangular. Our next activity looks at the relationship between the circulation along a rectangular curve with sides parallel to the coordinate axes and the double integral of the circulation density over the rectangular region inside the curve.

Activity 4.6.2.

For this activity, we will consider the region \(R\) satisfying \(a\leq x\leq c\) and \(b\leq y\leq d\text{.}\) The rectangular curve \(C\) will be the boundary of \(R\text{,}\) oriented counterclockwise. Let \(\vF = \langle F_1(x,y),F_2(x,y)\rangle\) be a vector field that is smooth on \(C\) and \(R\text{.}\) We can write \(C = C_1 + C_2 + C_3 + C_4\text{,}\) where \(C_1\) is the top edge of the rectangle and the edges are numbered counterclockwise around the rectangle.
(a)
Sketch the region \(R\) and curves \(C_1,C_2,C_3,C_4\text{.}\) Label the four corners of \(R\) with the coordinates of the vertices, and be sure to indicate the proper orientation on \(C\text{.}\)
(b)
Since we know that
\begin{equation*} \oint_C\vF\cdot d\vr = \int_{C_1}\vF\cdot d\vr + \int_{C_2}\vF\cdot d\vr + \int_{C_3}\vF\cdot d\vr + \int_{C_4}\vF\cdot d\vr\text{,} \end{equation*}
it makes sense to parametrize each curve \(C_i\) and write \(\int_{C_i}\vF\cdot d\vr\) as a one-variable integral. For example, \(C_1\) can be parametrized as \(\vr(x) = \langle c+a-x,d\rangle\) with \(a\leq x\leq c\text{.}\) Thus, we have
\begin{equation*} \int_{C_1}\vF\cdot d\vr = \int_a^c \langle F_1(x,d),F_2(x,d)\rangle\cdot \langle -1,0\rangle\, dx\text{.} \end{equation*}
Complete this process to write the integrals for the remaining three sides as integrals with respect to \(x\) or \(y\text{,}\) as appropriate.
Answer.
\begin{align*} \int_{C_2}\vF\cdot d\vr &= \int_b^d \langle F_1(a,y),F_2(a,y)\rangle\cdot \langle 0,-1\rangle\, dy\\ \int_{C_3}\vF\cdot d\vr &= \int_a^c \langle F_1(x,b),F_2(x,b)\rangle\cdot \langle 1,0\rangle\, dy\\ \int_{C_4}\vF\cdot d\vr &= \int_b^d \langle F_1(c,y),F_2(c,y)\rangle\cdot \langle 0,1\rangle\, dy \end{align*}
Solution.
The parametrizations rely on fixing \(x=a\) for \(C_2\text{,}\) \(y=b\) for \(C_3\text{,}\) and \(x=d\) for \(C_4\text{.}\)
\begin{align*} \int_{C_2}\vF\cdot d\vr &= \int_b^d \langle F_1(a,y),F_2(a,y)\rangle\cdot \langle 0,-11\rangle\, dy\\ \int_{C_3}\vF\cdot d\vr &= \int_a^c \langle F_1(x,b),F_2(x,b)\rangle\cdot \langle 1,0\rangle\, dy\\ \int_{C_4}\vF\cdot d\vr &= \int_b^d \langle F_1(c,y),F_2(c,y)\rangle\cdot \langle 0,1\rangle\, dy \end{align*}
(c)
Write \(\oint_C\vF\cdot d\vr\) as the sum of two one-variable integrals (one with respect to \(x\) and one with respect to \(y\)).
Answer.
\(\displaystyle\oint_C\vF\cdot d\vr = \int_b^d \left(F_2(c,y) - F_2(a,y)\right)\, dy + \int_a^c\left(F_1(x,b)-F_1(x,d)\right)\, dx\)
Solution.
Because the limits of integration for the \(dx\) integrals match and the limits of integration for the \(dy\) integrals match, we can add and obtain \(\displaystyle\oint_C\vF\cdot d\vr = \int_b^d \left(F_2(c,y) - F_2(a,y)\right)\, dy + \int_a^c\left(F_1(x,b)-F_1(x,d)\right)\, dx\text{.}\)
(d)
We would like to see if the equality of the line integral and double integral in Preview Activity 4.6.1 was a coincidence or holds for more general vector fields. (Soon we will expand our curves beyond rectangles, but for the moment, we focus on the curve \(C\) and region \(R\text{.}\)) To do this, we have to find some way to take our answer to the previous part of this activity, which includes two single integrals and turn them into a double integral over the rectangle \(R\text{.}\)
We have assumed that \(\vF\) is smooth on \(C\) and \(R\text{,}\) so the partial derivatives \(\partial F_2/\partial x\) and \(\partial F_1/\partial y\) are continuous functions on \(C\) and \(R\text{.}\) Rewrite the expression you have for \(\oint_C\vF\cdot d\vr\) in the previous step as the sum of two double integrals.
Hint.
With \(x\) fixed, the integrand in your \(dx\) integral matches what you see on the left-hand side of the Total Change Theorem
 1 
activecalculus.org/single/sec-4-4-FTC.html#wtR
, but now instead of prime notation for the derivative, you will have a partial derivative with respect to \(y\text{.}\)
Answer.
\(\displaystyle \oint_C\vF\cdot d\vr = \int\limits_b^d\!\! \int\limits_{a\mathstrut}^c \frac{\partial F_2}{\partial x}(x,y)\, dx\, dy + \int\limits_{a\mathstrut}^c\!\! \int\limits_d^b \frac{\partial F_1}{\partial y}(x,y)\, dy\, dx\)
Solution.
By the Total Change Theorem
 2 
activecalculus.org/single/sec-4-4-FTC.html#wtR
, with \(y\) fixed, we have that
\begin{align*} F_2(c,y) - F_2(a,y) &= \int_a^c \frac{\partial F_2}{\partial x}(x,y)\, dx\text{.}\\ \end{align*}

Similarly,

\begin{align*} F_1(x,b) - F_1(x,d) &= \int_d^b \frac{\partial F_1}{\partial y}(x,y)\, dy\text{.} \end{align*}
Therefore, we have
\begin{equation*} \oint_C\vF\cdot d\vr = \int\limits_b^d\!\! \int\limits_{a\mathstrut}^c \frac{\partial F_2}{\partial x}(x,y)\, dx\, dy + \int\limits_{a\mathstrut}^c\!\! \int\limits_d^b \frac{\partial F_1}{\partial y}(x,y)\, dy\, dx\text{.} \end{equation*}
(e)
Rewrite your sum of two double integrals in the previous step above as a double integral over the rectangle \(R\text{,}\) thereby showing that the line integral of \(\vF\) along the curve \(C\) is equal to a double integral of the region inside \(C\text{.}\)
Answer.
\(\displaystyle\oint_C\vF\cdot d\vr = \iint_R \left(\frac{\partial F_2}{\partial x}- \frac{\partial F_1}{\partial y}\right)\, dA\)
Solution.
After interchanging the order of integration so that both integrals are done \(dx\, dy\text{,}\) the other step we must do to add them is to make the integral of \(\partial F_1/\partial y\) have lower limit \(b\) and upper limit \(d\) (with respect to \(y\)). We do this by multiplying by \(-1\) and then can add the integrals to obtain
\begin{equation*} \displaystyle\oint_C\vF\cdot d\vr = \iint_R \left(\frac{\partial F_2}{\partial x}- \frac{\partial F_1}{\partial y}\right)\, dA\text{.} \end{equation*}
Now we can see that our result in Preview Activity 4.6.1 was not a coincidence and integrating the circulation density of a smooth vector field on a rectangular region gives the same result as calculating the circulation along the boundary of the region.
The phenomenon we verified for rectangular regions in Activity 4.6.2 is a more general result that will be the subject of our next section.

Subsection 4.6.2 Green’s Theorem

So far in this section, we have restricted ourselves to relatively nice closed curves when thinking about circulation. While the main theorem of this section will not allow us to consider arbitrary closed curves, it does cover more varied curves than we have discussed so far. A simple closed curve is a closed curve that does not cross itself, and these are the curves to which our next theorem applies.
In Activity 4.6.2, we showed that Green’s Theorem holds when the region \(R\) is a rectangle with sides parallel to the coordinate axes. The discussion that introduced the previous subsection may have you convinced that Green’s Theorem holds in its full form. Although our argument there was not a rigorous proof, we will not say more about why this theorem is true in full generality.
At first glance, it may seem that Green’s Theorem is of purely intellectual interest. However, we have already encountered situations where parametrization a curve can be complicated. This is particularly true when the curve has “corners” that require us to give separate parametrizations for several pieces of the curve, such as with the rectangular curve pictured in Figure 4.6.4. However, as with this rectangle, is often the case that a curve that is difficult to parametrize bounds a region that is not too complex to describe using rectangular or polar coordinates. For instance, the rectangular region in Figure 4.6.4 can be described as \(1\leq x\leq 4\) and \(2\leq y\leq 4\text{,}\) rather than parametrization each of the four sides of the rectangle separately. Additionally, the integrand of the double integral in Green’s Theorem involves partial derivatives that can sometimes result in an integrand that is very easy to work with. The purpose of Green’s Theorem is, at its core, to allow you to exchange one type of integration problem (a line integral) for another type of integration problem (a double integral). This will be a recurring theme as this chapter continues.
described in detail following the image
A rectangular curve in the \(xy\)-plane. The lower-left corner is labeled \((1,2)\text{.}\) The upper-right corner is labeled \((4,4)\text{.}\) The curve is oriented counterclockwise.
Figure 4.6.4. A rectangular curve
The next activity asks you to apply Green’s Theorem.

Activity 4.6.3.

For each of the curves described below, find the circulation of the given vector field around the curve. Do this both by calculating the line integral directly as well as by calculating the double integral from Green’s Theorem.
(a)
The curve \(C_1\) is the circle of radius \(3\) centered at the point \((2,1)\) (oriented counterclockwise) and the vector field is \(\vF = \langle y^2, 5x+2xy\rangle\text{.}\)
Answer.
\(45\pi\)
Solution.
Here we find that
\begin{equation*} \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y} = (5+2y)-2y = 5\text{.} \end{equation*}
Therefore, the integrand in Green’s Theorem is the constant \(5\text{.}\) Thus, the circulation of \(\vF\) around \(C\) is \(5\) times the area of the region inside \(C\text{.}\) That area is \(\pi\cdot 3^2 = 9\pi\text{,}\) which means that the circulation is \(45\pi\text{.}\)
(b)
The curve \(C_2\) is the triangle with vertices \((0,0)\text{,}\) \((3,0)\text{,}\) and \((3,3)\) (oriented counterclockwise) and the vector field is \(\vG = \langle y^2 , 3xy\rangle\text{.}\)
Answer.
\(27/6\)
Solution.
Here we find that
\begin{equation*} \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y} = 3y-2y = y\text{.} \end{equation*}
The region inside \(C_2\) can be described by \(0\leq x\leq 3\) and \(0\leq y\leq x\text{.}\) Therefore, by Green’s Theorem, we have
\begin{equation*} \oint_{C_2} \vG\cdot d\vr = \int_0^3\!\int_0^x y\, dy\, dx = \int_0^3 \left.\frac{y^2}{2}\right|_0^x\, dx = \int_0^3 \frac{x^2}{2}\, dx = \left.\frac{x^3}{6}\right|_0^3 = \frac{27}{6}\text{.} \end{equation*}
Generally, Green’s Theorem is useful for allowing us to calculate line integrals by instead calculating a double integral. Going the other direction is harder, since finding a vector field \(\langle F_1,F_2\rangle\) so that \(\partial F_2/\partial x - \partial F_1/\partial y\) is equal to the integrand in the double integral can be difficult. However, the exercises will explore some situations where calculating a suitable line integral is a viable alternative to the double integral.

Subsection 4.6.3 What happens when vector fields are not smooth?

Notice that the assumptions in Green’s Theorem require that the region \(R\) be bounded by a simple closed curve \(C\) and that the vector field have continuous partial derivatives on \(R\) and \(C\text{.}\) In Exercise 2, we will explore what happens when the region \(R\) cannot be bounded by a simple closed curve, as sometimes multiple applications of Green’s Theorem can be used in those circumstances. Now, however, we will take a look at what happens in some cases where the vector field \(\vF\) is not smooth.

Activity 4.6.4.

Consider the vector field \(\vF = \displaystyle\frac{-y}{{x^2+y^2}}\vi + \frac{x}{{x^2+y^2}}\vj\text{.}\) Notice that \(\vF\) is smooth everywhere in the plane other than at the point \((0,0)\text{.}\) This vector field is plotted in Figure 4.6.5, but we have not plotted the vectors close to the origin as their magnitudes get so large that they make it hard to interpret the figure. Here Green’s Theorem applies to any simple closed curve \(C\) that neither passes through \((0,0)\) nor bounds a region containing \((0,0)\text{.}\)
described in detail following the image
A vector field with vectors tangent to circles oriented counterclockwise and centered at the origin as well as the unit circle centered at the origin and oriented counterclockwise. The vectors are longer close to the origin and shorter farther from the origin.
Figure 4.6.5. The vector field \(\vF\)
(a)
Find the circulation density of \(\vF\) (i.e., the integrand of the double integral in Green’s Theorem).
Hint.
Make sure to note any exceptional points at which your formula is not valid.
Answer.
\(\displaystyle\frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y} = 0\) for \((x,y)\neq(0,0)\)
Solution.
For \((x,y)\neq (0,0)\text{,}\) we have
\begin{align*} \frac{\partial F_2}{\partial x} & = -\frac{2x^{2}}{{\left(x^{2} + y^{2}\right)}^{2}} + \frac{1}{{x^{2} + y^{2}}} = \frac{y^2-x^2}{{\left(x^{2} + y^{2}\right)}^{2}}\\ \frac{\partial F_1}{\partial y} & = \frac{2y^{2}}{{\left(x^{2} + y^{2}\right)}^{2}} - \frac{1}{{x^{2} + y^{2}}} = \frac{y^2-x^2}{{\left(x^{2} + y^{2}\right)}^{2}}\text{.} \end{align*}
Thus,
\begin{equation*} \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y} = 0\text{.} \end{equation*}
(b)
Suppose that \(C\) is the unit circle centered at the origin. Without doing any calculations, what can you say about \(\oint_C \vF\cdot d \vr\text{?}\) what does this tell you about if \(\vF\) is path-independent?
Answer.
Positive; \(\vF\) is not path-independent
Solution.
Notice that \(\vF\) is tangent to \(C\text{.}\) Thus, \(\oint_C \vF\cdot d \vr\) must be positive. (In fact, the circulation of \(\vF\) along \(C\) is the length of \(C\) or \(2\pi\) because \(\vF\) consists of unit vectors and the vector field is is tangent to \(C\text{.}\)) Since there is a closed curve with nonzero circulation, \(\vF\) is not path independent.
(c)
What would you get if you integrated the circulation density of \(\vF\) over the region bounded by \(C\text{?}\)
Answer.
\(0\)
Solution.
The circulation density is \(0\) everywhere inside \(C\) except at the origin, where it is undefined. That single point does not impact the value of the double integral of the circulation density, which is thus \(0\text{.}\)
(d)
Do the previous two parts contradict Green’s Theorem? Explain your reasoning.
Answer.
No.
Solution.
There is no contradiction because Green’s Theorem does not apply here. The region \(R\) inside the curve \(C\) includes the point \((0,0)\text{,}\) and \(\vF\) is not defined at that point. Green’s Theorem requires that the vector field be smooth at every point of \(R\text{.}\)
(e)
Is the vector field \(\vG = \displaystyle\frac{x}{x^2+y^2}\vi +\frac{y}{x^2+y^2}\vj\text{,}\) which is shown in Figure 4.6.6, path-independent? Why or why not?
described in detail following the image
A vector field with vectors radiating from the origin and orthogonal to circles centered at the origin. Vectors are longest near the origin and get shorter farther from the origin.
Figure 4.6.6. The vector field \(\vG\)
Answer.
Yes
Solution.
Let \(g(x,y) = \frac{1}{2}\ln(x^2+y^2)\text{.}\) Note that \(\grad g = \vG\text{,}\) and the domain of \(g\) is the same as the domain of \(\vG\) (all of \(\R^2\) except \((0,0)\)). Thus, \(\vG\) is a gradient vector field and is therefore path-independent.
(f)
Suppose that \(C\) is the unit circle centered at the origin. Find \(\oint_C \vG\cdot d \vr\text{.}\) Can you do this using Green’s Theorem?
Answer.
\(0\text{;}\) No
Solution.
We have already shown that \(\vG\) is path-independent, and therefore \(\oint_C \vG\cdot d\vr = 0\) since \(C\) is a closed curve. Because \(\vG\) is not defined at the origin, Green’s Theorem does not apply to this line integral.
We can now see that Green’s Theorem is a powerful tool. However it cannot be used for all line integrals. In particular, the restriction that the vector field be smooth on the entire region bounded by a simple closed curve \(C\) creates limitations. This can occur even in cases where holes in the vector field’s domain or points where the vector field is not continuously differentiable lie away from \(C\text{.}\)

Subsection 4.6.4 Computing area as a line integral

Recall that the area of a region \(D\) in the plane can be computed by integrating \(1\) over \(D\text{.}\) So if we can find a vector field \(\vF\) such that \(\frac{\partial F_2}{\partial x} - \frac {\partial F_1}{\partial y} = 1\text{,}\) then the double integral in Green’s Theorem becomes the area of \(D\text{.}\) Luckily, there are many such vector fields \(\vF\text{:}\)
\begin{equation*} \text{e.g.} \quad \vF = \langle 0,x \rangle, \quad \vF = \langle -y, 0 \rangle, \quad \vF = \langle -y/2, x / 2 \rangle \end{equation*}

Using Green’s Theorem to compute area.

If \(D\) is a region in the plane with boundary \(\partial D\text{,}\) then the area of \(D\) can be computed as any of the following line integrals around \(\partial D\text{.}\)
\begin{align*} \operatorname{area}(D) \amp = \iint_D 1\,dA \\ \amp = \oint_{\partial D} x\,dy\\ \amp = \oint_{\partial D} -y\,dx\\ \amp = \frac12 \oint_{\partial D} x\,dy - y\,dx \end{align*}

Example 4.6.7.

Compute the area of the ellipse
\begin{equation*} \left( \frac{x}{a} \right)^2 + \left( \frac{y}{b} \right)^2 = 1. \end{equation*}
Solution.
The ellipse is parametrized as follows.
\begin{equation*} \begin {array} {c} x(t) = a \cos(t) \\ y(t) = b \sin(t) \end {array}, 0\le t \le 2\pi \end{equation*}
Using Green’s Theorem, we compute the area of the ellipse as a line integral of \(\vF = \langle 0,x \rangle\) around the boundary of the ellipse. (Any of the other vector fields listed above would work just as well.)
\begin{align*} \iint_E 1\,dA \amp = \oint_{\partial E} x\,dy\\ \amp = \int_{t=0}^{2\pi} (a\cos(t))(b\sin(t))\,dt\\ \amp = ab \int_0^{2\pi} \cos^2(t)\,dt\\ \amp = \frac{ab}{2} \int_0^{2t} 1 + \cos(2t)\,dt\\ \amp = \frac{ab}{2} \left[ t + \frac{1}{2} \sin(2t) \right]_0^{2\pi}\\ \amp = \frac{ab}{2} (2\pi) = ab\pi \end{align*}

Subsection 4.6.5 Summary

  • Green’s Theorem tells us that we can calculate the circulation of a smooth vector field along a simple closed curve that bounds a region in the plane on which the vector field is also smooth by calculating the double integral of the circulation density instead of the line integral.
  • Integrating the circulation density of a smooth vector field on a region bounded by a simple closed curve gives the same value as calculating the circulation of the vector field around the boundary of the region with suitable orientation.

Exercises 4.6.6 Exercises

1.

Sometimes, the value in Green’s Theorem is in converting a double integral into a line integral. Recall that the area of a region \(R\) in the plane can be found by calculating \(\iint_R 1\, dA\text{.}\)
(a)
Find a smooth vector field \(\vF\) such that the circulation density of \(\vF\) is \(1\) everywhere.
Hint.
There are several ways to choose \(\vF = \langle F_1,F_2\rangle\) so that \(\partial F_2/\partial x - \partial F_1/\partial y = 1\text{.}\)
(c)
For positive real numbers \(a,b\text{,}\) the ellipse \(x^2/a^2 + y^2/b^2 = 1\) can be parametrized as \(\vr(t) = \langle a\cos(t),b\sin(t)\rangle\) with \(0\leq t\leq 2\pi\text{.}\) Find the area of this ellipse by calculating a line integral.
(d)
Find another vector field \(\vG\) (different from the one you found in part a) that has circulation density \(1\) everywhere, and calculate the area of the ellipse in part c by calculating the circulation of \(\vG\) along the ellipse.
(e)
How many possible vector fields are there with circulation density \(1\) everywhere? Why?

2.

The types of regions in \(\R^2\) to which Green’s Theorem applies are formally called simply connected regions. To be precise, a simply connected region \(D\) in \(\R^2\) is a set so that there is a path between every pair of points in \(D\) that stays inside \(D\) and any simple closed curve in \(D\) can be shrunk to a point while remaining inside \(D\text{.}\) We can think of a simply connected region as being a region that does not have any “holes”. There are many instances where we can find a way to apply Green’s Theorem multiple times to work with regions that are not simply connected, however.
Figure 4.6.8. A region that is not simply connected
Figure 4.6.9. Breaking up a region
  1. Explain why the region \(R\) in Figure 4.6.8 is not simply connected.
  2. If we let \(C = C_1 + C_2\text{,}\) then the orientation of \(C\) is as required by Green’s Theorem—when walking along the pieces of \(C\text{,}\) the region \(R\) is always on the left-hand side. Assume that \(\vF\) is a vector field that is smooth on \(R\) and \(C\text{.}\) Using Figure 4.6.9 as a guide, write \(\oint_C \vF\cdot d\vr\) as a sum of double integrals, one over \(R_1\) and the other over \(R_2\text{.}\)
  3. In Activity 4.6.4, we considered the vector field \(\vF = \displaystyle\frac{-y}{{x^2+y^2}}\vi + \frac{x}{{x^2+y^2}}\vj\) and found that the circulation of \(\vF\) around the unit circle centered at the origin (oriented counterclockwise) was \(2\pi\text{.}\) Show that for every simple closed curve \(C\) in the plane that bounds a region containing the origin, \(\oint_C\vF\cdot d\vr = 2\pi\text{.}\)
Hint.
For part c, create a region \(R\) as in Figure 4.6.8 using the curve \(C\) in place of \(C_1\text{.}\)

3.

This exercise presents another occasion where Green’s Theorem can be used to convert a double integral into a line integral.
  1. Recall from Section 3.4 that the centroid of a lamina \(D\) of area \(A\) is given by
    \begin{equation*} \overline{x} = \frac{1}{A}\iint_D x\, dA \qquad \text{and}\qquad \overline{y} = \frac{1}{A}\iint_D y\, dA\text{.} \end{equation*}
    Find vector fields \(\vF\) and \(\vG\) so that
    \begin{equation*} \overline{x} = \oint_{C} \vF\cdot d\vr \qquad \text{and}\qquad \overline{y} = \oint_C \vG\cdot d\vr\text{,} \end{equation*}
    where \(C\) is the boundary of the lamina \(D\text{.}\)
  2. Find the centroid of the triangle with vertices \((0,0)\text{,}\) \((a,0)\text{,}\) and \((0,b)\) for real numbers \(a,b>0\text{.}\)