Skip to content

Vector Calculus

3.1 Vector Fields

A vector field on Rn\mathbb{R}^n is a function F:DRnRn\mathbf{F} : D \subseteq \mathbb{R}^n \to \mathbb{R}^n.

A vector field F=(P,Q,R)\mathbf{F} = (P, Q, R) on R3\mathbb{R}^3 is conservative if there exists a scalar Potential ϕ\phi such that F=ϕ\mathbf{F} = \nabla \phi.

Theorem 3.1. F\mathbf{F} is conservative (on a connected domain) if and only if ×F=0\nabla \times \mathbf{F} = \mathbf{0}.

Proof. (\Rightarrow) If F=ϕ\mathbf{F} = \nabla \phi with ϕC2\phi \in C^2Then by Clairaut”s theorem fxy=fyxf_{xy} = f_{yx}Etc., which directly gives ×(ϕ)=0\nabla \times (\nabla \phi) = \mathbf{0}.

(\Leftarrow) If ×F=0\nabla \times \mathbf{F} = \mathbf{0} on a connected domain DDThen for any Closed curve CC in DDStokes’ theorem gives CFdr=S(×F)dS=0\oint_C \mathbf{F} \cdot d\mathbf{r} = \iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S} = 0. This means line integrals are path-independent, so we can define ϕ(x)=x0xFdr\phi(\mathbf{x}) = \int_{\mathbf{x}_0}^{\mathbf{x}} \mathbf{F} \cdot d\mathbf{r} (independent of path), And one verifies that ϕ=F\nabla \phi = \mathbf{F}. \blacksquare

3.2 Line Integrals

Definition. The line integral of a vector field F\mathbf{F} along a curve CC parameterised by r(t)\mathbf{r}(t) for atba \leq t \leq b is

CFdr=abF(r(t))r(t)dt\int_C \mathbf{F} \cdot d\mathbf{r} = \int_a^b \mathbf{F}(\mathbf{r}(t)) \cdot \mathbf{r}'(t)\, dt

Theorem 3.2 (Fundamental Theorem for Line Integrals). If F=ϕ\mathbf{F} = \nabla \phi and CC is a Piecewise smooth curve from AA to BBThen

CFdr=ϕ(B)ϕ(A)\int_C \mathbf{F} \cdot d\mathbf{r} = \phi(B) - \phi(A)

Proof. Parameterise CC by r(t)\mathbf{r}(t) for t[a,b]t \in [a,b] with r(a)=A\mathbf{r}(a) = A r(b)=B\mathbf{r}(b) = B.

CFdr=abϕ(r(t))r(t)dt=abddt[ϕ(r(t))]dt=ϕ(r(b))ϕ(r(a))=ϕ(B)ϕ(A)\int_C \mathbf{F} \cdot d\mathbf{r} = \int_a^b \nabla \phi(\mathbf{r}(t)) \cdot \mathbf{r}'(t)\, dt = \int_a^b \frac{d}{dt}\left[\phi(\mathbf{r}(t))\right]\, dt = \phi(\mathbf{r}(b)) - \phi(\mathbf{r}(a)) = \phi(B) - \phi(A)

By the chain rule. \blacksquare

Corollary 3.3. The line integral of a conservative field around any closed curve is zero.

Problem. Evaluate CFdr\int_C \mathbf{F} \cdot d\mathbf{r} where F=(y,x+ey,z+1)\mathbf{F} = (y,\, x + e^y,\, z + 1) and CC is the Curve r(t)=(t,t2,t3)\mathbf{r}(t) = (t,\, t^2,\, t^3) for 0t10 \leq t \leq 1.

Solution

First check if F\mathbf{F} is conservative. Compute the curl:

(×F)x=(z+1)y(x+ey)z=00=0(\nabla \times \mathbf{F})_x = \frac{\partial (z+1)}{\partial y} - \frac{\partial (x + e^y)}{\partial z} = 0 - 0 = 0

(×F)y=yz(z+1)x=00=0(\nabla \times \mathbf{F})_y = \frac{\partial y}{\partial z} - \frac{\partial (z+1)}{\partial x} = 0 - 0 = 0

