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 \(\mathbb{R}^2\text{:}\) 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 \(\mathbb{R}\text{.}\) Let \(f\) be a zero-form on \(U \subset \mathbb{R}\text{.}\) Let \([a,b] \in U\) with the canonical orientation of increasing numbers. The boundary of the interval, with its induced orientation, is \(\partial [a,b] = \{(b,+), (a,-) \}\text{.}\) Then the integral of the exact one-form \(\omega = df\) along \([a,b]\) can be evaluated as:

\begin{equation*} \int_{[a,b]} d f = \int_{\partial [a,b]} f = f(b) - f(a). \end{equation*}

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 \(\mathbb{R}^2\text{.}\)

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 \(\eta = f\ dx + g\ dy\) and \(d \eta = \left( \frac{\partial g}{\partial x} - \frac{\partial f}{\partial y} \right) dx \wedge dy\text{.}\)

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

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2 ~|~ x \in [a,b], u(x) \leq y \leq v(x) \}, \end{equation*}

for some continuous functions \(u,v: \mathbb{R} \to \mathbb{R}\text{.}\) The boundary \(\partial D\) is a closed simple curve, with counterclockwise orientation. It can be split into four curves:

  1. \(C_1\) is the vertical line \(x=b\) between \(y=u(b)\) and \(y=v(b)\text{;}\)

  2. \(C_2\) is the curve \(y = v(x)\) between \(x=b\) and \(x=a\text{;}\)

  3. \(C_3\) is the vertical line \(x=a\) between \(y=v(a)\) and \(y=u(a)\text{;}\)

  4. \(C_4\) is the curve \(y=u(x)\) between \(x=a\) and \(x=b\text{.}\)

In fact, for \(C_2\) and \(C_3\text{,}\) 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 (\(C_2\) and \(C_3\) have opposite orientations):

  1. \(\alpha_1: [0,1] \to \mathbb{R}^2\) with \(\alpha_1(t) = (b, u(b)+(v(b)-u(b))t)\text{;}\)

  2. \(\alpha_2: [a,b] \to \mathbb{R}^2\) with \(\alpha_2(t) = (t, v(t) )\text{;}\)

  3. \(\alpha_3: [0,1] \to \mathbb{R}^2\) with \(\alpha_3(t) = (a, u(a)+(v(a)-u(a))t)\text{;}\)

  4. \(\alpha_4: [a,b] \to \mathbb{R}^2\) with \(\alpha_4(t) = (t, u(t) )\text{;}\)

Calculating the pullbacks, we evaluate the line integrals:

\begin{align*} \int_{\alpha_1} \eta =\amp \int_{0}^{1} g(b,u(b)+(v(b)-u(b))t) (v(b)-u(b))\ dt,\\ \int_{\alpha_2} \eta =\amp \int_a^b \left( f(t,v(t) ) + g(t,v(t) ) v'(t) \right)\ dt, \\ \int_{\alpha_3} \eta =\amp \int_{0}^{1} g(a,u(a) + (v(a)-u(a))t ) (v(a)-u(a)) \ dt, \\ \int_{\alpha_4} \eta =\amp \int_a^b \left( f(t,u(t) ) + g(t,u(t) ) u'(t) \right)\ dt. \end{align*}

Putting those together, remembering that we need to add a negative sign for \(C_2\) and \(C_3\text{,}\) we get:

\begin{align*} \int_{\partial D} \eta =\amp \int_{\alpha_1} \eta - \int_{\alpha_2} \eta-\int_{\alpha_3} \eta+\int_{\alpha_4} \eta \\ =\amp \int_0^1 \left( 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)) \right) dt \\ \amp+ \int_a^b \left( f(t,u(t) ) - f(t,v(t)) + g(t,u(t) ) u'(t) - g(t,v(t) ) v'(t) \right)\ dt . \end{align*}

Next, we need to evaluate the integral \(\int_D d \eta\text{.}\) 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 \(\phi:D_2 \to D\text{,}\) where

\begin{equation*} D_2 = \{ (s,t) \in \mathbb{R}^2~|~ s \in [a,b], t \in [0,1] \}, \end{equation*}

