A Comprehensive Guide To Vector Analysis
A Comprehensive Guide To Vector Analysis
WW
bat of
= we etna roa
_ r
, as a
ies y:
- ipod!
“Sd 4 ‘i
— =a
ghee
i “=
org
KA. Stroud
Cale
x
with additions by
Dexter Booth
K. A. Stroud
Formerly Principal Lecturer
Department of Mathematics
Coventry University
Dexter J. Booth
Principal Lecturer
School of Computing and Engineering
University of Huddersfield
JUN 09 2006
7 wwbibrary of Congress Cataloging-in-Publication Data
Stroud, K. A.
Vector analysis / K.A. Stroud, Dexter Booth.
Ducky
ISBN 0-8311-3208-6
1. Vector analysis--Textbooks. I. Booth, Dexter J. II. Title.
QA433.S77 2005
515’.63--dc22
2005043301
Contents
Triple integrals 86
Applications 88
Review exercise 92
Alternative notation 93
Determination of volumes by multiple integrals 99
Review summary 102
Can You? Checklist 4 105
Test exercise 4 105
Further problems 4 106
Learning outcomes
Double integrals
Surface integrals
Space coordinate systems 163
Cartesian coordinates 163
Cylindrical coordinates 164
Spherical coordinates 165
Element of volume in space in the three coordinate systems 167
Volume integrals 168
Alternative method ive
Change of variables in multiple integrals Fs?
Curvilinear coordinates 179
Transformation in three dimensions 187
Review summary 189
Contents
Learning outcomes
Introduction: scalar and vector quantities
Vector representation
Components of a given vector
Components of a vector in terms of unit vectors
Vectors in space
Direction cosines
Scalar product of two vectors
Vector product of two vectors
Angle between two vectors
Direction ratios
Triple products
Properties of scalar triple products
Coplanar vectors
Vector triple products of three vectors
Review summary
Can You? Checklist 7
Test exercise 7
Further problems 7
Answers 340
Index 345
Hints on using
this book
This book contains lessons called Programs. Each Program has been written in
such a way as to make learning more effective and more interesting. It is like
having a personal tutor because you proceed at your own rate of learning and
any difficulties you may have are cleared before you have the chance to
practise incorrect ideas or techniques.
You will find that each Program is divided into numbered sections called
frames. When you start a Program, begin at Frame 1. Read each frame carefully
and carry out any instructions or exercise that you are asked to do. In almost
every frame, you are required to make a response of some kind, testing your
understanding of the information in the frame, and you can immediately
compare your answer with the correct answer given in the next frame. To
obtain the greatest benefit, you are strongly advised to cover up the following
frame until you have made your response. When a series of dots occurs, you
are expected to supply the missing word, phrase, number or mathematical
expression. At every stage you will be guided along the right path. There is no
need to hurry: read the frames carefully and follow the directions exactly. In
this way, you must learn.
Each Program opens with a list of Learning outcomes which specify
exactly what you will learn by studying the contents of the Program. The
Program ends with a matching checklist of Can You? questions that enables
you to rate your success in having achieved the Learning outcomes. If you
feel sufficiently confident then tackle the short Test exercise which follows.
This is set directly on what you have learned in the Program: the questions are
straightforward and contain no tricks. To provide you with the necessary
practice, a set of Further problems is also included: do as many of these
problems as you can. Remember, that in mathematics, as in many other
situations, practice makes perfect —- or more nearly so.
vil
Digitized by the Internet Archive
in 2024
https ://archive.org/details/vectoranalysisOO0Ostro
Useful background
information
Symbols used in the text
= is equal to — tends to
& is approximately equal to a is not equal to
> is greater than = is identical to
= is greater than or equal to < is less than
n! factorial w=1% 2x3 X24. x7 < is less than or equal to
|k| modulus of k, i.e. size of k 00 infinity
irrespective of sign Lim limiting value as n — co
% summation ie
ix
Useful background information
atan$
tan @ = x
1 7 tan
GED C18
(e) sinC +sinD
= 2sin 5 cos 5)
C+D_,.G—p)
SinG = sia /D) = Zeos 5 sin 5)
D E = ID)
Ak Resist cos 5)
D C=D
cos D — cosC = 2sin + sin 5)
3 Standard curves
(a) Straight line
Slope, NE
m =—— Cayean
Dae dx XN,
ie. xy=C"
sls
where ¢c =
4 Laws of mathematics
on
(a) Associative laws — for addition and multiplicati
a+(b+c)=(a+b)+c
a(bc) = (ab)c
on
(b) Commutative laws — for addition and multiplicati
a+b=b-+a
ab = ba
division
(c) Distributive laws — for multiplication and
a(b+c) =ab+ac
BEG DLC
265= (provided a # 0)
a a a
i uit CRED Mt wifardabion
Jitiye=. ). dias
[gh ie
et
~ ya
— ay
| = ea) ted
oF ce Oe ET: Sas
* wi ee *
Pe chkiigteaae
Ha FS ae aay ieee
ivi nage, - enainaps
yur a . rer ae Rie :
ws
, a5 ¥
oak
eu)
a | ;
ai Ae a ag ge
bi. o= a | ap eel ta Bins
rn ra >
tg Sef
1
=. y “Xves® vt Heftmae ws
It is now over 35 years since Ken Stroud first developed his approach to
personalized learning with his classic text Engineering Mathematics, now in its
fifth edition. That unique and hugely successful programmed learning style is
exemplified in this text and I am delighted to have been asked to contribute to
it. I have endeavored to retain the very essence of his style that has
contributed to so many students’ mathematical abilities over the years,
particularly the time-tested Stroud format with its close attention to
technique development throughout.
The first two Programs deal with the elements and application of partial
differentiation and present a standard approach with a thorough coverage.
The next Program introduces polar coordinates to acclimatize the student to
the possibility of different coordinate systems. Program 4 extends the
students’ awareness of the integration of functions of a single variable to
integrals of two and three variables. Program 5 then discusses the concept of a
differential and applies this to line integrals, culminating in Green’s theorem
in the plane. Program 6 looks at surface and volume integrals including the
change of variables; the Program ends with an introduction to curvilinear
coordinates. The next three Programs deal with the essence of the vector
calculus, covering the algebra of vectors, vector fields, their calculus, and the
integral theorems. The final Program deals with curvilinear coordinates as
applied to all that was learned in the previous Programs dealing with the
vector calculus.
To give the student as much assistance as possible in organizing their study
there are specific Learning outcomes at the beginning and Can You?
checklists at the end of each Program. In this way, the learning experience is
made more explicit and the student is given greater confidence in what has
been learnt. Test exercises and Further problems follow, in which the student
can consolidate their newly-found knowledge.
This is the fourth opportunity that I have had to work on the Stroud books,
having made additions to both the Engineering Mathematics and Advanced
Engineering Mathematics texts. It is as ever a challenge and an honor to be able
to work with Ken Stroud’s material. Ken had an understanding of his students
and their learning and thinking processes which was second to none, and this
is reflected in every page of this book. As always my thanks go to the Stroud
family for their continuing support for and encouragement of new projects
and ideas which are allowing Ken’s work an ever wider public.
xiii
— |
Partial ee
differentiation
Learning outcomes
When you have completed this Program you will be able to:
e Find the first partial derivatives of a function of two real variables
e Find second-order partial derivatives of a function of two real variables
e Calculate errors using partial differentiation
2 Vector Analysis
Partial differentiation
If we keep r constant and increase the height h, the volume V will increase. In
these circumstances, we can consider the derivative of V with respect to h —
but only if r is kept constant.
Ee
Len |= | is written aM
dh |, constant Ox
6 d
Notice the new type of ‘delta’. We already know the meaning of = and =.
OV OVE, ; ee ;
Now we have a new one, Dh oh called the partial derivative of V with
respect to h and implies that for our present purpose, the value of r is
considered as: Dello ke tien. eetcner
constant
OV
V =ar7h. To find ahr WE differentiate the given expression, taking all symbols
ee ey = 72h = 2a
Or
In the statement V = ar2h, V is expressed as a function of two variables, r and
h. It therefore has two partial derivatives, one with respect to
one with respect to
es
Partial differentiation
3
Another example:
OA OA
ar = 2rh and Oh = 2ur
Of course, we are not restricted to the mensuration of the cylinder. The same
will happen with any function which is a function of two independent
variables. For example, consider z = x’y°. }
Oz *
Here z is a function of x and y. We can therefore find = and ay
GZ =X 25.2
—— ape)
Oy x 3"oy = 3X y
constant
Example 1
u=x?+xy+y?
Ou
(a) To find Dx , we regard y as being constant.
Example 2
z=
4+y? — 2x*y
OL VE
ame +0 — 4xy = 3x° — 4xy
OZ _
or ee 2 KV
Dene) x 2
Example 3
Z = (2x — y)(x + 3y)
This is a product, and the usual product rule applies except that we keep
hast
y constant when findingtga (Ey
ee and x constant when finding oe
Ox Oy
OZ
Aion (2x — y)(1 + 0) + (x + 3y)(2 — 0) = 2x —y + 2x + 6y = 4x + Sy
OZ
Sp ae OE SE VO ea 3p =X =3y = 5x — 6p
Partial differentiation 5
Oz OZ
ay ak + 14y; ie 14x — 20y
Example 4
Example 5
OZ Oz
HZ ==)!sin(sx- 2y)2y) find
iim —
ayant and —
By
OZ
(b) By’ we treat x as constant.
constant
2 Z=(3x+2y)(4x
m x — Sy)
a oar
ayo a7 y os
6y— 7xX —— 20 2Vy
OZ Oz
3 z == t tan(3x + 4 4y) x
—_ = 3 sec* 2 (3x + 4y) ay
— =
4 sec“ 2 (3x + 4y)
ee sin(3x + 2y)
xy
Oz _ 3xcos(3x + 2y) — sin(3x + 2y) Oz 2ycos(3x + 2y) — sin(3x + 2y)
Ox x2y Oy xy2
If you have got all the answers correct, turn straight on to Frame 15. If you
have not got all these answers, or are at all uncertain, move to Frame 12.
Powe axes
Bxy 4 Sy"
To find Meregard y as a constant:
OZ OZ
ay ee Le, OX 4-OV eg ay,
Similarly, regarding x as constant:
3 z=tan(3x + 4y)
OZ
Now have another go at finding By in the same way.
Then check it with Frame 14
Here it is:
Pp sin(3x + 2y)
xy
OZ xy cos(3x + 2y) - (2) — sin(3x + 2y) - (x)
iy 0 xe
* 2y cos(3x + 2y) —sin(3x4-2y)
= oy?
That should have cleared up any troubles. This business of partial differentia-
tion is perfectly straightforward. All you have to remember is that for the time
being, all the independent variables except the one you are using are kept
constant — and behave like constant factors or constant terms according to
their positions.
Ny Og ae
The expression as 6x + 4y is itself a function of x and y. We could therefore
find its partial derivatives with respect to x or to y.
Partial differentiation
Oo {oz O a2
>- 4=~? and this is written - (much like an ordinary second derivative,
Ox | Ox Ox2
but with the partial 0)
coz 2
xz = ay OX + 4Y) =6
So we have this:
Z = 3x 4xy — Sy”
Oz OZ
~Ox = 6x + 4y pee
ay x 1) y
Cpa sl
Ox2
O92 d Kia
Oy - OX
OZ
Of course, we could carry out similar steps with the expression for Dy on the
right. This would give us:
2
Deeks t0
Oy?
2
fe ie
OX - Oy
N ote hat
tha ay.ax™
of eans 5 (5 ay'e means
10 Vector Analysis
Oz OF
Ox? ‘ Oy?
2
OW ay 2
| Ow yi
Oy - OX Ox - OY
Oz oO Zz
We see in this case, that By ak ae ay There are then, two first derivatives
and four second derivatives, though the last two seem to have the same value.
Here is one for you to do.
Oz Oz Oz Oz Oz OZ
If z=53
5x? + 3x*y y 3 find
+ 4y”, £3 Ax’ Dy’ Ox! By! Ox.Oyand By Ox
Oz Oz
Oz OZ
= —— = 6x
Oy - OX Ox - OY
a ; Oz Oz
Again in this example also, we see that —_ - = ———_. Now do this one.
Oy-OX Ox- Oy
It looks more complicated, but it is done in just the same way. Do not rush at
it; take your time and all will be well. Here it is. Find all the first and second
partial derivatives of z= xcosy —ycosx.
Then to Frame 19
—— ee ee ee
Partial differentiation
11
SO we get:
Oz ; OZ
5x 7 COSY ty: sinx jy = 7X Siny — cosx
a2 Q2
LEE Beers Maat
Ox? Oy?
ByozOx —siny yv+ sinx Aira
sis —siny + sinx
i
Pz Oz
And
feet eet Oy - Ox Ox: Oy
In fact this will always be so for the functions you are likely to meet, so that
there are really three different second partial derivatives (and not four). In
a2 2
practice, if you have found 5 Z it is a useful check to find Z separately.
Ox Ox - OY
They should give the same result, of course.
OV ae V.
If V = In(x* a + ae y*), prove that Boat = aye =3 0
This merely entails finding the two second partial derivatives and substituting
them in the left-hand side of the statement. So here goes:
V = In(x?
+ y”)
Ass Tae
Ox (x2
+ y2)
2x
Te
OV 2y* — 2x?
We had found that
Ox (2 + y2y?
So making a fresh start from V = In(x + y?), we get:
OVOy x (x7+y")
a1 alysis bay
Ae
OV (x? +y?)2 = 2y-2y
oy (x2 + y2)?
_ 2x? + 2y? —4y? 2x? — Dy? (b)
aye ee
Substituting now the two results in the identity, gives:
OV MOV ye ey
Ox2 ' Oy2 (x2 + y2)? (x2 + y2)?
_ 2y? — 2x% 42x? — 2? _
0
pe pete ea in ee
Now on to Frame 22
Example 1
OVe + HOV
If V = f(x?f(x? +y")
+ y2), show that
x— — y— = On
Here we are told that V is a function of (x* + y*) but the precise nature of the
function is not given. However, we can treat this as a ‘function of a function’
and write f’(x? + y) to represent the derivative of the function w.r.t. its own
combined variable (x? + y”).
Oe 8 —
By Hh? +) x 0? +97) =P 02 +9") 22
OV | ifn? 2 7 O. 2 DNA ene? 7EN
yr thie + y) op EV" Vie (Kore Peay
oe ei tat a eal al ee
< OVem ov Wifa 52, 2 Ve eager? 2
Example 2
= f'(ax+by)-a=a-f'(ax+by) (a)
OZ ; O
i
= f'(ax+by)-b=b-f'(ax+by) (0)
J pe - i = ab - f'(ax + by) — ab - f'(ax + by)
= ()
Move on to Frame 25
a
14 Vector Analysis
So to sum up so far.
dz
Ox
ot
Oy
On oy?
Oz OZ Oz Oz
Oy-Ox Ox Oy ee Oy- Ox Ox- OY
Now for a review exercise
F| Review exercise
1 Find all first and second partial derivatives for each of the following
functions:
(a) 2= 3x? + 2xy + 4y?
(D) 2 = sinxy
(c) Z= Ay
rca £
2 If z=In(e* +e”), show that cae os=< ik
Ox Oy
OZ OZ
Seliz=xa7 f (xy), — — Yay
(xy), express xa y— in
in its
its simplest
si form.
O77 =— 6
Oz =
Ox2 Oy?
Pz O*z
eae = Z =
Oy - Ox Ox - OY
Partial differentiation 15
(b) z= sinxy
Oz
Pe ORIOP Ss
Oz
py weeay
Oz (x—y)1—(a+y)1 a 2y
OX (x-y ss (x-y?
Oz (x—y)1—(x+y)(-1) 2x
oy (x-y)? (x—y)?
Oz (—2) Ay
2 Sy aaeene
(x=p) 7 (X=y)
&z ‘ee ra 4x
a ays (xy)
Oz — (x—y)?(-2)
— (—2y)2(x -—y)(-1)
SAE
LY ay eae ey
— (x-yP (x=?
Pz _ (x—y) (2)— 2x:2x-y)l
Ox- Oy (x —y)*
_ x-y) - 4x0-y)
(x-y)*
Zs 4x
(x-y)’ (x-y) 3
| Bx = 2y — 4x ex = Zy
~~ (x-yy ay)?
16 Vector Analysis
242 ain(e 2)
oz_ 1 Oia!achwan at
Ox e+e dy ete
OA OI A eY
Ox Oye+e e+e
BN =
ler they B
On,Oe_
Ox soy
3 Z=x-f (xy)
ax FO) “y+foy)
aePON 3
OZaa Oo,
oe ee es = x’ yf"(xy) + xf (xy) — x°yf" (xy)
OZ OZ
te aye ey) eC
That was a pretty good review test. Do not be unduly worried if you made a
slip or two in your working. Try to avoid doing so, of course, but you are doing
fine. Now on to the next part of the Program.
So far we have been concerned with the technique of partial differentiation.
Now let us look at one of its applications.
So move on to Frame 28
Small increments
OF =
= SRO
2ath; Oh = TI
If r becomes r + ér, and h becomes h + 6h, let V become V + 6V. Then the new
volume is given by:
V + 6V = n(r + 6r)*(h + 6h)
= n(r* + 2rér + [5r]*)(h + 6n)
= n(r7h + 2rhér + hior)’ + r26h + 2réréh + [5r]?6n)
Subtract V = xr*h from each side, giving:
6V = x(2rhér + hl6r)? + r6h + 2r6r6h + [5r|°5h)
~ (2rhér + 126h) since ér and 6h are small and all the remaining
terms are of a higher degree of smallness.
Therefore
On to Frame 29
6V = 1007(0.2) + 252(—0.1)
=A = Aon = Wien
*. 6V = 54.98 cm?
i.e. the volume increases by 54.98 cm?
Just like that!
aS
18 Vector Analysis
, but to any
This kind of result applies not only to the volume of the cylinder
function of two independent variables. Here is an example :
by small
If z is a function of x and y, i.e. z = f(x,y) and if x and y increase
y small. If we expand 6z
amounts 6x and 6y, the increase 6z will also be relativel
in powers of 6x and dy, we get:
6z = Aéx + Béy + higher powers of 6x and éy,
where A and B are functions of x and y.
If y remains constant, so that 6y = O, then:
6z = Aéx + higher powers of 6x
This is the key to all the forthcoming applications and will be quoted over and
Over again.
The result is quite general and a similar result applies for a function of three
independent variables. For example:
If z= f(x,y,w)
OZ Oz Oz
then en oz6z = Bye a ay!
- I aye
Example 1
V
hi R and V = 250 volts and R = 50 ohms, find the change in J resulting
from an increase of 1 volt in V and an increase of 0.5 ohm in R.
Ol Ol
f(V,R)
T=FVLR 2 tlhe BPN tae
—
Example 2
ws3
ii y= = , find the percentage increase in y when w increases by 2 per cent,
s decreases by 3 per cent and d increases by 1 per cent.
Notice that, in this case, y is a function of three variables, w, s and d. The
formula therefore becomes:
a Ow + ons
by = aL ay. 4 =28 6d
55
We have
2 : ; 2 _ 2w
Because éw is not 700 of a unit, but 2 per cent of w, i.e. dw = 700 of w= 00
= =35 d
Similarly,
imi 100 of s = 100 and 6 d — ——.
6s = —— 100 Now that we have cleared that
SWS We 2 9 4
me |100" 100 100
i
= 4See
waaSs SS = ML percent.oty
f
Move on to Frame 35
: fi » 2 1 1 4
oD
ciiergatis hdd+50sik raa ‘s+P (=55)
‘*. Maximum possible error in P is 4 per cent of P
Finally, here is one last example for you to do. Work right through it and then
check your results with those in Frame 37.
The two sides forming the right-angle of a right-angled triangle are denoted by
a and b. The hypotenuse is h. If there are possible errors of +0.5 per cent in
measuring a and b, find the maximum possible error in calculating (a) the area
of the triangle and (b) the length of h.
Ooh 1 2 2,4 a
a ae b 2(2, =
Oa 3 4 <a Va+b
ea
Oh 1 2
OE ed) ee 2,—4 b
h =a
<__(4+57-)
200/
+Va? : +b? (+555)
\ 200
z 1 ae b*
~~ 200 Va? + b2
toa}
=t559 V4 Bh +b pes= +50
oh |
iy,)
That brings us to the end of this particular Program. We shall meet partial
differentiation again in the next Program when we shall consider some more
of its applications. But for the time being, there remain only the Can You?
checklist and the Test exercise.
So on now to Frames 38 and 39
(4 Can You?
Checklist 1
Check this list before and after you try the end of Program test.
Test exercise 1
b) If z=y f(x
z=f(x++e ay)+F(x—ay),
= : One
2 a?
Oo”
(b) If find x2 and a2 = and hence prove that
Cee OZ
az ° x2"
Partial differentiation 23
2
3 The power P dissipated in a resistor is given by P = 7
If E = 200 volts and R = 8 ohms, find the change in P resulting from a drop of 5
volts in E and an increase of 0.2 ohm in R.
If p= KHL Va, where k is a constant, and there are possible errors of +1 per cent
in measuring H, L and V, find the maximum possible error in the calculated
value of 6.
That’s it
2 2
Prove that, if V = In(x? + y?), then td - a
2
If z = sin(3x + 2y), verify that a me =rOZ,
O22» Or,
If z=e*(xcosy—ysiny), ae oe =
92 2
If u = (1+ x) sinh(5x — 2y), verify that po +20 gH o 250 a0
ae By. Oy Oy?
50°Z Oo On
ifz=f(2), show that 755 + 2x ae esa
Rooke OZ
one oy
3
10 In the formula D = aaa h is given as 0.1+0.002 and v as 0.3 + 0.02.
12(1 — v7)’
Express the approximate maximum error in D in terms of E.
2
11 The formula z = ae ee is used to calculate z from observed values of
y"
x and y. If x and y have the same percentage error p, show that the percentage
error in z is approximately —2p(1 + 2). >
24 Vector Analysis
s
12 In a balanced bridge circuit, Rj = R2R3/Ra. If R>, R3,R4 have known tolerance
of +x per cent, +y per cent, +z per cent respectively, determine the maximum
percentage error in R;, expressed in terms of x, y and Z.
13 The deflection y at the centre of a circular plate suspended at the edge and
kwd 4 :
uniformly loaded is given by y= *; , where w = total load, d= diameter
By _ 2 By
Ot2 —— x2
Show that one solution of this equation is y =Asin2. sin(pt +a), where
€
A, p, c and a are constants.
18 If y= Asin(px + a) cos(qt + b), find the error in y due to small errors 6x and 6t in
x and t respectively.
19 Show that ¢ = Ae /? sin pt cos qx, satisfies the equation
Po 1 {a3 a k?
D2 2) ope + Kog fr Provided that p? = c*q* — mah
OV OV ORV '
20 Show that (a) the equation x2 + 5y2 + al O is satisfied by
1 02V 02V
= Wiiaapiaee? , and and th that (b) the equation
ionDaz + +52 =— (0
25
26 Vector Analysis
Introduction
a
In the first part of the Program on partial differentiation, we established
result which, we said, would be the foundation of most of the applications of
partial differentiation to follow.
You surely remember it: it went like this:
If zis a function of two independent variables, x and y, i.e. if z = f(x, y), then
OZ Oz
6z =——6x +—6
Ox Kioy y
We were able to use it, just as it stands, to work out certain problems on small
increments, errors and tolerances. It is also the key to much of the work of this
Program, so copy it down into your record book, thus:
OZ OZ
liz = (x,y) thentoz — Evie ceehO4
OZ OZ
Za 7 Cey)tbel o7— ayo cig
ner CZ, OZ in :
In this expression, ax and By are the partial derivatives of z with respect to x
and y respectively, and you will remember that to find:
OZ
(a) ay WE differentiate the function z, with respect to x, keeping all
OZ (
rs = sec? (x? = y’); == —2y sec? (x? = y’)
Because Z=tanie 9)
Hi OZ 2(-p2 2 6) 2 2
eae = SEC (AY SN
Ox
= sec2(x” — y”)(2x) = 2x sec”(x? — y*)
and Bp= CP) «5 -7)
0. 9 (
O72
pee ey 2=4-e
Ds A pet ey
Axe e
Pz ae ae 2x—3y
ya mes &(Gan) 2x—3y
aa Wein
2
0 Z te. e2x—3y a6 e2x-3y
OX - Oy
All correct?
You remember, too, that in the ‘mixed’ second partial derivative, the order
of differentiating does not matter. So in this case, since
Oz 2x—3yJoa Oz =RO
& 6-e
Gi. , itt) then By O he CES Cie
Ox- oy
Pz s. Oz a Gis e2x3y
Ox-Oy Oy- OX
Well now, before we move on to new work, see what you make of these.
Find all the first and second partial derivatives of the following:
(a) 2 — <olny
(b) z= (x + y) In(xy)
When you have found all the derivatives,
check your work with the solutions in the next frame
28 Vector Analysis
; OU sity Oz _ xcosy
BOX oy
ie =U) em = —xsiny
Ox? Oy?
oi = COS) a = Cosy
Oy - OX Ox - Oy
ea al jad 1
Oy-ax x i were ae
Vamiae
xy
(ee a er
Ox: dy y xy’ y xX
ee
SF
Well now, that was just by way of warming up with work you have done
before. Let us now move on to the next section of this Program.
Applications of partial differentiation 29
Rate-of-change problems
Let us consider a cylinder of radius r and height h as Sy
before. Then the volume is given by
V=arrh
ae OV 2
We, = 27 and a
dv
_av dr |av ah
di
~ dr dt’ Anstdé
This result is really the key to problems of the kind we are about to consider. If
we know the rate at which r and h are changing, we can now find the
corresponding rate of change of V. Like this:
Example 1
The radius of a cylinder increases at the rate of 0.2 cm/s while the height
decreases at the rate of 0.5 cm/s. Find the rate at which the volume is changing
at the instant when r = 8 cm and h= 12 cm.
Warning: The first inclination is to draw a diagram and to put in the given
values for its dimensions, i.e. r= 8 cm, h = 12 cm. This we must NOT do, for
the radius and height are changing and the given values are instantaneous
values only. Therefore on the diagram we keep the symbols r and h to indicate
that they are variables.
30 Vector Analysis
Here it is then:
V=arh
OV OV
dv ov dr AV dh
dt Ordt oh dt
OV ecOVe Rees,
Or — 2 th; Oh =e
dV dr 2 dh
dh = 2ath = + Tr dt
So you can now substitute these values in the last statement and finish off the
calculation, giving:
dV 5
af 20.1 cm?”/s
dv dr dh
B ecause dE
——_ = 2ath
Al aK:+r
—— ai
2 —
Example 2
B= Ve Bae) Nie
; OZ Oz.
4s OVA = By Ox + ape (The key to the whole business)
dz_ oz dx Oz dy
dt Ox dt’ dy dt
: az A et x
n
In th this case an
——— Stas a (x7Baz yoy2 2 2%) = ee
A Oe ae y
ay a ia ) 2(—2y) a D y
dz dx y dy
dz
dence
dz
apt 4.75 cm/s
Ze 2 2 3 3
Because we have
ieee ) aa )
re a? ete) e 1LOde Quite
—. A i 474 7475 cm/s
Example 3
The total surface area S of a cone of base radius r and perpendicular height h is
given by
S=ar +arvr2 +h
If r and h are each increasing at the rate of 0.25 cm/s, find the rate at which S is
increasing at the instant when r = 3 cm and h= 4 cm.
Do that one entirely on your own. Take your time: there is no need to hurry.
Be quite sure that each step you write down is correct.
Then move to Frame 14 and check your result
32 Vector Analysis
os Las) aes mh
—=1r= 22) — Vre + he
© oh ieee ey
eee 7125
Slee
di
Also we are given that Ss= 0:29 and =,1045)
GS) O47) Lael
dip ese5
167 omy LOT
eS ee
=-3.89 = 11.94 cm?2/s
So there we are. Rate-of-change problems are all very much the same. What
you must remember is simply this:
(a) The basic statement
OZ
LZ T= = ((xXey)
fy) tha then:oz ee = éy
ax On tay (a)
(b) Divide this result by 6t and make 6t — 0. This converts the result into the
form for rate-of-change problems:
dz_dzdx|dz dy
dt dx dt’ dy dt
:
tb)
The second result follows directly from the first. Make a note of both of
these in your record book for future reference.
Then for the next part of the work, move on to Frame 16
ee
Applications of partial differentiation 33
Partial differentiation can also be used with advantage in finding derivatives of | 16-
implicit functions.
d
For example, suppose we are required to find an expression for = when we
x
are given that x? + 2xy+y? =0.
We can set about it in this way:
Let z stand for the function of x and y, i.e. z= x? + 2xy + y3. Again we use
the basic relationship 6z = = ox + ~ by.
If we divide both sides by 6x, we get:
6z Oz : Oz dy
dx Ox Oy 6x
: Om. 02. Og uy
N f 6x —- 0, —=—4—-=
ea Oe cay = oy oly
Oz
If we now find expressions for < and we shall be quite a way towards
7) oy’
finding = (which you see at the end of the expression).
If only we knew ES we could rearrange this result and obtain an expression for
dx
dz
ee So where can we find out something about ane
b
Refer back to the beginning of the problem. We have used z to stand for
x? + 2xy + y? and we were told initially that x? + 2xy + y? = 0. Therefore z = 0,
Z
i.e. z is a constant (in this case zero) and hence = 0.
dy
, O= (2x4 2y) + (2x + 3495,
On to Frame 18
eS aL ei. Seana
34 Vector Analysis
dy 24-4 2y
dx 2x+3y?
dz_d2
02 dy[dz _,
abe aeee Wa abe”
abe
Make a note of this result.
Example 1
d
Ife’ +x+y=1, evaluate at (O, 0), The function can be written
Cre X= Vy — 10.
But zeeds a0
dx
yee
dx
= 8bral
x-eVY41]
dy il d
Atx=0,) y =0, ES == oe
i 1 ee
naa 1
All very easy so long as you can find partial derivatives correctly.
On to Frame 20
Applications of partial differentiation 35
Example2
' rey si
If ay -+ sin “s 2, find
y= =a
ay
Let z=xy+siny—2=0
aut Oz Oz
6z=— Zz be se
ay
UlanOn, Oe dy
dx "Ox “ayp® dx
Bzorden
te edz,
tear eed
oe dy
ee TREO) ae
dz
B Th
ae =a arr
== 0
Wis rey
" dx x+cosy
Here is one for you to do:
Example 3 2
Find an expression for — when x tany = ysin x. Do it all on your own.
Then check your working with that in Frame 21
dy _=tany —ycosx
dx —-xsec2y — sinx
Did you get that? If so, go straight on to Frame 22. If not, here is the working
below. Follow it through and see where you have gone astray!
xtany = y sinx Y, xtany—ysinx=0
Let Z=xtany—ysinx=0
Oz OZ
6Z are eG
dz_ 02 0% dy
dx Ox dy dx
Oz d
OZ Oz 2 f
—=tany—ycosx; =—=xsec*y—sinx
Ox Oy
C2
: ay any —V cosx) + (xsec 2 y —sinx)
in x)a
dz
Dit a Oe ax
dy _—itany
— ycosx
dx xsec?y
— sinx
On now to Frame 22
i
36 Vector Analysis
the
Right. Now here is just one more for you to do. They are really very much
same.
Example 4
d
If e+” = xy”, find an expression for rr
dz oz -0z dy
dx 0x, Oy dx
So continue with the good work and finish it off, finally getting that
oe.
Aye ae
dy 2xy?- et”
dx ett” — 2x2y
Because
z = e*+’— x7y =
= 7 aaky Sa
a ORY
TAS
Apps (Ory
es
4 (EY
ete dy
( — 2y2y)hae<7
dz
Buty Os
we = Nala cmt
dy (e#” = 29?)
OT ee (7 eee)
beady. -2xy7 Os
“ dx (etty
— 2x2y)
That is how they are all done. But now there is one more process that you
must know how to tackle.
So on to Frame 24
Applications of partial differentiation 37
Change of variables
If z is a function of x and y, ie. z=f(x,y), and x and y are themselves
functions of two other variables u and v, then z is also a function of u and v.
_ Oz Oz.
= Boe Bs Pm GBA Y
If v is kept constant for the time being, then a when éu — 0 becomes 2 and
éy oy
54 PeComes =.
—b A
OZ OZ 0% OB» OV
du Ox Ou dy Ou
Note these
Ger ve ee ey
and
Ov Ox OV Oy OV
Next frame
Example 1
0
If z= x? +y*, where x = rcosé and y = rsin 26, find ands
Oz OZ Ox OZ dy
Or Ox Or Oy Or
Of OZ O% £0Z. Oy
and 96
dx 06" ay 00
OZ
a
OZ
late AAG)
as fr |
rec sé Fe = sin 28]
OZ
as= 2x cos 9 + 2ysin 20
And in these two results, the symbols x and y can be replaced by rcosé@ and
rsin 26 respectively. Ee
Sr
38 Vector Analysis
=) y x-cos(u — v) |
u+v
OL OZ. OX" OZ, Oy . —— 5G o BG — COS(U — V
)
.
(
—
=e" le = x-cos(u— »)
u+vV
Now move on to Frame 27
Here is one for you to do on your own. All that it entails is to find the various
partial derivatives and to substitute them in the established results:
OZ) OZ OXe FOZ OY OZ OZ OX, FOZ, OY
du Ox Ou’ dy Ou a ey ae dy av
So you do this one:
When you have completed the work, check with the result in Frame 28
Z = Sin Ve ex = 1b Ve y= Duy,
ae
a ae(x+y); : PE= cos(x
By + y)
Ox _ oy Me 9%
Ou u
EDS Oe oy
du Ox du’ dy du
= COS(X + y) -2u + cos(x + y) + 2v
= 2(u+ v) cos(x + y)
OZ) OZ OX ez aay.
Al : ;
°° Ov Ox Ov Oy OV
Ox ad Oy
ae Dp oe 2u
OZ
ae cos(x + y) -2v+cos(x+y)-2u
= 2(u+Vv) cos(x + y)
aa
a
ee
Applications of partial differentiation 39
Example 2
If z= x? — y? and
x =rcos@ and y=rsin#, then
Oz OZ Ox Oz Oy
Or @x Or Ody Or
Of 02 DX Oz Oy
d = Z ‘
ee = 86 Ba. OB “By oe
We now need the various partial derivatives
Ox
eae OR Bp = 6089 —= sind
OZ Ox y
By Yi ae rsind rcosé
; OZ
— = 2x cos 8 — 2ysin9;
Or
e = —(2xr sin @ + 2yr cos @)
OZ
= 2r (cos? 6 — sin? @); = = —4r’ sin @cos6 |
Or
e OZ
Similarly ae Pe ena and panipo oe te
40 Vector Analysis
Oz
peek
S Oz ;
| Ap 22) Ar cos 26;ee
a0 = r Si Ip
Because
OTe
Ape ORNOZ )
( )= age sin 20) = —4r* cos 20
2
and nas = - (3) ae (2) sin.20) = —4rsin20
Example 3
If z= f(x, y) and x =}(u? — v) and y = wy, show that
BOZ BOLT, wea Oz
Ov Ou Oy Y Ox
Although this is much the same as the previous example, there is, at least, one
difference. In this case, we are not told the precise nature of f(x, y). We must
remember that z is a function of x and y and, therefore, of u and v. With that
in mind, we set off with the usual two equations.
Applications of partial differentiation 41
Whereupon <=
If we now multiply the first of these by (—v) and the second by u and add the
two equations, we get the desired result.
Oz -- uve 2 OZ
Ou Ox Oy
BAM taeBRS 2 OE
OV Ox Oy
"Ove AGL Oy Y Ox
With the same given data, i.e.
and
42 Vector Analysis
Because
O O a) O fa) >)
o = (ugt+v5) and 5 ( va Op
Likewise a ea pe pe age
’ Av2-Av \Av) Ov Ox Oy
Because
a yz 4 0% s dz Oz 0 (Oz
Ove Ov\ ax May) Ox dv \ax) | Ya (=)
=e v( een OF ae Ae oye
Ox Ox Oy ] OX Ox : "Oy Oy
ee oer CZ a Oz a OZ 1 gee
ax ax2 axdy ’ bxdy 7 By?
Wah AYA Sibel CL Were OrF
— Bv2 Ox” ae ae OxOy Oy?
Adding together results (1) and (2), we get
Inverse functions
Example1
éu _ a
du a
du . (1)
1
OV. OV
év as a
oe =a
ta (2)
2
Ghfe) By oc Ptie
; eet
aye Ov
ax aii (ees es 20, re , ay ea C0 2s Ce 661658 e
g 9
ou=e cosy; By = esiny
areca ; EW Bematosy
ane i oy
See —xX
CBROSY e"sinyLY
cos 2y cos 2y
Because
(yee = COSY: y 6u = cos” y 6x — siny cosy éy
e~* cos
(4) x e*sin y: e“siny 6v = — sin’ y 6x + siny cosy dy
Adding: e-*cosy 6u+e*siny év = (cos? y — sin” y) 6x
—x
% ee cos y ersiny
cos 2y cos 2y
Ox Ox
But Stoel ane d
Oke 4 COSY, Ox e*siny
” Ou cos 2y : dv cos2y
which are, of course, two of the expressions we have to find.
Starting again with equations (3) and (4), we can obtain
A Pee
Se
ft ELS ay : ercosy
cos 2y cos 2y
Because
oy ov
éy = oPbu + ay
oy _e*siny aes Oy _ e*cosy
du cos2y i av cos 2y
General case
If z=f(x, y) with u = g(x, y) and v= h(x, y), then we have
OuOu
ou = ey— (1)
il
OV OV
év acdd
_ 6x + odd
ay (2)
2
es
3
ou OX OV
(ie
OVBUate
2
OY OX» OY
Starting afresh from (1) and (2) and eliminating 6x, we have
46 Vector Analysis
by =
Ou Ov Ov Ou
OKeOy Ok Oy
ou ov
ive j= bu 2 and is often written as sa
Oy Oy
Ou OV
_ Ou,v) _ |x ax
2 = aG,y)— fou av
dy Ooy|
Applications of partial differentiation 47
du dv an. ay
Ou Ov
By ay
Ou ov
by by
We can now get a number of useful relationships.
wa Oy Ou
Similar]
imilarly Eveyes
ay =i
Summary
lig = firpj)ecand »u= 9(4;-y). and, v= h(x, y), then
ax_
du Oy
av, ax__
Ov
auf
sy
oy __ avf,ay_au
Ou — OX Ov OX
where, in each case
Ou OV
— O(u,v) _ |Ox Ox
J O(x,y) |Ou ov
dy oy
Let us put this into practice by doing again the same example that we
started with (Example 1 on page 43), but by the new method. First of all,
however, make a note of the important summary listed above for future
reference.
enn
eee
Vector Analysis
Example 1A
Ifz=f(x, vy) and u=ecosy and v=e“*siny, find the derivatives
OX © ORR OVER OY.
du’ Ov’ Ou’ OV
ue =e Cosy v=esiny
Example 2
Ox
If z=f(x,y) with u=x?—y? and v=xy, find expressions for
Ou’
oxOv’ ayOu’ ayOv’
First we need
Ou _ Ou
By : By
a du_
Ox
Because
du 2024p) bv x24 y2
Ove emieVaeners! ge.” -Xx
Ou 2(x2+y?)’ avs x2 + 2
Example 3
If z=2x7+3xy+4y? and u=x*+y? and v=x+2y, determine
CX OX OY le OY OZ OZ
(a)
we LE OY LF
Ou’ Ov’ Ou’ dv ©)
Bare
ou “MN!
Ee
By
Section (a) is just like the previous example. Complete that on your own.
50 Vector Analysis
ax 1 Ox =e BOY = ee ee
av = oxy! du 2(2x—y)’ Ov 2x-y
Ox Ou =)
a | y |2(2% V map ay
Oy Ou ‘- (eee:
oh a /I= 2x/2(2x 1 are
Ox Lea ox ae OF, SL Oy ee
~ Ou Py AMG Ix ey Mon TO ee
Now for part (b).
; : Oz OZ
Since z is also a function of u and v, the expressions for au and =p are
u
OZ 19 OZ OX. 02,09
du Ox Ou ‘ Oy “Ou
OZ 02 COX 02 SOV
dv
Ox Ov dy av
Because
n>)GZ. OR Oz. OY
du
Ox Ou Oy Ou
Be ynGeaee* allay Gy
Ov Ox OV” dy OV
: Of = x
atta b+ ax+ ey {A
Review summary
You have now reached the end of this Program and know quite a bit about
partial differentiation. We have established some important results during the
work, so let us list them once more.
1 Small increments
2 ey eager. (a)
2 Rates of change
2
epee eZescl5 oe
OZ dy
cA b
dt dx dt’oy dt (b)
3 Implicit functions
gz Oz oz dy
ie
= ote Oy Os
(c)
4 Change of variables
me Ore OX FOZ" OY
du ax ou dy Ou
OL 02 0x 02 OY
(d)
dv Ox ov dy Ov
52 Vector Analysis
5 Inverse functions
=f(%y) u=sy) v=hy)
LCM hes ot 3 fh
a uw I;
is Ou
a = /J; oe = /}
Olu, Vv) _ “ Ox
where /= By)
y) On av
dy oy
All that now remains is the Can You? checklist and the Test exercise, so
move to Frames 57 and 58 and work through them carefully at your own
speed.
4 Can You?
Checklist 2
Check this list before and after you try the end of Program test.
On a scale of 1 to 5 how confident are you that you can: Frames
e Derive the first- and second-order partial derivatives of a
function of two real variables? JioCz)
Yes No
e Apply partial differentiation to rate-of-change problems? (Cs) (23)
Yes No
e Apply partial differentiation to change-of-variable problems? to (39)
Yes No
e Use the Jacobian to obtain the derivatives of inverse functions
of two real variables? (40) +» Gs)
Yes | | a No
Applications of partial differentiation 53
Oz Oz Oz
(<) 7 Oxdy ~ “ax dy’
6 Two sides of a triangular plate are measured as 125 mm and 160 mm, each to
the nearest millimeter. The included angle is quoted as 60° + 1°. Calculate the
length of the remaining side and the maximum possible error in the result.
7 Vz (xe = yy? and x is increasing at 3.5 m/s, determine at what rate y must
change in order that z shall be neither increasing nor decreasing at the instant
when x¥=5 mand y=3 m.
8 x? + +44xy+ 3y2
If 2x7 3y? = 1,1, obtain expressions
p oP and ay
for ay ae:
Further problems 2
If F = f(x,y) where x = e“ cos v and y = e“sinv, show that
OF OF OF OF OF _OF
ae soe OF ae ov Yay tay’
: Th
Given that z=x>+y? and x?+y?=1, determine an expression for a, im
terms of x and y.
dy OZ OZ
If z= f(x,y) = 0, show that re ax! Dy’ The curves 2y* +3x—8=0 and
x3 + 2xy3 + 3y — 1 =0 intersect at the point (2, —1). Find the tangent of the
angle between the tangents to the curves at this point.
If u = (x? — y*)f(t) where t= xy and f denotes an arbitrary function, prove
2
that A (x? — y’){t- f"(t) + 3f'(O)}. [Note: f(t) is the second derivative
of f(t) w.r.t. t.]
If V = xy/(x2 +y2)* and x =rcosé, y=rsin6, show that
PV 10V 18V_
Ore + Or 1? 062
If u = f(x,y) where x = r? — s? and y = 2rs, prove that
Ou OH es Ou
Te gt
If f = F(x,y) and x = re’ and y = re~®, prove that
wor Wor A0F of of _of
“Ox ee a aed Or 00)
ik 7= sell 2 2s —1 Ng c Oz OZ
Z=xIn(x* + y*) — 2ytan (5) eo Oe ere
en es BP prey
By means of partial differentiation, determine ay in each of the following cases:
x
(a) xy+ 2y—x=4
(b) x°y? — 2x2y
+ 3xy” — 8xy=5
4
2x
() 2453
10 Lf OXVeoy, yea te verify that:
(a) 2
Oz Ev)
Oz O?z O2z - (5)
Ox- Oy dy-Ox at) Ox2 dy?
fe
Ox2 aye’ rn
dx-dy' ee ind the val of aes
tie valu and b.
— Pye,
Pz
1S if z= xy/(x? + y")’, verify that Ox? Tay
d
16 If sin*x —5Ssinxcosy +tany =0, find osby using partial differentiation.
d
17 ‘Find . by partial differentiation, when x tany = ysin x.
OV OV PV @V
(a) x5 + Vay =0 (bd) xb + aye =0
OZ 1 OZ OZ
(b) sap = ae
miie {" = vel.
56
Polar coordinates 57
Introduction
We already know that there are two main ways in which the position of a
point in a plane can be represented:
(a) by Cartesian coordinates, i.e. (x,y)
(b) by polar coordinates, i.e. (r, 4).
The relationship between the two systems can be seen from a diagram:
Y
Dest For instance, x and y can be ex-
'g
ay (xc, y) pressed in terms of r and 6.
pa & 6)
LST ee y= sire
Or, working in the reverse direction, the coordinates r and @ can be found if we
know the values of x and y.
This is just by way of review. You should already be familiar with polar
coordinates from past studies. In this Program, we are going to focus attention
more to the polar coordinates system and its applications.
First of all, some easy examples to warm up.
Example 1
Express in polar coordinates the position (—S, 2).
Important hint: always draw a diagram; it will enable you to see which
quadrant you are dealing with and prevent your making an initial slip.
(5.385, 158°12’)
Because
(ayer ae 252 = A 0
- r= 729 = 5.385.
2
(b) tanE=2=0.4
»)
Pe B= 21°48
SO =ASS TZ:
Example 2
Express (4, —3) in polar coordinates. Draw a sketch and you cannot go wrong!
When you are ready, move to Frame 5
5) 323°3"
Here it is:
ry Seve
(4,=3) = (5,323°8)
Example 3
Express in polar coordinates (—2, —3).
|3.606, 236°19'
(a) r? = 2? 4 3*=449=13
r=¥V 13'= 3.606
3
(b) tan B= >= LS le £= 56719"
ten 23619"
(=2, 3) = (3.606, 236°19))
Polar coordinates
59
Of course conversion in the opposite direction is just a matter of evaluating
xX =rcos@ and y = rsin#. Here is an example.
Example 4
Express (5, 124°) in Cartesian coordinates.
Do that, and then move on to Frame 7
|(-2.796, 4.145)|
Working:
(ajex= sicosl24° = —=5:c0s 56"
= —$(0.5592) = —2.7960
(b) y=Ssin 124° = 5sin 56°
= 5(0.8290) = 4.1450
Polar curves
In Cartesian coordinates, the equation of a curve is given as the general
relationship between x and jy, i.e. y = f(x).
Similarly, in the polar coordinate system, the equation of a curve is given in
the form r=f(6#). We can then take spot values for 6, calculate the
corresponding values of r, plot r against #, and join the points up with a
smooth curve to obtain the graph of r = f(6).
Example 1
To plot the polar graph of r = 2sin@ between 6 = O and 6 = 2r.
We take values of 6 at convenient intervals and build up a table of values
giving the corresponding values of r:
0 30 60 90 1207 “130
sin 6 on SU.OO 1 0.866 0.5
re 250g Ne) 1732 2 ies 1.0
(a) We choose a linear scale for r and indicate it along the initial line.
(b) The value of r is then laid off along each direction in turn, points plotted,
and finally joined up with a smooth curve. The resulting graph is as shown
above.
Note: When we are dealing with the 210° direction, the value of r is negative
(—1) and this distance is therefore laid off in the reverse direction which once
again brings us to the point A. So for values of 6 between 6= 180° and
§ = 360°, r is negative and the first circle is retraced exactly. The graph,
therefore, looks like one circle, but consists, in fact, of two circles, one on top
of the other.
Now here is an example for you to do.
Example 2
This time, r is always positive and so there are, in fact, two distinct loops.
So on to Frame 12
i ce ead oeots a A a
62 Vector Analysis
3) 7 =acos¢
§ r=asin26 6 r=asin36
150°s,
=-—_-— —
WH P= OCR LE
8 1r=acos36
120°
N
"Qu
I /
270° 240°
9 r=a(1+cosé6) 10 r=a(1+2cos@)
\
'
Cant
(Py,
a
|
y
!
11 12 =a’cos26 12 r=a0
wee ' as !
Xx
X
1 =
'
XN ' 4
No =a
1
4 r NS t
¢
7
i]
‘x
~
!
7 t N '
7 xX
7 i ‘
Sketch these 12 standard curves in your record book. They are quite common
in use and worth remembering.
'
i]
'
So sketch the graphs of the following. Do not compile tables of values.
(a) r=2+2cos0 (c) r=1+2cosé
(b) r=5+3cos0 (d) r=2+cos0
64 Vector Analysis
2 5
(cj) r— T-.cos¢ (a < b) (d) r=2+cosé (a > b)
al
a 3
2
If you have slipped up with any of them, it would be worthwhile to plot a few
points to confirm how the curve goes.
On to Frame 16
To find the area of the plane figure bounded by the polar curve r = f(@) and the
radius vectors at @ = 0; and 0 = @>.
i of ea enes
ee) 50
6A dA sin 60
If 660, —=5
sO —
qa’ ér —. O, 7
Next frame
Polar coordinates 65
ue otk i4
ae irl =o
2
A=| ~1°d0
AT
Example1
To find the area enclosed by the curve r= Ssiné@ and the radius vectors at
6=0 and @=7/3.
a pare
=| at 7]
1
—(1 —cos 26) dé
aoe.
— ee Seer Finish it off.
Because
25 (nm sin2n/3
eres: 2
Sa
B(g-@) -3
2 [nxn
pM)
V3
Wort el:
Example 2
Find the area enclosed by the curve r = 1 + cos @ and the radius vectors at 6 = 0
and @= 7/2.
First of all, what does the curve look like?
jE
66 Vector Analysis
+
i!
cow
‘|
[4 3m
ee
rean
Because
n/2
A=3[- rao = 5 (1 + 2cos 6 + cos” 6)d0
0)
26 1/2
=; |@+2sine+5 i I,
_-
D = (ae ye)
5 eS 3 +1=2.178
ie
1s
Example 3
Find the total area enclosed by the curve r = 2. cos 38. Notice that no limits are
given, so we had better sketch the curve to see what is implied.
This was in fact one of the standard polar curves that we listed earlier in this
Program. Do you remember how it goes? If not, refer to your notes: it should
be there.
Then on to Frame 22
SS
Polar coordinates 67
Since we are dealing with r = 2 cos 38, r will become zero when cos 36 = 0, i.e.
when 30 = 7/2, i.e. when 0 = 77/6.
We see that the figure consists of 3 equal loops, so that the total area, A, is
given by:
A = 3 (area of one loop)
= 6 (area between 6 = 0 and @ = 7/6)
om/6 4 Y 7/6
7 units”
Because
7/6
Example 4
Find the area enclosed by one loop of the curve r = asin 20.
First sketch the graph.
r=asin26@
So the integral denoting the area of the loop in the first quadrant will be
68 Vector Analysis
A = ra*/8 units”
Now on to Frame 27
To find the volume generated when the plane figure bounded by r =f (@) and the
radius vectors at 6 = 0, and 0 = 42, rotates about the initial line.
fo) D x
If we regard the elementary sector OPQ as approximately equal to the AOPQ,
} Poe DHE
then the centroid C is distance — from O.
i) 2 es 2, ey —
=9tlr + or) sin 60. 2m Srsiné = > ar'(r + 6r) sin 60 - sin
Ne ee sinéé
9 3 + Or) 7 -sin@
dv
Then when 60
> 0,—=............
; do
—_—_—
esd
Polar coordinates
69
and therefore V =
Oo p} .
V= | = xr sin
6dé
ihe
Example 1
Find the volume of the solid formed when the plane figure bounded by
r=2sin@ and the radius vectors at 6=0O and 6=7/2, rotates about the
initial line.
om /2
Well now, V= | amr sin
6dé
0
m/2 x m/2 16
=| 3/7: (2sind)* -sinede = | —rsin* 6d0
Jo 3 0 3
Since the limits are between O and 7/2, we can use Wallis’s formula for this.
(Remember?)
V =7 units?
Because
n/2
V= | sin* 6 d6
Sedo
Example 2
Find the volume of the solid formed when the plane figure bounded by
r = 2acosé@ and the radius vectors at 9=0 and 6=7/2, rotates about the
initial line.
Do that one entirely on your own.
When you have finished it, move on to the next frame
——————————————
LL
70 Vector Analysis
Because
m/29
v=| 3 mr sin ddd and r = 2acosé@
1/2
=| oT en costae de
o. 3
3 pr/2
ee | cos? @(— sin @) dé
0
7, Lona cos ¢ vez l6na? [ 1
i“ 3 Aga 3 4
3
oe ee
3
So far then, we have had:
i)
(a) A= | 77d
0, 2 Check that you have noted these results in
292 your record book.
(Dave | =ar sin
édé
0,
r=f(@)
Q
&r
AS 8s
r38
Pp
: y 5.8)
With the usual tieure'ds-¢s 7 = ob.
© 9 y) 2 6s2 9 6r2
on meee roi ee
60 662
2 2 ie
if 80 0,(55) = 1+ cs 63} is r2 4 ae
dé dé dé dé
Polar coordinates
71
Example 1
Find the length of arc of the spiral r = ae*’ from 6 = 0 to 6 = 2r.
dr
Now, r=ae*® |. —=
7 3ae?
ae
. we"
SS dr\? 2 60 2 60 = 10ge
aa = ge. + 9g-e"
2 60
20 dr 2 27
3 =| r2 4 (Si) do = | V10
- ae*’ da
0
Because
2m / 2 /
| V10-a-e"°d@ = = | ae — te = 1}
) 0
As you can see, the method is very much the same every time. It is merely a
question of substituting in the standard result, and, as usual, a knowledge of
the shape of the polar curves is a very great help.
Here is our last result again:
8 dr\2
s=[
= n+24{(—]|
(35 dé
Example 2
Find the length of the cardioid r = a(1 + cos@) between @ = 0 and @= z.
Finish it completely, and then check with the next frame
72 Vector Analysis
2
+ (5) = a*{1+2cos6
+ cos? 6+ sin’ 6}
= a’{2
+ 2cos 6} = 2a*(1 + cos8)
0
Now cosé@ can be rewritten as (2cos” 5 _ 1)
7 6 T
e | 2acos5 dé
= 2a}2sin5|
0 0
= 4a(1 — 0) = 4a units
Next frame
Let us pause for a moment and think back. So far we have established three
useful results relating to polar curves. Without looking back in this Program,
or at your notes, complete the following.
If r= f(0) (Dh ee
: _
From our previous work, we know that =~ r? +
ér\2
SO that —© & 2arsin 64/12 +7
} dry
And now, if 66 — 0, — = 2arsin64/712-+ |—
0
A dr
This looks a little more involved, but the method of attack is much the same.
An example will show.
Example 1
Find the surface area generated when the arc of the curve r= S(1 + cos@)
between 6 = 0 and 6 = 7, rotates completely about the initial line.
dr ;
Now, f= 3(1-+1cos@) -.’. ha —5 sind
74 Vector Analysis
50(1 + cos @)
Because
2
To dr\* = 25(14+ 2cos6 + cos“2 6 + sin*
ae 9 8)
= 25(2
+ 2.cos 8)
= 50(1
+ .cos.@)
We would like to express this as a square, since we have to take its root, so we
now write cos @ in terms of its half angle.
2}
6
Dimas (ss) = so(1 + 2.cos* 5 — 1)
ae
0
Died
=sLOOEOS 5
dr\~ 7]
r 2 + (3) - 10 cos eh5
What do we get?
i 6
S$ = 4007 |,cos? sins dé
Because
d 0 0 6g 0
= 1007] | 2cos?—:
cos 5 2sin—cos=-cos—
sill, COSs cos, dé
6
= 4007 |,cos4 A sins dé
Polar coordinates 75
0 inte
Now the derivative of COs; 1S hee ka
ne
LA eS 8007 |cost 5 = reez dé
=e a ee Finish it off.
Because
Tv @ sin —
5 6 .
II |
COs" 5) — —8007
OOj=)=) pS) = —
¢)
S§ = 160z units?
And finally, here is one for you to do.
Example 2
Find the area of the surface generated when the arc of the curve r= ae’
between @ = 0 and 6 = 7/2, rotates about the initial line.
Finish it completely and then check with the next frame
2/2
§ =——=- na’(2e" + 1)
Because we have:
m/2 dr 2
s=|. 2nr sin 0 r+ (5) dé
a, dr 6
r=ae —=ae|.
dé
2
“ley? i (5) — gre + a2e = 2a7e8
76 Vector Analysis
1/2
S= | 2nae sin 6- V2ae’ dé
0
m/2
= 2\/2ra* | e7" sind dé
0
Let p= |e?singaa = 2" cos) +2 |cosa” a8
028
1= {2sind — cosa}
29 m/2
At 5 =2v2-n-a*|& fasing — cosa |
0
- OD eae
{e"(2-0) =41(0— »}
* 5
a 24/2 mae
S (2e" + 1) units”
Sars
We are almost at the end, but before we finish the Program, let us collect our
results together.
So move to Frame 47
Review summary
ES Polar curves — applications.
A i
1 Area A= | —1* dé
0, 2
A 2
2 Volume V= 5 ar sin 6 dé
. 6;
Oo dr\?
3 Length of arc s= | r+ (—) dé
4 dé
02 | 2
4 Surface of revolution S= 2ar sin 64/12 + (55) dé
d 6,
(4 Can You?
Checklist 3
Check this list before and after you try the end of Program test.
$ Test exercise 3
All the questions are quite straightforward: there are no tricks. But take your time
and work carefully.
1 Calculate the area enclosed by the curve ré? = 4 and the radius vectors at
06=2/2 and 0=7n.
2 Sketch the polar curves:
(a) r=2sin0 (b) r= 5 cos’6
(c) r=sin26 (d) r=1+cosé
(e) r=1+3cosé (f) r=3+cosé
78 Vector Analysis
3 The plane figure bounded by the curve r= 2 + cos @ and the radius vectors at
.
@=0 and @=7, rotates about the initial line through a complete revolution
Determine the volume of the solid generated.
]
Find the length of the polar curve r = 4 sin’ 5 between 6 = 0 and 6=7.
Find the area of the surface generated when the arc of the curve r = a(1 — cos @)
between 6 = 0 and @ = 7, rotates about the initial line.
That completes the work on polar curves. You are now ready for the next Program
Show that sin* 6 = 5= 5005 20F 5005 46. Hence find the area bounded by
the curve r= 4sin? 6 and the radius vectors at 6 = 0 and 6=7.
6
Find the area of the plane figure enclosed by the curve r= asec? (5) and
the radius vectors at @=0 and 6=7/2.
Determine the area bounded by the curve r = 2sin@+3cosé@ and the radius
vectors at 6 = 0 and = 7/2.
2
Find the area enclosed by the curve r = —————— and the radius vectors at
1+ cos 20
6=0 and 0=7/4.
Plot the graph of r = 1 + 2cos@ at intervals of 30 and show that it consists of a
small loop within a larger loop. The area between the two loops is rotated about
the initial line through two right-angles. Find the volume generated.
Find the volume generated when the plane figure enclosed by the curve
: 0
r = 2asin” (5)between @ = O and 6 = 7, rotates around the initial line.
The plane figure bounded by the cardioid r = 2a(1+cos6@) and the parabola
r(1+cos@) = 2a, rotates around the initial line. Show that the volume
generated is 187a°.
Find the length of the arc of the curve r=acos® (5) between 6=0
and @ = 37.
10 Find the length of the arc of the curve r = 3sin@ + 4cos@ between 6 = 0 and
Cte) ar
Polar coordinates
ae
<= 79
11 Find the length of the spiral r = a0 between
6 = 0 and 6 = 2r.
12 Sketch the curve r = asin3 (5) and calculate
its total length.
Learning outcomes
When you have completed this Program you will be able to:
e Determine the area of a rectangle using a double integral
e Evaluate double integrals over general areas
e Evaluate triple integrals over general volumes
e Apply double integrals to find areas and second moments of area
e Apply triple integrals to find volumes
80
Multiple integrals 81
If we add together all the elements of area, like 6a, to form the vertical strip
PQ, then 6A, the area of the strip, can be expressed as 6A =............
y=m
6A = i? dy - 6x
y=k
0] r s x
30 |
To evaluate this expression, we start from the inside and work outwards.
0 foto OT
jag
xX=r
X=S
|
y=m
y=k
=
:
=)
pnt
y=m
V=K
= | (m — k)dx
Because
xX=S
A= lm= kx
x=r
A =(m—k)-(s—r)
which we know is correct, for it is merely A = length x breadth.
That may seem a tedious way to find the area of a rectangle, but we have done
it to introduce the method we are going to use.
First we define an element of area dy - 6x.
Then we sum in the y-direction to obtain the area of a............
Finally, we sum the result in the x-direction to obtain the area of
Te ata Pie
y=m ve
Y=MN ¢xX=S
y=k =F
Aj II QuBs Qu=
a
———— Sear
Se
Because
y=mn 5 m m
So the order in which we carry out our two summations appears not to matter.
Remember:
(a) We work from the inside integral.
(b) We integrate with respect to x when the limits are values of x.
(c) We integrate with respect to y when the limits are values of y.
Move to the next frame
OT
84 Vector Analysis
Double integrals
v2 x2
The expression | | f(x,y) dxdy is called a double integral (for obvious
v1 U1
reasons!) and indicates that:
Example 1
2 (4
Evaluate l= ||(x + 2y) dx dy
wJ2
So (x + 2y) is first integrated with respect to x between x = 2 and x = 4, with y
regarded as constant for the time being.
2 cae S “4
| | (x + 2y) dx dy
er A
Finish it off
Because
Example 2
2 3
Evaluate [= |x*ydx dy
J1 JO
3
Do this one on your own. Remember to start with |x’y dx with y constant.
0
Finish the double integral completely and then move on to Frame 12
Multiple integrals 85
= 18 — 4.5. 13:5
Example 3
2 pu
Evaluate = | |(3 + sin @) dé dr
1 0
Here it is:
2 pt
=| |(3 + sin 6) dé dr
1 Jo
2 1
=| [30— cost dr
1 0
2
= {(3x+ 1) -(-1) Sar
Triple integrals
Sometimes we have to deal with expressions such as
but the rules are as before. Start with the innermost integral and work
outwards.
eS Se Se a Se SS SS
All symbols are regarded as constant for the time being, except the one
variable with respect to which stage of integration is taking place. So try this
one on your own straight away.
Example 1
3 pl (2
Evaluate r=| | |(x + 2y —z) dx dy dz
=i
Because
1
ine ae
—_ie)= | st
\
"| ios) iN}
= [30°-4] — 192-—4=188
0)
Applications
Example 1
4
Find the area bounded by y = = the x-axis and the ordinate at x =5.
Area of element = dy - 6x
ve=yal
.. Area of strip S~ éy- 6x
y=0
Because
Example 2
x
Find the area under the curve y = 4sin; between x =5 and x=, by the
double integral method.
Steps as before:
Area of element = 6y - 6x
Area of vertical strip
yey
as dy - 6x
y=0
Aw ¥ oral
x=r/3 | y=0
A= | | a 0b eee
ae Wee
1/3 JO
Pom):
Complete it, remembering that y; = 4sin 3
A = 4V3 unit?
Because you get:
T V) oT val pill
A=| | dy dx = | y dx=| y; dx
1/3 JO n/3 0 T/ 3
7 i ' x17
= fas4 sins ax = |-8cos 5 #
Example 3
Find the area enclosed by the curves
2
yi? = 9x and y2 = >
As usual:
Area of element = dy - 6x
‘. Area of strip PQ
VA
Soy: x
VV 2
9 V1
If by — O and 6x — 0, A= | [dydx
Here it is.
9 Pi 9 YM
A=| | dydx = |y dx
0 Jy2 0 iy
9
= |(vi — Y2)dx
9 1 xX
)
= {32 - a}dx
a7 \5
= 54— 27 = 27 unit?
Now for a different one. So move on to the next frame
Multiple integrals 91
Double integrals can conveniently be used for finding other values besides
areas.
Example 4
Find the second moment of area of a rectangle 6 cm x 4 cm about an axis
through one corner perpendicular to the plane of the figure.
z
I = 416 cm*
Because
6 74 6 y? 4 6 64
1=| |(2 +P) dydx= | ima dx =| {axe Shax
0 JO 0 a] (6) 0 3
3 6
Example 5
Find the second moment of area of a rectangle 5 cm x 3 cm about one 5 cm
side as axis.
Complete it and then on to Frame 26
EEE
92 Vector Analysis
Area of element= 6a = dy - 6x
Second moment of area of 6a about OX
=6a-y”
= y" -dy-6x
y=3
x=5 y=3
Second moment of whole figure ~ y~ - dy - &x
x=0 y=0
If 5y — O and 6x — 0:
5 73
i |y~ dy dx
0 Jo
5137? 5 5
4 | | dx =| hake = ox
(0) 3 0) 0 0
On to Frame 27
Review exercise
Now a short review exercise. Finish both integrals, before moving on to the
next frame. Here they are.
Evaluate the following:
(a) if
i(y? — xy) dy dx
a2
(b) || G+) ayau.
OJ1
Le-3)-G-dje
0J1 0 1
[{ee)-(e-9)}
Now on to Frame 29
Alternative notation
Example 1
2 1/2
f=] dx | 5 cos @ dé
0 0
2 n/2
=| dx ssino]
0 0)
2 2 2
=] a5] =| Silex 5s]
0 0 0
= iC
Example 2
6 m/2
Evaluate I = |dy | 4sin 3x dx.
3 0
Here it is:
6 a |2,
| dy | 4 sin 3x dx
3 0
Example 3
iax (a— x) dy
0 10)
Example 4
idy im yyy
y
When you have finished them both, move on to the next frame
Multiple integrals 95
ew:
Example3 I= |dx (x — x”) dy
- fo
0
1
Example 4 |—
1
2 vs PA
9 Pe
al
oeed | ee oY
Fieie i LF
tel & 6. 6
Next frame
4 vx
(b) |dx| (2y— 5x)dy
1 0
When you have completed both of them, move on to Frame 34
i
96 Vector Analysis
Working:
4 p2y
(a) | | (2x + 3y) dx dy
0 Jy
4 X—Ly
= | 2 + 339] dy
0) xXx=y
=|0meaVay
6y3]*
- cae
4
= ay" 7128
0
4 vx
(b) | dx| (2y— 5x) dy
ie e)o
4 yavx
= |dx [Y= sx]
1 y=0
b
Lf
6=6>
Total area = SS (all such thin sectors)
A=,
0=n T=
Ax) > ror: 60
d=0 r=0
Le ial
A=| | r dr dé
0 Jo
2.I5 ue ay 2Piacos @)
Because
If we now sum the columns between y = y; and y; = y2, we obtain the volume
of the slice:
Y=V2 2=2)
OV; = iS S| 6x: dy - 6z
y=y1 z=0
Then, summing all slices between x = x; and x = X2, we have the total volume:
X=X2 VV 2 27h
Vrs) meen Sy BE
X=X1 Y=)1 Z=Zo
The result this time is a triple integral, but the development is very much the
same as in our previous examples.
Let us see this in operation in the following examples.
Next frame
100 Vector Analysis
Example 1
A solid is enclosed by the plane z = 0, the planes x = 1,x= 4, y= 2,y =Sand
the surface z= x+y. Find the volume of the solid.
First of all, what does the figure look like? The plane z = 0 is the x-y plane
and the plane x = 1 is positioned thus:
This is just preparing the problem, so that we can see how to develop
the
integral.
Volume of column ~* 6x - éy 6
z=0
y=5 Z=x+y
Volume of slice + 6x} éy 6Z
y=2 z=0
x=4 y=5 Z=X+y
=a ~
Volume of total solid = SY bx. dy 6Z
Kel y=2 z=0
And this you can now finish off without any trouble. (With this form of
notation, start at the right-hand end. Remember?)
4 5 xX+Y 4 5
dx| dy| dz =| dx| dy(x+y)
1 2 Jo 1 2
4 95 “4 y °
dx| (x+y) dy = |ax|ay +5]
1 1 2 Jo
4 Dik 4
dx|5x +7 — 2x—2| | (3x43 )an
1 2
{132- 24 — 54 unit?
Example 2
of the solid bounded by the planes z=0, x=1, X= 2,
Find the volume
— —1, y=1 and the surface z=x?+y’.
of
In the light of the previous example, can you conjure up a mental picture
what this solid looks like? As before it will give rise to a triple integral.
2 1 x2 + y2
Wee ax| dy | dz
ail =il 0
Evaluate this and so find V. SaaS oe
Because we have:
2 1 x? +? 2 1
v=| dx | ay | az =| dx| dy(x? + y”)
Hie #0 eet
= eles]
= [{(@+3)- Cra}
Next frame
Review summary
Introductory work on double and triple integrals was covered in detail by this
text’s authors in Programme 23 of Engineering Mathematics (Fifth Edition).
Another look at the main points is well worth while.
You will no doubt recognize the following.
1 Double integrals
v2 x2
| | f(x, y) dx dy
|
Multiple integrals 103
ifdy fef(x, y) dx
v1 x1
[a aeForen?
Y1 | x) |
ee |
——--—--—--—--—--F
Meee
ei ei ee
2. Triple integrals
Triple integrals follow the same procedure.
Z2 p¥2 pXr2
| | | f(x, y, z) dx dy dz is evaluated in the order
Z) YV¥) YX
3. Applications
(a) Areas of plane figures
ob v2
If 6x + 0 and dy +0, A II | dy dx
a
104 Vector Analysis
(b) Areas of plane figures bounded by a polar curve r = f(@) and radius vectors
at 0 = 0) and 6 = 02
ro
A pt=f (A)
> If or Oand 8-0, A=| | rdrdé
6, Jo
(c) Volume of solids
Volume of element 6V = 6x dy 6z
z=f (x, y)
Volume of column = S > dx6éyéz
z=0
y=y2 (Z=f(x, y)
Volume of slice ~ > dx by |
yan (220
.. Total volume V = sum of all such slices
ines ~ oe : : 6x dy 6Z
Multiple integrals 105
(4 Can You?
Checklist 4
Check this list before and after you try the end of Program test.
1 Evaluate:
2 Determine:
V3+2 pr/3
(a) | (2cos 6 — 3 sin 36) dé dr
0 0
(b) ifixy(z + 2) dx dy dz
2J1 Jo
i! Pe Xx
a pb pc
4 Evaluate |||(x* + y)dx dy dz
m pu/2 pr
5 pvaluate | x’ sind
dx dé dé
JoJo Jo
6 Find the area bounded by the curve y = x? and the line y = x +2.
7 Find the area of the polar figure enclosed by the circle r = 2 and the
cardioid
r=2(1+cos@).
ae ae be
8 bvaluate |dx dy |xy"zdz
eer ers
Multiple integrals 107
Qoaihers
9 Evaluate |dx |(x? + y”) dy
0 1
1 m/4
10 Evaluate |dr| rcos? 6dé
0 Jo
11 Determine the area bounded by the curves x = y? and x = 2y — yr)
12 Express as a double integral, the area contained by one loop of the curve
r = 2cos 36, and evaluate the integral.
m/2 ptan-1(2) 4
13 Evaluate | | |X sin y dx dy dz
0 1/4 0
m p4cosz py/16—y2
14 +Evaluate || | y dx dy dz
0 10) 10)
1S__A plane figure is bounded by the polar curve r = a(1 + cos @) between @ = 0 and
@ = 7, and the initial line OA. Express as a double integral the first moment of
area of the figure about OA, and evaluate the integral. If the area of the figure is
known to be =mn unit*, find the distance (h) of the centroid of the figure
from OA.
16 Using double integrals, find (a) the area and (b) the second moment about OX
of the plane figure bounded by the x-axis and that part of the ellipse
xy
a + ae 1 which lies above OX. Find also the position of the centroid.
17 The base of a solid is the plane figure in the x-y plane bounded by x = 0, x = 2,
y =xand y = x* + 1. The sides are vertical and the top is the surface z = x? + y?.
Calculate the volume of the solid so formed.
18 A solid consists of vertical sides standing on the plane figure enclosed by x = 0,
x=b, Y=aand y=c. The top is the surface z=xy. Find the volume of the
solid so defined.
19 Show that the area outside the circle r = a and inside the circle r = 2acos@ is
given by
a/3 p2acos@
AS 2 | | rdr dé
0 a
Differentials and
line integrals
Learning outcomes
When you have completed this Program you will be able to:
Understand the role of the differential of a function of two or more real
variables
Determine exact differentials in two real variables and their integrals
Evaluate the area enclosed by a closed curve by contour integration
Evaluate line integrals and appreciate their properties
Evaluate line integrals around closed curves within a simply connected
region
Link line integrals to integrals along the x-axis
Link line integrals to integrals along a contour given in parametric form
Discuss the dependence of a line integral between two points on the path
of integration
Determine exact differentials in three real variables and their integrals
Demonstrate the validity and use of Green’s theorem
108
Differentials and line integrals 109
Differentials
y y=f(x)
Xo Xo+Ox xX
y=f(x)
Note that the differentials dy and dx are finite quantities — not necessarily zero
— and can therefore exist alone.
Note too that dx = 6x.
110 Vector Analysis
That was easy enough. Let us now consider a function of two independent
vatiables, Z= f(xy).
If z= f(x, y) then z+ 6z = f(x + 6x,y + dy)
". 62 =f (x + 6x, y + by) — f(x, y)
Expanding 6z in terms of 6x and dy, gives
6z = Adx + Béy + higher powers of 6x and 6éy,
where A and B are functions of x and y.
If y remains constant, i.e. 5y = 0, then
. heeOz Oz
OZ
Z x 6X + By oy
Exercise
1 dz = 2(xdx + ydy)
2 dz=x"(3sin 2y dx + 2x cos 2y dy)
3 dz=e’{2dx
+ (6x — 3)dy}
4 dz=2(xdx
+ 2ydy + 3wdw)
5 dz=x’y(3ywdx + 2xwdy + xy dw)
Now move on
————
aaa
112 Vector Analysis
Exact differential
We have just established that if z = f(x, y)
OZ OZ
=—dx+—d
Fe a cae
We now work in reverse.
Any expression dz = P dx + Qdy, where P and Q are functions of x and y, is an
exact differential if it can be integrated to determine Z.
OZ OZ
OP 0-7 00 Oz 7, O*z
N —_ = — = t ——— = 3
ou Oy OyOx oy OX OX OV aE ee gue ee OyOx OxOdy
Therefore, for dz to be an exact differential 5= = and this is the test
we apply.
Example 1
dz = (3x? + 4y”) dx + 8xy dy.
If we compare the right-hand side with Pdx + Qdy, then
Pa x- 3s
+ 4y ae
ay 8y
JQ
Q = 8xy
= 8
ape)
SS
Oe
cio ., dz is an exact differential
Oy Ox
Example 2
dz = (1+ 8xy) dx + 5x? dy.
Frome his Wel tind paeeene eee
oe Ome ay 8x
0
Q = 5x? OY
OX
Differentials and line integrals 113
Exercise
We have just tested whether certain expressions are, in fact, exact differentials
— and we said previously that, by definition, an exact differential can be
integrated. But how exactly do we go about it? The following examples will
show.
Example 1
dz = (2xy + 6x) dx + (x* + 2y3) dy.
4
c= x2y 4) 5 + F(x) where F(x) is an arbitrary function of x only
Example 2
Integrate dz = (8e* + 2xy”) dx + (4cos 4y + 2xy) dy.
Argue through the working in just the same way, from which we obtain
Exercise
Integrate the following exact differentials to obtain the function z.
1 dz = (6x? + 8xy3) dx + (12x7y? + 12y3) dy
2 dz = (3x? + 2xy + y”) dx + (x? + 2xy + 3y?) dy
3 dz=2(y + 1)e* dx + (e* — 2y) dy
4 dz= (3y? cos 3x — 3sin 3x) dx + (2y sin 3x + 4) dy
5S dz = (sinhy + y sinh x)dx + (x cosh y + cosh x) dy
Finish all five before checking with the next frame.
a
a
Differentials and line integrals 115
zZ=x°+x*y+xy4y3
Z=e*(l+y)-y"
z=y” sin 3x + cos 3x + 4y
=
NO
WwW
nm z=xsinhy + ycoshx
In the last one, of course, we find that the two expressions for z agree without
any further addition of f(y) or F(x).
We shall be meeting exact differentials again later on, but for the moment let us deal
with something different. On then to the next frame
x4 X2 x
For example
X2
| F(x)dx gives Az as before, but integrating from B to A along c2 with
x1
x2
.. The result A = A; — Az = |
xX}
A ~ pyds
x1 x2
y=f(x) B
Xx} X2 Xx
Example 1
Determine the area enclosed by the graphs of y = x° and y = 4x for x > 0.
First we need to know the points of intersection. These are
se =O) Sigel se
A = 4 square units
Because
A= —bydx = {0 Sollee
Another example.
=the)
Example 2
Find the area of the triangle with vertices (0, 0), (5, 3) and (2, 6).
The equation of
B (2, 6)
COUARISK Ge aera
BATS S.-,. tee ee
OBAS*. 2 heeee
3 A (5, 3)
© 2 5 i
The limits chosen must progress the integration round the boundar
y of the
figure in an anticlockwise manner. Finishing off the integration,
we have
ae Differentials and line integrals 119
A = 12 square units
Line integrals
This is normally written | F,ds where A and B are the end points of the
AB
curve, or as |F, ds where the curve c connecting A and B is defined.
JSC
where c is the curve y = f(x) between A (x1, v1) and B (x2, yz).
There is in fact an alternative form of the integral which is often useful,
so let us also consider that
a a eee
120 Vector Analysis
If F, has a component
P in the x-direction
Q in the y-direction
then the work done from K to L can
be stated as P6x+ Qéy.
t= [F ds = |(wax +Qdy)
Example 1
Evaluate |(x + 3y)dx from A (0, 1) to B (2, 5) along the curve y= 14 x”.
(e
|(Pax +Qay)
Cc
x :
- bene =16
10)
Example 2
P= 2x.
Because
3
r=| {(x2 4.x +2) dx— (x2 43x44) dx}
0
3 a
Here is another.
Example 3
Evaluate I = | {(x?
+ y 2y)dx+xydy}
ydy from O (0,0) to B (1,4) along g the
fe
curve y = 4x’.
[e294
Because
Example 4
i= OF and Tae, = 84
Notice therefore that integration along two distinct paths joining the same
two end points does not necessarily give the same results.
Let us pause here a moment and list the main properties of line integrals.
i.e. the sign of a line integral is reversed when the direction of the
integration along the path is reversed.
3 (a) For a path of integration parallel to the y-axis, i.e. x =k,
y = fi(x) from A to K
y = f2(x) from K to B
(x1, vy)
Example
Rey sO 0.
For any value of x, y = +v1— x’. Therefore, we divide c into two parts
particular casey. ey
1 10)
" j=| Par=| (x+ VI=x) dx+ |(x= VI —3) ax
elite JO Jl
1 =— 1
= | + Vine x4 VIR9) dx = 2 [Via
é JO
Because
1
1=2| V1—x2dx x=sin@ .. dx=cosédé
0
V1—x2=cosé
- onze"
2 (@)
— A
The closed curves involved in problems in this Program all relate to simply
connected regions, so no difficulties will arise.
126 Vector Analysis
y
Positive direction (anticlockwise) line integral
denoted by f
O x
y Cc
(1) Use y= f(x) for ALB (2) Use y = fo(x) for BMA
Unless specially required otherwise, we always proceed round the closed curve
iN aT ok ners
anticlockwise direction
Example 1
three sides of the triangle joining O (0, 0), A (1, 0) and B (0, 1).
First draw the diagram and mark in c1, co and c3, the proposed directions of
integration. Do just that.
Differentials and line integrals
127
Cy 1 x
ae Pe Ea— i| 1
T 1 — Jaee dx eiF 3 See
‘ ls == 3
Example 2
Ms yaad oo V4 =
y is thus not single-valued. Therefore use
y=vV4-—x2 for ALB between x=2 and
x = —2 and y= —v4 — x? for BMA between
x=—2 and — 2.
ed =| ve ars [ {V4
— x? }dx
2 =2
=y) 2
22 WS Reb? = -2| \/ 4 92d
2 =2
2
=-4| V4 —x2dx
0
To evaluate this integral, substitute x = 2 sin 6 and finish it off.
Because
= -8|
7/2
(1 + cos 20) dé = -8|04
sin 267°
7) >) ne awh
0 0
Now for one more
Example 3
Evaluate l= |{xy dx + (1+ y”) dy} where c is the boundary of the rectangle
c
Now evaluate J; for AB; Iz for BC; [3 for CD; I, for DA; and finally I.
Complete the working and then check with the next frame
(Cc) CD es ise 2 dy
=0
1 1
a b= | 2x dx = [| 8 kk =-8
3 3
Finally
T=h+bh+h+h
=0+4%-8-42=-8 wel eB
Remember that, unless we are directed otherwise, we always proceed
round the closed boundary in an anticlockwise manner.
On now to the next piece of work
LL
130 Vector Analysis
| F, ds =| {P dx + Qdy}
AB AB
where F; denoted the tangential force along the curve c at the sample point
K (x, y).
The same kind of integral can, of course, relate to any function f(x, y) which
is a function of the position of a point on the stated curve, so that
(p= | egalaeks
Example
Evaluate n= |(4x
+ 3xy)ds where c is the straight line joining O (0, 0)
Cc
COCA 2s
Se ae
Use dy
2
=
he
ati xX=1
=| (4x + Say)(VS)dx. But y = 2x
|, (4x + 3ay)ds 1
Differentials and line integrals
131
T=4VJ5
Because
“1 1
r= | (4x + 6x?)(V/5)dx = 2v5 | (2x + 3x7)dx =4V5
0 0
Try another.
The path length of the parabola defined by y = x2 betwen the values x = 0 and
X = 2 is given by the integral
t=| d= Sor ee oa to 3 dp
Because
2 dy =
1=| as=| 1+ (3) ar
ne be 01 dx
2
== | V¥1+2xdx
x=0
u=1 oi
5
a 1 2 A 2
at)
= 5 (1254?1/2 =1)
os 1
= 3.393 to 3 dp
Parametric equations
When x and y are expressed in parametric form, e.g. x = f(t), y= g(t), then
ds ae
faa dx\*
ae /dy\?
bee : = dx\*
eet /dy\’
mass t
dt (2) +(#) “ie /(2) +(2) g
and result (1) above becomes
; re dx 2 dy 2
Example
Because
D
.= (2) +(2) = Vcos?t+sin?t =1
t 2 7/4
ef 1=| Foy) ae +(H) at = | 4sint costdt
ty dt JO
1/4 m/4
= | sin 2t dt = -2| = Sf
0 0
Example
EE
Differentials and line integrals 133
The path c, is y = x? ay = xr
2 z
Ty =| {3x2x* dx + 2x3x?2x dx} = | (3x°s0/4x*)ydx
Jo 0
2
= >| = 128 “, I, = 128
leet28
1 2 x
From (0, 0) to (0,4) x=0 .. dx=0 .. ba=0
2
From (0, 4) to (2,4) y=4 .. dv=0 ~-. In =48 \. x?dx = 128
fale co
On to the next frame
In the example we have just worked through, we took three different paths
and in each case, the line integral produced the same result. It appears,
therefore, that in this case, the value of the integral is independent of the path
of integration taken.
Let us investigate
Cy Toy = Ic,
Oo
xX
y A
Co If we reverse the direction of c2, then
To, “a ley
Leo idee, = 0
Cy 2
O
x
Hence, if (Pdx + Q dy) is an exact differential, then the integration taken round a
closed curve is zero.
Example 1
fyy=73 FX) =0
Hence z = 3xy + y”
(3, 5) |
. 1=| (3y dx + (3x + 2y) dy} = | d(3xy
+y?)
(e (1, 2)
(3, 5)
= xy + a
(1, 2)
= (45 +25) —(6+4)
= 60
Example 2
Evaluate J = |Hee + ye*)dx + (e* + y) dy} between A (0, 1) and B (1, 2).
Cc
pa aw tye! i A ih Senios
az
6 ea aa Sy rd 7
oy
Continue the working and complete the evaluation.
When you have finished, check the result with the next frame
x ;
Be rey oe Hh)
y
Z=ye +5 + FX)
y 3
For these expressions to agree, Ay) = 2) F(x) =—
x3 y? (1, 2)
Then i Fre + a
° 2}, 1)
5
Example
Verify that dw = (3x?yz + 6x)dx + (x3z — 8y)dy + (x3y + 1)dz is an exact
differential and hence evaluate |dw from A.(1, 2, 4) to B (2, 1, 3).
EE tO OD
’ OZ OZ _ OW
Now to find w. cae Cay’ rece
= xyz +Z+8(X, y)
For these three expressions for z to agree
POV, Zee gata as tos: eR re ees We rshtnminssSEG UN GV) ele hats weescaatatess
138 Vector Analysis
' waxyz
+ 3x? —4y? +z
(2,1,3)
= yz Bye Ay +2
(1,2,4)
Because
(2, 1,3)
I= |x°yz+ 3x? — 4y? +2
(1, 2,4)
= (244 12—4+43)-—(8+3-16+4) = 36
If the first partial derivatives are continuous within the region and on the
boundary, then Green’s theorem states that
Example 1
oP
P=2x-
x—y =H
By 1
AEB) a0
= a[.fia —1)dxdy
=2| Jaray
x*4+9y7 =16 Le x p=
4
a=4, b=,
A= nab =~"
Sys
I =2A ==
To demonstrate the advantage of Green’s theorem, let us work through the
next example (a) by the method of line integrals, and (b) by applying Green’s
theorem.
140 Vector Analysis
Example 2
round the triangle with vertices at O (0, 0), A (1, 0) and B (1, 2).
O
is
e
itd
iant
i
acl
ailia
There are clearly three stages with cj, cz, cz. Work through the complete
evaluation to determine the value of J. It will be good revision.
When you have finished, check the result with the solution in the next frame
(ae@)ecrisy
=0 F dy2-0
1 1
head =| 2x dx = | ea ee |
0 0
(2N5CoMS Naa 1a xO)
v2 z
i =| (3 — 2y) dy = sy- | = hema ez
0 (0)
ne | fer 0Q ;
\,|(3 =) a
BIniSHr it Otte t= eees oe
Because
ioe -| fa = 3)dxdy
R
=2| |axay=24
R
= 2x the area of the triangle
=o ag i fZ
Application of Green’s theorem is not always the quickest method. It is
useful, however, to have both methods available. If you have not already done
so, make a note of Green’s theorem.
Example 3
Evaluate the line integral J = {{xydx + (2x — y) dy} round the region bounded
Sree es
ie<a ial
--[j@-ap]
x? eieh
~ Gn
Before we finally leave this section of the work, there is one more result to
note.
= {way yan
Je
A= | axay = xp (x dy — ydx)
Note this result in your record book. Then let us see an example
SS
eet
Differentials and line integrals
143
Example 1
Determine the area of the figure enclosed by y = 3x? and y = 6x.
y Points of intersection:
A (2, 12)
3x*=6x .. x=Oor2
Area A =i (x dy — ydx)
JC
2A = | (x dy — ydx) +| (xdy—ydx)=h+h
c; (along OA) c2 (along AO)
Because
cZ2isy=6x .. dy=6dx
0
a h= | (6xdx —6xdx)=0 .. h=0
22
5 taht
b =84+0=8. ~.. A=4. square units
Example 2
Determine the area bounded by the curves y = 2x3, y= x3 +1 and the axis
a Oi TOne es 0:
Complete the working and see if you agree with the working in the next frame
144 Vector Analysis
im
nae
pa
©
9p
= 14 Vaxv=al «— #21
1
Area A =54 (x dy — ydx)
She! =4 (x dy — y dx)
aa = G-N=2 b=}
(Cc) BOoy Ceis
x OF a dye)
y=0
2A =T=h+h+h=1+54+0=13
And that brings this Program to an end. We have covered some important
topics, so check down the Review summary and the Can You? checklist
that follow and revise any part of the text if necessary, before working through
the Test exercise. The Further problems provide an opportunity for
additional practice.
ns
Differentials and line integrals
145
=>|Review summary
1 Differentials dy and dx
(a)
OZ OZ OZ
f(x,y, w) Z Ax cri agress
r= | Fas =| (P dx + Qdy)
AB AB
integration
e
al
e
iii
(c) | (P dx + Qdy + Rdz) is zero when c is a closed curve.
(
9 Green’s theorem
Check this list before and after you try the end of Program test.
|dw where
€
Further problems 5
Show that I = |.{xy*w* dx + x*yw? dy + x°y?w dw} is independent of the path of
integration c and evaluate the investal from A (1,3, 2) to B (2, 4,1).
Determine whether dz = 3x*(x? + y”) dx + 2y(x3 + y*) dy is an exact differential.
If so, determine z and hence evaluate f.dz from A (1, 2) to B (2, 1).
xdy — ydx
Evaluate the line integral r= $4 } where c is the boundary
x2 +y244
of the segment formed by the arc of the circle x? + y2 = 4 and the chord
y=2-xforx>0.
Show that
I = {.{(3x? siny + 2 sin 2x + y?) dx + (x3 cosy + 3xy?) dy}
is independent of the path of integration and evaluate it from A (0,0) to
where c is the boundary of the trapezium with vertices A (0, 1), B (S, 1), C (3, 3)
and On)
10 Evaluate I = {_{(3x2y? + 2.cos 2x — 2xy) dx + (2x3y + 8y — x”) dy}
(a) along the curve y = x* — x from A (0, 0) to B (2, 2)
(b) round the boundary of the quadrilateral joining the points (1, 0), (3, 1),
(2, 3) and (0, 3)
150 Vector Analysis
|aw
where c is the straight line joining (0, 0, 1) to (1, 2, 3) for either region z > 0 or
Lee)
-__—
eeeeSS
OoSSSSSSSeeSSSFSFSFMMMSSMMmmmeFeFeFsFseF
Program6
Surface and an
volume integrals
Learning outcomes
When you have completed this Program you will be able to:
Evaluate double integrals and surface integrals
Relate three-dimensional Cartesian coordinates to cylindrical and sphe-
rical polar forms
Evaluate volume integrals in Cartesian coordinates and in cylindrical and
spherical polar coordinates
Use the Jacobian to convert integrals given in Cartesian coordinates into
general curvilinear coordinates in two and three dimensions
151
152 Vector Analysis
Double integrals
Let us start off with an example with which we are already familiar.
Example 1
A solid is enclosed by the planes z=0, y=1, y=2, X=0, X= 3 and the
surface z =x + y?. We have to determin e the volume of the solid so formed.
r
—
(7)
In the plane y= 1, z=x+1, ie. a straight line joining (0, 1, 1) and (3, 1, 4)
In the plane y = 2, z=x+4, ie. a straight line joining (0, 2, 4) and (3, 2, 7)
2
In the plane x = 0, z=y’, i.e. a parabola joining (0, 1, 1) and (0, 2, 4)
In the plane x = 3, z=3 +”, i.e. a parabola joining (3, 1, 4) and (3, 2, 7)
Consideration like this helps us to visualize the problem and the time
involved is well spent.
Now we can proceed.
The element of volume 6év = 6x 6y 6z
We could also have said that the element of area on the z= 0 plane
6a = by 6x
and that the volume of the column =év,. = zda = z 6x by
Then, since z= x+y”, this becomes dv¢ = (x + y”) 6x 6y
Summing in the usual way then gives
Y= |za
=| [xt y)dvay
R
where R is the region bounded in the x-y plane.
Now we insert the appropriate limits and complete the integration
Because
Leal = | [re y) da
= \s[ro y) dy dx
In this particular case, R is the region in the x-y plane bounded by x= 0, x =3,
y=1,y=2.
154 Vector Analysis
Example 2
A triangular thin plate has the dimensions shown and a variable density
o Where p= 1k wey:
We have to determine
(a) the mass of the plate
(b) the position of its center of
gravity G.
da = 6x dy
The mass 6m of the element is
then
6m = p 6x dy
Because we have
I
X=2 py=2x
xy) dy dx
[hoot
M= [odxay =| | (1+x+
x=0 Jy=0
Ps 27 V=2x
0 oo
v2
= |{2x + 2x? + 2x3}dx
3 472
Se cm eeap. ey
Surface and volume integrals 155
58m xy) 6x dy
= x(1+x+
O
pe
‘. Sum of first moments = | Jo+x? + xy) dx dy
R
Kad pyaly
Because sum of first moments = | | (x + x? + xy) dy dx
x=0 Jy=0
2 24,27 Y=2x
=| prety] dx
y=0
2
= ee 2x9 4 2x" hx
iL Nes 2
We found previously that M = 17, ie (175): 20
33
which gives X = le =l:508
(2) To find y we proceed in just the same way, this time taking moments
about OX. Work right through it on your own.
156 Vector Analysis
(vy + xy + xy”) dx dy
X=2 py=2x
| (vy + xy + xy”) dy dx
0 Jy=0
bi 2
yay 7]ey Reape
Ay
2, WA es MR
2 4
=| {2x2 +248 4 dx
0 3
_ es , a
eee Es
3
: Se
“. My = 305
So we finally have:
Note that this again referred to a plane figure in the x-y plane.
Surface integrals
When the area over which we integrate is not restricted to the x-y plane,
matters become rather more involved, but also more interesting.
If S is a two-sided surface in
space and R is its projection on
the x-y plane, then the equa-
tion of S is of the form
z=f(x, y) where f is a single-
valued function and continu-
ous throughout R.
Let 6A denote an element of R
and 6S the corresponding ele-
ment of area of S at the point
Pts, ¥,2) 07a:
Let also @(x, y,z) be a function of position on S (e.g. potential) and let 7 denote
the angle between the outward normal PN to the surface at P and the positive
Z-axis.
6
Then 6A = 6Scos7 i.e. 6S = Re = 6Asecy and
cos7+
ye 0(x,y,z)d6S is the total value of ¢(x,y,z) taken over the surface S.
As 6S — 0, this sum becomes the integral
il Q(x, y; z) dS
s
and, since 6S ~ 6Asecy, the result can be written
| |oo. y, Z)secy
dx dy (1<5)
R
Notice that cos = n-k, where k is the unit vector in the z-direction and n is
the unit normal to the surface at P.
With limits inserted for x and y, the integral seems straightforward, except for
the factor sec, which naturally varies over the surface S.
We can, in eerie
fact, show ere
that secy = ee
Ay eee
ay \\
(a) I= Rc y, 2) ds (1)
pee
or (b) 1=| [os y, ayfr+ (2) +(Z) dx dy (2)
where z= f(x, y)
158 Vector Analysis
Note that, when ¢(x,y,z) = 1, then I = |dS gives the area of the surface S.
S
Example 1
Find the area of the surface z= \/x2+y* over the region bounded by
Ney = Ne
2 2
So we now find = and = and determine i1+ (=) +(=) which is
Because
za (pyle 9 Senet
Z|
Bh tlyt) MP 2x = irae.
aye x
OZ lt 2 Oni? y
a =-(xX + “2y = ———
fs BH) 22
$a
—_——_ el cal att soioay =
Surface and volume integrals
159
Example 2
Find the area of the surface S$ of the paraboloid z = x” + y? cut off by the cone
Z=2,/xX2 + y2.
We can find the point of intersection A by considering the y-z plane, i.e. put
yee).
Coordinates of"A‘are*:......4;..
A (2, 4)
|thecircle x7 +y? =4 |
/ mix 2 2
OZ OZ
= — — | dxd
: [Jy2+ Ge) +) ae
For this we use the equation of the
surface S$. The information from the
projection R on the x-y plane will
later provide the limits of the two
stages of integration.
160 Vector Analysis
X=2 py=V4—x?
Leese 4| | [1 + 4x2 + 4y2 dy dx
x=0 Jy=0
Gale pr=2
| | V 1+ 4r2 rdrdé
0=0 Jr=0
er COSY, y=rsin@g
x*4+y2=r* dxdy=rdrdé
(refer to Frame 67)
6=27r pr=2
Finish it off.
Surface and volume integrals 161
Because
O=2n pr=2 2a i oh
Now on to Example 3.
Example 3 Vv
To determine the moment of inertia of a thin spherical shell of radius a about
a diameter as axis. The mass per unit area of shell is p.
Equation of sphere
P+y4+7 =a"
Mass of element = m = péS
I~ Sr? & Dp bSr?
cee =| pr? dS
S
~~ no az\*_
ee (az\*
Be:
=|,fo \) @ ats aay
Now determine the partial derivatives and simplify the integral as far as
possible in Cartesian coordinates.
3 i Fi oa tt eS
dx dy
H R af /a2 — x2 — y2
In this particular example, R is, of course, the region bounded by the circle
x2 +y? =a’ in the x-y plane.
Converting to polar coordinates
x=rcosd#; y=rsind; dxdy=rdrdé
the integral becomes I =.....--.+---
ee
162 Vector Analysis
es 0/016
i Vaz — r2
=| OVa7—1
er
a r
Because
a r3 a r2
=| | et
0 Va? — 1 r=0 Va2 — r2
|
Oa du
u=a? Ju 2
SS a ue 0 ibe
414. 3/2 W
2) |‘ ci 2; Fce |
7a
3
var. ¢3
=¢2
2a?
i)
Now, to complete I; we have
27 3
ie pa| calls)
0
Surface and volume integrals 163
4rpat
Because
In = pa |
2a
2 a 3 qo — 24 4 ag 2n _ 4na*p
4
0 3 3 0 3
Therefore, the moment of inertia for the complete spherical shell is
8rat
I, = 3 P
av 2Ma?
The total mass of the shell M=4na2p .. I
3
Now let us turn our attention towards volume integrals and in preparation
review systems of space coordinates.
y-z plane
(x=0)
x-y plane
x (z=0)
The three coordinate planes, x =0, y=0, z=O, divide the space into
eight sections called octants. The section containing x > 0, y = 0, z = O is
called the first octant.
P (x,
y, 2)
x=rsindcos@ +2?
r= x+y?
y=rsinésin@g § = arccos(z/T)
Z= COS) ¢@= arctan(y/Xx)
Because
=o)
) — 3) oe
Because
év = r60 or 6z
OV = Tonogioz
rsinOdo
It is important to make a note of these results, since they are required when we
change the variables in various types of integrals. We shall meet them again
before long, so be sure of them now.
168 Vector Analysis
Volume integrals
22= F(x, y)
v= | | | dz dy dx
X=X% Jy=y1 JZ=Zy
Example 1
Find the volume of the solid bounded by the planes z=0, x=0, y=O,
x* 4+? = 4 andz=6-=—xy
forx>0, y>0,z>0.
First sketch the figure, so that we can see what we are doing. Take your time
over it.
SE
Surface and volume integrals 169
dv = 6x by 6z
zZ=6—xy
Volume of column = S~ dx éy6z
z=0
V4—x2 (6—xy
Volume of slice= S~ is 6x by |
y=0
2
Total volume ~ 5~ 6x by 62
x=0 y=0" “z=0
If 6x — 0, dy > O, 6z — O, then
2 pV4—x? p6—xy
v= | | | dz dy dx
0 JO 0
| i=) |
Because
Gy 4a (4 x?)
2,
Then finally V = |{o(4= iad1/2 - 2143}
x
dx
‘ 2
Now we are faced with | — pi le dx. You may remember that this is a
x
standard form |Va? — x2 dx =5 {xva? — x2 + d@ arcsin a
2
If not, to evaluate | V4 — x2dx, put x = 2sin 6 and proceed from there.
0
Finish off the main integral, so that we have
Because we had
2 x3
Vis |{o(4Sy 19x + hax
10)
2 an2
= 3 [v4 — x2 + 4aresin | ae[35
0 8 0
= 3{4 arcsin 1 — 4 arcsin0} —4+42
= 32m}
= 2 = 67 =2
=~ 16.8
Surface and volume integrals 171
Alternative method
év =r ér 60 6z
xX=rcosé; y=rsiné
*. Z2=6—xy
= 6-—r’sinécosé
r2
= 6 —— 5 sin
i 20
2 n/2 p6—(r2/2)
sin 20
% v= | r dr dé dz
r=0 J6=0 Jz=0
m/2 p2 6—(r?/2)
sin 20
| | | dzrdrdé
6=0' Jr=0 Jz=0
Finish it
V = 6n — 2 (as before)
n/2 2 ea
vel | (6-5 sin20)rdrao
6=0 Jr=0
m/2 2 r3
= | | (<r——sin 28)
dr dé
6=0 Jr=0 2
n/2 4 r=2
= | [37—" sin 24 do
8 r=0
a /2
= | (12 — 2sin 26) dé
m/2
= 1264-'COS 26|
0
1 V = On = Z
In this case, the use of cylindrical coordinates facilitates the evaluation.
Let us consider another example.
EE
172 Vector Analysis
Example 2
sphere
To find the moment of inertia and radius of gyration of a thick hollow
of material
about a diameter as axis. Outer radius = a; inner radius = b; density
= Gs
It is convenient to deal with one-eighth of the sphere in the first octant.
1
-. Total mass of the solid M; = at
dv = 17 sin
06r 66 56
= mp? = m(rsind)
=cr* sin 6 6r60 66 r sin? 6
=cr'sin? 6 6r 60 6¢
.. Total second moment for the solid
2, Te a
EE
a,
Because
1/2 pr/2 rs a
i =| | er sin? 6d6 dé
b
m/2 pr/2 C
ae
:
wf Be
hg
2
aeA x8
ie Fob . 9 (qd bs
65a — Bb = = (=)
Example 3
Find the total mass of a solid sphere of radius a, enclosed by the surface
x2 + y2 +z? =a? and having variable density c where c = 1 +7r/|z| and ris the
distance of any point from the origin.
This is a case where spherical coordinates can clearly be used with
advantage.
My. = J ocenees
Write out the integral and insert the limits.
es
Surface and volume integrals 175
20 m/2 wa
= I, =. Ip
2n pr/2 pa
i= | | |r’ sin 6 drdé dd gives ............
0 Jo Jo
Do not work it out. You can doubtless recognize what the result would
represent.
Because
2m pm/2 a
i= | 5 Sin @cos 6 dé de
0 0
m/2
a sin? 6
de
Bisse olen 0
da 27
= — 1d¢d
a
5
| na°
TG 0)
0)
5
Ta
Loe
So now finish it off. For the complete sphere
176 Vector Analysis
Because
Example 4
Determine the volume of the solid bounded by the planes x = 0, y=0, z=x,
z=2 and y =4 — x? in the first quadrant.
First we sketch the diagram.
Cartesian
2,
Ve 65 cubic units
2 4—x* 92
Aa | dz dy dx
x=0 Jy x
C (x3, y3)
B
A | (Xz, Y2)
(1, V1) | |
| |
| |
M | |N
x X3 X2 i
Curvilinear coordinates
eo)
These two sets of curves will therefore cover the region R and form a network,
and to any point P (xo, yo) there will be a pair of curves u = Uo (constant) and
Vv = Vo (constant) that intersect at that point.
Up Ug + Ou
The u- and v-values relating to any particular point are known as its
curvilinear coordinates and x =f(u, v) and y=g(u, v) are the transformation
equations between the two systems.
In the Cartesian coordinates (x, y) system, the element of area 6A = éxdéy
and is the area bounded by the lines x=X9, x=X0+6xX, y=yo, and
y=Yo + oy.
In the new system of curvilinear coordinates (u, v) the element of area 6A; can
be taken as that of the figure P, Q, R, S, i.e. the area bounded by the curves
U= Uo, U = Up + OU, V = Vo and V = Vo + OV.
Since 6A; is small, PQRS may be regarded as a parallelogram
i.e. 6A; © 2 x area of triangle PQS
and this is where we make use of the result previously revised that the area of a
triangle ABC with vertices (x1, yi), (X2, y2), (x3, y3) can be expressed in
determinant form as
Surface and volume integrals
181
1
Area == Xp x2 X3
Vi “V2 V3
Before we can apply this, we must find the Cartesian coordinates of P, Qand
S in the diagram in the previous frame where we omit the subscript 9 on the
coordinates.
If x = f(u, v), then a small increase 6x in x is given by
Now
(a) P is the point (x, y)
(b) Q corresponds to smal! changes from P.
Ox Ox oy oy
Ox ati Ou
u + bv and éy 5s gee Ou
1 1 i
Ox OX
pe —=()
Mreg oe fa rie sara,
2
Then, taking out the factor éu from the first column and the factor é6v from the
second column, this becomes
The area of the approximate parallelogram is twice the area of the triangle.
Ox Ox
; i : OW NOU
. Area of parallelogram
= 6A; = _ |Oudv
Oe OY
Ou Ov
Expressing this in differentials
QA = ok ay du dv
Ou OV
and, for convenience, this is often written
dA Bey) dudv
O(U, V)
Surface and volume integrals
183
Ox y)
A(u, v) is called the Jacobian of the transformation from the Cartesian
coordinates (x, y) to the curvilinear coordinates
(u, Vv).
ox Ox
J(u,v) = uv(4u 4+ v)
Because
- = Nea is =
OE a) ae
OY» Oy
i 1) ae 2uyv
2u+v u
J(u,v) = = 4u’v + 2uv? — uv?
y2 2uv
Sometimes the transformation equations are given the other way round. That
is, where uw and v are given as expressions in x and y. In such a case J(u, v) can
be found using the fact that
O(x, y) 1
O(Uu, V) O(U, V)
O(x, y)
For example, if the transformation equations are given as u = x*+y? and
v = 2xy then
184 Vector Analysis
Because
Ou Ou
OUD) NOx OV 2x may, Pee:
a(x, y) OV Ov Vox Feit)
Ox Oy
and so
- Ones iy x 1
De ay ce ») ~ 4(x2 — 2)
O(x, y)
Now u—v = x2 —2xy +y2 = (x—y)*
and u+v=x2 +2xy+y2 =(x+y)?
and so x2 —y2 = (x—y)(x+y) = Ju—v/ut+v = Vi? — v? giving
Co
Ae ye
Example 1
Stas
aa ae
ap 7 i
Sing
ae ne
ar
oy
ap 1 COSe
cos@ —rsing
Tiss?) = =rcos*6+rsin?26=r
sin @ rcosé
Example 2
Because
me Dawe Ox
anon Ox
eget
i ae roi au Ov i
Oy Oy
VU De ay 2v ame
Ox Oy
PC mals ae) Ri
" Olu,v) Jax ay
Ov OV
Also x2 +y? = (u? — vv’)? + (2uv)* = ut — Qu’ + v* + 40?Vv?
a lr ey Ge)
[= 4| [uw + v*)3dudv
R
One more.
Example 3
By substituting x = 2uv and y=u(1—v) where u > O and v > O, express
of Ox
Ox
X = 2uv Ox
ie 2u
Ou 2v
Oy ay ts
A
UV
=u—
1—v av u
Ou
,
Ox Oy
du
Ou = yi WEEE:
J(u,v) = ae y)_ |Ou
du
AuY Ox Oy 2u —u
Ov OV
y 1)
yl Ti Sos
= Ou
=— 2u
=4 l 0|
if
A(x, y)| _
en vy)|~ 2u
O(X, y, Z) 2 5
Then Ade O(u, Vv, W) lav av. av
ORE OV OZ
Ow OW Ow
and the element of volume dV = dx dy dz becomes
dV =|J (u,v,w) |dudvdw
[Jfoc.v O(u,
O(x, ¥;2Z) dudvdw
Vv, W)
Now for an example, so move on
eaee
eee eee eee ener ree nee... ——————_—_—_—_—_—_—
188 Vector Analysis
Example 4
OX Opals : OZ
OF
=—_ = sin 9 cos @ Ar
—— —— sin sin @ Be
— = cos Q
Ox Vv ; Oz
— 70 = rcosécos @ 20
——_— = rcos@sin @ 70
—_— SS rsin @
Ox F : oy 1. 4 OZ)
En ee ee Ag 7 SINE COse a6 >
r sin @
Because
} cos@ sing
J(r, 6, ¢) =r? cos? @siné|
—sing cos¢
; cos@ sing
Sar Sa)
—sing cos@
: : 9 cos@ sin
= (r* sin? 6 +r’ sin @ cos? 6) | ?
—sing cosd
= r* sin (sin? 6 + cos” @)(cos? ¢ + sin? 6) = r sing
Do eats: .
And that is about it. Check carefully down the Review summary and the
Can You? checklist that now follow, before working through the Test
exercise. The Further problems give additional practice.
Review summary
1 Surface integrals
1=| JR
f(x, y)da= | ‘F(x, y) dy dx
Ris
2 Surface in space
=|Josxay]t+ (Go) a
: faz\? (az\?
Pe (X,VitZ)
First octant:
A Oe 7 0)
x ye
Z
P (x,y, Z)
8 aT COS r= /fx2+y2
z y= tsin ? 6 = arctan(y/x)
o : L=Z Laz
9 y
(x, Y, 2)
(r, 8, )
190 Vector Analysis
4 Elements of volume
(a) Cartesian coordinates
dv = 6x by 6Z
Z roe
te 780 aa r)Zz
halt iz
fou | Ls,
< ag
Cela?= |
é6v =r 6r 60 6z
rsin6 do
x a) pdo
dv =? sin @ 6r 6056
5 Volume integrals
v=||fazayax
t= J fre: y, z)dzdydx
Surface and volume integrals
191
J(u; v; w) ==
OX, Y, Z) |ney ein?
FOX “Oy
ae
Oz
ie
Ou, v,w) |dv av av
ax
Ow
dy
OW
oz
Ow
= [| fou V, wee
(4 Can You?
Checklist 6
Check this list before and after you try the end of Program test.
Test exercise 6
Determine the area of the surface z = \/x2 + y? over the region bounded by
xv +y2 =4,
1
2 Evaluate the surface integral J = |é dS where ¢ = ————— over the surface of
s [x2 4 y?
Further problems 6
5 Determine the area of the paraboloid z= 2(x? +y) cut off by the cone
Z= f/x? +y2.
6 Find the area of the cone z* = 4(x?+y?) which is inside the paraboloid
2 == 2 (x72 7),
Cylinders x* + y? = a* and x? + z? = a? intersect. Determine the total external
surface area of the common portion.
Determine the surface area of the sphere x” + y? +z* =a? cut off by the
cylinder x? + y? = ax.
A cylinder of radius b, with the z-axis as its axis of symmetry, is removed from a
sphere of radius a, a > b, with center at the origin. Calculate the total curved
surface area of the ring so formed, including the inner cylindrical surface.
Find the volume enclosed by the cylinder x? + y* = 9 and the planes z = 0 and
Z=S-xX.
Determine the volume of the solid bounded by the surfaces y = x, x = y*, z=2
and x+y+z=4.
Find the volume of the solid bounded by the plane z = 0, the cylinder
x? + y? = a* and the surface z = x? + y?.
A solid is bounded by the planes x =0, y=0, z=2, z=x and the surface
x? + y* = 4. Determine the volume of the solid.
Find the position of the center of gravity of the part of the solid sphere
x? + y* + z* = a” in the first octant.
A solid is bounded by the cone z=2,/x2+y?, z>0, and the sphere
x? + y2 + (z—a)* = 2a”. Determine the volume of the solid so formed.
Be deg
Determine the volume enclosed by the ellipsoid 5S7 - ig“521h
Find the volume of the solid in the first octant bounded by the planes
x=0, y=0,z=0, z=x+y and the surface x? + y” =a’.
18 Express the integral ||e?+») axay in terms of u and v, using the
transformations u=x+y,v=x—y.
Lo Determine an expression for the element of volume dx dy dz in terms of u, v, w
using the transformations x = u(1 —v), y= uv, Z = uvw.
20 A solid sphere of radius a has variable density c at any point (x, y, Z) given by
c =k(a—z) where k is a constant. Determine the position of the center of
gravity of the sphere.
21 Calculate |[#y? axay over the triangular region in the x-y plane with
2 pr/4—y2 y
22 Evaluate the integral | dxdy by transforming to polar
0 ipso Reve
coordinates.
1 py
xy?
23 Evaluate [ = || dx dy.
0) 0 ye wae ale
24 Find the volume bounded by the cylinder x? + y* = a”, the plane z = 0 and the
surface z = x* + y”. Convert to polar coordinates and show that
na‘
25 By changing the order of integration in the integral
fen I[ We dy dx
Gunn) eye
show that I = 4a? In(1 + v2).
Vectors aie
Learning outcomes
When you have completed this Program you will be able to:
e Define a vector
e Represent a vector by a directed straight line
e Add vectors
e Write a vector in terms of component vectors
Write a vector in terms of component unit vectors
e Set up a coordinate system for representing vectors
e Obtain the direction cosines of a vector
e Calculate the scalar product of two vectors
e Calculate the vector product of two vectors
e Determine the angle between two vectors
e Evaluate the direction ratios of a vector
e Obtain the scalar and vector triple products and appreciate their geometric
significance
195
196 Vector Analysis
Ce) Physical quantities can be divided into two main groups, scalar quantities and
vector quantities.
(a) A scalar quantity is one that is defined completely by a single number with
appropriate units, e.g. length, area, volume, mass, time, etc. Once the
units are stated, the quantity is denoted entirely by its size or magnitude.
(b) A vector quantity is defined completely when we know not only its
magnitude (with units) but also the direction in which it operates, e.g.
force, velocity, acceleration. A vector quantity necessarily involves
direction as well as magnitude.
So (a) a speed of 10 km/h is a scalar quantity, but
(b) a-velocity of 10.Kmi)nvdite noth ised 4). .- a.2-ee quantity.
vector
So that:
(a) gAatemperature ol 100° Cus a. eee = quantity.
(b) An acceleration of 9.8 m/s” vertically downwards isa............ quantity.
(Cg lhe weishitrored /ke massasea eee quantity.
(dy The sum One SUO IS ain. er quantity.
(e) A north-easterly wind of 20 knots isa............ quantity.
|(a) scalar (b) vector (c) vector (d) scalar (e) vector
Since, in (b), (c) and (e) the complete description of the quantity includes not
only its magnitude, but also its ............
Move on to Frame 6
Vectors
197
Vector representation
A vector quantity can be represented graphically by a line, drawn so that:
(a) the length of the line denotes the magnitude of the quantity, according to
some stated vector scale
(b) the direction of the line denotes the direction in which the vector quantity
acts. The sense of the direction is indicated by an arrowhead.
e.g. A horizontal force of 35 N acting to the right, would be indicated by a line
and if the chosen vector scale were 1 cm = 10 N, the line would be
Be i cm long.
Note that BA would represent a vector quantity of the same magnitude but
with opposite sense.
Woe
o m
A AB=a A BA=-AB
On to Frame 8
So on now to Frame 11
To find the sum of two vectors a and b then, we draw them as a chain, starting
the second where the first ends: the sum ¢ is given by the single vector joining
the start of the first to the end of the second.
e.g. if p =a force of 40 N, acting in the direction due east
q = 4 force of 30 N, acting in the direction due north
then the magnitude of the vector sum 7 of these will forces will be............
aD Because
ma p 4 ra
f=V 2500250 N
Vectors 199
a+b=AC
AC+c=AD
a+b+c=AD
AD +d = AE
a+b+c+d=AE
i.e. the sum of all vectors, a,b, c,d, is given by the single vector joining the
start of the first to the end of the last —- in this case, AE. This follows directly
from our previous definition of the sum of two vectors.
R
Similarly:
Q
it PO 4,OR 4. RS-7S
lias arte ite
Now suppose that in another case, we draw the vector diagram to find the sum
of a,b,c,d,e, and discover that the resulting diagram is, in fact, a closed
figure.
What is the sum of the vectors
a,b,c,d,e in this case?
Because we said in the previous case, that the vector sum was given by the
single equivalent vector joining the beginning of the first vector to the end of
the last.
But, if the vector diagram is a closed figure, the end of the last vector
coincides with the beginning of the first, so that the resultant sum is a vector
with no magnitude.
200 Vector Analysis
Without drawing a diagram, we can see that the vectors are arranged in a
chain, each beginning where the previous one left off. The sum is therefore
given by the vector joining the beginning of the first vector to the end of the
last.
* Sum = AF
In the same way:
AK + Ki LPs PO = «rn csey See
Because
AB + BC — DC —- AD = AB+BC+CD+DA
and the lettering indicates that the end of the last vector coincides with the
beginning of the first. The vector diagram is thus a closed figure and therefore
the sum of the vectors is 0.
Now here are some for you to do:
(a) POOR
SRS sos Looe. aeena
(eA Ce a I ice eee
ee ae
(c) GH+HJ4JK+KL4+LG=
(d) AB+ BC+ CD + DB =
When you have finished all four, check with the results in the next frame
———
eee
Vectors
201
(a) PQ+QR+
RS + ST = PT
(b) AC + CL— ML =AC+CL+LM = AM
(c) GH+ HJ) +JK+ KL+LG=0
[Since the end of the last vector coincides
with the beginning of the first.]
(d) AB+BC+CD
+DB = AB
The last three vectors form a closed D
figure and therefore the sum of
these three vectors is zero, leaving Cc
only AB to be considered.
A B
Now on to Frame 18
Example 1
ABCD is a quadrilateral, with G and H the mid-points of DA and BC
respectively. Show that AB + DC = 2GH. cS
We can replace vector AB by any chain of
A
B vectors so long as they start at A and end at B
e.g. we could say
AB = AG-+ GH + HB
202 Vector Analysis
[Dc =DG+GH + AC
So we have:
AB = AG+ GH + HB
DE =]DG + Gi + HC
Example 2
Points L, M, N are mid-points of the sides AB, BC, CA of the triangle ABC.
Show that:
(a) AB+BC+CA=0
(b) 2AB + 3BC
+ CA = 2LC
(c) AM+BN+CL=0
Because
Ch= CAVE AL
So AM+BN+CL=AB+ BM+BC+CN+CA+ AL
Finish it off
[ AM +BN+Cr =0|
: em peat
Because AM + BN + CL = (AB + BC + CA) + 5(BC + CA + AB)
Now AB+BC+CA is a closed figure .. Vector sum =90
and BC + CA + AB is a closed figure... Vector sum = 0
-, AM+BN+CL=90
Here is another.
204 Vector Analysis
Example 3
ABCD is a quadrilateral in which P and Q are the mid-points of the diagonals
AC and BD respectively.
Show that AB+ AD+CB+CD = 4PQ
A D
To prove that AB + AD + CB + CD = 4PQ
Taking the vectors on the left-hand side, one at a time, we can write:
IN 2140) (Oe
IND) eau Noe 12085 2.0)By)
Che CPSPO--
Ob CD Cr Oo
AP +CP=0
Example 4
Prove by vectors that the line joining the mid-points of two sides of a triangle
is parallel to the third side and half its length.
B Cc
Now express DA and AE in terms of BA and AC respectively and see if you can
get the required results.
Then on to Frame 30
Z1 +z 5)it
(442)j
= OB=(2+ = 7i+ 6j
Z2-Z, =i+j
Because we have
OA + AB = OB (from diagram)
*. AB=OB-—OA
= (Si — 2j) — (31+ Sj) = 21 - 7j
On to Frame 36
Vectors in space
The axes of reference are defined by
the ‘right-hand’ rule.
OX, OY, OZ form a right-handed set
if rotation from OX to OY takes a
right-handed corkscrew action
along the positive direction of OZ.
x
Similarly, rotation from OY to OZ gives right-hand corkscrew action along the
positive direction of ............
208 Vector Analysis
x
Let i = unit vector in OX direction
j = unit vector in OY direction
k = unit vector in OZ direction
Then OP = ai+bj+ck
Also OL? = a? + b* and OP? = OL? + 7
OP =a +b? +c?
So, ifr=ai+bj+ck, thenr= V/a2+b?+¢c?
This gives us an easy way of finding the magnitude of a vector expressed in
terms of the unit vectors.
Now you can do this one:
If PQ = 4i+ 3j+ 2k, then |PQ|=............
Because
PQ = 41+ 3j + 2k
PQ| = Ve 43? 422
=V16+9+4
= n/29
= ONaSO
Direction cosines
Let OP =¥ = ai + bj-+'ck
Then
a
Pi ike typ aT COS
b
7 = Cos 6 b=rcosB
c
7 = Cosy ae COS
X
Also a? +b?+C=r 2
r = 3i— 2j+6k
.ad=5, Pe —2, c=6, r=v9 +44 36
tea
AD sey
Tie ee 7
poe = 2. “eo
Because we have:
B OA - OB = OA: OB-cos 6
? = 97 1COS45-
re ssve2
5 A —
90°
Vectors
211
[9]
Because in this case a- b = abcos 90° = ab0 = 0. So the scalar product of any
two vectors at right-angles to each other is always zero.
And in this case now, with two vectors in the same direction, 0 = 0°
Now suppose our two vectors are expressed in terms of the unit vectors i, j
and k.
Let a= d)i+ a2j + a3k
and b = bi+ boj + b3k
Because
a-b =a,b,-1+a,b2-0+4,b3 -0+4 azby -0 + dob2-1+d2b3-0
+ a3b, -0 + a3b2 -0+43b3-1
a- b=a,b; + aob2 + a3b3
i.e. we just sum the products of the coefficients of the unit vectors along the
corresponding axes.
212 Vector Analysis
For example:
If a= 2i+ 3j + 5k and b = 4i + 1j+ 6k
thena-b=2x4+3x1+5x6
= 843 +/30
— wl - a -b=41
for you: If p = 3i— 2+ 1k;
One q = 2i+ 3j-— 4k
Because
p-q=3x2+4+(-2) x3+1x(-4)
=6-6-4
= .. p-q=-4
Now on to Frame 47
6=0°, |axb|/=0
§ = 90°, |a x b| = ab
Then:
a x b = a, bi x i + ay boi x j + ay b3i x kK + azbyj x i+ agboj x j
+ a2b3j x k + a3b1k x i+ a3bok x j + a3b3k x k
Now with the results of (a) and (b), and this last reminder, you can simplify
the expression for a x b.
Remove the zero terms and tidy up what is left.
Because
ax b=a4,b,0+ aibok 35 a, b3(—j) ai dzb,(—k) + dob20 + dzb3i
and that is the easiest way to write out the vector product of two vectors.
Notes: (a) The top row consists of the unit vectors in order i, j, K.
(b) The second row consists of the coefficients of a.
(c) The third row consists of the coefficients of b.
For example, if p= 2i+4j+3k and q=i+ 5j-— 2k, first write down the
determinant that represents the vector product p x q.
Unit vectors
Coefficients of p
Coefficients of q
Because
ij k
pxq=|2
—s 4 3 — le - if; ; +k
J ;|
te > 5 2 2, I ats
Because
1 j ik
Pxq=/3 -4 2
2 5 -l
Rel i|-i), +k; a
oo =a) 2 —-l 2 S
= i(4 — 10) — j(-3 —4) + k(15 +8)
= —6i+ 7j + 23k
Remember that the order in which the vectors appear in the vector product is
important. It is a simple matter to verify that:
q x p = 61 — 7j — 23k = -(p x q)
On to Frame 53
Nu
+ '
Then
(PP’)? = (1-1)?
+ (m—m) 2
ein-n)
=P 9.04 +n? —2n-n +n?
48 =—2-m-m-+m
= (2 4m? +n?) + (7 +m? +n) — 2( + mm! + nn’)
But (2 + m2 +n?) = 1 and (/? +m? +n”) = 1 as was proved earlier.
- (PP)? =2—2(I +mm'+nn') (a) >
216 Vector Analysis
That is, just sum the products of the corresponding direction cosines of the
two given vectors.
So, if [1,m,n] = (0.54, 0.83, —0.14]
and [{I',m’,n’| = [0.25, 0.60, 0.76]
the-ansle berweens the Vectors 15 @ 2 e.
C=58 TH
Because, we have:
Cosi — + mm’ +. nn!
= (0.54)(0.25) + (0.83)(0.60) + (—0.14)(0.76)
= O15 + 0.4980 — 0.1064
= (OER) = OL! = 0.5266
() = Sor
56)
Because
Because
Because
eT re ety bed Bnd
COS? = 759 a9 * 29 29° V29 29
bbaacierds
Dolo 29
7
=—-=0,2414 .. 0@=76°2'
29
Now on to Frame 59
218 Vector Analysis
Direction ratios
Note: The direction ratios can be converted into the direction cosines by
dividing each of them by r (the magnitude of the vector).
Move on to Frame 60
Triple products
We now deal with the various products that we form with three vectors.
Loe jae
then B x C=\|b, by b;
Cy Cy Cz
OES fm
". A- (Bx C) = (a,i+ ayj+azk)- |b, by b,
Cy Cy Cz
Multiplying the top row by the external bracket and remembering that
i-j=j-k=k-i=0 and i-i=j-j=k-k=1
ay ay az
we have A- (Bx C) =|b, by bz
Vectors 219
Example
If A= 2i- 3j+ 4k; B=i-2j—3k; C= 2i+j+
2k;
2° 33 4|
Cie AthCy Sette a2 — 9
2 1 2
Because
2 -3 4
AP BoC) <1«:—2inee3
2 1 2
= 2(=4 3) £3024 6)44(1 4) = 429
As simple as that.
Cy Cy Cz Cy Cy Cz
’. B-(Ax
C) =—-A- (Bx C)
i.e. a change of vectors not in cyclic order, changes the sign of the scalar
triple product.
ay ay az
Example
If A=i+2j+3k; B= 2i- 3j+k; C = 3i+j —2k
Because
| a, 3
A-(Bx©)=|2 -3 +3(2+9) =52
1)=1(6—1)—2(-4—3)
3 1 -—2
Coplanar vectors
The magnitude of the scalar triple product |A - (B x C) |is equal to the volume
of the parallelepiped with three adjacent sides defined by A, B and C.
Example 1
Show that A=i+ 2j-3k; B= 2i—j+ 2k; and C = 3i+j—kare coplanar.
We just evaluate A> (Bi C= ae and apply the test.
Vectors 221
Because
i >
AS (BSC)(2., =—1 2|= 1(1 — 2) — 2(-2 — 6) — 3(2 + 3) =0.
3) |
Example 2
If A= 2i—-j+3k; B=3i+2j+k; C=i-+pj+
4k are coplanar, find the
value of p.
The method is clear enough. We merely set up and evaluate the determinant
and solve the equation A - (B x C) =0.
Because
2 =1; 3
il p 4
Example 3
Determine whether the three vectors A= 3i+2j—k; B= 2i—j+ 3k;
C =i — 2j + 2k are coplanar.
Work through it on your own. The result shows that
Because
Sa pe" eh
in this case A“(Bx C)=\2. =1 i= 13
1 —2 2
A x (BxC)
are called the vector triple products.
and (A x B) x C
1
; b, bz| .|by bz| . |bx By
Cy Cz Cy Cz Cx Cy
Cy Cy Cz
i j k
ax ay az
Then Ax(BxC)=
by b, by b; by b,
Gy G GanG, &X G
i j k
ay ay az
by Del |Ds Da Ds by
Gy Gane Cx Cy
Example 1
If A= 2i-3j+k,; B=i+2j—k; C=3i+j+3k; determine the vector
triple product A x (B x C).
We Start off with Bx CG=............
Vectors
223
he B x C= 7i- 6j —5k i
Because
ta es ae 5
B GL it Sy = i(6+ 1) -j(3+3)+k(1-6)
0 3 = 7i — 6j — 5k
Because
a od
Ax(BxC)=|2 -3 1
7 =6 5
= i(15 + 6) —j(—10 — 7) + k(—12 + 21)
= 21i+ 17j+ 9k
That is fundamental enough. There is, however, an even easier way of
determining a vector triple product. It can be proved that
A x (B x C) = (A: C)B—-
(A- B)C
and (Ax B)xC=(C-A)B-—(C-B)A (11)
Make a careful note of the expressions: then we will apply the method to the
example we have just compieted.
Example 2
If A= 3i+2j—2k; B=4i-—j+3k; x B) x C
C= 2i-—3j+k determine (A
using the relationship (A x B) x C= (C-A)B—(C-B)A.
(A Bix, Coe reste
—_————————
LL
224 Vector Analysis
Because
(A x B) x C=(C-A)B-(C-B)A
= (6 —6 —2)(4i -j + 3k) — (8 +3 + 3)(3i + 2j — 2k)
— —2(4i — j + 3k) — 14(3i + 2j — 2k)
= —50i— 26f 422k
Now one more.
Example3
IfA=i+3j+2k; B=2i+5j—k; 3k
C=i+2j+
Ax (BSC) =o cage aos
(ASC Bc Ca) feo eo
Example 4
If A= Si— 2j+ 3k; B=3i+j-
2k; C =i-—3j+
4k; determine
(a) the scalar triple product A - (B x C)
(b) the vector triple products (1) A x (B x C)
(2) (Ax B) x C.
Finish all these and then check with the next frame
(a) A-(BxC)=_12 |
(b) (1) A x (Bx C) = 624 + 44j — 74k
(2) (Ax B) x C = 1094 + 7j — 22k
5)Review summary
1 A scalar quantity has magnitude only; a vector quantity has both
magnitude and direction.
2 The axes of reference, OX, OY, OZ, are chosen so that they form a right-
handed set. The symbols i, j, k denote unit vectors in the directions OX,
OY, OZ, respectively.
=r = Va? +b? +2
If OP = ai + bj + ck, then |OP|
3 The direction cosines |l,m,n| are the cosines of the angles between the
vector and the axes OX, OY, OZ respectively.
a b Cc 5
For any vector: !=—,m=—,n=—; and 12 +m*+n*=1
r
226 Vector Analysis
So off you go
f4 Can You?
Checklist 7
Check this list before and after you try the end of Program test.
On a scale of 1 to 5 how confident are you that you can: Frames
e Define a vector? (7) oC)
Yes al | i LJ No
e Represent a vector by a directed straight line? Ce)» Go)
Yes te a No
e Add vectors? to (47)
Yes | Gi & tial No
e Write a vector in terms of component vectors? (48) (Bo)
Yes O LJ J LJ | No
e Write a vector in terms of component unit vectors? (B17)
>(35)
Yes * | i No
e Set up a coordinate system for representing vectors? (Ge) (Bs)
Yes [ it) | LJ a No
e Obtain the direction cosines of a vector? (39)
to(40)
Yes [ L] ha No
e Calculate the scalar product of two vectors? a
Yes | | bal No
e Calculate the vector product of two vectors? » GZ
(47)
Yes ie ] No
e Determine the angle between two vectors? > Ga)
(3)
Yes | 4 eh No
8 Find the unit vector perpendicular to each of the vectors 2i— j+k and
3i + 4j — k, where i, j, k are the mutually perpendicular unit vectors. Calculate
the sine of the angle between the two vectors.
9 IfA is the point (1, —1, 2), Bis the point (—1, 2, 2) and C is the point (4, 3, 0),
find the direction cosines of BA and BC, and hence show that the angle
ABC = 69°14’,
10 If a=3i—j-+ 2k, b =i+ 3j—2k, determine the magnitude and direction
cosines of the product vector (a x b) and show that it is perpendicular to a
vector c = 91 + 2j + 2k.
11 aand b are vectors defined by a = 8i + 2j — 3k and b = 3i — 6j + 4k, where i,
j, k are mutually perpendicular unit vectors.
(a) Calculate a-b and show that a and b are perpendicular to each other.
(b) Find the magnitude and the direction cosines of the product vector a x b.
12 If theposition vectors of P and Q are i + 3j — 7k and Si — 2j + 4k respectively,
find PQ and determine its direction cosines.
13 If position vectors, OA, OB, OC, are defined by OA=2i-—j+3k,
OB = 3i+ 2j — 4k, OC = —i + 3j — 2k, determine:
(a) the vector AB
(b) the vector BC
(c) the vector product AB x BC
(d) the unit vector perpendicular to the plane ABC.
14 lfA=2i+
3j —4k; B= 31+ 5j+ 2k; C=i-—2j+
3k; determine
A- (Bx C).
158 IfA=2i+j-—3k, B=i-2j+
2k; C=3i+2j—k; find Ax (Bx C).
16 IfA=i-2j+3k; B=2i+j-2k; C=3i+2j+k; find
(a) Ax(BxC); (b) (Ax B)xC.
17 Find the scalar triple product of
(a) A=i+2j-—3k; B=2i-—j+4k; C=3i+j-
2k.
(ob) A=2i-3j+k, B=3i+j+2k,; C=i+
4j- 2k.
(c) A=—2i+3j-2k; B=3i-—j+3k; C=2i-—5j+KkK.
18 Find the vector triple product A x (B x C) of the following.
(a) A=3i+j-—2k; B=2i1+4j+3k, C=i-2j+k.
(b) A=2i—-j4+3k; B=i+4j-—5k; C=3i-2j+Kk.
(c) A=4i+2j-3k; B=2i-3j+2k; C=3i—3j+Kk.
19 Determine the value of p such that the three vectors A, B, C are coplanar when
A= 2i+j+4k; B=3i+2j+pk; 2k.
C=i+4j+
20 IfA=pi-6j—3k; B=4i+3j—-k, C=i-3j+2k
(a) find the values of p for which
(1) A and B are perpendicular to each other
(2) A, B and C are coplanar.
(b) determine a unit vector perpendicular to both A and B when p = 2.
Frames
Vector jto} 68 |
differentiation
Learning outcomes
When you have completed this Program you will be able to:
Differentiate a vector field and derive a unit vector tangential to the vector
field at a point
Integrate a vector field
Obtain the gradient of a scalar field, the directional derivative and a unit
normal to a surface
Obtain the divergence of a vector field and recognise a solenoidal vector
field
Obtain the curl of a vector field
Obtain combinations of div, grad and curl acting on scalar and vector
fields as appropriate
230
Vector differentiation
231
Differentiation of vectors
In many practical problems, we often deal with vectors that change with time,
e.g. velocity, acceleration, etc. If a vector A depends on a scalar variable t, then
A can be represented as A(t) and A is then said to be a function of t.
eA =a: ayj+a,k then ay, ay, a, will also be dependent on the
parameter tf.
i.e. A(t) = a,(t)i+ a,(t)j + a,(t)k
A(u)
Example 1
If A = (3u2 + 4)i + (2u — 5)j + 4u°*k, then
dA
232 Vector Analysis
te = 6ui + 2j + 12u’k
du
2
GRAS yes
If we differentiate this again, we get aiden 6i + 24uk
Nap
When u = 2, oe 12i+2j+48k and —, = 6i+ 48k
du du
Because
dA
aula f12e AMY:
2 AGA21/2 ee IZAS2 Lza= 49 52
and d°A
ata {6°2 4+ 48°)
211/2* = (2340) 1/2 = 46.37
Example 2
If F = isin 2t + je* + k(t? — 4t), then when t = 1
dF
Gp Gp
iad re ce
dF e 36
Gp = 200824 + 3e j—k
d°F
= —4sin 2i + 9e°j + 6k
dt2
ar
TA ae ;ae
q@|a cc
rr
Vector differentiation 233
Because
Example 3
If A= (u+3)i—(2+u7)j + 2u°k, determine
dA dA se aA
du (b) du2 (©) du du2 ike
Work through all sections and then check with the next frame
Example 4
A particle moves in space so that at time ft its position is stated as
X= 2t+3, y=t?4+3t, z=t? + 2t?. We are required to find the components
of its velocity and acceleration in the direction of the vector 2i + 3j+4k
when f= 1.
First we can write the position as a vector r
r= (2t+ 3)i+ (7 4+30)j+ (0 + 227)k
Then, at ¢=1
234 Vector Analysis
dr dr :
ae
—— 21i + Sji + 7k;i df
a 2 ir 10k
Because
dr : 3 2
dr
ST3 Atg= = 1) at
—_ = 23 Hse Sylar
SI 7k
d’r i
and dz 2j+ (6t+ 4)k
d’r
*‘ Atta = 1: qe
a 21 \
bee 10k
21+ 33+ 4k 1 ; :
— 2i+ 3j+4k
/4+9+16 39 | }
Because
il } il
—— (21 + 5j + 7k) - (21+
3) + 4k) ==—(44 154 28
Mie 4
V29
se thf!
2
Because
d’r dr
1
= —= (2j + 10k)- (2i+ 3j + 4k)
= ——(6 +40)
denoting <{A(u)} is
fe) A(u)
lA)
< caw} Seg”
In simpler notation, this becomes:
If r = d,i + a,j + ak then the unit tangent vector T is given by
_ dr/du
~ |dr/du|
Example 1
Determine the unit tangent vector at the point (2, 4, 7) for the curve with
parametric equations x = 2u; y = u2 + 3; z= 2u? +5.
First we see that the point (2, 4, 7) corresponds to u = 1.
The vector equation of the curve is
ee :
ie= 21+ 2uj + 4uk
and at u= 1, m= 24 +2) + 4k
Because
dr
= {444416}?
= 2417 ~2V6
du
dr
T— du 2i + 2j+ 4k a
ar] 3v6 +j+ 2k}
= — = 1
du
Let us do another.
Example 2
Find the unit tangent vector at the point (2, 0, 7) for the curve with
parametric equations x = 2sin@; y= 3cos@; z= 20.
We see that the point (2, 0, 7) corresponds to 6 = 7/2.
Writing the curve in vectorform r= 6-5. <+.:.
r= 2sin0i+3cos6j+
26k
Finish it off
238 Vector Analysis
And now
Example 3
Determine the unit tangent vector for the curve
x=3t; y=2t?; z=t+t
at the point (6, 8, 6).
ODS VOUL OWS Las ire canes
_-!1 (3 4 914 5k
c fag ee )
Example
If F = 2uvi + (u* — 2v)j + (u+v?)k
OF OF
< = 2vi+ 2uj + k; —— = 2ui — 2j + 2vk
Ou Ov
OF . O°F
ou2 25; auav :
Example 1
If F= (3t? + 4t)i+ (2t — 5)j+ 4t°k, then
3 3 3 3
|Fat =i| Gr +4e)de+j| (2t~5)dr+k| af dt =... ©.clter
1 ] 1 1
42i — 2) + 80k
Because
3 3
|F dt = ite’PI ith sh) + Ke]
1 1
= (45i — 6j + 81k) — (3i— 4 +k) = 42i — 2j + 80k
Here is a slightly different one.
Example 2
If F = 3ui+ uj + (u+ 2)k
and V = 2ui — 3uj + (u— 2)k
2
evaluate |(F x V)du.
0
First we must determine F x V in terms of u.
240 Vector Analysis
Because
i j k
ExV=|3u w (u+2)
2u —3u (u—-2)
which gives the result above.
2
Then |(Be VU =o eee
JO
Because
See 3 us yt ra P
[Rx v)du = (F+5+30)t
|— == |(5 si) ee, (5
ee +30 k
2
|(x V)du = (4+$+ 12) G - 20)j- (8+ 24)k
0)
=4{14i+
13j— 24k}
Example 3
If F= A x (B x C) where
A—8i-{(2i— 3)j44tk
B= 21+ 4tj+ 3(1-tH)k
Caen S37) = rk
1
determine |Fdt.
0)
First we need to find A x (B x C). The simplest way to do this is to use the
relationship
Vector differentiation
241
Because
Examples of vector fields are force, velocity, acceleration, etc. F(x, y, z) can be
defined in terms of its components parallel to the coordinate axes, OX, OY,
OZ.
That is, FQ, y, z)=HKi+kj+ 2K.
Note that, while ¢ is a scalar function, grad ¢ is a vector function. For example,
if ¢ depends upon the position of P and is defined by ¢ = 2x*yz’, then
8 rad ¢ = 4xyz7i
Ve + 2x?z>¥J + 6x’yz’k
Me
Notation
aos 0 : 0 0
aha lintig tks te
(Ohat O Or :
where (ii ky) is called a vector differential operator and is
oy
denoted by the symbol V (pronounced ‘del’ or sometimes ‘nabla’)
0 O O
Woe V=(i j k
(ix +i + =)
Beware! V cannot exist alone: it is an operator and must operate on a stated
scalar function ¢(x, y, Z).
If F is a vector function, VF has no meaning.
So we have:
O O O
Vop= = grad¢
grad@ = ( ala =)¢
(i—+j—+k—}¢
Example 1
If 6 = xyz + xy*z?, determine grad ¢ at the point Py E32):
4(214
+ 8j+ 18k)
Because
Example 2
If A=x?zi+xyj+y?zk
and B=yz*i+ xzj+x?zk
determine an expression for grad (A - B).
This we can soon do since we know that A - B is a scalar function of x, y and z.
Rirstetijetigy oA 0h tee ore:
4 xyzey)
LA. B = xyz? + xyz
Because
if 6 = A-B= (x2zi+ xyj + y?zk) - (yz°i + xzj + x°zk)
= XP EX ZEN
Oo
= 3x7 yz? + x2y + 2x7y2z = x?y(3z? + 1 + 2yz)
VW (A-B) = 2xyz(z’ +14 yz)it+ x?z(2? + 1 + 2yz)j
x y(3z Sb eyzk
Now let us obtain another useful relationship.
ee If OP is a position vector r where
a r=xi+yj+zk and dr is a small
displacement corresponding to
3 changes dx, dy, dz in x, y, z respec-
tively, then
——
Oh M dr = dxi+dyj+dzk
Oe elo) 0@
d@¢@-dr @-dr =—eee
grad ets ie
Because
ce Od ‘ Odie Od '
grad@-dr (seis mh+ 7k) (dxi+ dyj+dzk)
_> Og To)
as dx iay© ay +5 7 az
= the total ae dé of é
That is
dé = dr - grad ¢
(14)
This will certainly be useful, so make a note of it
Vector differentiation
245
Directional derivatives
We have just established that z Q (r+ dr)
dé = dr- grad ¢ P(r) dr dz
If ds is the small element of arc between P (r) d S Wax
and Q (r+ dr) then ds = |dr|
dr dr
ds_ |dr|
dr
and a is thus a unit vector in the direction of dr.
eed¢ dr ,
ds ds ad?
d ‘
If we denote the unit vector = by a then the result becomes
d
5. =a: grad ¢
ae is thus the projection of grad ¢ on the unit vector a and is called the
directional derivative ot ¢ in the direction of a. It gives the rate of change of ¢
Example 1
Find the directional derivative of the function ¢ = x?z
+ 2xy* + yz at the
point (1, 2, -1) in the direction of the vector A = 2i + 3j — 4k.
Because
Ob _ 2. Og |
By 7 ke + V's Emer iy Se)
Panis + 2 LZ
Then, at (1, 2, — 1)
Because
A=2i+3j-4k ©. |A|=V4+9+16=v29
J A
a ees dk)
|A| 29
So we have V¢ = 61 + 9j — 3k eal (21 + 3j — 4k)
V29
d¢
ae ee
Because
Lee
ds =pa 4 3ji — 4k) )-- (64+
(61 + 9]ie— 3k )
il
Seed (1D 27) ee
V29 J29
That is all there is to it.
(a) From the given scalar function ¢, determine V¢.
(b) Find the unit vector a in the direction of the given vector A.
(c) Then a
do. a- Vo.
Vector differentiation 247
Example 2
Find the directional derivative of ¢ = xy + y2z + z2x at the point (1, — 1, 2) in
the direction of the vector A = 4i + 2j — Sk.
Same as before. Work through it and check the result with the next frame
3.43
ts
Because
op=x*y+y*2-z"x
", Vb = (2xy + 27)i + (x? + 2yz)j + (y? + 22zx)k
At (1, —1, 2), Vd=2i—3j+5k
A=4i+2j-5k .. |A|=V164+4+425
= V45.4.3V5
; 1
egies 3V5
ee |ay ne
j Sik)
de 1 a0 fg i
eto ae ae ee (2i — 3j + Sk)
1 iG anyCEE —23
a 43
7 BV ie 3V/5
Example 3
Find the direction from the point CL, de 0) pee gives the greatest rate of
increase of the function ¢@= (x + 3y)? + (2y — ie
This appears to be different, but it rests on the fact that the greatest rate of
increase of ¢ with respect to distarice is in
I the direction of vo |
Vo = 4(2i+8j —k)
Because
Vo
N=
| Vo |
(1S)
Vector differentiation
249
Example 1
Find the unit normal vector to the surface xy + 4x2? +
xy?z +2 =0 at the
point (1, 3,—1).
Vector normal = V¢ =
1
SNe
Tas | Ee Se
j+k)
Because
| Vo |= V164 254+1=
V42
Vo 1
and N= = 4i —5j+k
|Vo | Ji A
One more.
Example 2
Determine the unit normal to the surface
Unit normal = N =
Because
o=xyz+x’y
— 5yz—5
Vb = (yz + 2xy)i + (xz + x” — 5z)j + (xy — Sy)k
At (3, 1,2), ‘Vé=8i+5j—2k; |Vo|= V64+ 25 +4 = V93
: _ Uniti norma Ne
Vo 1
v6] = 753 |8i + 5jJ — 2k)
250 Vector Analysis
dy + Fez
¢ where d¢= soax+ Se
(b) dé = dr- grad
d A
(c) directional derivative < =a-grad¢
lg OA OA. GA {2 OB.
|OBy.
={Rit SS hy Ox Toy} Oz \
~ V(A+B)=VA+VB
(b) V(AB) = if5-(4B)} + if (4B) ofkf5.aB)}
Al Ox
i
Ox
rifareey
Oy Oy
beklare
OZ
+B}
OZ
= 458 a5[Akt
OB
- hy {Bae@A. OA. bP_OA
rit Pol asK|
OB, 2B.
=e) ea j+ OB 0A, OA, @A
+ oy) az 2k \43) Si + Fitz}
‘. V(AB) = A(VB) + B(VA)
Remember that in these results A and B are scalars. The operator V acting on a
VCCLO Lan sataycee
has no meaning
Example
If A = x?yz + xz? and B = xy*z — z?, evaluate V(AB) at the point (2, 1, 3).
We know that V(AB) = A(VB) + B(VA)
Ata iea
Vector differentiation
251
oe ee B, OB ae :
Nalin he ay) +a K = y’zit 2xyzj + (xy — 327)k
= 3i+12j—25k at (2, 1, 3)
OA,
MA tre ilar
DAcwBA
ay it okie
pee:
(2xyz
+ 27)i + x?zj + (x2y+ 2xz)k
=211+12j+16k at (2, 1, 3)
Now V(AB) = A(VB) + B(VA) =
Finish it
Because
If A= ai+ a,j + ak
io. © O : :
dvA=V-A= (12 +15 +key) (ad+ai+ak)
Example 1
If A = x2yi — xyzj + yz7k then
Example 2
If A = 2x*yi — 2(xy? + y?z)j + 3y?z*k, determine V - A, ie. div A.
Complete it. V- Asset
Because
Such a vector A for which V - A = 0 at all points, i.e. for all values of x, y, z, is
called a solenoidal vector. It is rather a special case.
ign ok
ahes -O* a
OX, OV A0Z
a) Ge
Example 1
i j k
O 6) O
cur rA=VVe A = —
Dx =
dy —_
a7
yt Pea: x2 4 2 ~x2yz
All that now remains is to obtain the partial derivatives and substitute the
values of x, y, z.
2i — 8j — 106k
Example 2
Determine curl F at the point (2, 0, 3) given that
i j k
a a a
Now expand the determinant and substitute the values of x, y and z, finally
obtaining curl F=...:........
Lanne
een eener ee ee
254 Vector Analysis
Because
V x F=if{2— 2x cosy} — j{1 — e”} + k{2z cosy— 2xze””}
- At (2,0,3) MoxFAi@=
4/2401 SIeek(6 —12)
ei ok = 201 3k)
Every one is done in the same way.
Check through that list, just to make sure. We shall need them all
Example 3
lf @ = xy? +x yz ye"
and F = xy*i — 2yzj+xyzk
determine for the point P (1, — 1, 2),
(a) V¢, (b) unit normal, (c) V-F, (d) VxE
Complete all four parts and then check the results with the next frame
Vector differentiation
255
Multiple operations
We can combine the operators grad, div and curl in multiple operations, as in
the examples that follow.
Example 1
If A = x7yi + yz3j — 2x°7k
O 0 ; a
then divA=V-A= (Zits! +o) C7 yi + yz>j — zx°k)
=27+274+x°=¢ say
Od, 9b, ,OG
Then _ grad (div A) = V(V -A)
i — . —
yi t ay! + Az k
Example 2
If ¢ = xyz — 2y?z + x22”, determine div grad ¢ at the point (2, 4, 1).
First find grad ¢ and then the div of the result.
At (2, 4, 1), div grad 6=V - (Vo) =. 2-2
a
06,
abe
96,
a
,Ob
Example 3
If F = x*yzi + xyz?j + y?zk determine curl curl F at the point (2, 1, 1).
Determine an expression for curl F in the usual way, which will be a vector,
and then the curl of the result. Finally substitute values.
CurleCurieh— ee eee
curl curl F= V x (V x F) = i+ 2j + 6k
Because
nS ON Pays Boz
xeyz ixyZ y2Z
= (2yz — 2xyz)i+ x*yj + (yz? — x*z)k
i j k
Remember that grad, div and curl are operators and that they must
act on a
scalar or vector as appropriate. They cannot exist alone and must be
followed
by a function.
One or two interesting general results appear.
(a) Curl grad @ where ¢ is a scalar
Od. Od. O¢
oe d¢g=ain ye aoe
Ines lioy K
dy Gy az
=
. (Oa, —
Oday
——_— ==
SV =
Cyie Ss
elee: k
Oay —
>)
(> a) 1(> 7m)+ (F oy
Then en div
di cur LA ==¥,-(Vi
Niel 4A)vial= jee
Ga lay ke
Dz -(V
x A)
Review summary
WN
1 Differentiation of vectors
If A, dy, dy, az are functions of u
dA a da Xx i+ da Ve mi da Z
dus du du > du s
2 Unit tangent vector T
dA
a iis
3 Integration of vectors
b b b b
|Adu=i| a.du+j] a, du +k |a, du
a a a
grad dd=
06,
d= Vo - ax
06, j O¢ k
@yit
Ae io, a A
(a) Directional derivative a =a-graddé=a-Vo where a is a unit
Oa " Sse at
divA=V-A=
aa as Fs
If V -A=0 for all points, A is a solenoidal vector.
6 Curl (curl of a vector function A)
a Lg beh
arlA = Win 'A ele a ao
Ox Op Oz
ay dy ay
Check this list before and after you try the end of Program test.
2 Determine the unit tangent vector at the point (2, 4, 3) for the curve with
parametric equations
x=2u- y=u+3; z=4u* —u.
3 IfF=2i+4uj+u2k and G=v7i-—2uj+ 4k, determine
2
|(F x G)du.
0
4 Find the directional derivative of the function ¢ = x’y — 2xz* + yz at the point
(1, 3, 2) in the direction of the vector A = 3i+ 2j — kK.
5 Find the unit normal to the surface ¢ = 2x°z + x*y? + xyz — 4 = 0 at the point
CRO:
6 IfA=x?yi+ (xy +yz)j+xz2k; B= yzi— 3xzj + 2xyk; and
b = 3x*y + xyz — 4y7z? — 3; determine, at the point (1, 2, 1)
(a) Vd; (b) V-A; (c) V xB; (d) grad div A; (e) curl curl A.
Vector differentiation
261
Further problems 8
1 If F = xi + (3x + 2)j+sinxk, find
dE dF dF d
(C) a
integration
Learning outcomes
When you have completed this Program you will be able to:
e Evaluate the line integral of a scalar and a vector field in Cartesian
coordinates
Evaluate the volume integral of a vector field
Evaluate the surface integral of a scalar and a vector field
Determine whether or not a vector field is a conservative vector field
Apply Gauss’ divergence theorem
Apply Stokes’ theorem
Determine the direction of unit normal vectors to a surface
Apply Green’s theorem in the plane
263
264 Vector Analysis
We dealt in some detail with line, surface and volume integrals in an earlier
Program, when we approached the subject analytically. In many practical
problems, it is more convenient to express these integrals in vector form and
the methods often lead to more concise working.
Line integrals
Let a point P on the curve c joining A
and B be denoted by the position
vector r with respect to a fixed origin
OF
If Q is a neighboring point on the
curve with position vector r+ dr,
then PQ = dr.
The curve c can be divided up into
many (n) such small arcs, approximat-
ing to drj, drz, dr3... dry)... so that
AB = ss dr,
p=1
where dry is a vector representing the element of arc in both magnitude and
direction.
Scalar field
n
If a scalar field V exists for all points on the curve, then Ss V dr, with dr — 0,
p=1
defines the line integral of V along the curve c from A to B,
dr = idx +jdy+kdz
Example 1
dr =i3du+j4udu+k3u?du |
Because
omarotlix NAGS ONG!
veur,- ~ dy=4udu
z=u, .. dz=3u? du
Limits: A (0; 0; @)» corresponds to v= 14 0.008
B(3, 2, 1) cormesponds to US... acs
Finish it off
Because
1
Example 2
|V=6°+20t?; dr=i2tdt+j2dt+k dt
Because
A.(0,0,
5).= tO) eBudeas7y=
1 S2
2
|Vdr=| (6t? + 20t2)(42t +j2 +k) dt
Oke 10)
8
75 (4441 + 290j + 145k)
2
|Vdr=2 |{(6t* + 2083)i + (6t3 + 2082)j + (3t3 + 10t2)k} dt
JC JO
The actual integration is simple enough and gives the result shown. All line
integrals in scalar fields are done in the same way.
nn nn nnn
SS SSS
Vector integration
267
Vector field
If a vector field F exists for all points of
the curve c, then for each element of
arc we can form the scalar product
F- dr. Summing these products for all
n
=| Far
Cc
In this case, since F- dr is a scalar product, then the line integral is a scalar.
To evaluate the line integral, F and dr are expressed in terms of x, y, z and
the curve in parametric form. We have
F= F,i+ F,j+ Fk
and dr=idx+jdy+kdz
Then F.- dr= (F,i+ F,j+ Fk) - (idx+jdy+kdz)
= F, dx + F, dy + F, dz
Example 1
Because
xy =4(1607)(2t7) = 32t- X= At -, dx =4dt
xz = (4t)(t?) = 4t* peel eS y= 4rd
2yz =(4P)(P) = 4P z= dz =3t2 dt
Then |F-dr= [sae +4t4j — 40k) -(i4dt +j4tdt + k3t* dt)
1
|F.dr= |CL2S eG =n at eer teiece
c 0)
If the vector field F is a force field, then the line integral |F-dr
Cc
represents the work done in moving a unit particle along the prescribed curve c
from A to B.
Now for another example.
Example 2
If F = x*yi + 2yzj + 3z?xk, evaluate | F- dr between A (O, O, O) and B (1, 2, 3)
(a) along the straight lines c, from (0, O, 0) to (1, 0, 0)
then > Trom(1, 0; 0) to (1, 2, 0)
and” cs froma pz, 0) ton 245)
(b) along the straight line cy joining (0, 0, 0) to (1, 2, 3).
Because
F- dr = (x*yi+ 2yzj+3z?xk)-(idx+j dy+k dz)
|F-dr = [Pyar + |2yzay 4 |32°x dz
| F-dr=07070=0
JC]
| F-dr=0+0+0=0
JC2
Because
3
| F-de=0+0+| Boge 27
Cr (0)
Because
F= 204+ 12t?j + 27°k
dr =idx+jdy+kdz=idt+j2dt+k3dt
a
| F-dr =| (2034 + 1207j + 27P°k) - (i+ 2 + 3k) dt
C4 0
1 1
= |(2t3 + 24t? + 81t?) dt = |(83t? + 24t?) dt
0 0
t4 a ae
= |83—4+8f|
va i=——=28.
= 8.75
So the value of the line integral depends on the path taken between the two
end points A and B
Example 3
If F=x’y?i+y3zj+ 22k, evaluate |F-dr along the curve x= 2v?,
Cc
.> 3u, z= uw’ between A (2, —3, — 1) and B (2, 3, 1). Proceed as before. You
will have no difficulty.
500
| FE: dr
r = —_
7] = ZO: 23. 7
point of V and on its boundary surface $, then FdV is the volume integral
JV
Example 1
Because
4 8 ves
3FdV = | yt +2924 — Fy] dz
ZO
oO
bh
= (91
+ 6zj — 8k)dz
ial
=!
4
= oi + 327} — Sz |
0
= 361 + 48j — 32k
= 4(9i + 12j — 8k)
Now another.
Vector integration
273
Example 2
So 2x +y+z=2 cuts the axes at A (1, 0, 0); B (0, 2, 0); C (0, 0, 2).
| FdV rR
Wl
=2(2i +k)
V
Because
= Gi x)*4 Fl nd =
Example 3
X == pcos o; y=psing
Rie dV = pdpdddz
FdV =
V
Vector integration
275
72n(i + 2))
Because
2m 3 3
I FdV =4 | i + 2p7§ + Fsin ok de
JO 10)
Dan
= 4| (91+ 18) + 9singk) dd
0
20
= 36 | (i+2j + sin¢gk)
dd
10]
2n
= 36 [oi
+ 2¢j — cos ok
0
= 36{(27i + 4nj —k) — (-k)}
= 72n(i+ 2j)
You will, of course, remember that in appropriate cases, the use of cylindrical
polar coordinates or spherical polar coordinates often simplifies the
subsequent calculations. So keep them in mind.
Now let us turn to surface integrals — in the next frame
Surface integrals
a 5
ABA If 9=>5) then |Ax B|/=AB in the
direction of the normal. Therefore, if n is
B a unit normai then
A Ax B=|A||Bjn=ABn
276 Vector Analysis
Scalar fields
Example 1
XyZ(xi + yj) dS
s
KcJase
yZi
Dee
+ xy“zj ) dS (1)
Vector integration
277
x=2c0s¢; y=2sind
Ihme fe qs =Zd@dz
nm/2 p3
|vas =| |(8cos? sin ozi + 8cos sin? 624) 2dzds
JS QO JO
m/2 3
33 |(cos? ¢sin di + cos ¢dsin? dj) 2zdzd¢
10)
n/2
= 4| (cos? ¢sin i+ cos sin” dj) 9d¢
0
and this eventually gives
278 Vector Analysis
|V dS = 12(i+))
S
Because
3 Aue} 1/2
Example 2
A scalar field V =x +y +z exists over the surface S defined by 2x + 2y+z=2
bounded by x = 0, y = 0, z= 0 in the first octant.
S$: 2x 2y FZ =z
bat PAZ 2y
cet
[Vél=V4444+1=V9=3
Therefore
~ Wo 2i+2j+k ‘ toss. Deke
Male 3 so that dS = nds = = (2i + 2j +k) ds
2
3 (21+ 2)+k)
Because
1 1—x
Vector fields
Example 1
A vector field F = yi + 2j +k exists over a surface S defined by x? + y* +z? =9
bounded by x = 0, y= 0, z = 0 in the first octant. Evaluate |F.- dS over the
S
surface indicated.
(PE dares
|F as = |Fonds =|(yi+2) +k) g(xd + yi + 2k) ds
a S &
=5| (xy + 2y + z) dS
3 Js
Before integrating over the surface, we convert to spherical polar coordinates.
emg See sie ; Pie ooo
280 Vector Analysis
Because
m/2
|F. as=9| (2sinocoso + Fsino +5 dé
§ 0 2 2
m/2
: 3r
= O [sin
keine p 7 COS poir
I. es 9(14 rab
7)
Example 2
4
—
\Vo|
45h 0 le
Vector integration
281
Because
x=2cos¢; y=2sind
oe. dS = 2d¢dz
Because
) “g
|Fas—4| | (1 — cos 2¢) dd dz
BS Z=0 J¢=0
5
= 4| rdz=4n|z| = 207
10)
Example 3
Evaluate FE. dS where F is the field x*i— yj +2zk and S is the surface
5
2x +y + 2z = 2 bounded by x = 0, y = 0, z = 0 in the first octant.
We can sketch the diagram by putting x =0, y=0, z=0 in turn in the
equation for S.
When 50)
F=x7i-yj+2zk; ¢: 2x+y+2z—2=0
0d. O06, O¢ comes .
Vo Bx ay! ae i+j+ |Vo|
Because
1
|.[52x —y + 4z)
dx dy
Because
|Fens =
S
2 3
[2x —y + 4z) 5 dxdy
R
[_(2x? - y + 42) axay
R
Limits: y=Otoy=2-—2x; x=Otox=1
r 1 1 p2—2x
|Faas —>] | (2x — y + 4z) dy dx
S 0 J0
In general, the value of the line integral |F- dr between two stated points
c
| F-dr= | F-dr
c, (AB) C2(AB)
es pF-dr=0
Note that this result holds good only for a closed curve and when the vector
field is a conservative field.
Now for an example.
Example
B (2, 4, 6)
Because
2
[Fear =| 15u4 du = [3u°],= 96 |Far =96
(e 0) c
(b) The diagram for (b) is as shown. We consider each straight line section in
turn.
; B (2,4,6)
C3
A
‘ yay
1 & /
x 2
| F-dr=0+0+0=0
Cy
In the same way, we evaluate the line integral along c2 and c3.
| F. dr =0; | F. dr = 96
C2 C3
| F-dr=0+0+0=0
C2
| F-dr=0
C2
C3), (4;,.4,0)
40 (2,4,6);
x=2, Pond, dk =0, op=
6 6
| F-dr=0+0+ | léedz= 162| = 96
(ey 0 0
| F. dr = 96
Two tests can be applied to establish that a given vector field is conservative.
If F is a conservative field
(a)Pcure kr [0
(b) F can be expressed as grad V where V is a scalar field to be determined.
For example, in the work we have just completed, we showed that
F = 2xyzi + x°zj + x’vk is a conservative field.
(a) If we determine curl F in this case, we have
Because
i j k
ie
ie ah)Aa oe
moll
Oy Oz
2XVE ORL MV
Sx = x7) (2x 2xy)j + (xz = 227)k =.0
| tel r= >
(b) We can attempt to express F as grad V where V is a scalar in x, y, Z.
lft Vf (x; y, 2)
7 QVeuoV.-
, , ov ie
grad V rial ay * Dz
Because
Exercise
Divergence theorem
(Gauss’ theorem)
For a closed surface S, enclosing a region
V in a vector field F,
In general, this means that the volume integral (triple integral) on the left-
hand side can be expressed as a surface integral (double integral) on the right-
hand side. Let us work through one or two examples.
Example 1
Verify the divergence theorem for the vector field F = x7i + zj + yk taken over
the tevion bounded bythe planessz— Oz = 2, X= 0, 1, v= 0, vy = 3.
Start off, as always, by sketching the relevant diagram, which is
dV = dxdydz
We have to show that
| div FdV =6
V
Because
1 p3 (2 1 73 2
| div Fav =| ||2xdzdydx = ||axe dy dx
V 0 JO JO JO JO 0
1 3 1 1
=| aay dx=| 12x dx = 6x2 =6
10) 6) 0 @)
F=17i+zj+yk
(1) S; (base): z=0; m=-—k (outwards and downwards)
FE
I x i+yk
Zz
dS, =dxdy
| F-nds =| (x°i + yk) - (—k) dy dx
Sy
(2) Sz (top): z= 2;
Vector integration
291
Because
I,F-ndS
n = |[02+
25 j
2+ ve) Cc) ay ax
3
9
Bydydx
Y PRS =—
JO
SO we go on.
Because
; i A
| F-nds= || (i + aj + yk) -(-j)dzdx =| |(-2azax
Sy IS« 0
1 7 2 1
=| 5 dx=| (2) dx = -2
(6) 2 (6) 10)
Check:
(S) os (nOnt)s ket, n=i dS; = dy dz
Now on to the next frame where we will collect our results together
[; B.ds=-34542-246+0=6
a3
and from our previous work in section (a) | div FdV = 6
V
Example 2
Verify the Gauss divergence theorem for the vector field F = xi + 2j + z7k
taken over the region bounded by the planes z = 0, z= 4, x = 0, y =O and the
surface x* + y” = 4 in the first octant.
Divergence theorem
oei+ A dhe
ok) Gd +2) 427k)
te s
Vector integration 293
Finish it
Because
/
n/2 2 4
div Fav =| ||(1 +22) dzpdpde
JV 0 0 JO
xf2 »2(2 4 /2 2
=| | e+ 2 papas = | |20p dp d@
0) 0 0 0 0
n/2 2 m/2
~ | 10/?| dé = | 40 d¢ = 20x (1)
0 0 0
F=xi+ 2j+2°k
294 Vector Analysis
=16() = 16"
Because we have
(3)53° y=0; Nesp Fox
2j4 7k
Finally we have
(5) Sg: Sx? sey? A ee
Because
47+ 16
Because we have
. 1
| F-nds—>| (x? +:2y) ds
Ss yi Ss
Stokes’ theorem
|curt Fd =} F- dr
Cc s c
This means that we can express a surface integral in terms of a line integral
.
round the boundary curve.
The proof of this theorem is rather lengthy and is to be found in the
Let us
Appendix of the authors’ Advanced Engineering Mathematics (4th edition).
bb
demonstrate its application in the following examples.
296 Vector Analysis
Example 1
A hemisphere S is defined by x7+y?+z?=4(z>0). A vector field
F = 2yi —xj+xzk exists over the surface and around its boundary c.
S: x+y? 4+27-4=0
F= 2yi — xj +xzk
c is the circle x7 + y* = 4.
i (2y dx — x dy + xz dz)
Because
20
}F-dr =| (4 sin 6[—2 sin 6 dé] — 2cos #2 cos 6 dé)
€ 0
2
=~4| (2 sin? 6 + cos? @) dé
10)
2r 27
- -4| (1 + sin?
6) dé = -2| (3 — cos 26) dé
0 0
sin 26]7"
5) bs 127 (1)
curl F = —zj — 3k
Because
ini ek
ip—|2 OY aj. .
cur ax by oz 1(0 — 0) —j(z—0) + k(-1- 2) = -gj: -3k
2y —X XZ
| curl F-ndS =
S
Because
| curl Fd =} F-.dr
Sj Cc
Before we proceed with another example, let us clarify a point relating to the
direction of unit normal vectors now that we are dealing with surfaces.
So on to the next frame
298 Vector Analysis
ni
c (é
adeCc A A
A unit normal n is drawn perpendicular to the surface S$ at any point in the
direction indicated by applying a right-handed screw sense to the direction of
integration round the boundary c.
Having noted that point, we can now deal with the next example.
Example 2
A surface consists of five sections formed by the planes x =0, x= 1, y=0O,
y = 3, z=2 in the first octant. If the vector field F = yi + z*j + xvk exists over
the surface and around its boundary, verify Stokes’ theorem.
| Far=[(0+0+0)=0
Cy ‘
(2) Along
cz: x=1; z=0; dx=0; dz=0
F-dr= [(0+040) =0
Because
(3)-Along32 y= 1 = 08 dre dz 20
0 0)
| F-dr =| (3dx+040) = [33] =-—3
C3 1 1
: } F-dr=0+0-3+0=-3
Cc
Because
i Oj k
O00
curl F= V x F=|7" ay dz
Menta: XY
=i(x —2z) -j(v—0)+k(O- 1) = (x- 2z)i —yj —k
Because
Then, likewise
=| (as
S2
aS,
Because
= | y ds
JS3
But y = 0 over S3
| curl F-ndS=0
S3
Because
(4)(S2.(iront): |n=
= | (x — 2z) dS
Ss
Buty = 1 oveto.
3 p2 3 Z
=| (—2) dy = |-2y| = —6
0 0
(5) Ss (back): m= —i with x = 0 over Ss
So, referring back to our result for section (a) we see that
Example 3
A surface S consists of that part of the cylinder x” + y? = 9 between z = 0 and
z=4 for y>0O and the two semicircles of radius 3 in the planes z=0
and z = 4. If F= zi+ xyj + xzk, evaluate | curl F- dS over the surface.
S
The surface S consists of three sections
(a) the curved surface of the cylinder
(b) the top and bottom semicircles.
We could therefore evaluate
| curl F- dS
§
-—
—
——
\
——
over each of these separately.
However, we know by Stokes’ theorem that
| curlih.dS= -.. oe
S
Now we can work through this easily enough, taking c), C2, ¢3, C4 in turn, and
summing the results, which gives
curl F-d8 =} F- dr =
Js JC
Vector integration
303
—24
| F-dr=| (0+0+0)=0
cy
3721"
| F-dr=| (0+0~32dz) = | =| amg
C2 C2 2 10)
322]°
| F-dr=| (040+ 32dz)= |] = —24
C4 C4 2 4
pF-dr=0-24+24~24=-24
i
This working is a good deal easier than calculating | curl F- dS over the three
S
separate surfaces direct.
So, if you have not already done so, make a note of Stokes’ theorem:
Green's theorem
[rasa
ff(2) aw
where the line integral is taken round c in an anticlockwise manner.
In vector terms, this becomes:
y S is a two-dimensional space
dS=dxdy enclosed by a simple closed curve c.
dS = dx dy
rn c dS = ndS =k dxdy
x
Because
LAE ER ees) JQ aP
=i( =) i(0 eye inf see |
O
OQ OP OQ OP
; curl F-d8 = |ie
(22-27). acyas =| { y
|
s s \Ox oy Ge JJs\Ox dy ean
0Q’ OP
.
ra
ipoghky =
tee = ou ;
(1)
Now by stokes’ theoreni.7.........
eee
Vector integration
305
curl F- dS =} F- dr
ls ; Cc
= |(P dx + Qdy)
OQ OP
Green’s theorem || (3 ~ er)dydy= iG dx + Qdy)
S
Example
py = O<x<2 ce
x=0- O<y<2
fe)
C4 Z “:
aQ aP |
Green’s, theorem: : IG——
=,4 Urn
x dy =o i(Paxr+Ood
Qdy)
ee Make any necessary substitutions and evaluate the line integral for Co.
306 Vector Analysis
nm—4
Because we have
Putting
x =2sin6, V4—x2=2cos@ dx=2cosédé
Limits: x = 2, = 5 20 C=O.
2 4sin?6
| wars aay) = [4 = 4siné
4sin d= wal 2 cos
cos 6 dé
1 sin 20) ]°
=4]2sing sin® 6 5 C 5) ihion
Finally
=s[-(@-1-9) =*-5 16
(3) cs: 4 = 0 dx 0
0 0
C3 J2 2
{[ (C22)axay= [fanaa
Q=x+2y
= 2 Fe
ae
——
1
OP x |
OQ
|J.(ae~ay)
w/2 p2
II |(1 — 2rsin 6)r drd@
Jo Jo
a Complete it
Vector integration
307
Here it is:
po 8, n/2 (2 is
IAF ~ Oy, dxdy = ii Kc — 2r° sin@) dr dé
1/2
= | {2— sina} dé
0 3
4.
|
0
=Tt-
uP3
—_—
(2)
So we have established once again that
And that brings us to the end of this particular Program. We have covered a
number of important sections, so check carefully down the Review
summary and the Can You? checklist, and then work through the Test
exercise that follows. The Further problems provide valuable additional
practice.
5|Review summary
1 Line integrals
F=F,i+ F,j+ 2k
dr =idx+jdy+kdz
2 Volume integrals
F is a vector field; V a closed region with boundary surface S.
X2 py2 ¢Z2
| Fav =| | | Fdzdy dx
V Xy Jy FZ)
308 Vector Analysis
A= (COS0, y.Tsing
fo ay Chee
mrelmele,
|
y=psing
Z Za
iy 7 dS = pdddz
o p dV = pdpdddz
Xx
P (r, 8, 0) x=rsinécos¢@
y=rsinésing
Za COSY
dS = 1? sin@déd¢
dV =r’siné drdéd¢
(b) curl F= 0
(c) F= grad V where V is a scalar
Vector integration
309
7 Stokes’ theorem
| curl Fd = fF-dr
S (e
xX
8 Green’s theorem
y Thecurve~c~is~a~simple=closed
curve enclosing a plane space S$ in
the x-y plane. P and Q are func-
tions of both x and y.
x
(4 Can You?
Checklist 9
Check this list before and after you try the end of Program test.
Test exercise 9
3 Evaluate to one decimal place the integral |(xyz + 4x’y) dr along the curve c
Cc
where c is the path with parametric equations x = 2u, y = Fang Wey Heh I
6 A unit particle is moved in an anticlockwise manner round a circle with center
(0,0, 4) and radius 2 in the plane z=4 in a force field defined as
F = (xy + z)i+ (2x+y)j + (x+y +2)K. Find the work done.
312 Vector Analysis
7 Evaluate | FEdV where F = i— yj + k and Vis the region bounded by the plane
V
z= 0 and the hemisphere x? + y? + z* = 4, for z > 0.
8 V is the region bounded by the planes x =0, y=0, z=0 and the surfaces
y=4-x? (z>0) andy=4-2? (y>0).
If F= 2i+y’j —k, evaluate | FdV throughout the region.
V
Curvilinear
coordinates
Learning outcomes
When you have completed this Program you will be able to:
e Derive the family of curves of constant coordinates for curvilinear
coordinates
e Derive unit base vectors and scale factors in orthogonal curvilinear
coordinates
e Obtain the element of arc ds and the element of volume dV in orthogonal
curvilinear coordinates
e Obtain expressions for the operators grad, div and curl in orthogonal
curvilinear coordinates
314
Curvilinear coordinates 315
This short Program is an extension of the two previous ones and may not be
required for all courses. It can well be bypassed without adversely affecting the
rest of the work.
Curvilinear coordinates
Let us consider two variables u u=a
and v, each of which is a func-
tion of x and y
ie. m= (x, y)
v= 8(x, y)
a, u=constant
by ~~ v=constant
O x x
A pair of curves u = a; and v = b, pass through each point in the plane. Hence,
or in
any point in the plane can be expressed in rectangular coordinates (x, y)
curvilinear coordinates (u, V).
Let us see how this works out in an example, so move on
316 Vector Analysis
Example 1
Let us consider the case where u = xy and v = x” —y.
4 :
(a) With u=xy, if we put u=4, then y= z and we can plot y against x
to obtain the relevant curve.
Similarly, putting u = 8, 16, 32, ... we can build up a family of curves,
all of the pattern u = xy.
3 4 |]
1 8
5 12
9 16
13 20
17 24
21 28
25 32
Note again that we label each graph with its own v-value.
This again is a family of curves with the common pattern v = x* — y, the
members being distinguished from each other by the value assigned to v in
each case.
Now we draw both sets of curves on a common set of x-y axes, taking
the range of x from x =0.tox=4
and the range of y from y = 0 to y = 20.
It is worthwhile taking a little time over it - and good practice!
When you have the complete picture, move on to the next frame
i
318 Vector Analysis
The position of any point in the plane can now be stated in two ways. For
example, the point P has Cartesian rectangular coordinates x = 2, y = 8. It can
also be stated in curvilinear coordinates u = 16, v = —4, for it is at the point of
intersection of the two curves corresponding to u = 16 and v= —4.
Likewise, for the point Q, the position in rectangular coordinates is
x = 2.65, y=5.0 and for its position in curvilinear coordinates we must
estimate it within the network. Approximate values are u = 13, v= 2.
Similarly, the curvilinear coordinates of R (x= 1.8, y= 14) are approxi-
mately
7.0 WV = I
Example 2
If u = x” + 2y and v= y— (x +1)’, these can be rewritten as y =4(u—x?) and
y=v+(x+1)*. We can now plot the family of curves, say between x = 0 and
x=4, with u=5(5)30 and v= —20(5)5, ie. values of u from 5 to 30 at
intervals of 5 units and values of v from —20 to 5 at intervals of 5 units.
The resulting network is easily obtained and appears as
Curvilinear coordinates
319
Example 3
Given the curvilinear coordinates u and v where u = xy and v = x* — y” then
u and v form a coordinate system that is ............
Ee
320 Vector Analysis
Because
) 0 OV
u = xy so Clas and ee v=x*%—y" so a = 2x and ees Then
OX oy
BAD ae Cuey = 2xy — 2xy = 0 and so uw and v form a coordinate system that
OxOx Oy doy
is orthogonal.
Example 4
Given the curvilinear coordinates u and v where u = x? + 2y and
vay—(x+1)? then
u and v form a coordinate system that is ............
not orthogonal
Because
ana eee he
oy
Then
Ou OV snOudov
4x(x+1)+2+#0 and so uw and v form a coordinate
Ox Ox Oydy —
system that is not orthogonal.
Let us extend these ideas to three dimensions. Move on
Situations can arise, however, where the directions of the unit vectors do
not remain fixed, but vary from point to point in space according to
prescribed conditions. Examples of this occur in cylindrical and spherical
polar coordinates, with which we are already familiar.
_43or ot
dp|
ee Or /\or
Oz
Note that I, J, K are mutually perpendicular and form a right-handed set. But
note also that, unlike the unit vectors i, j, kK in the Cartesian system, the unit
vectors I, J, K, or base vectors as they are called, are not fixed in directions, but
change as the position of P changes.
So we have, for cylindrical polar coordinates
OF /|or
~ Op/ \dp
Or /\or
Je)
Oo / |0d
Or /\or
ayes
If F is a vector associated with P, then F(r) = F,I + F,J + F,K where F,, F,, F,
are the components of F in the directions of the unit base vectors I, J, K.
Now let us attend to spherical coordinates in the same way.
Curvilinear coordinates
323
I- Or /\or
~ Or / \6r
(b)wa If r and # remain constant and 6
increases, P will move along the
‘meridian’ through P, i.e. = is a
tangent vector to this circle at P
and the unit vectorJ is given by
Pesor
Once again, the three unit vectors at P (base vectors) are mutually
perpendicular and form a right-handed set. Their directions in space, however,
change as the position of P changes.
A vector F associated with P can therefore be expressed as
F=F,1+ FoJ + FjK where F,, Fg, Fy are the components of F in the directions
of the base vectors I, J, K
Both cylindrical and spherical polar coordinate systems are
orthogonal
Scale factors
Collecting the recent results together, we have:
pe 1 Or h apes Or
-5 dz| hy oz Sens ee
2 For spherical polar coordinates, the unit base vectors are
Or /\or 1 Or or
. =f Or|. Ty OF TLL aide
Or /\or 1 Or or
= By 56 = hy 00 where he = 30
_ or Or| 1 Or Or
Od a¢| hy 06 where he — ad
r=xi+yj+
zk
Ls) COS.
y=psing
ZZ
oi, 1; he So ig |
P (r,
8, 0)
r= xi+yj+
zk
x=rsinécos¢
Z y=rsiné@sin¢
Z y Z=rcos0
‘ r=rsinédcos¢i+rsinésingj+rcosék
Then working as before
326 Vector Analysis
heals =e hg =rsin@
Because
r=rsinédcos¢i+rsinésing j+rcosék
Or or “4 ior
r= / Or ih, Or
h, =|F|=1 sin@cos¢i+sind@sin
¢ j + cos é k |
2% /|ar|_ Lar
~ 060/ |00| he 00
or * : 3 :
hg = 36| = |r cosé cos¢i+r cosé singj—r sind k|
K-22 Or; 1 or
— O¢/ |0¢| hg Ad
or , : Fit abe A
hg = a6 =|-rsinésingi+rsinécos
¢j |
or 10r 1 Or
=—: = : ‘ K = - se
r 00 rsin@0@
SS
Curvilinear coordinates 327
orthogonal
Exercise
Determine the unit base vectors in the directions of the following vectors and
determine whether the vectors are orthogonal.
1 i-2j+ 4k 2 2i — 3j + 2k
2i+3j+ k i 2] 2k
—2i+ j+ k — 10i — 2j+ 7k
3 4i1+2j- k 4 314+2j+ k
3i — Sj+ 2k 13] faak
i+ 2j+ 6k 6i+ j-— k
1 bars :
2 GGn
1=——(2i-— 3j+2k);
J J==(i+
3 2j+ 2k);
1
K = —— (-—10i + 2j + 7k)
Vis3 |
I.J=0; J-K=0; K-I=0_ .. orthogonal »
328 Vector Analysis
it 1
3
wo
Il=—=~(4i+ 2j-k
ee aR —— (3i — Sj + 2k);
it
K = —= (i+ 2j + 6k
Jal | ;
Ij=0-— J 2K 0 “. not orthogonal
Transformation equations
where the functions f, g, h are continuous and single-valued, and whose partial
derivatives are continuous.
Then
r=xi+yj+zk=f(u, v, w)i+g(u, v, w)j+h(u, v, w)k
and coordinate curves can be formed by keeping two of the three variables
constant.
A : ; _ Or or ? or
Now r = xi+ yj + zk 4 dr =~ du+— di +2 dw (1)
Or. ;
a, 84 tangent vector to the u-coordinate curve at P
Or ,
am is a tangent vector to the v-coordinate curve at P
Ore
aw 8a tangent vector to the w-coordinate curve at P
1 ,0r<7\0E or
micy, au au = ul where h, lai
Or /\or or
=z / my a hy J where hy=
or Or Or he
aw /
K = — | |—
law ties
Wise hy hy K where h, ,=|—
fe
|
O(x, y, Z)
where au, v, w) is the Jacobian of the transformation.
———$§
er
eee
330 Vector Analysis
Let a scalar field V exist in space and let dV be the change in V from P to Q.
If the position vector of P is r then that of Q is r+dr.
OV OV OV
Then dV =A bt dw
or Or (|or
B t a =| |= By, L v ,
* : Ov== hy J
ne
and awn atIK= hy K.
oo = (WV yh Y (WV y=
c= (VV yh (VV), = an
iv Jo KoO
i.e. grad operator V =
h,0u h,Ov h,ydOw
Example 1
Show that the curvilinear expression for divF agrees with the earlier
definition in Cartesian coordinates.
In Cartesian coordinates x, y, z we have h, =h,=h,=............ so that
Civ: Fy te eas
Example 2
Show that the curvilinear expression for curl F agrees with the earlier
definition in Cartesian coordinates.
In Cartesian coordinates x, y, z we have hy = hy, =h, =............ and I, J,
Kien ee oO tial
332 Vector Analysis
Because in Cartesians
h,=h,=h,=1 and Ij, K=i,j,k so that
hl) hh, ak
i) a) O O
Vex FE
hyhtyhy| Ou Woy, Ow
HE pe Ee Wel a
Lye iiak
Ee ees
Ox Oy OZ
F, Fy F,
OF, OF,\ .(OF, OF; OF, =»
(5 a) +i(F x) + (F Oy
= 5 | SS Saeee —— ———— k ———
Example 3
Show that the curvilinear expression for V7V agrees with the earlier definition
in Cartesian coordinates.
In Cartesian coordinates x, y, z we have h, =h, =h,=......... so that
VeVi. ok.
h, = h, = h, = 1 so that
2 V= OV OV ley
¥ x2 * Oy? + O22
Let’s try another example, this time in coordinates other than Cartesians.
Example 4
If V(u, v,w) =u + v* + w® with scale factors h, = 2, h, = 1, hy = 1, find V2V at
the point (5, 3, 4).
There is very little to it. All we have to do is to determine the various partial
derivatives and substitute in the expression above with relevant values.
div grad V=
ER eee
Curvilinear coordinates
333
Because
vy - 1 y hy ONE ne hy, pe geX n Omi fig OV
Kiri ou My Ou Ov\ hy, Ov Ow\ hy Ow
OV OV OV
In this case, V=u+v24+w> 5
Ou Teas © ce
Also’ h, =2, hy =1, hy
Wf ve at O O
WV= :—(6w"
2 ‘a (3)ape ayo” )}
=${0+4+12w}
. Atw=4, V2V=26
That is all there is to it. Here is another.
Example 5
If V = (uv? + v*)w> with h, = 3, hy = 1, hy = 2, find div grad V at the point
(2, - 2, 1).
Because
OV V
Ve == US 2 EV 2 aw 3 oeeS Du uw”;fees
a = 2vw?; ——
ISOM oe ae ae
‘i 1(a/20V\
2 =a/.aV\. a (3aV paces es
3 v=2\= (sor) tay (Sor) tam (aan) =
cal O (4
3 3 0
=o han(gum) + 5p (1207) + a5 (5 2 2
5 (ul + v*)w a
at (2, —2, 1)
prs F;,
Bea lee) ina ene alA
Op pdp |page az
OV Ae Lav InoV
tad V = —I
5 Or TT 50) rind Od
diviv F= Zand t 6 F,)+ taal
ae sin
oO
alt sin 6 Fa) + 5p (tho) }
I rr rsindK
1 OO 0
curl R= :
es rsin6\Or 06 Oo
F, rFg rsindF,
way VY ,2OV 1 eV cotaV 1 &v
Or2 = r Or i r2 062 iu r 00 i:r2 sin 6 Od2
The Review summary follows as usual. Make any further notes as necessary:
then you can work through the Can You? checklist and the Test exercise
without difficulty. The Program ends with the usual Further problems.
Review summary
1. Curvilinear coordinates in two dimensions
=f (x,y); v=s(%,y)
2 Orthogonal coordinate system in space
(a) Cartesian rectangular coordinates (x, y, Z)
F=F,i+F,j+F,k Scale factors hy =h, =h,; = 1
=> %
ig ah
Op
a
or or
J= ad o h,d =|=—
Aad p)
Or /\or or
= — = h, =|—|=
ie Oz / \Oz erty
F=F,1+F,J+F,K
336 Vector Analysis
r= /- hy =|>|=1
8/8) tele
K-= . »=[5,[=rsing
Or or
or or
ane hyK where hy, = =|
SOX V2)
= alu, v, w) dudvdw
hl hy hyK
(c) Curl F= : a a 2
hyhyhy| du Ov Ow
fabalaly Iabsliby (ph dP.
dive = {te} gi [= | ob
p Op p | 0d OZ
I pj K
Cot: Ks a a !
pi\Oop Od Oz
Py ers: FP,
AT Vy
ae 19V
SA 92ot as By)aa
Ns
Bee oe ae ae
(b) Spherical coordinates (r, 6, @)
OV 10V 1aey.
rad V = — [+—— =
6 Or Hr ap Fin he
div: 10
Bae {7° F) siO lo
ne @
CA areas aint
pe 2 ale reineoes”
I rr rsinéK
1 Bs ¥2 O
DF=
Age rsin@\Or 06 Od
F, Fg rsinéF¢
yaad ak 20V 10@V_ cotd9dV 1 Oty
Or csr Or Te 0e2 r2 00 Loe sin?
@ O¢2
Check this list before and after you try the end of Program test
6 If Visa scalar field such that V = u?vw? and scale factors are h, = 1, h, = 2,
hy =4, determine VV at the point (2, 3, — 1).
Curvilinear coordinates 339
F, = F, cos @ — F,sin@
5 If V isa scalar field, determine an expression for V?V
(a) in cylindrical polar coordinates
(b) in spherical polar coordinates
6 Transformation equations from rectangular coordinates (x, y, z) to parabolic
cylindrical coordinates (u, v, w) are
tem ye
bi : voy, 2= Ww
2
V is a scalar field and F a vector field.
(a) Prove that the (u, v, w) system is orthogonal
(b) Determine the scale factors
(c) Find div F
(d) Obtain an expression for V*V
a
Answers
1 ue 12x" saga
(a) nate — Sy ey —10xy + 9y pe iees24x
oe Fe
ay =p x pets y
wes = —10y ie =-10y (b) oe = —2sin(2x + 3y) = = —3 sin(2x + 3y)
Oy - OX Ox - OY Ox oy
Oz Oz Pz
—=-
aya 4 cos(2x + 3y) —=-
ay 9 cos(2x + 3y) ——
By - Ox = cos(—6c0s(2x+3y)
dySots
Ox —6x? sin(2x + 3y) + 6x cos(2x + 3y)
s ”
facie cos(2x + 3y)
shad
aS 2 si
sin(2x OAT)
+ 3y) oe arldu a.
rare x sin(2x + Sys3y)+ Gace
6x cos(2x + 3y)
340
Program title 341
O-
(d) r=1+cosé
D 2
Rose aes —s
37a
9 "10 PYG
41 21.250 12 37a
14 *a {vb a Fi bo; Se
erryAa, | )
15 ra2(2— V2)
15 M=| = 2 2 he
Zandteds dG ee
14 $4 32-4) 0 Jo 3 Or 2
9
b? 2 2 [>
(b) sabe: centroid ji 17 19.56 18 —(c—a‘*) 19 elt + 3V3)
8 37 4 6
20 232
3 4nab 3
15 “ {4v2 - 3} 16 — 175, 18 z|| +) duav
19 wvdudvdw 20 z=—2 21 23 *2 670 1)
(c) 13(i
+2j+k) @) a+ 2+8) 14 61 15 291-10)
6 A A
fick
16 (a) 22i+14j+2k (b) —2i+14j— 22k
17 (a) 15° (b) -33 ().7 IB (2) -6i+4j-7k "(b) G2110) — 38k
(c) 18i-21j+10k 19 p=6 20 (a) (1) p=15/4 (2) p= —33
(b) 5(3i — 2j + 6k)
=
(oil Ok
mig (121iwat!
eG —— + 4j — 21k) (a SiS oe —o. 165
Vise J ) V601
344 Vector Analysis
5 (a) V*V=
Aan OV 1 &V &vV
Op? p Op 3 p2 Od? * Oz?
b) VV
1107 30V i) aV 1 &V
(b) 2 Ss: ar («
4
) * sind a0 (sino) *sin26 Oe
6 (b) hy, — hy — VV u2 =F v2; hy, = 1
:
(OuGivare= aan tau (Ve ee 0 iF.
ma) tal ie aye a) + OF
Ou Ow
217 1 OV OV O2V
oD) ag u2 + v2 | du2 = Ov2 t Ow2
Se ee
Index
345
344 346 Index
Distinctive Features
© Uses a unique and proven “programmed” approach that takes you through
each topic in a step-by-step fashion with a wealth of worked examples and exercises.
e Serves as a "personal tutor", allowing you to proceed at your own rate of learning —
and master each topic before going on to the next.
® Provides many opportunities for self testing through numerous Quizzes, Learning
Outcomes, and Can You? checklists.
© Equally effective as a classroom textbook or a self-study tutorial.
q 90000
— ae , _— i q }
9"780831'132088 |