(×F)z=(x+ey)xyy=11=0(\nabla \times \mathbf{F})_z = \frac{\partial (x + e^y)}{\partial x} - \frac{\partial y}{\partial y} = 1 - 1 = 0

Since ×F=0\nabla \times \mathbf{F} = \mathbf{0}, F\mathbf{F} is conservative. Find ϕ\phi:

ϕx=y    ϕ=xy+g(y,z)\frac{\partial \phi}{\partial x} = y \implies \phi = xy + g(y,z)

ϕy=x+gy=x+ey    gy=ey    g=ey+h(z)\frac{\partial \phi}{\partial y} = x + g_y = x + e^y \implies g_y = e^y \implies g = e^y + h(z)

ϕz=h(z)=z+1    h(z)=z22+z+C\frac{\partial \phi}{\partial z} = h'(z) = z + 1 \implies h(z) = \frac{z^2}{2} + z + C

ϕ(x,y,z)=xy+ey+z22+z\phi(x,y,z) = xy + e^y + \frac{z^2}{2} + z

Now apply the fundamental theorem:

CFdr=ϕ(1,1,1)ϕ(0,0,0)=(1+e+12+1)(1+1)=e+12\int_C \mathbf{F} \cdot d\mathbf{r} = \phi(1, 1, 1) - \phi(0, 0, 0) = \left(1 + e + \frac{1}{2} + 1\right) - (1 + 1) = e + \frac{1}{2}

\blacksquare

3.3 Green’s Theorem

Theorem 3.4 (Green’s Theorem). Let CC be a positively oriented, piecewise smooth, simple closed Curve bounding a region DD. If PP and QQ have continuous partial derivatives on an open set Containing DDThen

CPdx+Qdy=D(QxPy)dA\oint_C P\, dx + Q\, dy = \iint_D \left(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right) dA

Proof (for a Type I region). Assume DD is a Type I region: D=(x,y):axb,g1(x)yg2(x)D = \\{(x,y) : a \leq x \leq b,\, g_1(x) \leq y \leq g_2(x)\\}. The boundary CC consists of Four pieces: bottom C1C_1Right C2C_2Top C3C_3And left C4C_4.

We first prove CPdx=DPydA\oint_C P\, dx = -\iint_D \frac{\partial P}{\partial y}\, dA.

On C1C_1: y=g1(x)y = g_1(x), xx goes from aa to bbSo C1Pdx=abP(x,g1(x))dx\int_{C_1} P\, dx = \int_a^b P(x, g_1(x))\, dx.

On C3C_3: y=g2(x)y = g_2(x), xx goes from bb to aaSo C3Pdx=baP(x,g2(x))dx=abP(x,g2(x))dx\int_{C_3} P\, dx = \int_b^a P(x, g_2(x))\, dx = -\int_a^b P(x, g_2(x))\, dx.

On C2C_2 and C4C_4: xx is constant, so dx=0dx = 0Hence C2Pdx=C4Pdx=0\int_{C_2} P\, dx = \int_{C_4} P\, dx = 0.

Therefore:

CPdx=abP(x,g1(x))dxabP(x,g2(x))dx\oint_C P\, dx = \int_a^b P(x, g_1(x))\, dx - \int_a^b P(x, g_2(x))\, dx

Meanwhile:

DPydA=abg1(x)g2(x)Pydydx=ab[P(x,g2(x))P(x,g1(x))]dx-\iint_D \frac{\partial P}{\partial y}\, dA = -\int_a^b \int_{g_1(x)}^{g_2(x)} \frac{\partial P}{\partial y}\, dy\, dx = -\int_a^b \left[P(x, g_2(x)) - P(x, g_1(x))\right]\, dx

=abP(x,g1(x))dxabP(x,g2(x))dx=CPdx= \int_a^b P(x, g_1(x))\, dx - \int_a^b P(x, g_2(x))\, dx = \oint_C P\, dx

An identical argument (using Type II regions) proves CQdy=DQxdA\oint_C Q\, dy = \iint_D \frac{\partial Q}{\partial x}\, dA. Adding the two equalities gives the result. For general regions, decompose DD into finitely many Type I and Type II regions and note that the line integrals along shared boundaries cancel. \blacksquare

