Skip to main content
Mathematics LibreTexts

5.5: Green’s Theorem

  • Page ID
    119741
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)
    Learning Objectives
    • Apply the circulation form of Green’s theorem.
    • Apply the flux form of Green’s theorem.
    • Calculate circulation and flux on more general regions.

    In this section, we examine Green’s theorem, which is an extension of the Fundamental Theorem of Calculus to two dimensions. Green’s theorem has two forms: a circulation form and a flux form, both of which require region \(D\) in the double integral to be simply connected. However, we will extend Green’s theorem to regions that are not simply connected.

    Put simply, Green’s theorem relates a line integral around a simply closed plane curve \(C\) and a double integral over the region enclosed by \(C\). The theorem is useful because it allows us to translate difficult line integrals into more simple double integrals, or difficult double integrals into more simple line integrals.

    Extending the Fundamental Theorem of Calculus

    Recall that the Fundamental Theorem of Calculus says that

    \[\int_a^b F′(x)\,dx=F(b)−F(a). \nonumber \]

    As a geometric statement, this equation says that the integral over the region below the graph of \(F′(x)\) and above the line segment \([a,b]\) depends only on the value of \(F\) at the endpoints \(a\) and \(b\) of that segment. Since the numbers \(a\) and \(b\) are the boundary of the line segment \([a,b]\), the theorem says we can calculate integral \(\int_a^b F′(x)\,dx\) based on information about the boundary of line segment \([a,b]\) (Figure \(\PageIndex{1}\)). The same idea is true of the Fundamental Theorem for Line Integrals:

    \[\int_C \vecs \nabla f·d\vecs r=f(\vecs r(b))−f(\vecs r(a)). \nonumber \]

    When we have a potential function (an “antiderivative”), we can calculate the line integral based solely on information about the boundary of curve \(C\).

    A graph in quadrant 1 of a generic function f(x). It is an increasing concave up function for the first quarter, an increasing concave down function for the second quarter, a decreasing concave down function for the third quarter, and an increasing concave down function for the last quarter. In the second quarter, a point a is marked on the x axis, and in the third quarter, a point b is marked on the x axis. The area under the curve and between a and b is shaded. This area is labeled the integral from a to b of f(x) dx.
    Figure \(\PageIndex{1}\): The Fundamental Theorem of Calculus says that the integral over line segment \([a,b]\) depends only on the values of the antiderivative at the endpoints of \([a,b]\).

    Green’s theorem takes this idea and extends it to calculating double integrals. Green’s theorem says that we can calculate a double integral over region \(D\) based solely on information about the boundary of \(D\). Green’s theorem also says we can calculate a line integral over a simple closed curve \(C\) based solely on information about the region that \(C\) encloses. In particular, Green’s theorem connects a double integral over region \(D\) to a line integral around the boundary of \(D\).

    Circulation Form of Green’s Theorem

    The first form of Green’s theorem that we examine is the circulation form. This form of the theorem relates the vector line integral over a simple, closed plane curve \(C\) to a double integral over the region enclosed by \(C\). Therefore, the circulation of a vector field along a simple closed curve can be transformed into a double integral and vice versa.

    GREEN’S THEOREM (CIRCULATION FORM)

    Let \(D\) be an open, simply connected region with a boundary curve \(C\) that is a piecewise smooth, simple closed curve oriented counterclockwise (Figure \(\PageIndex{2}\)). Let \(\vecs F=⟨P,Q⟩\) be a vector field with component functions that have continuous partial derivatives on \(D\). Then,

    \[ \begin{align} \oint_C \vecs F·d\vecs r =\oint_C P\,dx+Q\,dy \\[4pt] =\iint_D (Q_x−P_y)\,dA. \end{align} \nonumber \]

    A vector field in two dimensions with all of the arrows pointing up and to the right. A curve C oriented counterclockwise sections off a region D around the origin. It is a simple, closed region.
    Figure \(\PageIndex{2}\): The circulation form of Green’s theorem relates a line integral over curve \(C\) to a double integral over region \(D\).

    Notice that Green’s theorem can be used only for a two-dimensional vector field \(\vecs F\). If \(\vecs F\) is a three-dimensional field, then Green’s theorem does not apply. Since

    \[\displaystyle \int_C P\,dx+Q\,dy=\int_C \vecs F·\vecs T\,ds \nonumber \]

    this version of Green’s theorem is sometimes referred to as the tangential form of Green’s theorem.

    The proof of Green’s theorem is rather technical, and beyond the scope of this text. Here we examine a proof of the theorem in the special case that \(D\) is a rectangle. For now, notice that we can quickly confirm that the theorem is true for the special case in which \(\vecs F=⟨P,Q⟩\) is conservative. In this case,

    \[\oint_C P\,dx+Q\,dy=0 \nonumber \]

    because the circulation is zero in conservative vector fields. \(\vecs F\) satisfies the cross-partial condition, so \(P_y=Q_x\). Therefore,

    \[\iint_D (Q_x−P_y)\,dA=\int_D 0\,dA=0=\oint_C P\,dx+Q\,dy \nonumber \]

    which confirms Green’s theorem in the case of conservative vector fields.

    Proof

    Let’s now prove that the circulation form of Green’s theorem is true when the region \(D\) is a rectangle. Let \(D\) be the rectangle \([a,b]×[c,d]\) oriented counterclockwise. Then, the boundary \(C\) of \(D\) consists of four piecewise smooth pieces \(C_1\), \(C_2\), \(C_3\), and \(C_4\) (Figure \(\PageIndex{3}\)). We parameterize each side of \(D\) as follows:

    \(C_1: \vecs r_1(t)=⟨t,c⟩\), \(a≤t≤b\)

    \(C_2: \vecs r_2(t)=⟨b,t⟩\), \(c≤t≤d\)

    \(−C_3: \vecs r_3(t)=⟨t,d⟩\), \(a≤t≤b\)

    \(−C_4: \vecs r_4(t)=⟨a,t⟩\), \(c≤t≤d\).

    A diagram in quadrant 1. Rectangle D is oriented counterclockwise. Points a and b are on the x axis, and points c and d are on the y axis with b > a and d > c. The sides of the rectangle are side c1 with endpoints at (a,c) and (b,c), side c2 with endpoints at (b,c) and (b,d), side c3 with endpoints at (b,d) and (a,d), and side c4 with endpoints at (a,d) and (a,c).
    Figure \(\PageIndex{3}\): Rectangle \(D\) is oriented counterclockwise.

    Then,

    \[\begin{align*} \int_C \vecs F·d \vecs r &=\int_{C_1} \vecs F·d \vecs r+\int_{C_2} \vecs F·d \vecs r+\int_{C_3} \vecs F·d \vecs r+\int_{C_4} \vecs F·d \vecs r \\[4pt] &=\int_{C_1} \vecs F·d \vecs r+\int_{C_2} \vecs F·d \vecs r−\int_{−C_3} \vecs F·d \vecs r−\int_{−C_4} \vecs F·d \vecs r \\[4pt] &=\int_a^b \vecs F( \vecs r_1(t))· \vecs r_1'(t)\,dt+\int_c^d \vecs F( \vecs r_2(t))· \vecs r_2'(t)\,dt−\int_a^b \vecs F( \vecs r_3(t))· \vecs r_3'(t)\,dt−\int_c^d \vecs F( \vecs r_4(t))·\vecs r_4'(t)\,dt\\[4pt] &=\int_a^b P(t,c)\,dt+\int_c^dQ(b,t)\,dt−\int_a^bP(t,d)\,dt−\int_c^dQ(a,t)\,dt \\[4pt] &=\int_a^b(P(t,c)−P(t,d))\,dt+\int_c^d(Q(b,t)−Q(a,t))\,dt\\[4pt] &=−\int_a^b(P(t,d)−P(t,c))\,dt+\int_c^d(Q(b,t)−Q(a,t))\,dt. \end{align*}\]

    By the Fundamental Theorem of Calculus,

    \[P(t,d)−P(t,c)=\int_c^d \dfrac{\partial}{\partial y}P(t,y)dy \nonumber \]

    and

    \[Q(b,t)−Q(a,t)=\int_a^b \dfrac{\partial}{\partial x} Q(x,t)\,dx. \nonumber \]

    Therefore,

    \[−\int_a^b(P(t,d)−P(t,c))\,dt+\int_c^d(Q(b,t)−Q(a,t))\,dt=−\int_a^b\int_c^d \dfrac{\partial}{\partial y} P(t,y)\,dy\,dt+\int_c^d\int_a^b \dfrac{\partial}{\partial x}Q(x,t)\,dx\,dt. \nonumber \]

    But,

    \[\begin{align*} −\int_a^b\int_c^d \dfrac{\partial}{\partial y}P(t,y)\,dy\,dt+\int_c^d\int_a^b \dfrac{\partial}{\partial x}Q(x,t)\,dx\,dt &=−\int_a^b\int_c^d \dfrac{\partial}{\partial y}P(x,y)\,dy\,dx+\int_c^d\int_a^b \dfrac{\partial}{\partial x}Q(x,y)\,dx\,dy \\[4pt] &=\int_a^b\int_c^d(Q_x−P_y)\,dy\,dx\\[4pt] &=\iint_D(Q_x−P_y)\,dA. \end{align*}\]

    Therefore, \(\displaystyle \int_C \vecs F\cdot d\vecs r=\iint_D(Q_x−P_y)\,dA\) and we have proved Green’s theorem in the case of a rectangle.

    \(\square\)

    To prove Green’s theorem over a general region \(D\), we can decompose \(D\) into many tiny rectangles and use the proof that the theorem works over rectangles. The details are technical, however, and beyond the scope of this text.

    Example \(\PageIndex{1}\): Applying Green’s Theorem over a Rectangle

    Calculate the line integral

    \[\oint_C x^2y\,dx+(y−3)\,dy, \nonumber \]

    where \(C\) is a rectangle with vertices \((1,1)\), \((4,1)\), \((4,5)\), and \((1,5)\) oriented counterclockwise.

    Solution

    Let \( \vecs F(x,y)=⟨P(x,y),Q(x,y)⟩=⟨x^2y,y−3⟩\). Then, \(Q_x(x,y)=0\) and \(P_y(x,y)=x^2\). Therefore, \(Q_x−P_y=−x^2\).

    Let \(D\) be the rectangular region enclosed by \(C\) (Figure \(\PageIndex{4}\)). By Green’s theorem,

    \[\begin{align*} \oint_C x^2y\,dx+(y−3)\,dy &=\iint_D (Q_x−P_y)\,dA \\[4pt] &=\iint_D −x^2 \,dA=\int_1^5\int_1^4−x^2\,dx\,dy \\[4pt] &= \int_1^5 -\frac{x^3}{3}\Bigg|_1^4\,dy\\[4pt] &= \int_1^5 \left( -\frac{64}{3} + \frac{1}{3} \right)  \,dy \\[4pt] &=\int_1^5−21\,dy=−84.\end{align*}\]

    A vector field in two dimensions with focus on quadrant 1. The arrows near the origin are short, and the arrows further away from the origin are longer. A rectangle has endpoints at (1,1), (4,1), (4,5), and (1,5). The arrows in quadrant 3 are pointing to the right. At the y axis, they split at y = 3. Arrows above that line curve up at the y axis and shift until they are horizontally pointing to the right in quadrant 1. Arrows below that line and above the x axis curve down at the y axis and shift until they are horizontally pointing to the right. Arrows below the x axis point to the left and down, pointing back to the y axis.
    Figure \(\PageIndex{4}\): The line integral over the boundary of the rectangle can be transformed into a double integral over the rectangle.
    Analysis

    If we were to evaluate this line integral without using Green’s theorem, we would need to parameterize each side of the rectangle, break the line integral into four separate line integrals, and use the methods from the section titled Line Integrals to evaluate each integral. Furthermore, since the vector field here is not conservative, we cannot apply the Fundamental Theorem for Line Integrals. Green’s theorem makes the calculation much simpler.

    Example \(\PageIndex{2}\): Applying Green’s Theorem to Calculate Work

    Calculate the work done on a particle by force field

    \[\vecs F(x,y)=⟨y+\sin x,e^y−x⟩ \nonumber \]

    as the particle traverses circle \(x^2+y^2=4\) exactly once in the counterclockwise direction, starting and ending at point \((2,0)\).

    Solution

    Let \(C\) denote the circle and let \(D\) be the disk enclosed by \(C\). The work done on the particle is

    \[W=\oint_C (y+\sin x)\,dx+(e^y−x)\,dy. \nonumber \]

    As with Example \(\PageIndex{1}\), this integral can be calculated using tools we have learned, but it is easier to use the double integral given by Green’s theorem (Figure \(\PageIndex{5}\)).

    Let \(\vecs F(x,y)=⟨P(x,y),Q(x,y)⟩=⟨y+\sin x,e^y−x⟩\). Then, \(Q_x=−1\) and \(P_y=1\). Therefore, \(Q_x−P_y=−2\).

    By Green’s theorem,

    \[\begin{align*} W &=\oint_C(y+\sin(x))dx+(e^y−x)\,dy \\[4pt] &=\iint_D (Q_x−P_y)\,dA \\[4pt] &=\iint_D−2\,dA \\[4pt] &=−2(area(D))=−2\pi (2^2)=−8\pi. \end{align*}\]

    A vector field in two dimensions. The arrows further away from the origin are much longer than those near the origin. The arrows curve out from about (.5,.5) in a clockwise spiral pattern.
    Figure \(\PageIndex{5}\): The line integral over the boundary circle can be transformed into a double integral over the disk enclosed by the circle.
    Exercise \(\PageIndex{2}\)

    Use Green’s theorem to calculate line integral

    \[\oint_C \sin(x^2)\,dx+(3x−y)\,dy. \nonumber \]

    where \(C\) is a right triangle with vertices \((−1,2)\), \((4,2)\), and \((4,5)\) oriented counterclockwise.

    Hint

    Transform the line integral into a double integral.

    Answer

    \(\dfrac{45}{2}\)

    In the preceding two examples, the double integral in Green’s theorem was easier to calculate than the line integral, so we used the theorem to calculate the line integral. In the next example, the double integral is more difficult to calculate than the line integral, so we use Green’s theorem to translate a double integral into a line integral.

    Example \(\PageIndex{3}\): Applying Green’s Theorem over an Ellipse

    Calculate the area enclosed by ellipse \(\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}=1\) (Figure \(\PageIndex{6}\)).

    A horizontal ellipse graphed in two dimensions. It has vertices at (-a, 0), (0, -b), (a, 0), and (0, b), where the absolute value of a is between 2.5 and 5 and the absolute value of b is between 0 and 2.5.
    Figure \(\PageIndex{6}\): Ellipse \(\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}=1\) is denoted by \(C\).
    Solution

    Let \(C\) denote the ellipse and let \(D\) be the region enclosed by \(C\). Recall that ellipse \(C\) can be parameterized by

    • \(x=a\cos t\),
    • \(y=b \sin t\),
    • \(0≤t≤2\pi\).

    Calculating the area of \(D\) is equivalent to computing double integral \(\iint_D \,dA\). To calculate this integral without Green’s theorem, we would need to divide \(D\) into two regions: the region above the x-axis and the region below. The area of the ellipse is

    \[\int_{−a}^a\int_0^{\sqrt{b^2−{(bx/a)}^2}} \,dy\,dx+\int_{−a}^{a} \int_{−\sqrt{b^2−{(bx/a)}^2}}^{0} \,dy\,dx. \nonumber \]

    These two integrals are not straightforward to calculate (although when we know the value of the first integral, we know the value of the second by symmetry). Instead of trying to calculate them, we use Green’s theorem to transform \(\iint_D \,dA\) into a line integral around the boundary \(C\).

    Consider vector field

    \[\vecs F(x,y)=⟨P,Q⟩=⟨−\dfrac{y}{2},\dfrac{x}{2}⟩. \nonumber \]

    Then, \(Q_x=\dfrac{1}{2}\) and \(P_y=−\dfrac{1}{2}\), and therefore \(Q_x−P_y=1\). Notice that \(\vecs F\) was chosen to have the property that \(Q_x−P_y=1\). Since this is the case, Green’s theorem transforms the line integral of \(\vecs F\) over \(C\) into the double integral of 1 over \(D\).

    By Green’s theorem,

    \[\begin{align*} \iint_D \,dA &=\iint_D (Q_x−P_y)\,dA \\[4pt] &=\int_C \vecs F\cdot d\vecs r=\dfrac{1}{2}\int_C −y\,dx+x\,dy \\[4pt] &=\dfrac{1}{2}\int_0^{2\pi}−b \sin t(−a\sin t)+a(\cos t)b\cos t\,dt \\[4pt] &=\dfrac{1}{2}\int_0^{2\pi} ab \cos^2 t+ab \sin^2 t\,dt \\[4pt] &=\dfrac{1}{2}\int_0^{2\pi} ab\,dt =\pi ab. \end{align*}\]

    Therefore, the area of the ellipse is \(\pi ab\;\text{units}^2\).

    In Example \(\PageIndex{3}\), we used vector field \(\vecs F(x,y)=⟨P,Q⟩=⟨−\dfrac{y}{2},\dfrac{x}{2}⟩\) to find the area of any ellipse. The logic of the previous example can be extended to derive a formula for the area of any region \(D\). Let \(D\) be any region with a boundary that is a simple closed curve \(C\) oriented counterclockwise. If \(\vecs F(x,y)=⟨P,Q⟩=⟨−\dfrac{y}{2},\dfrac{x}{2}⟩\), then \(Q_x−P_y=1\). Therefore, by the same logic as in Example \(\PageIndex{3}\),

    \[ \text{area of} \; D=\iint_D \,dA=\dfrac{1}{2}\oint_C−y\,dx+x\,dy. \label{greenarea} \]

    It’s worth noting that if \(\vecs F=⟨P,Q⟩\) is any vector field with \(Q_x−P_y=1\), then the logic of the previous paragraph works. So. Equation \ref{greenarea} is not the only equation that uses a vector field’s mixed partials to get the area of a region.

    Exercise \(\PageIndex{3}\)

    Find the area of the region enclosed by the curve with parameterization \(\vecs r(t)=⟨\sin t\cos t,\sin t⟩\), \(0≤t≤\pi\).

    An image of a curve in quadrants 1 and 2. The curve begins at the origin, curves up and to the right until about (.5, .8), curves to the left nearly horizontally, goes through (0,1), continues until about (-1, .7), and then curves down and to the right until it hits the origin again.

    Hint

    Use Equation \ref{greenarea}.

    Answer

    \(\dfrac{2}{3}\;\text{units}^2\).

    Solution

    For convenience, let's use the same vector field we used in Example \(\PageIndex{3}\), where \(Q_x−P_y=1\).
    \[\vecs F(x,y)=⟨P,Q⟩=⟨−\dfrac{y}{2},\dfrac{x}{2}⟩. \nonumber \]
    For the curve \(r(t)=⟨\sin t\cos t,\sin t⟩\) with \(0≤t≤\pi\), we have
    \[\begin{align*} \vecs r'(t) &= \langle \cos^2t - \sin^2t, \cos t \rangle \text{ and} \\[5pt]
    \vecs F(t) &= \left\langle -\frac{\sin t}{2}, \frac{\sin t \cos t}{2} \right\rangle \end{align*}\]
    By Green’s theorem,

    \[\begin{align*} \text{Area of the region }D =\iint_D \,dA &=\iint_D (Q_x−P_y)\,dA \\[4pt] &=\int_C \vecs F\cdot d\vecs r=\dfrac{1}{2}\int_C −y\,dx+x\,dy \\[4pt] &=\dfrac{1}{2}\int_0^{\pi}\big(−\sin t(\cos^2t - \sin^2t)+\sin t \cos^2 t \big)\,dt \\[4pt]
    &=\dfrac{1}{2}\int_0^{\pi} \big(−\sin t\cos^2t + \sin t\sin^2 t + \sin t\cos^2 t \big) \,dt \\[4pt]
    &=\dfrac{1}{2}\int_0^{\pi} \big(\cancel{−\sin t\cos^2t} + \sin t(1 - \cos^2 t) + \cancel{\sin t\cos^2 t} \big) \,dt \\[4pt]
    &=\dfrac{1}{2}\int_0^{\pi} \big( \sin t - \cos^2 t\sin t \big) \,dt \\[4pt]
    &= \frac{1}{2} \left[  -\cos t + \frac{\cos^3 t}{3} \right]\Bigg |_0^{\pi}\\[4pt]
    &= \frac{1}{2} \left[  -\cos \pi + \frac{\cos^3 \pi}{3} + \cos 0 - \frac{\cos^3 0}{3}  \right]\\[4pt]
    &= \frac{1}{2} \left[  1 - \frac{1}{3} + 1 -  \frac{1}{3} \right]\\[4pt]
    &= \frac{1}{2} \left[  2 - \frac{2}{3} \right]\\[4pt]
    &= \frac{1}{2} \left[ \frac{6}{3} - \frac{2}{3} \right]\\[4pt]
    &= \frac{1}{2} \left[ \frac{4}{3} \right]\\[4pt]
    &= \frac{2}{3} \;\text{units}^2 \end{align*}\]

    Therefore the area of this region is \(\dfrac{2}{3}\;\text{units}^2\).

    Flux Form of Green’s Theorem

    The circulation form of Green’s theorem relates a double integral over region \(D\) to line integral \(\displaystyle \oint_C \vecs F·\vecs Tds\), where \(C\) is the boundary of \(D\). The flux form of Green’s theorem relates a double integral over region \(D\) to the flux across boundary \(C\). The flux of a fluid across a curve can be difficult to calculate using the flux line integral. This form of Green’s theorem allows us to translate a difficult flux integral into a double integral that is often easier to calculate.

    GREEN’S THEOREM (FLUX FORM)

    Let \(D\) be an open, simply connected region with a boundary curve \(C\) that is a piecewise smooth, simple closed curve that is oriented counterclockwise (Figure \(\PageIndex{7}\)). Let \(\vecs F=⟨P,Q⟩\) be a vector field with component functions that have continuous partial derivatives on an open region containing \(D\). Then,

    \[\oint_C \vecs F·\vecs N\,ds=\iint_D P_x+Q_y\,dA. \label{GreenN} \]

    A vector field in two dimensions. A generic curve C encloses a simple region D around the origin oriented counterclockwise. Normal vectors N point out and away from the curve into quadrants 1, 3, and 4.
    Figure \(\PageIndex{7}\): The flux form of Green’s theorem relates a double integral over region \(D\) to the flux across curve \(C\).

    Because this form of Green’s theorem contains unit normal vector \(\vecs N\), it is sometimes referred to as the normal form of Green’s theorem.

    Proof

    Recall that \(\displaystyle \oint_C \vecs F·\vecs N\,ds=\oint_C −Q\,dx+P\,dy\). Let \(M=−Q\) and \(N=P\). By the circulation form of Green’s theorem,

    \[\begin{align*} \oint_C−Q\,dx+P\,dy &=\oint_C M\,dx+N\,dy\\[4pt] &=\iint_D N_x−M_y \,dA\\[4pt] &=\iint_D P_x−{(−Q)}_y \,dA\\[4pt] &=\iint_D P_x+Q_y \,dA. \end{align*}\]

    \(\square\)

    Example \(\PageIndex{4A}\): Applying Green’s Theorem for Flux across a Circle

    Let \(C\) be a circle of radius \(r\) centered at the origin (Figure \(\PageIndex{8}\)) and let \(\vecs F(x,y)=⟨x,y⟩\). Calculate the flux across \(C\).

    A vector field in two dimensions. The arrows point away from the origin in a radial pattern. They are shorter near the origin and much longer further away. A circle with radius 2 and center at the origin is drawn.
    Figure \(\PageIndex{8}\): Curve \(C\) is a circle of radius \(r\) centered at the origin.
    Solution

    Let \(D\) be the disk enclosed by \(C\). The flux across \(C\) is \(\displaystyle \oint_C \vecs F·\vecs N\,ds\). We could evaluate this integral using tools we have learned, but Green’s theorem makes the calculation much more simple. Let \(P(x,y)=x\) and \(Q(x,y)=y\) so that \(\vecs F=⟨P,Q⟩\). Note that \(P_x=1=Q_y\), and therefore \(P_x+Q_y=2\). By Green’s theorem,

    \[\int_C \vecs F\cdot\vecs N\,ds=\iint_D 2\,dA=2\iint_D \,dA. \nonumber \]

    Since \(\displaystyle \iint_D \,dA\) is the area of the circle, \(\displaystyle \iint_D \,dA=\pi r^2\). Therefore, the flux across \(C\) is \(2\pi r^2\).

    Example \(\PageIndex{4B}\): Applying Green’s Theorem for Flux across a Triangle

    Let \(S\) be the triangle with vertices \((0,0)\), \((1,0)\), and \((0,3)\) oriented clockwise (Figure \(\PageIndex{9}\)). Calculate the flux of \(\vecs F(x,y)=⟨P(x,y),Q(x,y)⟩=⟨x^2+e^y,x+y⟩\) across \(S\).

    A vector field in two dimensions. A triangle is drawn oriented clockwise with vertices at (0,0), (1,0), and (0,3). The arrows in the field point to the right and up slightly. The angle is greater the closer they are to the axis.
    Figure \(\PageIndex{9}\): Curve \(S\) is a triangle with vertices \((0,0)\), \((1,0)\), and \((0,3)\) oriented clockwise.
    Solution

    To calculate the flux without Green’s theorem, we would need to break the flux integral into three line integrals, one integral for each side of the triangle. Using Green’s theorem to translate the flux line integral into a single double integral is much more simple.

    Let \(D\) be the region enclosed by \(S\). Note that \(P_x=2x\) and \(Q_y=1\); therefore, \(P_x+Q_y=2x+1\). Green’s theorem applies only to simple closed curves oriented counterclockwise, but we can still apply the theorem because \(\displaystyle \oint_C \vecs F·\vecs N\,ds=−\oint_{−S} \vecs F·\vecs N\,ds\) and \(−S\) is oriented counterclockwise. By Green’s theorem, the flux is

    \[\begin{align*} \oint_C \vecs F·\vecs N\,ds &= \oint_{−S} \vecs F·\vecs N\,ds\\[4pt] &=−\iint_D (P_x+Q_y)\,dA \\[4pt] &=−\iint_D (2x+1)\,dA.\end{align*}\]

    Notice that the top edge of the triangle is the line \(y=−3x+3\). Therefore, in the iterated double integral, the \(y\)-values run from \(y=0\) to \(y=−3x+3\), and we have

    \[\begin{align*} −\iint_D (2x+1)\,dA &= −\int_0^1\int_0^{−3x+3}(2x+1)\,dy\,dx \\[4pt]
    &=−\int_0^1(2x+1)y\Big|_0^{-3x+3}\,dx \\[4pt]
    &=−\int_0^1(2x+1)(−3x+3)\,dx \\[4pt]
    &=−\int_0^1(−6x^2+3x+3)\,dx\\[4pt]
    &=−{\left[−2x^3+\dfrac{3x^2}{2}+3x\right]}\Bigg|_0^1 \\[4pt]
    &=−\dfrac{5}{2}. \end{align*}\]

    Exercise \(\PageIndex{4}\)

    Calculate the flux of \(\vecs F(x,y)=⟨x^3,y^3⟩\) across a unit circle oriented counterclockwise.

    Hint

    Apply Green’s theorem and use polar coordinates.

    Answer

    \(\dfrac{3\pi}{2}\)

    Example \(\PageIndex{5}\): Applying Green’s Theorem for Water Flow across a Rectangle

    Water flows from a spring located at the origin. The velocity of the water is modeled by vector field \(\vecs v(x,y)=⟨5x+y,x+3y⟩\) m/sec. Find the amount of water per second that flows across the rectangle with vertices \((−1,−2)\), \((1,−2)\), \((1,3)\),and \((−1,3)\), oriented counterclockwise (Figure \(\PageIndex{10}\)).

    A vector field in two dimensions. A rectangle is drawn oriented counterclockwise with vertices at (-1,3), (1,3), (-1,-2), and (1,-2). The arrows point out and away from the origin in a radial pattern. However, the arrows in quadrants 2 and 4 curve slightly towards the y axis instead of directly out. The arrows near the origin are short, and those further away from the origin are much longer.
    Figure \(\PageIndex{10}\): Water flows across the rectangle with vertices \((−1,−2)\), \((1,−2)\), \((1,3)\),and \((−1,3)\), oriented counterclockwise.
    Solution

    Let \(C\) represent the given rectangle and let \(D\) be the rectangular region enclosed by \(C\). To find the amount of water flowing across \(C\), we calculate flux \(\displaystyle \int_C \vecs v\cdot d\vecs r\). Let \(P(x,y)=5x+y\) and \(Q(x,y)=x+3y\) so that \(\vecs v=⟨P,Q⟩\). Then, \(P_x=5\) and \(Q_y=3\). By Green’s theorem,

    \[\begin{align*} \int_C \vecs v\cdot d\vecs r &=\iint_D (P_x+Q_y)\,dA \\ &=\iint_D 8\,dA \\ &=8(area\space of\space D)=80. \end{align*}\]

    Therefore, the water flux is 80 m2/sec.

    Recall that if vector field \(\vecs F\) is conservative, then \(\vecs F\) does no work around closed curves—that is, the circulation of \(\vecs F\) around a closed curve is zero. In fact, if the domain of \(\vecs F\) is simply connected, then \(\vecs F\) is conservative if and only if the circulation of \(\vecs F\) around any closed curve is zero. If we replace “circulation of \(\vecs F\)” with “flux of \(\vecs F\),” then we get a definition of a source-free vector field. The following statements are all equivalent ways of defining a source-free field \(\vecs F=⟨P,Q⟩\) on a simply connected domain (note the similarities with properties of conservative vector fields):

    1. The flux \( \displaystyle \oint_C \vecs F·\vecs N\,ds\) across any closed curve \(C\) is zero.
    2. If \(C_1\) and \(C_2\) are curves in the domain of \(\vecs F\) with the same starting points and endpoints, then \(\displaystyle \int_{C_1} \vecs F·\vecs N\,ds=\int_{C_2} \vecs F·\vecs N\,ds\). In other words, flux is independent of path.
    3. There is a stream function \(g(x,y)\) for \(\vecs F\). A stream function for \(\vecs F=⟨P,Q⟩\) is a function g such that \(P=g_y\) and \(Q=−g_x\).Geometrically, \(\vecs F=\langle a,b\rangle\) is tangential to the level curve of \(g\) at \((a,b)\). Since the gradient of \(g\) is perpendicular to the level curve of \(g\) at \((a,b)\), stream function \(g\) has the property \(\vecs F(a,b)\cdot\vecs\nabla g(a,b)=0\) for any point \((a,b)\) in the domain of \(g\). (Stream functions play the same role for source-free fields that potential functions play for conservative fields.)
    4. \(P_x+Q_y=0\)
    Example \(\PageIndex{6}\): Finding a Stream Function

    Verify that rotation vector field \(\vecs F(x,y)=⟨y,−x⟩\) is source free, and find a stream function for \(\vecs F\).

    Solution

    Note that the domain of \(\vecs F\) is all of \(ℝ^2\), which is simply connected. Therefore, to show that \(\vecs F\) is source free, we can show any of items 1 through 4 from the previous list to be true. In this example, we show that item 4 is true. Let \(P(x,y)=y\) and \(Q(x,y)=−x\). Then \(P_x=0=Q_y\), and therefore \(P_x+Q_y=0\). Thus, \(\vecs F\) is source free.

    To find a stream function for \(\vecs F\), proceed in the same manner as finding a potential function for a conservative field. Let \(g\) be a stream function for \(\vecs F\). Then \(g_y=y\), which implies that

    \(g(x,y)=\dfrac{y^2}{2}+h(x)\).

    Since \(−g_x=Q=−x\), we have \(h′(x)=x\). Therefore,

    \(h(x)=\dfrac{x^2}{2}+C\).

    Letting \(C=0\) gives stream function

    \(g(x,y)=\dfrac{x^2}{2}+\dfrac{y^2}{2}\).

    To confirm that \(g\) is a stream function for \(\vecs F\), note that \(g_y=y=P\) and \(−g_x=−x=Q\).

    Notice that source-free rotation vector field \(\vecs F(x,y)=⟨y,−x⟩\) is perpendicular to conservative radial vector field \(\vecs \nabla g=⟨x,y⟩\) (Figure \(\PageIndex{11}\)).

    Two vector fields in two dimensions. The first has arrows surrounding the origin in a clockwise circular pattern. The second has arrows pointing out and away from the origin in a radial manner. Circles with radii 1.5, 2, and 2.5 and centers at the origin are drawn in both. The arrows near the origin are shorter than those much further away. The first is labeled F(x,y) = <y, -x> and the second is labeled for the gradient, delta g = <x, -y>.
    Figure \(\PageIndex{11}\): (a) In this image, we see the three-level curves of \(g\) and vector field \(\vecs F\). Note that the \(\vecs F\) vectors on a given level curve are tangent to the level curve. (b) In this image, we see the three-level curves of \(g\) and vector field \(\vecs\nabla g\). The gradient vectors are perpendicular to the corresponding level curve. Therefore, \(\vecs F(a,b)\cdot\vecs\nabla g(a,b)=0\) for any point in the domain of \(g\).
    Exercise \(\PageIndex{6}\)

    Find a stream function for vector field \(\vecs F(x,y)=⟨x \sin y,\cos y⟩\).

    Hint

    Follow the outline provided in the previous example.

    Answer

    \(g(x,y)=−x\cos y\)

    Vector fields that are both conservative and source free are important vector fields. One important feature of conservative and source-free vector fields on a simply connected domain is that any potential function \(f\) of such a field satisfies Laplace’s equation \(f_{xx}+f_{yy}=0\). Laplace’s equation is foundational in the field of partial differential equations because it models such phenomena as gravitational and magnetic potentials in space, and the velocity potential of an ideal fluid. A function that satisfies Laplace’s equation is called a harmonic function. Therefore any potential function of a conservative and source-free vector field is harmonic.

    To see that any potential function of a conservative and source-free vector field on a simply connected domain is harmonic, let \(f\) be such a potential function of vector field \(\vecs F=⟨P,Q⟩\). Then, \(f_x=P\) and \(f_x=Q\) because \(\vecs \nabla f=\vecs F\). Therefore, \(f_{xx}=P_x\) and \(f_{yy}=Q_y\). Since \(\vecs F\) is source free, \(f_{xx}+f_{yy}=P_x+Q_y=0\), and we have that \(f\) is harmonic.

    Example \(\PageIndex{7}\): Satisfying Laplace’s Equation

    For vector field \(\vecs F(x,y)=⟨e^x\sin y,e^x\cos y⟩\), verify that the field is both conservative and source free, find a potential function for \(\vecs F\), and verify that the potential function is harmonic.

    Solution

    Let \(P(x,y)=e^x\sin y\) and \(Q(x,y)=e^x \cos y\). Notice that the domain of \(\vecs F\) is all of two-space, which is simply connected. Therefore, we can check the cross-partials of \(\vecs F\) to determine whether \(\vecs F\) is conservative. Note that \(P_y=e^x \cos y=Q_x\), so \(\vecs F\) is conservative. Since \(P_x=e^x \sin y\) and \(Q_y=e^x \sin y\), \(P_x+Q_y=0\) and the field is source free.

    To find a potential function for \(\vecs F\), let \(f\) be a potential function. Then, \(\vecs \nabla f=\vecs F\), so \(f_x(x,y)=e^x \sin y\). Integrating this equation with respect to x gives \(f(x,y)=e^x \sin y+h(y)\). Since \(f_y(x,y)=e^x \cos y\), differentiating \(f\) with respect to y gives \(e^x\cos y=e^x\cos y+h′(y)\). Therefore, we can take \(h(y)=0\), and \(f(x,y)=e^x\sin y\) is a potential function for \(f\).

    To verify that \(f\) is a harmonic function, note that \(f_{xx}(x,y)=\dfrac{\partial}{\partial x}(e^x\sin y)=e^x \sin y\) and

    \(f_{yy}(x,y)=\dfrac{\partial}{\partial x}(e^x\cos y)=−e^x\sin y\). Therefore, \(f_{xx}+f_{yy}=0\), and \(f\) satisfies Laplace’s equation.

    Exercise \(\PageIndex{7}\)

    Is the function \(f(x,y)=e^{x+5y}\) harmonic?

    Hint

    Determine whether the function satisfies Laplace’s equation.

    Answer

    No

    Green’s Theorem on General Regions

    Green’s theorem, as stated, applies only to regions that are simply connected—that is, Green’s theorem as stated so far cannot handle regions with holes. Here, we extend Green’s theorem so that it does work on regions with finitely many holes (Figure \(\PageIndex{12}\)).

    A nonsimply connected, oval-shaped region with three circular holes.
    Figure \(\PageIndex{12}\): Green’s theorem, as stated, does not apply to a nonsimply connected region with three holes like this one.

    Before discussing extensions of Green’s theorem, we need to go over some terminology regarding the boundary of a region. Let \(D\) be a region and let \(C\) be a component of the boundary of \(D\). We say that \(C\) is positively oriented if, as we walk along \(C\) in the direction of orientation, region \(D\) is always on our left. Therefore, the counterclockwise orientation of the boundary of a disk is a positive orientation, for example. Curve \(C\) is negatively oriented if, as we walk along \(C\) in the direction of orientation, region \(D\) is always on our right. The clockwise orientation of the boundary of a disk is a negative orientation, for example.

    Let \(D\) be a region with finitely many holes (so that \(D\) has finitely many boundary curves), and denote the boundary of \(D\) by \(\partial D\) (Figure \(\PageIndex{13}\)). To extend Green’s theorem so it can handle \(D\), we divide region \(D\) into two regions, \(D_1\) and \(D_2\) (with respective boundaries \(\partial D_1\) and \(\partial D_2\)), in such a way that \(D=D_1\cup D_2\) and neither \(D_1\) nor \(D_2\) has any holes (Figure \(\PageIndex{13}\)).

    Two regions. The first region D is oval-shaped with three circular holes in it. Its oriented boundary is counterclockwise. The second region is region D split horizontally down the middle into two simply connected regions with no holes. It still has a boundary oriented counterclockwise.
    Figure \(\PageIndex{13}\): (a) Region \(D\) with an oriented boundary has three holes. (b) Region \(D\) split into two simply connected regions has no holes.

    Assume the boundary of \(D\) is oriented as in the figure, with the inner holes given a negative orientation and the outer boundary given a positive orientation. The boundary of each simply connected region \(D_1\) and \(D_2\) is positively oriented. If \(\vecs F\) is a vector field defined on \(D\), then Green’s theorem says that

    \[\begin{align} \oint_{\partial D} \vecs F·d\vecs{r} &=\oint_{\partial D_1}\vecs F·d\vecs{r}+\oint_{\partial D_2}\vecs F·d\vecs{r} \\ &=\iint_{D_1}Q_x−P_y\,dA+\iint_{D_2}Q_x−P_y\,dA \\ &=\iint_D (Q_x−P_y)\,dA.\end{align} \nonumber \]

    Therefore, Green’s theorem still works on a region with holes.

    To see how this works in practice, consider annulus \(D\) in Figure \(\PageIndex{14}\) and suppose that \(F=⟨P,Q⟩\) is a vector field defined on this annulus. Region \(D\) has a hole, so it is not simply connected. Orient the outer circle of the annulus counterclockwise and the inner circle clockwise (Figure \(\PageIndex{14}\)) so that, when we divide the region into \(D_1\) and \(D_2\), we are able to keep the region on our left as we walk along a path that traverses the boundary. Let \(D_1\) be the upper half of the annulus and \(D_2\) be the lower half. Neither of these regions has holes, so we have divided \(D\) into two simply connected regions.

    We label each piece of these new boundaries as \(P_i\) for some \(i\), as in Figure \(\PageIndex{14}\). If we begin at \(P\) and travel along the oriented boundary, the first segment is \(P_1\), then \(P_2\), \(P_3\), and \(P_4\). Now we have traversed \(D_1\) and returned to \(P\). Next, we start at \(P\) again and traverse \(D_2\). Since the first piece of the boundary is the same as \(P_4\) in \(D_1\), but oriented in the opposite direction, the first piece of \(D_2\) is \(−P_4\). Next, we have \(P_5\), then \(−P_2\), and finally \(P_6\).

    A diagram of an annulus – a circular region with a hole in in like a donut. Its boundary is oriented counterclockwise. One point P on the outer boundary is labeled. It is the right endpoint of the horizontal diameter. The annulus is split horizontally down the middle into two separate regions that are each simply connected. Point P is labeled on both of these regions, D1 and D2. Each region has boundaries oriented counterclockwise. The upper curve of D1 is labeled P1, the left flat side is P2, the lower curve is P3, and the right flat side is P4. The lower curve of D2 is P6, the left flat side is –P2, the upper curve is P5, and the right flat side is –P4.
    Figure \(\PageIndex{14}\): Breaking the annulus into two separate regions gives us two simply connected regions. The line integrals over the common boundaries cancel out.

    Figure \(\PageIndex{14}\) shows a path that traverses the boundary of \(D\). Notice that this path traverses the boundary of region \(D_1\), returns to the starting point, and then traverses the boundary of region \(D_2\). Furthermore, as we walk along the path, the region is always on our left. Notice that this traversal of the \(P_i\) paths covers the entire boundary of region \(D\). If we had only traversed one portion of the boundary of \(D\), then we cannot apply Green’s theorem to \(D\).

    The boundary of the upper half of the annulus, therefore, is \(P_1\cup P_2\cup P_3\cup P_4\) and the boundary of the lower half of the annulus is \(−P_4\cup P_5\cup −P_2\cup P_6\). Then, Green’s theorem implies

    \[\begin{align} \oint_{\partial D}\vecs F·d\vecs{r} &=\int_{P_1}\vecs F·d\vecs{r}+\int_{P_2}\vecs F·d\vecs{r}+\int_{P_3}\vecs F·d\vecs{r}+\int_{P_4}\vecs F·d\vecs{r}+\int_{−P_4}\vecs F·d\vecs{r}+\int_{P_5}\vecs F·d\vecs{r}+\int_{−P_2}\vecs F·d\vecs{r}+\int_{P_6}\vecs F·d\vecs{r} \\ &=\int_{P_1}\vecs F·d\vecs{r}+\int_{P_2}\vecs F·d\vecs{r}+\int_{P_3}\vecs F·d\vecs{r}+\int_{P_4}\vecs F·d\vecs{r}+\int_{P_4}\vecs F·d\vecs{r}+\int_{P_5}\vecs F·d\vecs{r}+\int_{−P_2}\vecs F·d\vecs{r}+\int_{P_6}\vecs F·d\vecs{r} \\ &=\int_{P_1}\vecs F·d\vecs{r}+\int_{P_3}\vecs F·d\vecs{r}+\int_{P_5}\vecs F·d\vecs{r}+\int_{P_6}\vecs F·d\vecs{r} \\ &=\oint_{\partial D_1}\vecs F·d\vecs{r}+\oint_{\partial D_2}\vecs F·d\vecs{r}\\ &=\iint_{D_1}(Q_x−P_y)\,dA+\iint_{D_2}(Q_x−P_y)\,dA \\ &=\iint_D(Q_x−P_y)\,dA. \end{align} \nonumber \]

    Therefore, we arrive at the equation found in Green’s theorem—namely,

    \[\oint_{\partial D}\vecs F·d\vecs{r}=\iint_D (Q_x−P_y)\,dA. \nonumber \]

    The same logic implies that the flux form of Green’s theorem can also be extended to a region with finitely many holes:

    \[\oint_C \vecs F·\vecs N\,ds=\iint_D (P_x+Q_y)\,dA. \nonumber \]

    Example \(\PageIndex{8A}\): ​​​​​​Using Green’s Theorem on a Region with Holes

    Calculate the integral

    \[\oint_{\partial D}(\sin x−\dfrac{y^3}{3})dx+(\dfrac{y^3}{3}+\sin y)dy, \nonumber \]

    where \(D\) is the annulus given by the polar inequalities \(1≤r≤2\), \(0≤\theta≤2\pi\).

    Solution

    Although \(D\) is not simply connected, we can use the extended form of Green’s theorem to calculate the integral. Since the integration occurs over an annulus, we convert to polar coordinates:

    \[\begin{align*} \oint_{\partial D}(\sin x−\dfrac{y^3}{3})\,dx+(\dfrac{x^3}{3}+\sin y)\,dy &=\iint_D (Q_x−P_y)\,dA \\ &=\iint_D (x^2+y^2)\,dA\\ &=\int_0^{2\pi}\int_1^2 r^3\,drd\theta=\int_0^{2\pi} \dfrac{15}{4}\,d\theta \\ &=\dfrac{15\pi}{2}. \end{align*}\]

    Example \(\PageIndex{8B}\): Using the Extended Form of Green’s Theorem

    Let \(\vecs F=⟨P,Q⟩=\left\langle\dfrac{y}{x^2+y^2},-\dfrac{x}{x^2+y^2}\right\rangle\) and let \(C\) be any simple closed curve in a plane oriented counterclockwise. What are the possible values of \(\displaystyle \oint_C \vecs F·d\vecs{r}\)?

    Solution

    We use the extended form of Green’s theorem to show that \(\displaystyle \oint_C \vecs F·d\vecs{r}\) is either \(0\) or \(−2\pi\)—that is, no matter how crazy curve \(C\) is, the line integral of \(\vecs F\) along \(C\) can have only one of two possible values. We consider two cases: the case when \(C\) encompasses the origin and the case when \(C\) does not encompass the origin.

    Case 1: C Does Not Encompass the Origin

    In this case, the region enclosed by \(C\) is simply connected because the only hole in the domain of \(\vecs F\) is at the origin. We showed in our discussion of cross-partials that \(\vecs F\) satisfies the cross-partial condition. If we restrict the domain of \(\vecs F\) just to \(C\) and the region it encloses, then \(\vecs F\) with this restricted domain is now defined on a simply connected domain. Since \(\vecs F\) satisfies the cross-partial property on its restricted domain, the field \(\vecs F\) is conservative on this simply connected region and hence the circulation \(\displaystyle \oint_C \vecs F·d\vecs{r}\) is zero.

    Case 2: C Does Encompass the Origin

    In this case, the region enclosed by \(C\) is not simply connected because this region contains a hole at the origin. Let \(C_1\) be a circle of radius a centered at the origin so that \(C_1\) is entirely inside the region enclosed by \(C\) (Figure \(\PageIndex{15}\)). Give \(C_1\) a clockwise orientation.

    A diagram in two dimensions. A circle C1 oriented clockwise is centered at the origin completely inside a generic curve C that is in all four quadrants. Curve C is oriented counterclockwise.
    Figure \(\PageIndex{15}\): Choose circle \(C_1\) centered at the origin that is contained entirely inside \(C\).

    Let \(D\) be the region between \(C_1\) and \(C\), and \(C\) is orientated counterclockwise. By the extended version of Green’s theorem,

    \[\begin{align*} \int_C \vecs F·d\vecs{r}+\int_{C_1}\vecs F·d\vecs{r} &=\iint_D Qx_−P_y \,dA \\[4pt] &=\iint_D−\dfrac{y^2−x^2}{{(x^2+y^2)}^2}+\dfrac{y^2−x^2}{{(x^2+y^2)}^2}dA \\[4pt] &=0, \end{align*}\]

    and therefore

    \[\int_C \vecs F·d\vecs{r}=−\int_{C_1} \vecs F·d\vecs{r}. \nonumber \]

    Since \(C_1\) is a specific curve, we can evaluate \(\displaystyle \int_{C_1}\vecs F·d\vecs{r}\). Let

    \[ x=a\cos t, \;\; y=a\sin t, \;\; 0≤t≤2\pi \nonumber \]

    be a parameterization of \(C_1\). Then,

    \[\begin{align*} \int_{C_1}\vecs F·d\vecs{r} &=\int_0^{2\pi} \vecs F(\vecs r(t))·\vecs r′(t)\,dt \\[4pt] &=\int_0^{2\pi} ⟨−\dfrac{\sin(t)}{a},−\dfrac{\cos(t)}{a}⟩·⟨−a\sin(t),−a\cos(t)⟩\,dt \\[4pt] &=\int_0^{2\pi}{\sin}^2(t)+{\cos}^2(t)\,dt \\[4pt] &=\int_0^{2\pi}\,dt=2\pi. \end{align*}\]

    Therefore, \(\displaystyle \int_C \vecs F·d\vecs r=−2\pi\).

    Exercise \(\PageIndex{8}\)

    Calculate integral \(\displaystyle \oint_{\partial D}\vecs F·d\vecs{r}\), where \(D\) is the annulus given by the polar inequalities \(2≤r≤5\), \(0≤\theta≤2\pi\), and \(\vecs F(x,y)=⟨x^3,5x+e^y\sin y⟩\).

    Hint

    Use the extended version of Green’s theorem.

    Answer

    \(105 \pi\)

    MEASURING AREA FROM A BOUNDARY: THE PLANIMETER

    Imagine you are a doctor who has just received a magnetic resonance image of your patient’s brain. The brain has a tumor (Figure \(\PageIndex{16}\)). How large is the tumor? To be precise, what is the area of the red region? The red cross-section of the tumor has an irregular shape, and therefore it is unlikely that you would be able to find a set of equations or inequalities for the region and then be able to calculate its area by conventional means. You could approximate the area by chopping the region into tiny squares (a Riemann sum approach), but this method always gives an answer with some error.

    An MRI image of a patient’s brain with a tumor highlighted in red.
    Figure \(\PageIndex{16}\): This magnetic resonance image of a patient’s brain shows a tumor, which is highlighted in red. (credit: modification of work by Christaras A, Wikimedia Commons)

    Instead of trying to measure the area of the region directly, we can use a device called a rolling planimeter to calculate the area of the region exactly, simply by measuring its boundary. In this project you investigate how a planimeter works, and you use Green’s theorem to show the device calculates area correctly.

    A rolling planimeter is a device that measures the area of a planar region by tracing out the boundary of that region (Figure \(\PageIndex{17}\)). To measure the area of a region, we simply run the tracer of the planimeter around the boundary of the region. The planimeter measures the number of turns through which the wheel rotates as we trace the boundary; the area of the shape is proportional to this number of wheel turns. We can derive the precise proportionality equation using Green’s theorem. As the tracer moves around the boundary of the region, the tracer arm rotates and the roller moves back and forth (but does not rotate).

    Two images. The first shows a rolling planimeter. A horizontal bar has a roller attached to it perpendicularly with a pivot. It does not rotate itself; it only moves back and forth. To the right of the roller is the tracer arm with a wheel and a tracer at the very end. The second shows an interior view of a rolling planimeter. The wheel cannot turn if the planimeter is moving back and forth with the tracer arm perpendicular to the roller.
    Figure \(\PageIndex{17}\): (a) A rolling planimeter. The pivot allows the tracer arm to rotate. The roller itself does not rotate; it only moves back and forth. (b) An interior view of a rolling planimeter. Notice that the wheel cannot turn if the planimeter is moving back and forth with the tracer arm perpendicular to the roller.

    Let \(C\) denote the boundary of region \(D\), the area to be calculated. As the tracer traverses curve \(C\), assume the roller moves along the y-axis (since the roller does not rotate, one can assume it moves along a straight line). Use the coordinates \((x,y)\) to represent points on boundary \(C\), and coordinates \((0,Y)\) to represent the position of the pivot. As the planimeter traces \(C\), the pivot moves along the y-axis while the tracer arm rotates on the pivot.

    Watch a short animation of a planimeter in action.

    Begin the analysis by considering the motion of the tracer as it moves from point \((x,y)\) counterclockwise to point \((x+dx,y+dy)\) that is close to \((x,y)\) (Figure \(\PageIndex{18}\)). The pivot also moves, from point \((0,Y)\) to nearby point \((0,Y+dY)\). How much does the wheel turn as a result of this motion? To answer this question, break the motion into two parts. First, roll the pivot along the y-axis from \((0,Y)\) to \((0,Y+dY)\) without rotating the tracer arm. The tracer arm then ends up at point \((x,y+dY)\) while maintaining a constant angle \(\phi\) with the x-axis. Second, rotate the tracer arm by an angle \(d\theta\) without moving the roller. Now the tracer is at point \((x+dx,y+dy)\). Let ll be the distance from the pivot to the wheel and let \(L\) be the distance from the pivot to the tracer (the length of the tracer arm).

    A diagram in quadrants 1 and 2 showing the motion of the planimeter. Two points are labeled on the y axis: (0, Y) and (0, Y + dY), where Y is less than Y + dY. The first point is the pivot. Three points are labeled further up and to the right in quadrant 1: (x, y), (x, y + dy), and (x + dx, y + dy). Note that the uppercase Y and the lowercase y are not the same; y is much larger. A line segment is drawn between (0,Y) and (x,y). About midway down this line is a mark labeled for the wheel, and the (x,y) endpoint is labeled for the tracer. Let l be the distance from the pivot to the wheel, and let L be the distance from the pivot to the tracer. Line segments are also drawn from (0, Y + dY) to each of the other points in quadrant 1. The angle between the line segment with (0,Y) as an endpoint and the y axis is labeled phi. The angle between the line segments with (0, Y+dY) as an endpoint is “d theta.” A curve is drawn going through the wheel, the tracer, and the three points in quadrant 1, up and across the y axis, down and back across the y axis at a smaller y value lose to the height of the tracer, and down across the line segments and back to the wheel.
    Figure \(\PageIndex{18}\): Mathematical analysis of the motion of the planimeter.
    1. Explain why the total distance through which the wheel rolls the small motion just described is \(\sin \phi dY+ld\theta=\dfrac{x}{L}dY+ld\theta\).
    2. Show that \(\oint_C d\theta=0\).
    3. Use step 2 to show that the total rolling distance of the wheel as the tracer traverses curve \(C\) is
      Total wheel roll \(=\dfrac{1}{L}\oint_C xdY\).
      Now that you have an equation for the total rolling distance of the wheel, connect this equation to Green’s theorem to calculate area \(D\) enclosed by \(C\).
    4. Show that \(x^2+(y−Y)^2=L^2\).
    5. Assume the orientation of the planimeter is as shown in Figure \(\PageIndex{18}\). Explain why \(Y≤y\), and use this inequality to show there is a unique value of \(Y\) for each point \((x,y)\): \(Y=y=\sqrt{L^2−x^2}\).
    6. Use step 5 to show that \(dY=dy+\dfrac{x}{L^2−x^2}dx.\)
    7. Use Green’s theorem to show that \(\displaystyle \oint_C \dfrac{x}{L^2−x^2}dx=0\).
    8. Use step 7 to show that the total wheel roll is

      \[\text{Total wheel roll}\quad =\quad 1L\oint_C x\,dy. \nonumber \]

      It took a bit of work, but this equation says that the variable of integration Y in step 3 can be replaced with y.

    9. Use Green’s theorem to show that the area of \(D\) is \(\displaystyle \oint_C x\,dy\). The logic is similar to the logic used to show that the area of \(\displaystyle D=12\oint_C −y\,dx+x\,dy\).
    10. Conclude that the area of \(D\) equals the length of the tracer arm multiplied by the total rolling distance of the wheel.

    You now know how a planimeter works and you have used Green’s theorem to justify that it works. To calculate the area of a planar region \(D\), use a planimeter to trace the boundary of the region. The area of the region is the length of the tracer arm multiplied by the distance the wheel rolled.

    Key Concepts

    • Green’s theorem relates the integral over a connected region to an integral over the boundary of the region. Green’s theorem is a version of the Fundamental Theorem of Calculus in one higher dimension.
    • Green’s Theorem comes in two forms: a circulation form and a flux form. In the circulation form, the integrand is \(\vecs F·\vecs T\). In the flux form, the integrand is \(\vecs F·\vecs N\).
    • Green’s theorem can be used to transform a difficult line integral into an easier double integral, or to transform a difficult double integral into an easier line integral.
    • A vector field is source free if it has a stream function. The flux of a source-free vector field across a closed curve is zero, just as the circulation of a conservative vector field across a closed curve is zero.

    Key Equations

    • Green’s theorem, circulation form
      \(\displaystyle ∮_C P\,dx+Q\,dy=∬_D Q_x−P_y\,dA\), where \(C\) is the boundary of \(D\)
    • Green’s theorem, flux form
      \(\displaystyle ∮_C\vecs F·\vecs N\,ds=∬_D P_x+Q_y\,dA\), where \(C\) is the boundary of \(D\)
    • Green’s theorem, extended version
      \(\displaystyle ∮_{\partial D}\vecs F·d\vecs{r}=∬_D Q_x−P_y\,dA\)

    Glossary

    Green’s theorem
    relates the integral over a connected region to an integral over the boundary of the region
    stream function
    if \(\vecs F=⟨P,Q⟩\) is a source-free vector field, then stream function \(g\) is a function such that \(P=g_y\) and \(Q=−g_x\)

    Contributions 

    Paul Seeburger created a corrected image and the worked out solution for Exercise \(\PageIndex{3}\).


    5.5: Green’s Theorem is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by LibreTexts.