with

\begin{equation*} \phi(s,t) = (s, u(s) + (v(s)-u(s))t ). \end{equation*}

We see that \(\phi(D_2) = D\text{,}\) and that the function is bijective. So by Lemma 5.3.7, we know that

\begin{equation*} \int_{D_2} \phi^* (d \eta) = \int_{D_1} d \eta, \end{equation*}

so we may as well calculate the left-hand-side. In fact, we also know that \(\phi^* ( d \eta) = d(\phi^* \eta)\text{.}\) So let us calculate this pullback. First,

\begin{equation*} \phi^* \eta = f(\phi(s,t))\ ds + g(\phi(s,t) ) (u'(s) + (v'(s)-u'(s))t)\ ds + g(\phi(s,t) ) (v(s)-u(s))\ dt. \end{equation*}

Then

\begin{align*} d(\phi^*\eta) =\amp \left( - \frac{\partial f(\phi(s,t))}{\partial t} - \frac{\partial}{\partial t} \left( g(\phi(s,t) ) (u'(s) + (v'(s)-u'(s))t) \right) \right.\\ \amp \left. + \frac{\partial}{\partial s} \left( g(\phi(s,t) ) (v(s)-u(s)) \right) \right) ds \wedge dt. \end{align*}

Thus

\begin{align*} \int_{D_2} \phi^*(d \eta)= \amp \int_a^b \int_0^1 \left( - \frac{\partial f(\phi(s,t))}{\partial t} - \frac{\partial}{\partial t} \left( g(\phi(s,t) ) (u'(s) + (v'(s)-u'(s))t) \right) \right) dt ds \\ \amp + \int_0^1 \int_a^b \left( \frac{\partial}{\partial s} \left( g(\phi(s,t) ) (v(s)-u(s)) \right) \right)ds dt, \end{align*}

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:

\begin{align*} \int_{D_2} \phi^*(d \eta)= \amp \int_a^b \left( f(\phi(s,0)) - f(\phi(s,1) ) + g(\phi(s,0) ) u'(s) - g(\phi(s,1)) v'(s) \right) ds\\ \amp + \int_0^1 \left( g(\phi(b,t) ) (v(b)-u(b) ) - g(\phi(a,t) ) (v(a)-u(a) ) \right)\ dt. \end{align*}

Substituting back the expression for \(\phi(s,t)\text{,}\) we get:

\begin{align*} \int_{D_2} \phi^*(d \eta)= \amp \int_a^b \left( f(s, u(s)) - f(s,v(s)) + g(s,u(s)) u'(s) - g(s,v(s) ) v'(s) \right)\ ds\\ \amp+ \int_0^1 \left( 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) ) \right)\ dt. \end{align*}

Magic: this is exactly the same as the expression that we found many lines above for the line integral \(\int_{\partial D} \eta\text{!}\) Woot woot! Therefore

\begin{equation*} \int_{D} d\eta = \int_{\partial D} \eta. \end{equation*}

Remark 5.7.2.

If \(D\) is simply connected, then \(\partial D\) is a simple closed curve, and \(D\) is the closed curve with its interior. In this case the line integral \(\int_{\partial D}\eta\) is simply the line integral of the one-form \(\eta\) 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 \(\partial D\) would have two components: the inner and outer circle. The line integral \(\int_{\partial D} \eta\) would then be the sum of the line integrals of the one-form \(\eta\) 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 \eta\) over the region \(D\subset \mathbb{R}^2\) with canonical orientation is equal to the integral of the one-form \(\eta\) over the boundary \(\partial D\) with the induced orientation.

  2. The integral of the one-form \(\eta\) over the simple closed curve \(\partial D \subset \mathbb{R}^2\) with canonical orientation is equal to the integral of its exterior derivative \(d \eta\) 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 \(\omega = x y\ dx +(x+y)\ dy\) over the rectangle with vertices \((-2,-1), (2,-1), (2,0), (-2, 0) \text{,}\) 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\text{.}\) It is the boundary of the region \(D\) consisting of the interior of the rectangle with its boundary:

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2~|~ x\in [-2,2], y\in [-1,0] \}. \end{equation*}

Thus we can use Green's theorem to rewrite the line integral along \(C\) as a surface integral over \(D\text{.}\)

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:

\begin{equation*} \int_{C_-} \omega = - \int_{C_+} \omega = - \int_D d \omega. \end{equation*}

So instead of calculating the line integral in the problem, we can calculate \(-\int_D d \omega\text{.}\)

The two-form \(d \omega\) is:

\begin{equation*} d \omega = x dy \wedge dx + dx \wedge dy = (1-x) dx \wedge dy. \end{equation*}

The integral can be evaluated:

\begin{align*} \int_D d \omega =\amp \int_{-1}^0 \int_{-2}^2 (1-x) \ dx dy \\ =\amp \int_{-1}^0 \left[x - \frac{x^2}{2} \right]_{-2}^2 \ dy\\ =\amp \int_{-1}^0 (2 -(-2) - 2 + 2 )\ dy \\ =\amp 4 (0-(-1))\\ =\amp 4. \end{align*}

Therefore, the line integral of \(\omega\) along \(C_-\) is

\begin{equation*} \int_{C_-} \omega = - \int_{D} d \omega = - 4. \end{equation*}

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

\begin{equation*} \frac{x^2}{a^2} + \frac{y^2}{b^2} = 1. \end{equation*}

The area is given by integrating the basic two-form \(\omega = dx \wedge dy\) 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 \(\omega\) is exact. But it certainly is: for instance, we can write \(\omega = d \eta\) with \(\eta = \frac{1}{2} (x \ dy - y\ dx)\) (we could choose other one-forms such that \(d \eta =\omega\) as well, but this one gives a nice and easy integral). Then Green's theorem states:

\begin{equation*} \int_D \omega = \int_C \eta, \end{equation*}

where on the right-hand-side we are evaluating the line integral of \(\eta = \frac{1}{2}(x \ dy-y \ dx)\) along the ellipse with counterclockwise orientation.

We can parametrize the ellipse as \(\alpha: [0,2 \pi] \to \mathbb{R}^2\) with \(\alpha(t) = (a \cos(t), b \sin(t) )\text{.}\) The orientation is counterclockwise as required. The pullback \(\alpha^* \eta\) is:

\begin{equation*} \alpha^* \eta = \left( \frac{a b}{2} \cos^2(t) - \frac{a b}{2} \sin(t) (- \sin(t) ) \right)\ dt= \frac{ab}{2} dt. \end{equation*}

The line integral is then

\begin{align*} \int_C \eta =\amp \int_0^{2 \pi} \frac{ab}{2}\ dt\\ =\amp \pi a b, \end{align*}

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 \(\mathbb{R}^3\) where we embed it trivially in the \(xy\)-plane. That is, \(\alpha: D \to \mathbb{R}^3\text{,}\) with \(\alpha(x,y) = (x,y,0)\text{.}\) The tangent vectors are \(\mathbf{T}_x = (1,0,0)\) and \(\mathbf{T}_y = (0,1,0)\text{,}\) and the normal vector is of course \(\mathbf{n} = (0,0,1) = \mathbf{e}_3\text{.}\)

To the one-form \(\eta = f\ dx + g\ dy\text{,}\) we associate the vector field \(\mathbf{F} = (f,g,0)\text{.}\) To the two-form \(\omega = d \eta\) is then associated the vector field \(\boldsymbol{\nabla} \times \mathbf{F}\text{.}\) Using Corollary 5.6.5, we then see that we can rewrite the surface integral as

\begin{equation*} \int_D d \eta = \iint_D \left( (\boldsymbol{\nabla} \times \mathbf{F}) \cdot \mathbf{e}_3 \right)\ dA. \end{equation*}

In other words, the integrand is the \(z\)-component of the curl \(\boldsymbol{\nabla} \times \mathbf{F}\text{.}\)

As for the other side of Green's theorem, it is a line integral of the one-form \(\eta\) over the curve \(\partial D\text{.}\) Assuming that we parametrize the curve with a position function \(\mathbf{r}\text{,}\) we can rewrite the line integral as

\begin{equation*} \int_{\partial D} \eta = \int_{\partial D} \mathbf{F} \cdot d \mathbf{r}. \end{equation*}

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

\begin{equation*} \iint_D \left( (\boldsymbol{\nabla} \times \mathbf{F}) \cdot \mathbf{e}_3 \right)\ dA = \int_{\partial D} \mathbf{F} \cdot d \mathbf{r}. \end{equation*}

Exercises 5.7.3 Exercises

1.

Evaluate the line integral of the one-form

\begin{equation*} \omega = e^y\ dx + e^x\ dy \end{equation*}

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

Solution.

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

The rectangle.
Figure 5.7.6. The path \(\partial D\) along the rectangle.
To evaluate the line integral directly, we need to parametrize separately the four line segments. We write the parametrizations as \(\alpha_i : [0,1] \to \mathbb{R}\text{,}\) \(i=1,\ldots,4\text{,}\) with

\begin{align*} \alpha_1(t) =\amp (3t,0), \\ \alpha_2(t) =\amp (3,2t),\\ \alpha_3(t)=\amp (3-3t,2),\\ \alpha_4(t)=\amp (0,2-2t). \end{align*}

We calculate the pullbacks:

\begin{align*} \alpha_1^* \omega =\amp 3\ dt,\\ \alpha_2^* \omega =\amp 2 e^3\ dt,\\ \alpha_3^* \omega =\amp - 3 e^2 \ dt,\\ \alpha_4^* \omega =\amp -2\ dt. \end{align*}

Putting all this together, the line integral becomes

\begin{align*} \int_{\partial D} \omega=\amp \int_{\alpha_1} \omega +\int_{\alpha_2} \omega +\int_{\alpha_3} \omega +\int_{\alpha_4} \omega \\ =\amp \int_0^1 \left( 3 + 2 e^3 - 3 e^2 - 2 \right)\ dt\\ =\amp 1 + 2 e^3 - 3 e^2. \end{align*}

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

\begin{equation*} \int_{\partial D} \omega = \int_D d \omega, \end{equation*}

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

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2\ |\ x \in [0,3], y \in[0,2] \}, \end{equation*}

with canonical orientation. We calculate the exterior derivative:

\begin{align*} d\omega =\amp e^y\ dy \wedge dx + e^x\ dx \wedge dy\\ =\amp (e^x-e^y) \ dx \wedge dy. \end{align*}

The integral of the two-form is:

\begin{align*} \int_D d\omega =\amp \int_D (e^x-e^y)\ dx \wedge dy\\ =\amp \int_0^2 \int_0^3 (e^x-e^y) dx dy\\ =\amp \int_0^2 \left(e^3-1-3 e^y \right)\ dy\\ =\amp 2 e^3 - 2 - 3 e^2 + 3\\ =\amp 1 + 2 e^3 - 3 e^2. \end{align*}

Therefore,

\begin{equation*} \int_{\partial D} \omega = 1 + 2 e^3 - 3 e^2, \end{equation*}

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 \(\mathbf{F}(x,y) = (x^{1/3}+y, x^{4/3}-y^{1/3}) \) along the triangle with vertices \((1,3), (2,2), (1,1)\text{,}\) clockwise.

Solution.

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

The triangle.
Figure 5.7.7. The path \(\partial D\) along the triangle.
We associate to the vector field \(\mathbf{F}\) the one-form \(\omega = (x^{1/3}+y)\ dx + (x^{4/3}-y^{1/3})\ dy\text{.}\) Since the path \(\partial D \) is oriented clockwise (negative orientation), by Green's theorem we know that

\begin{equation*} \int_{\partial D} \omega = - \int_D d \omega, \end{equation*}

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\text{;}\) the green line, \(y=x\text{;}\) the orange line, \(y=-x+4\text{.}\) Using these equations, we can describe \(D\) as an \(x\)-supported region:

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2\ |\ x \in [1,2], x \leq y \leq -x+4 \}. \end{equation*}