Worked Example. Evaluate C(x2y)dx+(y2+x)dy\oint_C (x^2 - y)\, dx + (y^2 + x)\, dy where CC is the unit circle Traversed counterclockwise.

Solution. By Green’s theorem with P=x2yP = x^2 - y and Q=y2+xQ = y^2 + x:

Qx=1,Py=1\frac{\partial Q}{\partial x} = 1, \quad \frac{\partial P}{\partial y} = -1

CPdx+Qdy=D(1(1))dA=2DdA=2π12=2π\oint_C P\, dx + Q\, dy = \iint_D (1 - (-1))\, dA = 2 \iint_D dA = 2 \cdot \pi \cdot 1^2 = 2\pi

\blacksquare

3.4 Curl and Divergence

Definition. Let F=(P,Q,R)\mathbf{F} = (P, Q, R) be a C1C^1 vector field on R3\mathbb{R}^3.

The curl of F\mathbf{F} is

×F=(RyQz,PzRx,QxPy)\nabla \times \mathbf{F} = \left(\frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z},\, \frac{\partial P}{\partial z} - \frac{\partial R}{\partial x},\, \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right)

The divergence of F\mathbf{F} is

F=Px+Qy+Rz\nabla \cdot \mathbf{F} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z}

Physical interpretation. If F\mathbf{F} represents the velocity field of a fluid:

  • Curl ×F\nabla \times \mathbf{F} measures the local rotational tendency (vorticity) of the fluid. At a point p\mathbf{p}The component (×F)n(\nabla \times \mathbf{F}) \cdot \mathbf{n} gives twice the angular velocity of a small paddle wheel placed at p\mathbf{p} with axis along n\mathbf{n}.

  • Divergence F\nabla \cdot \mathbf{F} measures the net rate of outward flux per unit volume at a point. If F>0\nabla \cdot \mathbf{F} \gt 0 at p\mathbf{p}There is a net source at p\mathbf{p}; if F<0\nabla \cdot \mathbf{F} \lt 0There is a net sink.

Proposition 3.5. For any C2C^2 vector field F\mathbf{F}:

(×F)=0(div of curl is zero)\nabla \cdot (\nabla \times \mathbf{F}) = 0 \quad \mathrm{(div\ of\ curl\ is\ zero)}

×(ϕ)=0(curl of gradient is zero)\nabla \times (\nabla \phi) = \mathbf{0} \quad \mathrm{(curl\ of\ gradient\ is\ zero)}

Proof. Both follow from Clairaut’s theorem on equality of mixed partials. For the first:

(×F)=x(RyQz)+y(PzRx)+z(QxPy)\nabla \cdot (\nabla \times \mathbf{F}) = \frac{\partial}{\partial x}\left(\frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z}\right) + \frac{\partial}{\partial y}\left(\frac{\partial P}{\partial z} - \frac{\partial R}{\partial x}\right) + \frac{\partial}{\partial z}\left(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right)

Each pair cancels by Clairaut: 2Rxy=2Ryx\frac{\partial^2 R}{\partial x\,\partial y} = \frac{\partial^2 R}{\partial y\,\partial x}Etc. \blacksquare

3.5 Stokes’ Theorem

Theorem 3.6 (Stokes’ Theorem). Let SS be an oriented surface with piecewise smooth boundary curve CC (positively oriented). If F\mathbf{F} has continuous partial derivatives on an open set containing SSThen

CFdr=S(×F)dS\oint_C \mathbf{F} \cdot d\mathbf{r} = \iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S}

Where dS=ndSd\mathbf{S} = \mathbf{n}\, dS is the vector surface element with unit normal n\mathbf{n}.

Proof (sketch). Parametrise SS by r(u,v)\mathbf{r}(u,v) over a region DD in the uvuv-plane. The boundary CC of SS corresponds to the boundary D\partial D of DD. The left-hand side becomes:

CFdr=DF(r(u,v))(rudu+rvdv)\oint_C \mathbf{F} \cdot d\mathbf{r} = \oint_{\partial D} \mathbf{F}(\mathbf{r}(u,v)) \cdot \left(\frac{\partial \mathbf{r}}{\partial u}\, du + \frac{\partial \mathbf{r}}{\partial v}\, dv\right)

Define P~(u,v)=F(r(u,v))ru\tilde{P}(u,v) = \mathbf{F}(\mathbf{r}(u,v)) \cdot \mathbf{r}_u and Q~(u,v)=F(r(u,v))rv\tilde{Q}(u,v) = \mathbf{F}(\mathbf{r}(u,v)) \cdot \mathbf{r}_v. Applying Green’s theorem in the uvuv-plane:

DP~du+Q~dv=D(Q~uP~v)dudv\oint_{\partial D} \tilde{P}\, du + \tilde{Q}\, dv = \iint_D \left(\frac{\partial \tilde{Q}}{\partial u} - \frac{\partial \tilde{P}}{\partial v}\right) du\, dv

Expanding the partial derivatives and using the identity ru×rv=nru×rv\mathbf{r}_u \times \mathbf{r}_v = \mathbf{n}\, \lVert \mathbf{r}_u \times \mathbf{r}_v \rVertOne Verifies that the integrand equals (×F)nru×rv(\nabla \times \mathbf{F}) \cdot \mathbf{n}\, \lVert \mathbf{r}_u \times \mathbf{r}_v \rVert Which gives S(×F)dS\iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S}. \blacksquare

Remark. Green’s theorem is the special case of Stokes’ theorem where SS is a planar region in R2\mathbb{R}^2.

Problem. Use Stokes’ theorem to evaluate CFdr\oint_C \mathbf{F} \cdot d\mathbf{r} where F=(y2,xz,x2)\mathbf{F} = (y^2,\, xz,\, x^2) and CC is the triangle with vertices (1,0,0)(1,0,0), (0,1,0)(0,1,0) (0,0,1)(0,0,1) traversed counterclockwise when viewed from above.

Solution

The triangle lies in the plane x+y+z=1x + y + z = 1. Compute ×F\nabla \times \mathbf{F}:

(×F)x=(x2)y(xz)z=0x=x(\nabla \times \mathbf{F})_x = \frac{\partial (x^2)}{\partial y} - \frac{\partial (xz)}{\partial z} = 0 - x = -x

(×F)y=(y2)z(x2)x=02x=2x(\nabla \times \mathbf{F})_y = \frac{\partial (y^2)}{\partial z} - \frac{\partial (x^2)}{\partial x} = 0 - 2x = -2x

(×F)z=(xz)x(y2)y=z2y(\nabla \times \mathbf{F})_z = \frac{\partial (xz)}{\partial x} - \frac{\partial (y^2)}{\partial y} = z - 2y

So ×F=(x,2x,z2y)\nabla \times \mathbf{F} = (-x,\, -2x,\, z - 2y).

Parametrise the triangle in the xyxy-plane: 0x10 \leq x \leq 1, 0y1x0 \leq y \leq 1 - x. On the plane z=1xyz = 1 - x - yThe surface element dS=3dxdydS = \sqrt{3}\, dx\, dy.

S(×F)dS=13S(x2x+z2y)dS\iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S} = \frac{1}{\sqrt{3}} \iint_S (-x - 2x + z - 2y)\, dS

On the plane: 3x2y+z=3x2y+1xy=4x3y+1-3x - 2y + z = -3x - 2y + 1 - x - y = -4x - 3y + 1.

=130101x(4x3y+1)3dydx=0101x(4x3y+1)dydx= \frac{1}{\sqrt{3}} \int_0^1 \int_0^{1-x} (-4x - 3y + 1)\, \sqrt{3}\, dy\, dx = \int_0^1 \int_0^{1-x} (-4x - 3y + 1)\, dy\, dx

=01[(4x+1)y3y22]01xdx=01(4x+1)(1x)3(1x)22dx= \int_0^1 \left[(-4x + 1)y - \frac{3y^2}{2}\right]_0^{1-x}\, dx = \int_0^1 (-4x + 1)(1 - x) - \frac{3(1-x)^2}{2}\, dx

