15.

1 Vector Fields (page 554)

CHAPTER 15 VECTOR CALCULUS

15.1 Vector Fields (page 554)

A vector field assigns avector to each point (x, y) or (x, y, 2). In two dimensions F(x, y) = M(x, ~ ) i + N ( xy)j.
,
+
An example is the position field R = x i y j (+x k). Its magnitude is lRl = r and its direction is o u t &om t h e
+
1 2 y2). The level curves are circles, and they are perpendicular
origin. It is the gradient field for f = z(x
to the vectors R.

Reversing this picture, the spin field is S = -y i + x j. Its magnitude is IS1 = r and its direction is a r o u n d
the origin. It is not a gradient field, beeause no function has a f / a x = -y and a f / a y = x. S is the velocity
field for flow going around t h e origin. The streamlines or field lines or integral curves are circles. The flow
field pV gives the rate at which mase is moved by the flow.

A gravity field from the origin is-proportional to F = a/rswhich has IF / = l/r2. This is Newton's inverse
square law. It is a gradient field, with potential f = l/r. The equipotential curves f (x, y) = c are circles. They
are perpendicular to the field lines which are rays. This illustrates that the gradient of a function f (x, y) is
perpendicular to its level curves.

+
The velocity field y i x j is the gradient of f = xy. Its streamlines are hyperbolas. The slope dy/dx of a
streamline equals the ratio N/M of velocity components. The field is tangent to the streamlines. Drop a leaf
onto the flow, and it goes along a streamline.

lf(x,y)=z+2y 5f(x,y)=sin(x+y) 5f(z,y)=1n(x2+d)=21nr

9
+
7 F = zyi + $j, f ( z , y) = 9 = 0 so f cannot depend on x; streamlines are vertical (y = constant)
11F=3i+j lSF=i+2yj lSF=2xi-2yj 17F=ez-'i-eZ-gj
19 % = - l ; y = - ~ + C 2 1 2 = - ~I/'
- ~ 2 + y 2 = C 23$= - I Y = Y~ ; X ~ + ~ ~2Sparallel = C
27 F = y i + yj
29 F = -:fG(xi + mMG
y j) - ((Z-l,a+va)a/r
((x- 1) i + ~ j )
3 1 ~ = $ ~ i - e ~35~$ = S d ~=,/=-
za = - I2-1& I-- 2
35 = gE = g:; = g ~ ;
f(r) = c gives circles

+ +
57 T; F (no equipotentials); T; F (not multiple of xi yj zk)

+
39 F and F i and 2F have the same streamlines (different velocities) and equipotentials (different potentials).

But if f is given, F must be grad f.

Answers 2 - 8 includes extra information about streamlines.

+ sx + +
2 xi j is the gradient of f (x, y) = 1 2
y , which has parabolas i x 2 y = c as equipotentials (they open
down). The streamlines solve dy/dx = l/x (this is NIM). So y = In z + C gives the streamlines.
4 i/y - xj/# is the gradient of f (x, y) = x/y, which has rays x/y = C as equipotentials (compare Figure 13.2;
the axis y = 0 is omitted). The streamlines solve dy/dx = N/M = -x/Y. So y dy = -x dx and the
+
streamlines are y2 x2 = constant (circles).
3 +
6 x2i+ y2j is the gradient of f (x, y) = (xs + ys), which has closed curves x3 y3 = constant as equipotentials.
The streamlines solve dy/dx = y2/x2 or dy/y2 = dx/x2 or y-' = x-' + constant.
8 The potential can be f (x, y) = xfi. Then the field is V f = f i i +;
x j / 4 . The equipotentials are curves
15.2 Line htegrsb (page 562)

- #x2 = c.
10 If *
z f i = C or y = @/z2. The streamlines solve dy/dz = N / M = z/2y so 2y dy = x dz or y2
= -y then f = -yx+ a n y function C(y). In this c w a, = -x+ dC which can't give = x.
12 % = 1and % = -3; F = i - 3j has parallel l i e s z - 3y = c as equipotentials.
14 % = %c - 2 and % = 2y;F = (22 - 2)i + 2gj leads to circles (z - 1)2+ y2 = c around the center (1,O).
16 = iX cos y and 5= -iX sin y ; F = ez (cosyi - sin yl) leads to curves e' cos y = c which stay inside
a strip like Iyl < g. (They come in along the top, turn near the y axis, and leave along the bottom.)
18 %=
X
3 and = 4; +
F = - 3 1 f j has the r a y s f = c as equipotentials (omit the axis z = 0).
20 $ = x gives y = j1x 2 + C (parabolas). 22 $ = - Y gives y2 + x2 = C (circles).
24 2 3
= gives y = 2x +
1 C (parallel lines).

26 f (x, Y) = f h ( x 2+ y2) = ~n dw. his comes from = +orf=J*.

28 The gradient 3z2i + 3y2j is perpendicular. For unit length take F (or V) as (x2i + y 2 j ) / d w .
SO The field is a multiple of i +j. To have speed 4 take F (or V ) as &(i +j).
Sf Fkom the gradient of y - z2, F must be -2xi +j (or this is -F).
34 The slope 2 is -fx/fy from the h t equation. The field is f,i + fJ so this slope is -M/N. The product
with the streamline slope N/M is -1, so level curves are perpendicular to streamlines.
ax + +
36 F is the gradient of f = 1 2 b q i c y 2 . The equipotentials are ellipses if a c > b2 and hyperbolas
if a c < b2. (If a c = b2 we get straight lines.)
+ + +
40 (a) R S = (x - y)i (y z)j has magnitude fir. (b) The magnitude is now fi (difference of
+
perpendicular unit vectors). (c) The direction stays parallel to i j (at 45').
(d) yi is a shear field, pointing in the x direction and growing in the y direction.

15.2 Line Integrals (page 562)

Work is the integral of F dR. Here F is the force and R is the position. The d o t product finds the
component of F in the direction of mwement dR = dx i + dy j. The straight path (z, y) = (t, 2t) goes from (0,0)
+
at t = 0 to (1,2) at t = 1with dR = dt i 2 d t j.

Another form of dR is Tds, where T is the unit tangent vector to the path and the arc length has ds =
\/(dx/dt)t + +a)/&
(dy/dt)2. For the path (t, 2 t ) , the unit vector T is (i and ds = fidt. For F = 3i j, +
+ +
F . T ds is still bdt. This F is the gradient of f = Sx y. The change in f = 3z y from (0,0) to (1,2)is 6.

+
When F = grad f, the dot product F dR is (af /ax)dx (af/*)dy = df. The work integral from P to Q
I
is df = f(Q) - f(P).In this case the work depends on the endpoints but not on the path. Around a closed
+ +
path the work is zero. The field is called conservative. F = (1 y) i x j is the gradient of f = x + xy. The
work from (0,O) to (1,2) is 3, the change in potential.

+
For the spin field S = -y i x j, the work does depend on the path. The path (x, y) = (3cos t, 3 sin t) is a
+
circle with S-dR = -y dx x d y = 9 dt. The work is 1 8 r around the complete circle. Formally / g(z, y)ds is
the limit of the sum g(q,yi)Asi.

The four equivalent properties of a conservative field F = M i + N j are A: zero work a r o u n d closed paths,
 15.2 Line Integrals (page 562)

B: work depends only on endpoints, C: gradient field, D: aM/ay = aN,/ax. Test D is passed by F =
(y + 1)i + 4.The work J F -*dRaround the circle (cost, sin t) is Bero. The work on the upper semicircle equals
the work on the lawer semicircle (clockwise). This field is the gradient of f = x + *y, so the work to (-1,O)
is -1 starting f t o m (0,0).

5 I': (-3 sin t)dt = 0 (gradient field); $ ": -9 sin2 t dt = -9r = - area
7 No, xy j is not a gradient field; take line z = t, y = t from (0,O) to (1,l) and / t2dt # $
9 No, for a ci& (2rr)' #02 +02 11f = X + ? d ; f ( O , l ) - f(1,o) = -;
l t f = ~ ~ ~ ~ ~ ; f ( O , l ) - f ( l , O ) = O 16f=r=dm;f(0,1)-f(l,0)=0
1 7 Gradient for n = 2; after calculation 8U - = ~3
' 'I
1 9 ~ = ~ ~ 0 ~ t , ~ = a ~ i n t , d ~ = a( ad +t a, sM
i n t=) aJd~t = 2 r a 2
$tff
at T= ,a +
2 1 x = a cos t, y = a sint, ds = a dt, M =

20t dt; P -*dR

= l+tl'
= (6ti
a3 cos2 t dt = r a g , (2, p) = (O,O) by symmetry
= 3 x i + 4 j = 6 t i + 4 j , d s = 2J
4). (2 dti + 2t dt j) = 20t dt; work = J:
G~~~,F
= (.6 tT
i +~
4 j )s- ( - $ & ) t t / m d t
20t dt = 30
=

25 If 9 9 = then M = ay + b, N = ax + c, constants a, b, c
27 F = 4xj (work = 4 from (1,0) up to (1,l)) 29 f = [X - (1 1) = -1 2 (1 1)
31 f = [xy ](o:o) = 1
+
53 Not conservative; /'(ti - tj) =(dt i dt j) = J O dt = 0; (t2i - tj) =(dt i + 2t dt j) = J ' -t2dt = -$
as! = ax,
56 iX +
= 2x b, so a = 2, b is arbitrary 37
as!
= 2ye-' = E; f = -y2e-'
~ gas! ~ = ~ = ~ ; f = r = ~ ~ = l ~ i + d l
4 1 ~ = ( x - ~ ) i + ( x + ~h) aj , = - l , E = l , n o f 4821r;O;O

2 ~ o t e d s = t / s i ~ t + c o s 2 t d t = d t .hen ~ x d s = ~ ~ ' ~ c o s t dand f.

t =~l x ~ d s = / ~ ' ~ s i n t c o s t d t =
4 Around the square 0 5 x, y 5 3, dx = -9 along the top (backwards) and :/
-x dy = -9 up the
right side. All other integrals are sero: answer -18. By Section 15.3 this integral is always -2 x area.
6 /$dt=/ds=arc length = 5.
8 Yes The field xi is the gradient of f = ?x2. Here M = x and N = 0 so we have Mdx+ Ndy = f (Q) - f (P).
More directly: up and down movement has no effect on x dx. I
10 N o t much. Certainly the limit of C ( A S ) ~is rero.
12 = 0 and = 1; n o t conservative, take straight path x = 1- t, y = t : /F .*dB= / y dx dy = +
/,1 t(-dt) + dt = f.
14 = a d F is the gradient of f = xeg. Then J P . dR = f (Q) - f ( P ) = -1.
16 # z;not conservative, choose straight path x = 1-t, y = t : / -y2dz+x2dy = /t2dt+(l-t)2dt = ).
18 5 has M = '
(o'+~,.,l and = -zny(x2 + y2)-(n/2)-1.
,2-n
This agrees with R is a
so 7
gradient field for all n. The potential is f = or f = In r when n = 2.
20 The semicircle has x = a cos t, y = a sin t, ds = adt, 0 5 t 5 x. The mass is M = pds = padt = par. I
The moment is M' = / py ds = / pa2 sint dt = 2pa2. Then it = 0 (by symmetry) and 9 = = $.
22 (a) For a gradient field F . dR = f (Q) - f (P).Here Q = (1,1,1)and P = (0,0,O) so f (Q) - f (P)= 2.
+ +
(b) J Mdx Ndy Pdz = /t2dt - t(2t dt) t3(3t2dt) = ). +
24 P = 0 means = 0. So f is f (x, y). So M = 2
and N = cannot depend on z.
26 (a) / y3dx + 3 z d d y = ~ , ' ( y t ) ~ (dt)
x + ~ z t ( y t ) ~ ( y d=t x
) ys. Then = y3 and as! = 3 x d (conservative).
(b) W = I,' ( x t ) ~ ( xdt) + 3(yt) (xt)' (y dt) = (x + tY2x2). But
1 4 # M (not conservative).
15.3 Green's Theorem (page 571)

(c) W = $01 2 (O dt) + s ( y dt) = k2 + f2 . But # M (not consewative).

(d) W = I,' cd+"(z dt + y dt) = eZ+Y- 1. Then = cc'+ and 6Y
= ee'+ (conservative).
28 F = x2j on the circle x = c ~ s t , ~ = s i has
n t J F - ~ R ~=~ ' c o s ~ t ( c odt)
s t= 0.
SO / x2dy = J
: t2dt = but /,' t2(2t dt) = 1. 1

82 # + $01 4
(not conservative): $ x2Ydx x#dy = 2t3dt = but :$ t2(t2)dt t(t2)2(2tdt) = 17 + s.
34 The potential is f = )ln(x2 + +
y2 1). Then f (1'1) - f (0'0) = 1x1 S. f
86 -t2(-2t dt) (1 +
- t2)(2t dt) = 1(as before). On the quarter-circle ending at t = :
J,"'~(- sin 2t) (-2 sin 2t dt) + (cos2t)(2cos2t dt) = 2q = as before. 8
18 - = -2ycZ - 2yC # 0. N o potential f (2, y).
4 0 ~ y'+za
= ~ h a s ~ # ~ .
4 2 BatE = - dy ifandonlyifb=c.Then f ( x , y ) = f a x 2 + b x y + ~ d d .
44 Puc because / F dl?,= $ y dx. False because F = yi is not conservative. (The area underneath
depends on the curve.) !the because the area is 7r (and $ y dx = sin t(sin t dt) = r.)c"

15.3 Green's Theorem (page 571)

+
The work integral f M dx N dy equals the double integral J J ( N ~-My)& d y by Green's Theo-
rem. For F = 3i +
4 the work is sero. For F = and -yi the work equals the area of R. When
M = a f /ax and N = a f l a y , the double integral is zero because fxy = fF. The line integral is zero because
+
f(Q) = f (P) when Q = P (closed curve). An example is F = y i x j. The direction on C is counterclock-
wise around the outside and clockwise around the boundary of a hole. If R is broken into very simple pieces
+
with crosscuts between them, the integrals of M dx N d y cancel along the crosscuts.

Test D for gradient fields is a M / a y = aN/&. A field that passes this test has . dR = 0. There is aIF
+
solution to f. = M and f, = N. Then df = M dx N dy is an exact differential. The spin field 8/13 passes test
D except at r = 0. Its potential f = 9 increases by 2 r going around the origin. The integral I$(N, - M,)dx dy
is not rero but 2r.

The flow form of Green's Theorem is fC M d y - N dx = $ h ( M x + Ny)& dy. The normal vector in F-nds
points o u t across C and In1 = 1and n ds equals dy i - dx j. The divergence of Mi Nj is Mx Ny. For F = + +
zi the double integral is Ill
dt = area. There is a source. For F = y i the divergence is zero. The divergence
of R/? is zero except at r = 0. This field has a point source.

A field with no source has properties E = aero flux through C, F = equal flux across all paths
from P t o Q , G = existence of s t r e a m function, H= zero divergence. The stream function g satisfies
+
the equations ag/* = M and a g / a x = -N. Then aM/ax aN/ay = 0 because aag/azay = aag/ay &.
The example F = yi has g = )y2. There is not a potential function. The example F = xi - yj has g =
xy and also f = )x2 - iy2.
This f satisfies Laplace's equation f= +
fyy = 0, because the field F is both
conservative and source-free. The functions f and g are connected by the Cauchy-Riemann equations
 15.3 Green 's Theorem (page 571)

af = ag/ay and af/& = -ag/ax.

1"$
: (a cos t) a cos t dt = aa2; Nz - My = 1,$1dz dy = area aa2
3$ ~ X ~ X + $ ~ X ~ ~ = O , N , - M , = O , $ $ O ~ X ~ ~ = O

5 $ z2y dz = $:"(a cos t)2(a sint)(-a sin t dt) = -$ $:"(sin 2t)2dt = -$;

Nz - My = -z2,$$(-z2)dz dy = So So -?cos2B r drd9 = -$

2" a

7 $ z d y - y d z = $ ~ ( c o s 2 t + s i n 2 t ) d t = r ; $ $ ( 1 + l ) d z d y = 2 (area) = r ; $ z 2 d y - z y d z = $ + l ;
$l +
:$ (22 z)dz dy =
9 $ $:"(3 cos4 t sin2 t + 3 sin4 t cos2 t)dt = $ $,a* 3 cos2 t sin2 tdt = 2 (see Answer 5)
11/ F . d R , = o around any loop; F = :i+ F j and$/-dR.=$tT[-sintcost+sintcostldt =0;
= gives $$ O dz dy
I S z = cos 2t, y = sin 2t, t from 0 to 21; 1":
-2 sin2 2t dt = -2x = -2 (area);
$tu-2dt = -4n = -2 times Example 7
15 $ ~ d ~ - ~ d z = $ ~ " 2 s i n t c o s t d t = 0 ; $ $ ( ~ ~ + ~ , ) d z d ~ = ~ $ 0 d z d ~ = 0
1 7 M = :,N = ~ , $ M d y - N d z = $ ~ " ( ~ o s ~ t + s i n ~ t ) d t = 2 a ; $ $ ( ~ ~ + ~ , ) d z d y = $ $ ( f - $ + f - f ) d z d y =
$$ !dz dy = $$ dr dB = 2 r
19 $ Mdy - Ndz = $ -zZy dz = $ : -z2(1 -
x)dz = &;
z2dz dy =
2 1 $$(M. +
Ny)dz dy = $/ div F dz dy = 0 between the circles
+
2 3 Work: $ a dz b dy = $$(g
- g ) d z dy; Flux: same integral
25 g = tan-'($) = 9 is undefined at (0,O) +
2 7 Test My = Nz : z2dz y2dy is exact = d($z3 + iy3)
29divF=2y-2y=O;g=zy2 31divF=2~+2y;nog SSdivF=O;g=eZsiny
35 div F = 0;g = $
3 7 Nz - My = -22, -6zy, 0,2z - 2y, 0, -2eZ+y; in 31 and SS f = $(z3 y3) and f = eZcosy +
3 9 F = ( 3 ~ ~ - 3 ~ ~ ) i - 6 z ~ j ; d i v F =4O1 f = z 4 - 6 z 2 y 2 + y 4 ; g = 4 z 3 y - 4 z y 3
43 F = eZcosy i - ezsiny j ; g = e2siny

45 N = f(z), $ ~ d + Ndy
z
' f (1)dy +
=$ f(0)dy = f(1) - f(0); $$(N, - M,)dz dy =

$$ E d z dy = I,' dz (Fundamental Theorem of Calculus)

2 f z a y d y = $ ~ ~ a 2 ~ o s 2 t ( a s i n t ) ( a= ~O~; M
~ t=dOt ), N = z 2 y , ~ / 2 z y d z d y =
"$: :$ 2r cos B(r sin 9)rdr dB = 0
4 fydx=$,'t(-dt)= -#;M=y,N=o,$$(-l)dzdy=- area= -- 1
2'

6 f z2y d2 = $:(I -rn;

- t ) ~ t ( - d t ) = 1 M = zay, N = 0, Jo1 Jo1-y -z2dz dy = -
W d y= -$.
8M =2 + + +
3 and N = z2y 22 SO f M d z Ndy = J$[(2xy 2) - 2xyldz dy = 2 times area.
+
10 M = by and N = cz : f M d z Ndy = $$(c - b)dz dy = (c - b ) t i m e s area;
b = 7 and c = 7 make the integral sero.
1 2 Let R be the square with base from a to b on the z axis. Set F = f (x)j so M = 0 and N = f (x). The
+
line integral Mdz Ndy is ( b - a)f (b) up the right side m i n u s ( b - a)f (a) down the left side. The
double integral is $$ g d z dy = ( b - a) $,b gdx.Green's Theorem gives equality; cancel b - a.
,
8 S - dB = f -y dx
14 $ + x dy since the integrals along the axes are zero. By Green's Theorem this is
$$2dx dy = 2 times area between path and axes.
16 f F . nds = I zy dy = 3 up the right side of the square where n = i (other sides give zero).
Also + 0)dz dy = i.
1 8 In the double integral M, = K
a x2+v
-Y () = * = a n~d N , = & (~
-~+)=~& '+!I
15.4 Surface Integrals (page 581)

so M2 + Nu = 0 : Double integral = 0. Along the bottom edge (where y = 0 and n = -j) the line integral
is F e n d s = 1,' = -1. Theright side ( x = 1a n d n = i) yields J: -
d i e= = 1- -dm]: 6.
a
,*

2 +O

Back across the top ( y = 1,n = j , notice ds = -dx!) = fi - l.-Down the left side (notice
ds = -dy !) gives +l. Adding the four sides I
n ds = 0.
20 F = grad r = F) has F . n = 0 along the + axis where n = -j and y = 0. On the unit circle
(2,
n is equal to F (unit vector pointing outward) so F .n = 1. Amund the semicircle
F . nds = j," 1 d8 = r . The double integral has M2 = $(:) = f - ja, - ; r-;
rr a = f . Similarly
zr

!.
N, = a6u( n+ ) = $ and Mz + Nu = f = The double integral is I,' i(r
dr d8) = T .
22 f F n ds is the same through a square and a circle because the difference is jJ(E+ +)dz dy =
div Fdx dy = 0 over the region in betwacn.

,*
- .

24 j(cos3 y dy - sin3 x dx) = J J ( 0 - O)dxdy = 0. A different example would be more revealing.

26 div 3=B (+~) +
82 2 +Y
= ,m = O. Integrating gives g = #ln (x2 +y2) = h r -
This is infinite at x = y = 0.
28 & = M and 2 = -N are compatible when M, N, = g y - + ~ gxy = 0. I f also N2 = My then
+
+ g y y = -N2 Mu = 0 and g solves Laplace's equation.
SOaM+snr= 3y2 - 3 3 = 0. Solve
Bz 6y
8
= 3xy2 for g = xys and check = y3.2
32 + BY = 0 + 0. Solve 9~
BY
= y2 for g = f Y3 + C ( x ) and add C(x) = $ x 3 to give = x2. 2
Then g = I(y1 3 +xS).

8% +
34 aM = ez+u - ez+r = 0. Solve a = ez+' for g = eX+Y and check 2= eZ+Y.
16 + + +
BY
= y x # 0 (no stream h c t i o n ) .
38 g(Q) = J ' F -ndr starting from g(P) = 0. Any two paths give the same integral because forward on one
and back on the other gives f F .n ds = 0, provided the tests E - H for a stream function are passed.
40 With M2 + N, = 0 we can solve ag/ay = M = 3x2 - 3y2 and ag/ax = - M = 6x9 to find
g = SYC2y -ys. Then f, = g, = A4 and f y = -g, = N.
42 Mdy - Ndx is an exact differential if =- z.
(Then there is a stream function g.)
4 4 ~ I ~ . d R = f - ~ d x + x d ~x =area#O.
2
46 Simply connected: 2, 3, 6 ( ? ) , 7. The other regions contain circles that can't shrink to points.

15.4 Surface Integrals (page 581)

A small piece of the surface a = f ( x , y) is nearly flat. When we go across by dz, we go up by (az/ak)dx. That
+
movement is A d z , where the vector A is i d z / d x k . The other side of the piece is Bdy, where B = j (az/ay)k.+
+
The cross product A x B is N = -az/ax i - as/* j k. The area of the piece is dS = INldx dy. For the
surface x = xy, the vectors are A = fl J-dx dy and N = -y i -x + k. The area integral is
IldS=i+yk.

With parameters u and u, a typical point on a 45' cone is x = u cos u, y = u sin v,z = u. A change
+
in u moves that point by A d u = (cos v i + sin v j k)du. The change in u moves the point by Bdu =
(-u sinv i + u cos v j)dv. The normal vector is N = A x B = -u cos v i - u sin v j u k. The area is +
dS = 6 udu dv. In this example A -B- 0 so the small piece is a a rectangle and dS = IAIIBldu du.

For flux we need ndS. The u n i t n o r m a l vector n is N = A x B divided by INI.For a surface z = f (z,Y ) ,
 15.4 Surface Integrals (page 581)

the product n d S is the vector N dz dy (to memorize from table). The particular surface z = zy has n d S =
+ + +
(-yi - qj k)dz dy. For F = zi yj zk the flux through z = zy is F .n d S = -xy dz dy.

On a 30' cone the points are z = 2ucos v, y = 2usin v, z = u. The tangent vectors are A = 2 cos v i
+ 2 s i n v j + k and B = - 2 u a i n v i + 2 u c o s v j . This cone has n d S = A x B du dv = ( - 2 u c o s v i - 2
+ + +
u sin v j 4 u k ) d u dv. For F = x i y j z k , the flux element through the cone is F . n d S = zero. The
reason for this answer is that F is along t h e cone. The reason we don't compute flux through a Miibius strip
is that N c a n n o t b e d e 5 e d (the s t r i p is n o t orientable).

1N = -221 - 2yj k; d S = + dl +
4z2 4 9 dz dy; + So2% so2 d w r dr d0 = f ( 1 7 ~ 1-~1)
3~=-i+j+k;d~=fidzd~;area&%

9 N = (y2 - z2)i - Zzyj + k; = dl + (y2 - z2)2 + 4z2y2dzdy = dl + (y2 + z2)2dz dy;

5;" 5; 6.- r dr d0 =
Jz +
+ +
11N = 2i 2j k; dS = 3dz dy; 3(area of triangle with 22 2y 5 1) = f +
13 r a J W 1 5 $:-'5; zy(& dz dy) = $

17 5:"~:~~ sin2)cos)sin0cos0(sin) d) d0) = 0 1 9 A = i + j + f k ; B = j+k; N = -i-j+k; d S = f i du dv

2 1 A = -sinu(cosv i + s i n v j ) + c o s u k ; B = - ( 3 + c o s u ) s i n v i + (3+cosu)cosvj;

N = - ( 3 + c o s u ) ( c o s u c o s v i + c o s u s i n v j + s i n u k ) ; d S = (3+cosu)dudv
23 !!(-ME - N% + +
P)dz dy = I/(-2z2 - 2y2 z)dx dy = $$ -?(r dr dB) = -8a
25 F . N = - z + y + z = O o n p l a n e
27N=-i-j+k,F=(v+u)i-uj,$5~.Nd~=$$-vd~dv=0
29 55dS = so so
2n 2n
+
(3 cos u)du dv = 12a2 3 1 Yes 3 3 No
3 5 A = i + flcosOj+ f 1 s i n 0 k ; B = - f s i n B j + f c o s B k ; N = f f ' i - f c o s 0 j - f s i n ~ k ; d ~ = l N l d z d # =
f ( z ) J W dz d0

2 ~ = - 2 x i - 2 ~ ~ + k a n d d ~ = ~ ~ d x d ~ . ~ 2nh Jzefind $s r$ dd r ~d =0 =
~ ~
g ( J J J / ~- 17312).

AN=-3i-qj+k a n d d s = m d x d y . Thenarea=$;$,'adzdy=&.

6N=-
4-
xi - L + k a n d d S =
,/- ,/&.d x d Then area = so2 I~l l1-=
,~
/-
[-2*JT-7'i]:,, = fir.

~ N r= r-+Ik - and
d S = ~ V z fi dx dy. Then area = J02n Jab a r d r dB = fia(b2 - a2).
d dy =
1 0 N = -i - j + k and dS = 4 dx dy. Then surface area = times base area = 2 h .
1 2 z = d= gives N = + k and dS = c* Then area = 4 J"

dx a dy

1 4 N = -2zi + k and dS = d a d x dy. Area = /z2

5: d-dz dy = 4 J: , / w d x =
=/*,

85: 4-dx = 8[:4- + + In 1% + Jwll~ = 1 2 6 E + In 13 +1- -

(In +) = 39.
16 On the sphere d S = sin 4 d4 d0 and g = x2 + y2 = sin2 4. Then I.
so
2n n/2
sin3 4 d4 d0 = 2 4 5 ) = k.
1 8 x = 2 cos v, y = 2 sin v, and d S = 2 d u dv. Then $5 g d S = so .fo 2 cos v(2du dv) = 0.
2lr 3
3

20~=vi+j+k,~=ui+j-k,~=~~~=-2i+(u+v)j+(v-~)k,d~=,/4+2~~+2v~dud~.
22 ~ = c o s v i + s i n u j =
, ~- u s i n v i + u c o s v j + k , =~ sinvi- cosvj + u k , d ~= f i d u d v .
24 $$ F . n d S = 5;" s2&
-r3dr d0 = -24a. 26 5sF ndS = .f5
0 d S = 0.
15.5 The Divergeace Thmrem (page 588)

+ + + + +
28 F .ndS = ((u u)i - uvj). (-21 (u u)j (u - u)k)du du = (2u 2v - u2v - v2u)du dv.
Integrate with u = r cos 0, u = r sin B :So 2% 1
+
(2r cos B 2r sin B - r3 cos2 B sin B - ? sin2 B cos B)rdrdB = 0.
30 A = cosBi+sin#j - 2rk,B = -rsinBi+ r c o s # j , N = A x B = 2?cosBi+2r2sinBj+rk,
2u a
$ $ k . n d ~ = / $ k . ~ d u d u = / ~I, r d r d B = r a 2 a s i n ~ x a m p l e 1 2 .
1 2 I think a 'triple Mijbius stripn is orientable.
+
1 4 The plane z = ax by has roof area = I,/- times base area. So choose for example a = 1 and b = a.

15.5 The Divergence Theorem (page 588)

In words, the basic balance law is flow in = flow out. The flux of F through a closed surface S is the
+ + +
double integral I/ F . ndS. The divergence of Mi Nj Pk is Mx + Ny Pg, and it measures t h e source
at t h e point. The total source is the triple integral /fl
div F dV. That equals the flux by the Divergence
Theorem.

For F = 5zk the divergence is 5. If V is a cube of side a then the triple integral equals 5aS. The top surface
where z = a has n = k and F n = 5a. The bottom and sides have F . n = gero. The integral // F .ndS equals
5a3.

The field F = has div F = 0 except a t t h e origin. // F ndS equals 4 r over any surface around the
+
origin. This illustrates Gauss's Law: flux = 4 r times source strength. The field F = x i y j - 22 k has div
F = 0 and / 'F ndS = 0. For this F, the flux out through a pyramid and in through its base an equal.

+ +
The symbol V stands for (a/&)i (a/ay)j (a/az)k. In this notation div F is V . F. The gradient of f is
Vf. The divergence of grad f is V Vf o r v 2 f . The equation div grad f = 0 is Laplace's equation.

+
The divergence of a product is div(uV) = u div V (grad u ) V. Integration by parts in 3D is
$11 u div Vdx dy dz = - 1's +
V . grad u dx d y d~ /$ u V n dS. In two dimensions this becomes
+ + +l
$lu(aM/& aN/t3y)dx d y = - / ( M au/& N au/t3y)dx d y u V .n ds. In one dimension it becomes
integration by parts. For steady fluid flow the continuity equation is div pV = +/at.

l d i v F = l , $ $ f d ~ = !f SdivF=2~+2y+2z,$$$divFdV=0
2u u/2 a
5divF=3,$$3d~= :=;
7 F N = p2, $Jp=. $dS = 4ra4 9 div F = 22, Jo Jo 2pcos lo
sin 4 dp d) dB) = Bra4

11:/ /:(2x + 1)dx dy dz = a' + a3; -2a2 + +a2 + O + a' + 0 + a3

1 3 d i v ~ =: , / ~ ~ : ~ V = O ; F . ~ = X , / ~ X ~ 1S6 d
= iO
vF= l;$llldv= f;$l/ldV=)
R div R
1 7 d i ~ ( ~ ) = ~ - + & - g r a dg $- S=R - g r a d p
19 Two spheres, n radial out, n radial in, n = k on top, n = -k on bottom, n = *
JG on side;
n = -i,-j,-k,i+a +
3k on 4 faces; n = k on top, n = l ( k i + f j - k ) on cone
Jz
+
2 1 V = cylinder, // / div F dV = l l ( z +)dz dy (a integral = 1); F ndS = I/
/ Mdy - Ndx, z integral = 1on side, F .n = 0 top and bottom; Green's flux theorem.
23 div F = -tyM
= -4rG; at the center; F = 2R inside, F = 2(:)3R outside
2 ~ d i v u , = ~ P, ~ = ~ , / ~ ~ - n d ~ = / / l d 2 7~F =( d4i vnF = O ) ; F ; T ( F - n s 1 ) ; F
 15.6 Stokes' Theorem and the Curl of F (page 595)

29 Plane circle; top half of sphere; div F = 0

2 /IF.ndS=$$/0dV = O
1 1 1

4 $$F.P~s=$,$ , $o(2z+2y+2z)dzdydz=1+l+1=~.
6 $1 F .ndS = (directly) $1 dS = 4ra2. By the Divergence Theorem: $ lo
, 2 s s a ;p2 2 sin 4 dp d4 de = 4ra2
8$$~.nds=$, I,
2s a a
.

$, 3 p 4 ~ i n 4 d p d 4 d t 9 12n
=Ta
b

1 0 div F = O + zersinz - z ~ s i n = z 0 SO $$F * n d S= 0.

+
1 2 An integral over a box with small side a is near ca3. Here div F = 22 1has integral a' + a3, which is
near a3 because a is small. Then c = 1,which equals div F on the plane x = 0.
1 4 R - n = ( z i + g j + z k ) . i = z = 1 onone faceofthe box. On the fiveotherfaces R-n=2,3,0,0,0.
The integral is :$ ldy dz :$ 1 + + + +
,' 2dz dz :$ f 3dz dy = 18.Also div R = 1 1 1= 3 and
# 3dx dy dz = 18.

16 The normal vectors to the cube are n = fi, fj,f k . Then $1 I,' 1:
F -ndS =
+
z dx dy
1 +
: $. (-z)dx dy I, + + + I,'
$', 'z dz dz I,' J,' (-z)dx dz $', $', Ody dz $. ldy dz = 1.

Also $$/ div FdY = $, So 1 1 1

$, ldz dy dz = 1.

18 grad f .n is the directional derivative in the normal direction n (also written 2).
The Divergence Theorem gives $11 $l
div (grad f ) dV = grad f . n d S = $$ g d ~ .
But we are given that div (grad f) = j. + + j,, f, is sero.
-
20 Suppose F is perpendicular to n on the surface; then $$ F ndS = 0. Example on the unit sphere:
+ .
F is any q( z, y, z) times the spin field -yi xj
22 The spin field F = -yi + z j has div F = 0 and F - n = 0 on the unit sphere.
24 The flux of F = RIP3 through an area A on a sphere of radius p is Alp2, because I F I= 113and F
is outward. The spherical box has A/$ = sin)d) dB on both faces (minus sign for face pointing in).
No flow through sides of box perpendicular to F. So net flow = rero.
26 When the density p is constant (incompressible Bow), the continuity equation becomes div V = 0. If the Bow
is irrotational then F = grad f and the continuity equation is div (p grad f ) = -dp/dt.
If also p = constant, then div grad f = 0: Laplace's equation for the 'potential."
28 Extend E-F-G-H in Section 15.3 to 3 dimensions: E The total flux $$ F - n d S through every closed
surface is rero F. Through all surfaces with the same boundary $IF - ndS is the same
G There is a stream W d g for which F = curl g H. The divergence of F is sero (this is the quick test).
SO The boundary of a solid ball is a sphere. A sphere has no boundary. Similarly for a cube or a cylinder - the
boundary is a closed surface and that s d a c e ' a boundary is empty. This is a crucial fact in topology.

15.6 Stokes' Theorem and the Curl of F (page 595)

+ + + +
The curl of Mi N j Pk is the vector (Py - Ng )i (MI - Px)j (Nx - My )k. It equals the 3 by 3
i j k
+
determinant a/& slay alas The curl of x2i z2k is xero. For S = yi - (z z)j yk the curl is + +
M N P
2i - 2k. This S is a spin field a x R = f ( curl F) x R, with axis vector a = i - k. For any gradient field
f.i+ +fuj fsk the curl is sero. That is the important identity curl grad f = rero. It is based on f., = f,,and
15.6 Stokes' Theorem and the Curl of F [page 595)

fu = frx and fys = fry. The twin identity is d i v curl F = 0.

The curl measures the spin (or turning) of a vector field. A paddlewheel in the field with its axis along
n has turning speed f n . c u r l F. The spin is greatest when n is in the direction of c u r l F. Then the angular
1
velocity is ZI c u r l FI.

Stokes' Theorem is fCF . dR = JJS(curl F ) n dS. The curve C is the boundary of the surface S. This is
Green's Theorem extended to t h r e e dimensions. Both aides are sero when F is a gradient field because t h e
curl is sero.

The four properties of a conservative field are A : f F dR = 0 and B : Q F dR depends only

Ip
o n P and Q and C : F is the gradient of a potential f'unction f (x,y,r ) and D : c u r l F = 0.
+
The field y2z2i 2 x d z k fails test D. This field is the gradient of n o f. The work J F dB from (O,O,O)to
(1,1,1) is along t h e straight path x = y = = t. For every field F, JJ curl F . n d S is the same out through
a pyramid and up through its base because t h e y have t h e same boundary, so f F dl?,is t h e same.

lcurlF=i+j+k ScurlF=O 6curlF=O 7f = ~ ( X + ~ + Z ) ~

9 c u r l x m i = O ; xnj h a s r e r o c u r l i f n = O 11c u r l F =2yi; n = j on circle so I $ F - n d s = ~
13curl~=2i+2j,n=i,~$curl~.nd~=~$2d~=2r
15 Both integrals equal / F 6&;Divergence Theorem, V = region between S and T, always div curl F = 0
1 7 Always div curl F = 0 +
1 9 f = xz y 2 1 f = eZ-' 25 F = yk
+ +
25 curl F = (a3b2 - a2b3)i (al b3 - a3bl )j (a2bl - al b2)k 27 curl F = 2wk; curl F =
6 = 2w/,/3
+ + +
29 F = x(a3z a2y)i y(alx a3z)j z(alx a2y)k + 2% %/a
+
31 curl F = -2k, JJ -2k R d S = JO Jo -2 cos d ( s h 4 d4 dB) = -2r; J y dx - x dy =
(- sin2 t - cos2t)dt = -2r
2rr rr/2
SS curl F = 2a,2JJ(alx+ a2y+ a3z)dS = 0 + O + 2a3JO Jo cos)sind d$ dB = 2ra3
55 c u r l F = - i , n = m4, j / ~ - n d ~ =
-&r?
57 g = f - f = stream function; sero divergence
19 div F = div ( V + W) = div V so y = div V so V = $-j (has rero curl). Then W = F - V = xyi - 2 J
4 1 curl (curl F) = curl (-2yk) = -2i; grad (div F ) = grad 2x = 21; F., + F,, + Fzz= 4i
-=
43 c u r l E = sB = a s i n t so E = ?(a ~ R ) s i n t

2 curl F = 0 because curl of gradient is always sero. +

4 curl F = -i j - k from equation (1).
6 curl F = 2i + 2j from Example 2: curl (a x R) = 2a.
+
8 f ( x , y, r ) = ~ + ' / 2 ( n 1) has grad f = pnR (so its curl is rero).
10 curl (alx + a2y + a3z)k = a2i - a jwhich is sero when a1 = 0 and a2 = 0.
1 2 curl (i x R ) = 21 directly (or by Example 2 with a = i). Then f F . dR = curl F . n d S = 0 since n = j is
perpendicular to i.
1 4 F = (x2 + y2)k so curl F = 2(yi - xj). (Surprise that this F = a x R has curl F = 2a even with nonconstant
a.) Then F dB = JJ curl F n d S = 0 since n = k is perpendicular to curl F.
1 6 C is the equator (the common boundary of S and T);V is the whole ball (the earth). Note that n doesn't
point out in the bottom half T, or the direction around C would be opposite.
 15.6 Stokes' Theorem and the Curl of F (page 595)

For F = R (position vector), $Js F ndS = - $JT F ads.

18 If curl F = 0 then F is the gradient of a potential: F = grad f. Then div F = 0 is div grad f = 0
which is Laplace's equation.
20 The potential is f = x2y. 22 The potential is f = xyz )z3.+
24 Start with one field that has the required curl. (Can take F = f i x R = -gj + gk). Then add any F with
curl s a o (particular solution plus homogeneous solution as ahays). The fields with curl F = 0 are
+
gradient fields F = grad f , since curl grad = 0. Answer: F = i i x R any grad f.
;
26F = yi-xk has curlP = j -k. (a) Angular velocity = curl F - n = i f n = j . i
(b) Angular velocity = 1;urlFI = 9 (c) Angular velocity = 0.
28 One possibility: F = *k has curl F = spin field 5. Other possibilities: F = v + k any grad f.
M) False (curl F = curl G meam curl (F - G ) = 0 but not F - G = 0). !hue (curl ( F - 6)= 0 makes F - G
a gradient field). False ( F = xi + gj + zk and G = 0 have the same curl (nero) but div F = 3).
i j k
32 curl R / =~ a~/ a ~ a/ay a/& h a i component z g p - 2 - yKp-2 8
= 0. Similarly for j and k:
xlp2 ylp2 .Id
t h u s c u r l F = ~ a n d ~ ~ c u r l F ~ n d S =(separately)
~and I F - d R = f ~ d z + ~ d y = $ x d xydy=O.
+
34 Based on Problem 47 of Section 11.3, the triple vector product (a x R ) x R is F = (a R ) R - ( R . R ) a =
(ax + by + cz)R - (x2 + y2 + z2)a. Then by Problem 42 b of this section, or directly, the curl is
grad ( a x + b y + c z ) x ~ - g r a d ( x 2 + y 2 + z 2 ) x a = a x ~ - 2 ~ x a = ~ a N x ~o. w ~ ~ c u r l ~ . n d ~ = ~
R
since n = igi is perpendicular to the cross product curl F = 3a x R.
Also, f F . ' d i L = ~ ( a . ~ ) ~ - d( R .-R ) a - d R = 0 because R . ~ R = oon the circle and R - R = 1.
i j k
36 curl F = +
a/ay a/a+ = i(xz) j(1- YZ)+ k(1) and n = xi + + zk. SOcurl F . n =
Z x syz
-
x2z + y y2z + z. By symmetry ll x2z dS = IS y2zdS on the half sphere and ydS = 0.
2u r/2
This leaves J J z d S = J o Jo cos#(sin#d#dB) = i(2a) = r .
+
38 (The expected method is trial and error) F = 5yzi 2xyk any grad f . +
-
40 Work = f B dR = jJ (curl B) ndx dy = - pJ .ndx dy. So work is p times current through C.
42 (a) curl vi = 3j - e k . Then curl (curl F) = (-$ - 9)+ i &j +&k. Also
grad (div F ) = ei+ +
&k. The difference is (v. + +
vy, v,.)i. Note: The same steps
for the j and k components give identity (a) for any F. My favorite is to square this matrix:
][ ] [
[2
(b) curl (fvi) = (f.v
curl
-div
+
grad
0
= curl curl - grad div

+
0 -divgrad
O I
= V21!!

fv,U - (fv, f,v)k. This is f curl F = f ( v j - v,k) added to (grad f ) x F =

f,v j - f,vk. Again the identity extends to any F.
+ +
44 P x G = (Np - Pn)i (Pm - Mp)j (Mn - Nm)k. Its divergence is the sum of z, y, and z derivatives:
1N.p + Np. - P.n - Pn.] + + + +
[P,m Pm, - Mvp - Mp,] [M,n Mn, - Nzm - Nm,]. Note that m
-
multiplies P, - Nz,the first component of curl F. This starts G *curl F - F curl G , as we want.
46 False. Certainly G x F would be perpendicular to F but V x F is something different. For example P = i yk +
h a s V x F = i s o ( V x F ) - F = 1.
15.6 Stokes' Tlreonm and the Curl of F (page 595)

48 S = roof, its shadow = ground floor, C = edge of roof, shadow of C = boundary of ground floor. Similarly
+ + i.
for spherical cap P y2 z2 = 1 above z = Note C is on the plane z = and its shadow is a circle
around the shadow of the cap, down on the plane z = 0.
60 curl V = curl (-xk) = j. A wheel in the zz plane has n = j so it spins at full speed. A wheel perpendicular
to j will not spin, if it is in the r y plane with n = k.
MIT OpenCourseWare
http://ocw.mit.edu

Resource: Calculus Online Textbook

Gilbert Strang

The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.