We calculate the exterior derivative \(d \omega\text{:}\)

\begin{align*} d \omega =\amp dy \wedge dx + \frac{4}{3} x^{1/3} dx \wedge dy\\ =\amp \left( \frac{4}{3} x^{1/3} - 1 \right) dx \wedge dy. \end{align*}

Its integral along \(D\) is:

\begin{align*} \int_D d \omega =\amp \int_1^2 \int_x^{-x+4} \left( \frac{4}{3} x^{1/3} - 1 \right)\ dy dx\\ =\amp \int_1^2 \left( \frac{4}{3} x^{1/3} - 1 \right) \Big[ y \Big]_{x}^{-x+4}\ dx\\ =\amp \int_1^2 \left( \frac{4}{3} x^{1/3} - 1 \right) (-x+4-x)\ dx\\ =\amp \int_1^2 \left( - \frac{8}{3} x^{4/3} + 2 x + \frac{16}{3} x^{1/3} - 4 \right)\ dx\\ =\amp \left( - \frac{8}{7} x^{7/3} + x^2 + 4 x^{4/3} - 4 x \right)_{x=1}^{x=2}\\ =\amp - \frac{32}{7} 2^{1/3} + 4 + 8 \cdot 2^{1/3} - 8 + \frac{8}{7} - 1 - 4 + 4\\ =\amp \frac{3}{7} \left( 8 \cdot 2^{1/3} - 9 \right). \end{align*}