=01[4x25x+132+3x3x22]dx=01[5x222x12]dx= \int_0^1 \left[4x^2 - 5x + 1 - \frac{3}{2} + 3x - \frac{3x^2}{2}\right]\, dx = \int_0^1 \left[\frac{5x^2}{2} - 2x - \frac{1}{2}\right]\, dx

=[5x36x2x2]01=56112=23= \left[\frac{5x^3}{6} - x^2 - \frac{x}{2}\right]_0^1 = \frac{5}{6} - 1 - \frac{1}{2} = -\frac{2}{3}

\blacksquare

3.6 Divergence Theorem

Theorem 3.7 (Divergence Theorem / Gauss’s Theorem). Let EE be a solid region bounded by a closed Surface SS with outward normal n\mathbf{n}. If F\mathbf{F} has continuous partial derivatives on an Open set containing EEThen

SFdS=EFdV\iint_S \mathbf{F} \cdot d\mathbf{S} = \iiint_E \nabla \cdot \mathbf{F}\, dV

Where F=Px+Qy+Rz\nabla \cdot \mathbf{F} = \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} Is the divergence of F\mathbf{F}.

Proof (sketch for a Type I region). Assume EE is a Type I region: E=(x,y,z):(x,y)D,g1(x,y)zg2(x,y)E = \\{(x,y,z) : (x,y) \in D,\, g_1(x,y) \leq z \leq g_2(x,y)\\}. The boundary consists of Bottom surface S1S_1 (normal pointing downward), top surface S2S_2 (normal pointing upward), And the lateral surface S3S_3 (where the normal is horizontal).

We prove the result for the RR-component, i.e., SRkdS=ERzdV\iint_S R\, \mathbf{k} \cdot d\mathbf{S} = \iiint_E \frac{\partial R}{\partial z}\, dV.

The right-hand side:

ERzdV=Dg1(x,y)g2(x,y)RzdzdA=D[R(x,y,g2)R(x,y,g1)]dA\iiint_E \frac{\partial R}{\partial z}\, dV = \iint_D \int_{g_1(x,y)}^{g_2(x,y)} \frac{\partial R}{\partial z}\, dz\, dA = \iint_D \left[R(x,y,g_2) - R(x,y,g_1)\right]\, dA

On S2S_2 (top): dS=(g2x,g2y,1)dAd\mathbf{S} = (-g_{2x}, -g_{2y}, 1)\, dA (upward), so RkdS=R(x,y,g2)dAR\, \mathbf{k} \cdot d\mathbf{S} = R(x,y,g_2)\, dA.

On S1S_1 (bottom): dS=(g1x,g1y,1)dAd\mathbf{S} = (g_{1x}, g_{1y}, -1)\, dA (downward), so RkdS=R(x,y,g1)dAR\, \mathbf{k} \cdot d\mathbf{S} = -R(x,y,g_1)\, dA.

On S3S_3: kn=0\mathbf{k} \cdot \mathbf{n} = 0 (the normal is horizontal), so RkdS=0R\, \mathbf{k} \cdot d\mathbf{S} = 0.

Therefore SRkdS=D[R(x,y,g2)R(x,y,g1)]dA\iint_S R\, \mathbf{k} \cdot d\mathbf{S} = \iint_D [R(x,y,g_2) - R(x,y,g_1)]\, dAMatching the Volume integral. The PP and QQ components follow by an identical argument for Type II and Type III Regions. For general regions, decompose into finitely many regions of each type. \blacksquare

Worked Example. Compute the flux of F=(x3,y3,z3)\mathbf{F} = (x^3, y^3, z^3) through the unit sphere SS.

Solution. By the divergence theorem:

F=3x2+3y2+3z2=3(x2+y2+z2)=3ρ2\nabla \cdot \mathbf{F} = 3x^2 + 3y^2 + 3z^2 = 3(x^2 + y^2 + z^2) = 3\rho^2

Using spherical coordinates:

E3ρ2ρ2sinϕdρdϕdθ=302π0π01ρ4sinϕdρdϕdθ\iiint_E 3\rho^2 \cdot \rho^2 \sin\phi\, d\rho\, d\phi\, d\theta = 3 \int_0^{2\pi} \int_0^{\pi} \int_0^1 \rho^4 \sin\phi\, d\rho\, d\phi\, d\theta

=32π215=12π5= 3 \cdot 2\pi \cdot 2 \cdot \frac{1}{5} = \frac{12\pi}{5}

\blacksquare

Problem. Compute the flux of F=(x2,y2,z2)\mathbf{F} = (x^2,\, y^2,\, z^2) outward through the surface of the Cylinder x2+y21x^2 + y^2 \leq 1, 0z20 \leq z \leq 2.

Solution

By the divergence theorem:

F=2x+2y+2z\nabla \cdot \mathbf{F} = 2x + 2y + 2z

Use cylindrical coordinates. The region EE' is 0r10 \leq r \leq 1, 0θ2π0 \leq \theta \leq 2\pi 0z20 \leq z \leq 2.

E(2x+2y+2z)dV=E2zdV\iiint_E (2x + 2y + 2z)\, dV = \iiint_E 2z\, dV

Since ExdV=EydV=0\iint_E x\, dV = \iint_E y\, dV = 0 by symmetry (odd functions over a symmetric domain).

=202π0102zrdzdrdθ=202π01r[z22]02drdθ= 2 \int_0^{2\pi} \int_0^1 \int_0^2 z \cdot r\, dz\, dr\, d\theta = 2 \int_0^{2\pi} \int_0^1 r\left[\frac{z^2}{2}\right]_0^2\, dr\, d\theta

=202π012rdrdθ=202π1dθ=22π=4π= 2 \int_0^{2\pi} \int_0^1 2r\, dr\, d\theta = 2 \int_0^{2\pi} 1\, d\theta = 2 \cdot 2\pi = 4\pi

\blacksquare

3.7 Conservative Fields and Potential Functions

Definition. A vector field F\mathbf{F} on a domain DRnD \subseteq \mathbb{R}^n is conservative if There exists a scalar function ϕ:DR\phi : D \to \mathbb{R} (called a potential function) such that F=ϕ\mathbf{F} = \nabla \phi.

Proposition 3.8 (Equivalent conditions for conservative fields). Let F=(P,Q)\mathbf{F} = (P, Q) be a C1C^1 vector field on a connected domain DR2D \subseteq \mathbb{R}^2. The following are equivalent:

  1. F\mathbf{F} is conservative: F=ϕ\mathbf{F} = \nabla \phi for some ϕ\phi.
  2. CFdr=0\oint_C \mathbf{F} \cdot d\mathbf{r} = 0 for every closed curve CC in DD.
  3. CFdr\int_C \mathbf{F} \cdot d\mathbf{r} is path-independent in DD.
  4. Py=Qx\frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x} everywhere in DD.

Procedure for finding a potential function. Given F=(P,Q,R)\mathbf{F} = (P, Q, R) with ×F=0\nabla \times \mathbf{F} = \mathbf{0}:

  1. Integrate PP with respect to xx: ϕ=Pdx+g(y,z)\phi = \int P\, dx + g(y, z).
  2. Differentiate with respect to yy and set equal to QQ to determine gyg_y.
  3. Integrate gyg_y with respect to yy: g=gydy+h(z)g = \int g_y\, dy + h(z).
  4. Differentiate with respect to zz and set equal to RR to determine h(z)h'(z).
  5. Integrate to find h(z)h(z) and assemble ϕ\phi.

Problem. Determine whether F=(2xy+z2,x2+2yz,2xz+y2)\mathbf{F} = (2xy + z^2,\, x^2 + 2yz,\, 2xz + y^2) is conservative, And if so, find a potential function.

Solution

Check the curl:

(×F)x=y(2xz+y2)z(x2+2yz)=2y2y=0(\nabla \times \mathbf{F})_x = \frac{\partial}{\partial y}(2xz + y^2) - \frac{\partial}{\partial z}(x^2 + 2yz) = 2y - 2y = 0