Therefore,

\begin{equation*} \int_{\partial D} \omega = - \int_D d \omega = - \frac{3}{7} \left( 8 \cdot 2^{1/3} - 9 \right). \end{equation*}

3.

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

\begin{equation*} \mathbf{F}(x,y) = (x^2 + y^2, x y) \end{equation*}

while moving an object first along the parabola \(y=x^2\) from the origin to the point \((2,4)\text{,}\) 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 \(\partial D\text{.}\) It is shown in the figure below.

The path.
Figure 5.7.8. The closed path \(\partial D\text{.}\)
From the description of the path we know that the object is moving counterclockwise along the path. Therefore, if we associate the one-form \(\omega = (x^2+y^2)\ dx + x y\ dy\) to the vector field \(\mathbf{F}\text{,}\) by Green's theorem we know that

\begin{equation*} \int_{\partial D} \omega = \int_D d \omega, \end{equation*}

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:

\begin{equation*} D = \{(x,y) \in \mathbb{R}^2\ |\ x \in [0,2], x^2 \leq y \leq 4 \}. \end{equation*}

The exterior derivative \(d \omega\) is:

\begin{align*} d \omega =\amp 2 y\ dy \wedge dx + y\ dx \wedge dy\\ =\amp - y\ dx \wedge dy \end{align*}

Its integral along \(D\) is:

\begin{align*} \int_D d \omega =\amp - \int_D y \ dx \wedge dy\\ =\amp - \int_0^2 \int_{x^2}^4 y\ dy dx\\ =\amp - \int_0^2 \left[ \frac{y^2}{2} \right]_{y=x^2}^{y=4}\ dx\\ =\amp - \int_0^2 \left( 8 - \frac{x^4}{2} \right)\ dx\\ =\amp - \left( 16 - \frac{16}{5} \right)\\ =\amp -\frac{64}{5}. \end{align*}

Therefore, by Green's theorem, the work done by the force \(\mathbf{F}\) while moving an object along the path \(\partial D\) is:

\begin{equation*} \int_{\partial D} \omega = -\frac{64}{5}. \end{equation*}

4.

Suppose that a polygon has vertices \((x_1,y_1), (x_2, y_2), \ldots, (x_n,y_n)\text{,}\) in counterclockwise order. The area \(A\) of the polygon is given by integrating the basic two-form \(\omega = dx \wedge dy\) over the polygon. Use Green's theorem to show that

\begin{equation*} A = \frac{1}{2} \left( (x_1 y_2 - x_2 y_1) + (x_2 y_3 - x_3 y_2) + \ldots + (x_{n-1} y_{n} - x_n y_{n-1}) + (x_n y_1 - x_1 y_n) \right). \end{equation*}

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 \(\partial D\) the closed path that goes around the polygon counterclockwise. It consists of \(n\) line segments \(L_i\text{,}\) \(i=1,\ldots, n\text{,}\) with \(L_i\) being the line segment between the points \((x_i,y_i)\) and \((x_{i+1},y_{i+1}) \) (for simplicity of notation, we define \((x_{n+1}, y_{n+1}) := (x_1, y_1)\)).