(×F)y=z(2xy+z2)x(2xz+y2)=2z2z=0(\nabla \times \mathbf{F})_y = \frac{\partial}{\partial z}(2xy + z^2) - \frac{\partial}{\partial x}(2xz + y^2) = 2z - 2z = 0

(×F)z=x(x2+2yz)y(2xy+z2)=2x2x=0(\nabla \times \mathbf{F})_z = \frac{\partial}{\partial x}(x^2 + 2yz) - \frac{\partial}{\partial y}(2xy + z^2) = 2x - 2x = 0

Since ×F=0\nabla \times \mathbf{F} = \mathbf{0}, F\mathbf{F} is conservative. Find ϕ\phi:

ϕx=2xy+z2    ϕ=x2y+xz2+g(y,z)\frac{\partial \phi}{\partial x} = 2xy + z^2 \implies \phi = x^2 y + xz^2 + g(y,z)

ϕy=x2+gy(y,z)=x2+2yz    gy(y,z)=2yz    g(y,z)=y2z+h(z)\frac{\partial \phi}{\partial y} = x^2 + g_y(y,z) = x^2 + 2yz \implies g_y(y,z) = 2yz \implies g(y,z) = y^2 z + h(z)

ϕz=2xz+y2+h(z)\frac{\partial \phi}{\partial z} = 2xz + y^2 + h'(z)

This must equal 2xz+y22xz + y^2So h(z)=0h'(z) = 0Giving h(z)=Ch(z) = C.

Therefore ϕ(x,y,z)=x2y+xz2+y2z+C\phi(x,y,z) = x^2 y + xz^2 + y^2 z + C. \blacksquare

3.8 Common Pitfalls

:::caution Common Pitfalls

  • Singularities. When applying Green’s, Stokes’, or the Divergence theorem, verify that the field has continuous partial derivatives on the region (including interior). If there are singularities inside the region, the theorems do not apply directly; the singularity must be handled separately.
  • ** connected domains.** The condition ×F=0\nabla \times \mathbf{F} = \mathbf{0} guarantees that F\mathbf{F} is conservative only on a connected domain. For example, F=(y,x)x2+y2\mathbf{F} = \frac{(-y, x)}{x^2 + y^2} has zero curl on R2(0,0)\mathbb{R}^2 \setminus \\{(0,0)\\} but is not conservative there (the domain is not connected).
  • Orientation. Green’s and Stokes’ theorems require positive orientation (counterclockwise for planar curves, right-hand rule for surfaces). The divergence theorem requires the outward normal. Reversing orientation changes the sign of the result.

3.9 Relationships Among the Fundamental Theorems

The three major integral theorems of vector calculus are deeply connected:

Remark. Green’s theorem is the planar special case of Stokes’ theorem. Stokes’ theorem relates the Circulation around a curve to the curl through the surface it bounds. The divergence theorem relates The flux through a closed surface to the divergence inside the volume it encloses. Together, these Form the higher-dimensional analogues of the Fundamental Theorem of Calculus:

abf(x)dx=f(b)f(a)(FTC)\int_a^b f'(x)\, dx = f(b) - f(a) \quad \mathrm{(FTC)}

Cϕdr=ϕ(B)ϕ(A)(FTLI)\int_C \nabla \phi \cdot d\mathbf{r} = \phi(B) - \phi(A) \quad \mathrm{(FTLI)}

CFdr=S(×F)dS(Stokes)\oint_C \mathbf{F} \cdot d\mathbf{r} = \iint_S (\nabla \times \mathbf{F}) \cdot d\mathbf{S} \quad \mathrm{(Stokes)}

SFdS=E(F)dV(Divergence)\iint_S \mathbf{F} \cdot d\mathbf{S} = \iiint_E (\nabla \cdot \mathbf{F})\, dV \quad \mathrm{(Divergence)}

In each case, the integral of a “derivative” over a region equals the integral of the original function Over the boundary of that region. This is the generalised Stokes’ theorem:

Ωω=Ωdω\int_{\partial \Omega} \omega = \int_{\Omega} d\omega

Where Ω\Omega is a kk-dimensional manifold with boundary Ω\partial \Omega, ω\omega is a (k1)(k-1)-form, And dωd\omega is its exterior derivative.

:::