Further, we notice that the two \(\omega = dx \wedge dy\) is exact. We write it as \(\omega = d \eta\) for the one-form

\begin{equation*} \eta = x \ dy. \end{equation*}

Note that this is certainly not the only choice of \(\eta\text{,}\) but any \(\eta\) such that \(d \eta = \omega\) will work.

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

\begin{align*} A =\amp \int_D dx \wedge dy \\ =\amp \int_{\partial D} \eta\\ =\amp \sum_{i=1}^n \int_{\alpha_i} \eta, \end{align*}

where in the last line \(\alpha_i\) is a parametrization of the line segment \(L_i\text{.}\) More explicitly, since the line segment \(L_i\) joins the points \((x_i, y_i)\) and \((x_{i+1}, y_{i+1})\text{,}\) we can parametrize it as \(\alpha_i: [0,1] \to \mathbb{R}^2\) with

\begin{equation*} \alpha_i(t) = (x_i + t (x_{i+1}-x_i), y_i + t (y_{i+1}-y_i) ). \end{equation*}

Then

\begin{equation*} \alpha_i^* \eta = (x_i + t(x_{i+1}-x_i) ) (y_{i+1}-y_i)\ dt. \end{equation*}

The line integral along \(\alpha_i\) is then:

\begin{align*} \int_{\alpha_i} \eta =\amp \int_0^1 (x_i + t(x_{i+1}-x_i) ) (y_{i+1}-y_i)\ dt\\ =\amp (y_{i+1}-y_i) \left( x_i + \frac{1}{2} (x_{i+1}-x_i) \right)\\ =\amp \frac{1}{2} (y_{i+1}-y_i) (x_i + x_{i+1}) \\ =\amp \frac{1}{2} ( x_i y_{i+1} - x_{i+1} y_i + x_{i+1} y_{i+1} - x_i y_i ). \end{align*}

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

\begin{align*} A =\amp \sum_{i=1}^n \int_{\alpha_i} \eta\\ =\amp \frac{1}{2} \sum_{i=1}^n ( x_i y_{i+1} - x_{i+1} y_i + x_{i+1} y_{i+1} - x_i y_i )\\ =\amp \frac{1}{2} ( (x_1 y_2 - x_2 y_1 + x_2 y_2 - x_1 y_1) + (x_2 y_3 - x_3 y_2 + x_3 y_3 - x_2 y_2) + \ldots \\ \amp + (x_n y_1 - x_1 y_n + x_1 y_1 - x_n y_n) ), \end{align*}

where in the last line we used the fact that \(x_{n+1} = x_1\) and \(y_{n+1} = y_1\text{.}\) We see that the terms \(x_1 y_1, x_2 y_2, \ldots\) all cancel out, and we are left with

\begin{equation*} A = \frac{1}{2} ((x_1 y_2 - x_2 y_1) + (x_2 y_3 - x_3 y_2) + \ldots + (x_n y_1 - x_1 y_n) ). \end{equation*}

Bingo!

5.

We already studied the one-form

\begin{equation*} \eta = -\frac{y}{x^2+y^2}\ dx + \frac{x}{x^2+y^2}\ dy, \end{equation*}

defined on \(U=\mathbb{R}^2 \setminus \{(0,0) \}\text{,}\) 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 \(\eta\) 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 \(C_0\) with canonical orientation that does not pass through or enclose the origin. Use Green's theorem to show that

    \begin{equation*} \int_{C_0} \eta = 0. \end{equation*}

  2. Let \(C_1\) 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 \(\omega\) along \(C_1\text{.}\)

  3. As in (b), let \(C_1\) 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 \(C_1\text{.}\) Give \(K\) a counterclockwise orientation. Use Green's theorem to show that

    \begin{equation*} \int_{C_1} \eta = \int_K \eta. \end{equation*}

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

    \begin{equation*} \int_{C_1} \eta = 2 \pi. \end{equation*}

You have showed that the line integral of \(\eta\) along an arbitrary simple closed curve that encloses the origin is non-vanishing (and is in fact equal to \(2 \pi\)), 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 \(\eta\text{.}\) The one-form is define on the open set \(U = \mathbb{R}^2 \setminus \{(0,0) \}\text{.}\)

As a starting point, recall that the one-form \(\eta\) is closed, that is, \(d \omega = 0\text{.}\) Indeed,

\begin{align*} d \eta =\amp \left(-\frac{x^2+y^2 - 2 y^2}{(x^2+y^2)^2} \right)dy \wedge dx + \left( \frac{x^2+y^2 - 2 x^2}{(x^2+y^2)^2} \right) dx \wedge dy\\ =\amp \frac{(x^2-y^2) + (y^2 - x^2)}{(x^2+y^2)^2}\ dx \wedge dy\\ =\amp 0. \end{align*}

This will be useful in this problem.

(a) We assume that \(C_0\) is a simple closed curve that does not pass through or enclose the origin, with canonical orientation. Therefore, \(C_0 \subset U\text{.}\) Moreover, if we denote by \(D_0\) the region consisting of the closed curve \(C_0\) and its interior, then \(D_0 \subset U\) as well, since the origin is not in \(D_0\text{.}\) Therefore, by Green's theorem,

\begin{align*} \int_{C_0} \eta =\amp \int_{D_0} d \eta\\ =\amp 0, \end{align*}

since \(d \eta = 0\text{.}\)

(b) If \(C_1\) is a simple closed curve that encloses the origin, we cannot apply the argument in (a). The reason is that, if we denote by \(D_1\) the region consistin of the closed curve \(C_1\) and its interior, then \(D_1\) is not a subset of \(U\text{,}\) since \(D_1\) contains the origin. In other words, the one-form \(\eta\) is not defined on \(D_1\text{,}\) 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 \(\eta\) along \(C_1\) vanishes, which is not true, as we will see.)

(c) We suppose that \(C_1\) 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 \(C_1\text{,}\) also oriented counterclockwise. For instance, \(C_1\) could be the orange curve in the figure below, and \(K\) the blue circle.

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

\begin{align*} \int_{D_1} d \eta =\amp \int_{\partial D_1} \eta\\ =\amp \int_{(C_1)_+} \eta + \int_{K_-} \eta\\ =\amp \int_{(C_1)_+} \eta - \int_{K_+} \eta, \end{align*}

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 \eta = 0\text{.}\) Therefore,

\begin{equation*} \int_{(C_1)_+} \eta = \int_{K_+} \eta, \end{equation*}

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 \(\eta\) along an arbitrary simple closed curve \(C_1\) that encloses the origin by the evaluation of the line integral of \(\eta\) 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\text{,}\) for some \(R \in \mathbb{R}_{>0}\text{,}\) then we can parametrize \(K\) by \(\alpha: [0,2\pi] \to \mathbb{R}^2\) with

\begin{equation*} \alpha(\theta) = (R \cos(\theta), R \sin(\theta) ). \end{equation*}

The pullback of \(\eta\) is

\begin{align*} \alpha^* \eta =\amp \frac{1}{R^2} \left( - R \sin(\theta) (- R \sin(\theta) ) + R \cos(\theta) (R \cos(\theta) ) \right)d \theta\\ =\amp d \theta. \end{align*}

Therefore,

\begin{align*} \int_{\alpha} \eta =\amp \int_{[0,2\pi]}\alpha^* \eta\\ =\amp \int_0^{2 \pi} d \theta\\ =\amp 2 \pi. \end{align*}

So the line integral of \(\eta\) along any circle \(K\) centered at the origin is always non-zero and equal to \(2 \pi\text{.}\) By (c), we thus conclude that

\begin{equation*} \int_{C_1} \eta = 2 \pi \end{equation*}

for arbitrary simple closed curves \(C_1\) that enclose the origin. Very powerful result!