Mathematical Methods for Integration and Analysis
Mathematical Methods for Integration and Analysis
S T E M L E R ,
L . S T O YA N O V
M AT H E M AT I C A L
METHODS 2
2 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
5 Vector fields 71
5.1 Path integrals of vector fields in two dimensions 72
5.2 Flux through a surface in R3 76
5.3 Work done by the gradient of a scalar field 77
7 Fourier series 91
7.1 Introduction 91
7.2 Calculation of the Fourier coefficients 94
7.3 Functions of an arbitrary period 96
7.4 Convergence of Fourier series 97
7.5 Functions defined over a finite interval 98
7.6 Even and odd functions 99
7.7 Fourier cosine series for even functions 100
7.8 Fourier sine series for odd functions 101
7.9 Half-range expansions 102
7.10 Parseval’s Theorem 103
7.11 Differentiation of Fourier series 105
7.12 Integration of Fourier series 107
11 Index 185
1
Integration and some applications
f f −1 ( y ) = y f −1 f ( x ) = x.
and
Later on, we might discard the variable y and also use x as the
variable in the domain of f −1 ( x ) on the understanding that here x
is in the range of f .
π/2
-1 0 1 y
−π/2
π/2
-10 0 10 y
−π/2
P( x )
Given a rational function f ( x ) = Q( x) with deg( P( x )) <
deg( Q( x )) find polynomials gi ( x ) and Ai ( x ) such that
P( x ) A1 Ak
f (x) = = +...+ ,
( x − a1 ) . . . ( x − a k ) x − a1 x − ak
where a1 , . . . , ak are pairwise distinct. Note that then
P ( x ) = A 1 ( x − a 2 ) · · · ( x − a k ) + · · · + A k ( x − a 1 ) · · · ( x − a k −1 ).
P ( a i ) = A i ( a i − a 1 ) · · · ( a i − a i −1 ) · ( a i − a i +1 ) · · · ( a i − a k )
2 1
Hence f ( x ) = x +1 + x −1 .
mathematical methods 2 11
P( x ) B1 Bc
f (x) = = +...+ .
( x − a)c x−a ( x − a)c
P( x ) = B1 ( x − a)c−1 + · · · + Bc−1 ( x − a)1 + Bc .
So choosing x = a we see that P( a) = Bc . If we then evalu-
ate the expression for further (simple) values of x, we can deduce
equations for the remaining Bi s.
Example 1.2.
P( x ) 2x + 3 B1 B2 B3
f (x) = = 3
= + 2
+ .
Q( x ) ( x − 1) ( x − 1) ( x − 1) ( x − 1)3
x=0 −3 = − B1 + B2 − 5 (1.1)
x=2 7 = B1 + B2 + 5 (1.2)
2 5
Hence B1 = 0 and B2 = 2. Thus f ( x ) = ( x −1)2
+ ( x −1)3
.
P( x ) A1 C x + C2
f (x) = = + 21 .
( x − a)( x2 + bx + c) x−a x + bx + c
P( x ) = A1 ( x2 + bx + c) + (C1 x + C2 )( x − a). Hence (as before)
P( a)
A1 = .
( a2 + ba + c)
Example 1.3.
P( x ) 2x + 3 A1 C x + C2
f (x) = = 2
= + 12 .
Q( x ) ( x − 1)( x + 1) ( x − 1) x +1
x=0 −3 = − 25 + C2 C2 = −1/2
x = −1 − 41 = − 54 − C21 + C2
2 C1 = −5/2
Z b b Z b
0
f ( x ) g( x )dx = f ( x ) g( x ) − f ( x ) g0 ( x )dx,
a a This
a R 0 rule is sometimes remembered as
u v = [uv] − uv0 .
R
Z
an antiderivative for f 0 ( x ) g( x ) is f ( x ) g( x ) − f ( x ) g0 ( x )dx.
Z 4
Example 1.4. Find x ln xdx.
2
Solution: Notice from the above formula that the f function is integrated
while the g function is differentiated, so we need to choose which one is f
and which one is g to make the integral do-able using simple algebra. Here
we are going to choose f ( x ) = x and g( x ) = ln x (in this case choosing
the reverse is possible too but it makes the algebra a little bit harder). Thus
F ( x ) = x2 /2 and g0 ( x ) = 1/x.
4
x2
Z 4 Z 4 2
x 1
x ln xdx = ln x − dx
2 2 2 2 x 2
2 4
x Recall that ln( ab) = ln( a) + ln(b) and
= 8 ln(4) − 2 ln(2) − = 14 ln(2) − 3
4 2 so ln( an ) = n ln( a).
Hence Z
x2 cos xdx = x2 sin x + 2x cos x − 2 sin x + c.
Z
f g( x ) g0 ( x )dx
F g( x ) =
Rr √
Example 1.6. Find −r r2 − x2 dx where r is a positive constant.
14 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Motivation : Area
Let f ( x ) ≥ 0 be a continuous function on some finite interval
[ a, b], and let A be the area of the region under the graph of f , i.e.
the area of {( x, y) ∈ R2 : x ∈ [ a, b], 0 ≤ y ≤ f ( x )}.
To calculate an approximation to A, choose n ∈ N and consider a
partition a b
A = A1 + A2 + ... + An
n
≈ f (c1 )∆x1 + f (c2 ) ∆x2 + ... + f (cn )∆xn = ∑ f (c j )∆x j .
j =1
where the limit is taken over all possible partitions and all possible
choices of c j s. These considerations lie beneath the following defini-
tion of the Riemann integral (there are other types of integrals).
a
f ( x )dx = lim
|| P||→0
∑ f (c j )∆x j ,
j =1
Notes:
• Finally, note that the definition of the definite integral has not
used derivatives or anti-derivatives.
• For the interested student: the exact meaning of the limit in the
Rb
definition is the following: ∃ L ∈ R (L = a f ( x )dx) such that
∀e > 0, ∃δ > 0 s.t. ∀ partitions P of [ a, b] and for any choice of
the points c j ∈ [ x j−1 , x j ]( j = 1, 2, ..., n), | L − ∑nj=1 f (c j )∆x j | <
e ∀ P s.t. || P|| < δ.F
So Z b
cdx = lim c(b − a) = c(b − a).
a n→∞
This is just the area of the rectangle of base length b − a and height c, as
would be expected.
n ( n +1)
1 + 2 + · · · + ( n − 1) + n 2 n+1
= = =
n2 n2 2n
1
which converges to 12 as n → ∞. Therefore 0 xdx = 1/2, which is
R
expected since this is the area of the triangle formed by the graph of f ( x ),
the x-axis and the line x = 1.
A+ = Area{( x, y) : x ∈ [ a, b], f ( x ) ≥ 0, 0 ≤ y ≤ f ( x )}
is the area of the region above the x-axis and under the graph of f ,
while
A− = Area{( x, y) : x ∈ [ a, b], f ( x ) ≤ 0, f ( x ) ≤ y ≤ 0}
is the (positive) area of the region under the x-axis and above the
graph of f .
It is clear from the definition of the integral that every Riemann
sum gives an approximation to the value of the integral. In general,
partitions of smaller size subintervals lead to better approximations.
A0 ( x ) = f ( x )
and therefore
A( x + h) − A( x )
f (u) ≤ ≤ f ( v ).
h
20 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Note: The FTC says that (for a continuous function and for x in
an appropriate interval) d a f (t)dt = f ( x ) and so d x f (t)dt =
Rx Ra
dx dx
− f ( x ).
x3 /3 + c, c ∈ R.
R x4
Example 1.20. Let g( x ) = 1 1+1ln t dt, for x > 1. Find g0 ( x ).
Rx
Solution: Consider the function A( x ) = 1 1+1ln t dt, so A0 ( x ) = 1+1ln x
by the FTC. Since g( x ) = A( x4 ), application of the Chain Rule gives
3 3
g0 ( x ) = A0 ( x4 )( x4 )0 = 1+4x
ln( x4 )
= 1+44xln(x) for any x > 1.
Generalising this example, we can prove the following theorem.
0
h20 ( x ) − f h1 ( x ) h10 ( x ).
g ( x ) = f h2 ( x )
the bucket is at level x, the length of the part of the rope that remains to be
lifted is 10 − x, so its weight is f ( x ) = (10 − x )0.8 N. Thus, the work
done to lift the rope is
Z 10
W2 = 0.8(10 − x )dx = 40 Joules.
0
Hence the total work done is W = W1 + W2 = 50 + 40 = 90 J.
Example 1.25. A tank has the shape of an inverted cone with height 10m
and base radius 4m. It is filled with water to a height of 8m. It is required
to pump all the water out of the tank by inserting a hose through the top
to which a pump is attached. Calculate the work done against gravity to
achieve this. (Assume that the gravity constant is g = 9.8m/s2 and the
density (mass per unit volume) of the water is 1000kg/m3 ). 4
r j /4 = c j /10. 8
The force required to raise this volume against the force of gravity is
then given by (1000Vj ) g = 10(9.8)16πc2j ∆x j .
The work done in raising the jth layer to the top of the tank is then
Wj ≈ (1568πc2j ∆x j )(10 − c j ), and so the total work done is
Z 8
2048
W= 1568πx2 (10 − x )dx = 1568π ≈ 3.36 × 106 J.
0 3
Note: In some problems of this kind the density will be given as
weight per unit volume (i.e. the gravity constant has already been
taken into account); or SI2 units (Newtons, Joules) are not being 2
International System of Units (ab-
used. breviated SI from French: Système
international d’unités
By the Mean Value Theorem there exists c j ∈ ( x j−1 , x j ) such that Recall the Mean Value Theorem:
f ( x j ) − f ( x j−1 ) = f 0 (c j )( x j − x j−1 ). Let f be continuous on an interval
q [ a, b] and differentiable on ( a, b),
Therefore L j ≈ 1 + ( f 0 (c j ))2 ∆x j . Hence for the total length of then there exists c ∈ ( a, b) such that
f (b)− f ( a)
q f 0 (c) = b− a .
the curve we obtain L ≈ ∑nj=1 (1 + ( f 0 (c j ))2 ∆x j .
p
This is a Riemann sum for the function 1 + ( f 0 ( x ))2 , so letting
Rbp
|| P|| → 0 gives L = a 1 + ( f 0 ( x ))2 dx.
mathematical methods 2 23
If the limit exists then the improper integral is called convergent. If the
limit does not exist then the improper integral is called divergent.
(b) Similarly
Z b Z b
f ( x )dx = lim f ( x )dx.
−∞ t→−∞ t
Remark 1.28. It can be shown that the choice of c is not important; i.e.
if for one particular choice of c the integrals on the right-hand-side are
convergent, then the same is true for any other choice of c and the sum of
the two integrals is always the same3 . 3
Try to prove this.
Z ∞
1
Example 1.29. Find the improper integral dx if it is convergent
1 x3
or show that it is divergent.
Solution: For any t > 1,
Z t t
1 1 1 1 1
3
dx = − 2 = − 2+ →
1 x 2x 1 2t 2 2
as t → ∞. The improper integral is convergent and its value is 1/2.
(This example shows that the area A under the graph of f ( x ) = 1/x3
and above the interval [1, ∞) is finite even though the ‘boundaries’ of the
area are infinitely long.)
24 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Z ∞
1
Example 1.30. Find the improper integral dx if it is convergent or
1 x
show that it is divergent.
Solution: For any t > 1 we have
Z t
dx
= [ln x ]1t = ln t → ∞
1 x
R∞
as t → ∞. Hence the integral 1 1/x dx is divergent. (The area under
the graph of f ( x ) = 1/x and above the interval [1, ∞) is unbounded.)
Example 1.31. Find all constants p ∈ R such that the improper integral
Z ∞
dx
p is convergent.
1 x
Solution: For p = 1 the previous example shows that the integral is
divergent. Suppose p 6= 1. Then
− p +1 t
t 1− p − 1
Z t
dx x
p = = .
1 x −p + 1 1 1− p
When 1 − p < 0, limt→∞ t1− p = 0 so the integral is convergent. When
1 − p > 0, limt→∞ t1− p = ∞ and therefore the integral is divergent.
Hence the integral is divergent for p ≤ 1 and otherwise convergent.
If the limit exists, then the improper integral is called convergent, other-
wise it is divergent.
(b) In a similar way, if f is continuous on ( a, b] however it has some
kind of singularity at a, e.g. f ( x ) → ∞ or −∞ as x → a+ . We define
Z b Z b
f ( x )dx = lim f ( x )dx.
a t→ a+ t
R1 √
Example 1.33. Consider the integral 0 1/ x dx. This is an improper
√ √
integral, since 1/ x is not defined at 0 and 1/ x → ∞ as x → 0.
For 0 < t < 1 we have
Z 1 √ √ √
1/ x dx = [2 x ]1t = 2 − 2 t → 2
t
as t → 0+ . So the improper integral is convergent and its value is 2.
Z 2
1
Example 1.34. Consider the integral dx. It is improper, since
0 x−1
f ( x ) = 1/( x − 1) is not defined at x = 1 and is unbounded near x = 1.
However, f ( x ) is continuous on [0, 1) ∪ (1, 2].
Therefore
Z 2 Z 1 Z 2
dx dx dx
= + ,
0 x−1 0 x−1 1 x−1
and the integral on the left-hand-side is convergent if and only if both
integrals on the right-hand-side are convergent.
Consider Z 1 Z t
1 1
dx = lim dx.
0 x − 1 t → 1 − 0 x − 1
Given 0 < t < 1, we have
Z t
1/( x − 1)dx = [ln | x − 1|]0t = ln |t − 1| − 0 = ln(1 − t) → −∞
0
as t → 1− .
Hence this integral is divergent and therefore the original integral is
also divergent4 . 4
Note that if one integral is shown
Z 1 to be divergent we do not need to
Exercise 1.7.1. For what values of p is x p dx improper, and in that check if the other one is convergent or
0 not, we already know the answer: the
case for what values of p is the integral divergent? Compare with Example original improper integral is divergent.
1.31.
1.8 Quadrature
R1
Example 1.35. Approximate 0 1+1x2 dx using the composite trapezoidal
rule for N = 4. Deduce an interval for the true value of the integral.
1
Solution: Let f ( x ) = 1+ x 2
. We have h = 1/4. So
Z 1
1 1
dx ≈ ( f (0) + 2 f (1/4) + 2 f (1/2) + 2 f (3/4) + f (1))
0 1 + x2 8
1
= (1 + 2 · 16/17 + 2 · 4/5 + 2 · 16/25 + 1/2)
8
= 5323/6800 ≈ 0.78279 . . .
(6x2 −2)
To estimate the error, we need f 00 ( x ) = (1+ x 2 )3
, which is an increasing
function, so for c ∈ (0, 1), −2 < f 00 (c) < 1/2. Since the error is
3
− Nh 00
12 f ( c ) for some c ∈ (0, 1), we have
1 1 1
− · < error < − · (−2).
16 · 12 2 16 · 12
Hence the true value is in (0.78279 − 0.02605, 0.78279 + 0.01042) =
(0.75674, 0.79321). Exercise: evaluate the integral directly and check the
value is in this interval.
Example 1.36. Show that Simpson’s rule on [0, h] gives exact values of
the integral for all cubic polynomials.
Solution: Let an arbitrary cubic be p( x ) = c0 + c1 x + c2 x2 + c3 x3 . Then
Rh
the definite integral is 0 p( x )dx = c0 h + c1 h2 /2 + c2 h3 /3 + c3 h4 /4.
The Simpson rule approximation is 6h ( p(0) + 4p(h/2) + p(h))
2 3 ! !
h h h h 2 3
= c0 + 4 c0 + c1 + c2 + c3 + c0 + c1 h + c2 h + c3 h
6 2 2 2
h 3
= 6c0 + 3c1 h + 2c2 h + c3 h3
2
= c0 h + c1 h2 /2 + c2 h3 /3 + c3 h4 /4.
6 2
The error term is actually proportional to the fourth derivative of the inte-
grand at some unknown point, so using the fact that the fourth derivative
of a cubic polynomial is zero everywhere, we also get the result.
Exercise 1.8.1. Use the first Taylor polynomial of f about h/2 to show
1 3 00
that the error for the mid-point rule on [0, h] is 24 h f (c), for some un-
known point c ∈ (0, h). Note that using the Taylor polynomial approxima-
tion assumes that certain derivatives exist everywhere in the interval.
2
Double and triple integrals
b
Px : a = x0 < x1 < x2 < . . . < x N −1 < x N = b. c d
x
Volume under z = f ( x, y).
Similarly, given another arbitrary integer M ≥ 1, consider a parti-
tion of the y–interval [c, d]
pick a point ( xi∗ , y∗j ) (with xi−1 ≤ xi∗ ≤ xi and y j−1 ≤ y∗j ≤ y j ), and yj
∆ yj = yj − yj−1
the value of z = f ( xi∗ , y∗j ) describes approximately the height of the
surface above the xy-plane at the point ( xi∗ , y∗j ) and hence f ( xi∗ , y∗j ) y0 = c
f(xi ,yj )
y
a
xi
b
c yj d
x
M N
M,N →∞
∑ ∑ f (xi∗ , y∗j )∆xi ∆y j
lim
∆x,∆y→0 j=1 i =1
M N
VPx ,Py = ∑ ∑ f (xi∗ , y∗j )(∆xi ∆y j ) height of column · base area
j =1 i =1
!
M N
= ∑ ∑ f ( xi∗ , y∗j )∆xi ∆y j Take limit as N → ∞
j =1 i =1
M Z b
The wall is the intersection of the solid
→ ∑ a
f ( x, y∗j ) dx ∆y j area of wall at y = y∗j · thickness with the plane y = y∗j
j =1
Z d Z b
→V= f ( x, y) dx dy Taking limit as M → ∞
c a
mathematical methods 2 29
Notice that
and that in the calculation above we could have first taken the limit
M → ∞ and then N → ∞. We formalise this idea with the fol-
lowing theorem, which allows us to evaluate double integrals over
rectangular regions using well-known methods of integration with
respect to a single variable and shows that we can change the order
of integration:
R = {( x, y) : 1 ≤ x ≤ 2, −1 ≤ y ≤ 1}.
Then
ZZ Z 1 Z 2
f ( x, y)dA = ( x2 + xy) dxdy
R y=−1 x =1
Z 1 Z 2
= ( x2 + xy) dx dy
−1 1
1 3 1 2 2
Z 1
= x + x y dy
−1 3 2 1
Z 1
7 3
= ( + y) dy
−1 3 2
1
7 3
= y + y2
3 4 −1
7 3 2 7 3 2 14
= (1) + (1) − (−1) + (−1) = .
3 4 3 4 3
Note that attention needs to be paid to the order of the dx, dy
relative to the order of the integral symbols!
R
f ( x, y)dA = lim
M,N →∞
∑ f ( xi∗ , y∗j )(∆xi ∆y j ),
∆x,∆y→0 Ri,j ⊂ R
Theorem 2.6. ZZ
1 dA = Area( R).
R
R = {( x, y) : a ≤ x ≤ b , g( x ) ≤ y ≤ h( x )},
y y = h( x )
y = g( x )
c
a b x
Let f ( x, y) be continuous on R.
Take an arbitrary rectangle R b = [ a, b] × [c, d] containing R in its
interior and consider the extension fˆ of f which is 0 outside R.
Then, using the Remark about Definition 2.4 made earlier and
Fubini’s Theorem for rectangular regions, we get
ZZ ZZ Z b Z d
f ( x, y)dxdy = ˆ
f ( x, y)dxdy = ˆ
f ( x, y)dy dx
R R
b a c
Z b Z g( x ) Z h( x )
= [ fˆ( x, y) dy + fˆ( x, y)dy
a c g( x )
Z d
+ fˆ( x, y)dy ] dx
h( x )
Z b Z h( x )
= f ( x, y)dy dx.
a g( x )
32 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
R = {( x, y) : c ≤ y ≤ d , G (y) ≤ x ≤ H (y)},
Example 2.7. Let R be the triangle bounded by the lines x = 0, y = 1 Since this triangle has half the area of a
and y = x. Compute the area of R using a double integral. unit–square, its area is A = 1/2.
y
Solution: The most important tool that we have is that we are able to y=1
R = {( x, y) | 0 ≤ x ≤ 1, x ≤ y ≤ 1}
and deduce:
ZZ Z 1 Z 1
A= 1 dA = 1 dy dx
R x =0 y = x
1
x2
Z 1 Z 1
1
= [y]1x dx = (1 − x ) dx = x − = .
0 0 2 0 2
2. If R is defined by
For each i, one can then evaluate the double integral of f over Ri
using Fubini’s Theorem.
and
√
R2 = {( x, y) | 1 ≤ x ≤ 4 and x − 2 ≤ y ≤ x }.
The definition of the double integral by Riemann sums allows us to split
the integral over R into two different regions. Therefore we get by Fubini’s
more general theorem:
ZZ ZZ ZZ
I= f ( x, y) dy dx = f ( x, y) dy dx + f ( x, y) dy dx
R R1 R2
Z 1 Z √x Z 4 Z √x
= √ f ( x, y) dy dx + f ( x, y) dy dx
x =0 y=− x x =1 y = x −2
This example illustrates the fact that we can avoid a lot of calcu-
lations by being smart when choosing our order of integration.
R
f ( x, y, z)dV =
L,M,N →∞
lim ∑ ∑ ∑ f (xi∗ , y∗j , z∗k ) ∆xi ∆y j ∆zk ,
∆x,∆y,∆z→0 k =1 j=1 i =1
ZZZ ZZZ
1. k f dV = k f dV for a constant, k
R R
ZZZ ZZZ ZZZ
2. ( f 1 + f 2 ) dV = f 1 dV + f 2 dV
R R R
ZZZ ZZZ ZZZ
3. f dV = f dV + f dV
R1 ∪ R2 R1 R2
for solids R1 and R2 with disjoint interiors.
where
S = {( x, y, z, t) ∈ R4 : ( x, y, z) ∈ R , 0 ≤ t ≤ f ( x, y, z)}.
RRR
In general if f is not necessarily non-negative, then R f dV rep-
resents an oriented 4D volume (which may be positive, negative or
zero).
ZZZ
Theorem 2.14. 1 dV = Vol( R).
R
B = [ a, b] × [c, d] × [ p, q]
= {( x, y, z) ∈ R3 : a ≤ x ≤ b, c ≤ y ≤ d, p ≤ z ≤ q}
Z b Z q Z d
= f ( x, y, z)dy dz dx, etc.
a p c
ZZZ
That is, all 6 possible iterated integrals are equal to f dV.
B
Important Remark:
Z bZ dZ q Z b Z d Z q
f ( x, y, z) dz dy dx = f ( x, y, z) dz dy dx.
a c p a c p
That is:
R = {( x, y, z) : ( x, y) ∈ D , g( x, y) ≤ z ≤ h( x, y)},
This can be derived using an argument similar to that for double in-
tegrals (see "Evaluations of double integrals over bounded regions"
above).
Similarly, if
R = {( x, y, z) : ( x, z) ∈ D , g( x, z) ≤ y ≤ h( x, z)},
etc.
−1 ≤ x ≤ 1
2
x −1 ≤ z ≤ 0
z/3 ≤ y ≤ −2z
The channel weir.
ZZZ
M= ρ( x, y, z) dV.
R
Myz = xρ( x, y, z) dV
Z Z ZR
Mxz = yρ( x, y, z) dV
Z Z ZR
Mxy = zρ( x, y, z) dV
R
ZZ
My = xρ( x, y) dV
Z ZR
Mx = yρ( x, y) dV
R
My M x
(Cx , Cy ) = ( , ).
M M
Example 2.18. Calculate the mass and the centre of mass for a half-disk
x2 + y2 ≤ r2 , y ≥ 0, where the mass density function is ρ( x, y) = 1. (so the mass is just the area)
ZZ Z r Z √r 2 − x 2
Mdisk = M = ρ( x, y) dA = dy dx
disk x =−r y =0
√
πr2
Z r Z r
2 2
[y]y=r 0− x dx =
p
= r2 − x2 dx = .
x =−r x =−r 2
1
E= Mv2 , p = Mv, if the body moves in a straight line
2
1
E = Iω 2 , L = Iω, if the body is rotating.
2
Notice that in the second equation, the role of the total mass M is
replaced by the moment of inertia I and the role of the velocity v
is replaced by the angular velocity ω. The symbol E represents the
kinetic energy in the first equation and the rotational kinetic energy
in the second one.
If A is the axis of rotation, it can be shown that the scalar I is
given by the expression
ZZZ
I= ( distance from ( x, y, z) to A)2 ρ( x, y, z) dV.
R
where Ixx , Iyy , Izz give the moment of inertia about the x, y or z–axis
respectively. Intuitively, as a mass is moved away from the x-axis
Ixx increases.
For the interested student:
If we want to know the value of E and I for an arbitrary axis A,
RRR 2
the integral R ( distance from ( x, y, z ) to A) ρ ( x, y, z ) dV is often
not easy to compute. Instead we can compute further moments:
ZZZ
Iyz = yzρ( x, y, z) dV
Z Z ZR
Ixz = xzρ( x, y, z) dV
ZZZ R
Ixy = xyρ( x, y, z) dV,
R
which defines the last moments in the moment of inertia matrix (or
sometimes called inertia tensor):
Ixx − Ixy − Ixz
I = − Ixy Iyy − Iyz .
− Ixz − Iyz Izz
x = g ( u ),
dx = g0 (u) du,
x=a =⇒ u = g−1 ( a) = c,
x=b =⇒ u = g−1 (b) = d.
where S is the new region of integration, i.e. the domain of the new
coordinates u,
v.
∂g(u,v)
The matrix ∂(u,v) is the Jacobian matrix of the transformation
g(u, v) at (u, v) (recall the definition from Section 7.5 of the Mathe-
matical Methods 1 notes), that is,
!
∂φ(u,v) ∂φ(u,v)
∂g(u, v)
= ∂ψ(u,v) ∂ψ(u,v) ,
∂u ∂v
∂(u, v)
∂u ∂v
and its determinant J (u, v) is called the Jacobian of g(u, v) at (u, v).
Note that the change of coordinates requires a transformation
function g that maps a two-dimensional region onto another two-
dimensional region. Moreover g can be differentiated on S. This
mapping has to be one–to–one, so that for each point ( x, y) in R
there is a unique point (u, v) in S with g(u, v) = ( x, y).
More precisely, we have the following.
Then
∂g(u, v)
ZZ ZZ
f ( x, y) dx dy = f (g(u, v)) det du dv.
R S ∂(u, v)
The above formula remains true if the conditions (i) and (ii) are sat-
isfied on S except possibly on a finite union of continuously differen-
tiable curves.
S R
u x
Proof. We see from the picture that the area of the parallelogram is
θ
| a|.H = | a||b| sin θ, where θ is the angle between a and b. a
Now we are considering both vectors as part of the plane z = 0
in 3 dimensions. We set a0 = ( a1 , a2 , 0), b0 = (b1 , b2 , 0). Notice that
the angle between a0 and b0 is also θ, and that | a0 | = | a|, |b0 | = |b|.
46 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
we get that
! !
0 0 a1 a2 a1 b1
| a||b| sin θ = | a ||b | sin θ = det = det
b1 b2 a2 b2
∂φ(u j , vk )
φ(u j + ∆u j , vk ) − φ(u j , vk ) ≈ ∆u j .
∂u
Using similar approximations for ∆vk and for ψ, we get an area
approximation in R that is:
∂φ(u j ,vk ) ∂φ(u j ,vk ) ∂φ(u j ,vk ) ∂φ(u j ,vk )
! !
∆u j ∆vk
det ∂ψ(u ,v )
∂u ∂v
∂ψ(u j ,vk ) = det ∂ψ(u j ,vk ) ∂ψ(u j ,vk ) ∆u j ∆vk
∂u ∂v
∂u
j k
∆u j ∂v ∆v k ∂u ∂v
∂g(u j , vk )
= det ∆u j ∆vk .
∂(u, v)
We can use Riemann sums over our partitioning to approximate the
double integral of f ( x, y) on R:
N M
∂g(u j , vk )
∑∑ f (φ(u j , vk ), ψ(u j , vk )) det
∂(u, v)
∆u j ∆vk
j =1 k =1
Solution: The difficulty with this problem is that the boundaries of the
region do not align conveniently with the coordinate axes. Guided by the
fact that a linear change of variables maps parallelograms in one coor-
dinate system to parallelograms in the other, it would seem that putting
u = 3x − 2y and v = 2x + y would be a good choice. Then the domain in
(u, v) space is defined by −1 ≤ u ≤ 2 and 0 ≤ v ≤ 14.
To calculate the Jacobian we first need to find x and y in terms of u
and v; it is easy to show that with u and v defined as above then x =
(u + 2v)/7 and y = (3v − 2u)/7. Then
! !
∂x ∂x 1 2
7 7 1
det ∂y ∂y = det
∂u ∂v
2 3 = .
∂u ∂v
−7 7 7
1 2 2
Z 14 Z 2 Z 14
3 14
1 1
ZZ Z
(3x − 2y) dA = (u) · dudv = u dv = dv = 3.
D 0 −1 7 7 0 2 −1 14 0
RR y4
Example 3.4. Evaluate I = D x dxdy over the region D contained
between the parabolas x = 1 − y2 and x = 4(1 − y2 ).
Solution: First we need to understand the geometry of the region D. It is
clear that the two parabolas are symmetric about the x-axis; they meet at
x = 0, y = ±1 and they cut the x-axis at x = 1 and x = 4 respectively.
If we therefore define new coordinates x = v(1 − u2 ) and y = u then
the two parabolas are given by v = 1 and v = 4 and they meet at u = ±1.
Notice that in these coordinates the appropriate region of integration is
again just a simple rectangle with sides parallel to the u and v axes; this is
much simpler than the corresponding region in ( x, y) space.
Once more we must compute the Jacobian
! !
∂x ∂x
−2uv 1 − u2
J = det ∂y ∂y = det
∂u ∂v = −1(1 − u2 ) = u2 − 1.
∂u ∂v
1 0
y4 u4
ZZ Z 4Z 1
I = dxdy = | J |dudv
D x 1 −1 v (1 − u2 )
u4
Z 4Z 1 Z 4Z 1 4
u
= (1 − u2 )dudv = dudv
1 v (1 − u2 )
−1 1 −1 v
Z 4 5 1 Z 4
u 2 1 2
= dv = dv = ln 4.
1 5v −1 5 1 v 5
48 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
( x, y, z) = g(u, v, w),
∂g(u, v, w)
dx dy dz = det du dv dw,
∂(u, v, w)
S = g −1 ( R ),
then
1 ( u,v,w ) ∂g1 (u,v,w) ∂g1 (u,v,w)
∂g
∂u ∂v ∂w
∂g(u, v, w) ∂g2 (u,v,w) ∂g2 (u,v,w) ∂g2 (u,v,w)
= .
∂(u, v, w)
∂u ∂v ∂w
∂g3 (u,v,w) ∂g3 (u,v,w) ∂g3 (u,v,w)
∂u ∂v ∂w
Then
ZZZ
f ( x, y, z) dx dy dz =
R
∂g(u, v, w)
ZZZ
f (g(u, v, w)) det du dv dw.
S ∂(u, v, w)
The above formula remains true if the conditions (i) and (ii) are sat-
isfied on S except possibly on a finite union of continuously differen-
tiable surfaces.
Example 3.6. A torus (doughnut) has major radius b and minor radius a
b
(0 < a < b). a
To evaluate this determinant the zero in the third row indicates that the
simplest calculation will be using an expansion by this third row and then
J = r cos α(b + r cos α) cos α(sin2 θ + cos2 θ ) + sin α(b + r cos α)r sin α(sin2 θ + cos2 θ )
The example of the torus above illustrates one of the main obsta-
cles to using coordinate transformations in three dimensions. Fre-
quently the calculation of the Jacobian is somewhat tedious and
lengthy but many practical problems involve circular, cylindrical
or spherical geometries. There are three sets of coordinates that
are specifically designed to assist with the solution of these classes
of problem and, rather that having to work from first principles
each time, it is convenient to be able to use the standard results for
these coordinates. Accordingly, in the remainder of this chapter,
we introduce these three special coordinate systems known as polar
coordinates, cylindrical coordinates and spherical coordinates.
The choice of substitution/new coordinate system for solving a
double or triple integral depends on the particular problem. Above
we discussed a problem with a geometry that was simple on S but
complex on R. In these cases it is useful to make the transformation
and evaluate the original problem in the new coordinate system. In
the applied sciences, like physics, the new coordinate systems are
sometimes called natural coordinate systems.
x = ρ cos θ,
y = ρ sin θ.
adapt results.
Using the proposed transformation:
ρ cos θ x
g(ρ, θ ) = = ,
ρ sin θ y
we see that:
!
−ρ sin θ
∂g(ρ, θ ) cos θ
det = det = ρ(cos2 θ + sin2 θ ) = ρ.
∂(ρ, θ ) sin θ ρ cos θ
θ0 = arctan ( a) ,
θ1 = arctan (b) ,
0≤ρ≤r
θ0 ≤ θ ≤ θ1 ,
4 ≤ x 2 + y2 ≤ 9
1
x ≤ y ≤ 2x
2
y θ
R S
r
x
Solution: The region is a ring segment. From the first inequality we get
bounds on the radius
2 ≤ ρ ≤ 3.
For the lower bound on θ note that we need to use y = x/2. Substitution
of g gives us θ0 = arctan 1/2 and similarly by using y = 2x we find
θ1 = arctan 2. Therefore the area is:
Z arctan 2 Z 3 Z arctan 2 2 3
ρ 5 1
A= ρ dρ dθ = dθ = arctan 2 − arctan .
arctan 1/2 2 arctan 1/2 2 2 2 2
x2 y2
+ ≤ 1.
a2 b2
(The change of coordinates you need to make here is similar to polar coor-
dinates but is slightly different.)
Example 3.14. Evaluate the volume of an ice cream cone given by:
π
0 ≤ r ≤ R, 0 ≤ θ ≤ 2π, 0≤φ≤ .
4
Solution: Given these inequalities, we can determine the volume Vcone of
the cone by:
ZZZ ZZZ
Vcone = 1 dx dy dz = r2 sin φ dr dφ dθ
cone S
Z 2π Z π/4 Z R Z 2π Z π/4 3 R
2 r
= r sin φ dr dφ dθ = sin φ dφ dθ
0 0 0 0 03 0
√ !
R3
Z 2π
R3 2 2πR3 √
= [− cos φ]0π/4 dθ = 1− [θ ]2π
0 = 2− 2 .
3 0 3 2 6
π 2π
0 ≤ r ≤ R, ≤θ≤ , 0 ≤ φ ≤ π.
4 3
Solution:
Given these inequalities, we can determine the volume Vwedge of the
wedge by:
ZZZ ZZZ
Vwedge = 1 dx dy dz = r2 sin φ dr dφ dθ
wedge S
2π R
r3
Z R
5π 3
Z πZ
3
= 2
r sin φ dr dθ dφ = [θ ]2π/3 π
π/4 [− cos φ ]0 = R .
0 π
4 0 3 0 18
4
Path and surface integrals
dr
x
v(t) = (t) = ṙ(t) a path in the xy–plane.
dt
58 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Example 4.1. Find a parametric form of the straight line from (1, 2) to
(3, 6).
Solution: Obviously we are trying to find functions x (t) and y(t) that
together form the components of the parametric form r(t) = ( x (t), y(t)).
Both of them will be linear and therefore of the form:
x ( t ) = x0 + m x t y(t) = y0 + my t.
3. the range of r(t) includes both the starting and end points:
r( a) is the starting point and r(b) is the end point.
Exercise 4.1.1. Find two parametric forms r(t) of the straight line from
y
(1, 2) to (3, 6) where in the first case t ∈ [0, 2] and in the second case
t ∈ [−1, 1].
Note that the inequality y ≥ 0 implies that we only recover the positive
square root solution of y2 = 1 − x2 . The domain of our explicit curve is
x ∈ [−1, 1] and the range y ∈ [0, 1]. One possible parametric form can be
found by simply choosing t = x and then we get:
x (t) √ t
r(t) = = , t ∈ [−1, 1].
y(t) 1 − t2
cos(π − t)
x (t)
r(t) = = , t ∈ [0, π ].
y(t) sin(π − t)
If we further require that a particle moves once around the circle starting
and ending at (1, 0) we can write it as the union of these rather messy
parametric forms:
1−t
r(t) = p , t ∈ [0, 2],
1 − (1 − t )2
p t−1
r(t) = , t ∈ [0, 2].
− 1 − ( t − 1)2
With the help of these examples we have now revised the meth-
ods for finding parametric forms of curves. We can summarise
these methods in a more general context in the following way.
Curve given as an explicit equation: Let C be the curve given by
y = f ( x ) on the domain x ∈ D . Let ( x0 , y0 ) be the starting point
and ( x1 , y1 ) the end point of the curve. A parametric form r(t) of C
can be found by choosing t = x, where t ∈ [ x0 , x1 ]. Then
t
r(t) = , t ∈ [ x0 , x1 ].
f (t)
1. The ellipse:
x2 y2
+ = 1.
a2 b2
Let ( x0 , y0 ) = ( x1 , y1 ) = ( a, 0).
2. The astroid:
x2/3 + y2/3 = 1.
3. A spiral whose distance from the origin increases by 1 each time the
path completes one revolution.
7−→
(u, v) ∈ R2 7−→ ( x, y, z) ∈ R3
( x, y, z) = S(u, v),
∂S
Su (u, v) = (u, v),
∂u
∂S
Sv (u, v) = (u, v).
∂v
Since the surface S depends on two parameters, we obviously have
two tangent vectors at each point (u, v). From MATH1001 we re-
member that it is often useful to work with the vector that is nor-
mal to the surface:
N = Su × Sv ;
in general this vector is a function of u and v.
As we have seen, it is not a difficult task to find a parametric
form when an explicit function is given for a path. The same holds
for surfaces:
x = u, y = v, z = u2 − v3 .
Therefore
S(u, v) = (u, v, u2 − v3 ), 0 ≤ u ≤ 1, 0 ≤ v ≤ 4.
Again, like in the previous section, the simplest way to find the
parametric form of an implicitly given surface is to transform it into
an explicit one.
r(t 3) r(t 4)
r(t2) r(t5)
r(t1)
r(t0)
As the number of points in the partition increases (that is, as n
increases and the points on the curve become closer together) the
approximation to the length of the curve increases – a result that is
essentially a consequence of the triangle inequality. Indeed, in the
limit as n → ∞ the approximation length will converge to the actual
length of the curve C. That is,
n
L(C ) = lim
n→∞
∑ |r(tk ) − r(tk−1 )|.
k =1
(Informally this is saying that distance along the path between two
points that are near to each other is approximately the same as the
straight-line distance between them.)
If we divide both sides of this equation by ∆t = tk − tk−1 and
take the limit as ∆t → 0 then
s ( t k ) − s ( t k −1 ) |r(tk ) − r(tk−1 )|
lim = lim
∆t→0 ∆t ∆t→0 ∆t
so that
ds dr
= .
dt dt
Using this result we can conclude that it is possible to write the line
integral as
Z b Z b
dr ds
L(C ) = dt = dt,
a dt a dt
Example 4.7. Determine the length of the spiral S given by the parame-
terisation r : [0, 2π ] 7→ R3 , where r(t) = (cos t, sin t, t).
Solution: We compute
Example 4.8. Find the length of the part of the astroid x2/3 + y2/3 = 1
which is contained in the second quadrant (x ≤ 0, y ≥ 0). y
cos3 t
r(t) = .
sin3 t
Since our curve is in the second quadrant, it goes from (0, 1) to (−1, 0)
which corresponds to π/2 ≤ t ≤ π.
Now
−3 cos2 t sin t
ṙ(t) =
3 sin2 t cos t
so
q p
|ṙ| = 9 cos2 t sin2 t(cos2 t + sin2 t) = 9 cos2 t sin2 t = −3 cos t sin t.
Notice that the quantity |ṙ| must be positive and that for the range
π/2 ≤ t ≤ π it is the case that cos t ≤ 0 and sin t ≥ 0 so the prod-
uct cos t sin t ≤ 0. Hence |ṙ| = −3 cos t sin t. Then
3h 3 3
Z π iπ
L= −3 cos t sin t dt = − sin2 t = − (0 − 1) = .
π/2 2 π/2 2 2
mathematical methods 2 65
and hence
Z Z π √
( x + y + z)ds = (cos t + sin t + t) 2 dt
C 0
π √ √
t2 π2
= sin t − cos t + 2 = 2 2+ .
2 0 2
is well-defined.
Consider some infinitesimal rectangle with sides of length ∆u
and ∆v, based at (u, v) in the domain D. This rectangle has an area
∆u∆v in D. Moreover, as the surface S is parameterised in terms of
u and v, the rectangle in D will correspond to some infinitesimal
surface ∆S in S:
v y
v+ v
S
v
S
u u+ u u x
z
Note that in general the element ∆S in S will not be rectangular
but will have some other shape. Unless the parameterisation S is
extremely special it is not possible to find the area of ∆S precisely
but we can derive a very good approximation for its value. We use
now an argument similar to one we used in Chapter 3. Consider
the vertex S(u, v) of ∆S. The shape of ∆S is approximately that of
the parallelogram defined by the vectors
S(u + ∆u, v) − S(u, v) and S(u, v + ∆v) − S(u, v).
The area of the parallelogram is the length of the cross-product
of these vectors. That is
∆S ≈ |(S(u + ∆u, v) − S(u, v)) × (S(u, v + ∆v) − S(u, v))|
∂S ∂S
≈ (u, v)∆u × (u, v)∆v
∂u ∂v
≈ |Su (u, v) × Sv (u, v)| ∆u∆v
≈ |N(u, v)| ∆u∆v.
Hence the total surface area of S is approximately the Riemann sum
∑ N(ui , v j ) ∆ui ∆v j . Taking the limit for size of partitions going to
0, we get the following theorem.
Example 4.11. Find the area of the part of the surface z = x + y2 that
lies above the triangle with vertices (0, 0), (1, 1) and (0, 1).
Solution: We use u = x and v = y to parametrize the surface:
S(u, v) = (u, v, u + v2 ) v
1
and D
D = {(u, v)|0 ≤ v ≤ 1, 0 ≤ u ≤ v}.
0 1 u
Then the tangent vectors are Su = (1, 0, 1) and Sv = (0, 1, 2v) and hence
e1 e2 e3 p
Su × Sv = det 1 0 1 = (−1, −2v, 1) ⇒ |Su × Sv | = 2 + 4v2
0 1 2v
Therefore we get,
Z 1 Z v p Z 1 p
A = 2 + 4v2 du dv = v 2 + 4v2 dv
v =0 u =0 v =0 We use a substitution w = 2 + 4v2 ,
i1 √ √ dw = 8vdv to solve the last integral.
1 h 6 2
= (2 + 4v2 )3/2 = − .
12 0 2 6
Note: Computing the integral in the other order of integration is much
harder.
Example 4.12. Find the surface area of the torus (doughnut)
Sα = (− a sin α cos θ, − a sin α sin θ, a cos α) = a(− sin α cos θ, − sin α sin θ, cos α),
Sθ = (−[b + a cos α] sin θ, [b + a cos α] cos θ, 0) = [b + a cos α](− sin θ, cos θ, 0).
Then
e1 e2 e3
N = Sα × Sθ = a[b + a cos α] det − sin α cos θ − sin α sin θ cos α
− sin θ cos θ 0
= − a[b + a cos α](cos α cos θ, cos α sin θ, sin α),
68 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
and so
p
|N| = a|b + a cos α| cos2 θ cos2 α + sin2 θ cos2 α + sin2 α = a[b + a cos α]. Note that b − a ≤ b + a cos α ≤ b + a,
and so 0 < b + a cos α.
Note: This is the product of the perimeters of two circles of radii a and b.
∂u ( u, v ) × ∂v ( u, v ).
∂S ∂S
where N(u, v) =
The points (1, 0, 0), (0, 2, 0), (0, 0, 2) correspond respectively to (u, v) =
(0, 0), (1, 0), (0, 1), so the domain is
We compute
∂S ∂S √
× = |(−1, 2, 0) × (−1, 0, 2)| = |(4, 2, 2)| = 2 6
∂u ∂v
and
f (S(u, v)) = (1 − u − v)2u.
mathematical methods 2 69
! p( x, y, z)
p( x, y)
F( x, y) = in R2 , F( x, y, z) = q( x, y, z) in R3 .
q( x, y)
r ( x, y, z)
We see that for each point v of R p the vector field gives us a vector
F(v), whose components are given by the corresponding scalar functions
p, q or r evaluated at v.
F = αT + βn
for some quantities α and β, then these can be found by taking the Note that α and β depend on t!
scalar product with T and n, respectively. If we take the first of
these then
C T(t)
r (t) n (t)
Z Z b Z b
F · T ds = F · T |ṙ| dt = F · ṙ dt.
C a a
Z Z b Z b
F · n ds = F · n |ṙ| dt = F · N dt.
C a a
Note that we used here the definition of path integral seen in Chap-
ter 4.
1.2
0.8
0.6
0.4
0.2
−0.2
−1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1
1
Example 5.6. Calculate the flux of F( x, y, z) = 3 ( x , y , z) through the
surface
S(u, v) = (v, (1 − v2 ) cos u, (1 − v2 ) sin u) 0 ≤ u < 2π, −1 < v < 1 .
Solution: We can compute
e1 e2 e3
∂S ∂S
N= × = det 0 −(1 − v2 ) sin u (1 − v2 ) cos u
∂u ∂v
1 −2v cos u −2v sin u
= (1 − v2 )(2v, cos u, sin u).
We can check that N is indeed pointing outwards. In this example it is
probably easiest to proceed by looking at the point ( x, y, z) = (0, 1, 0)
which corresponds to the parameter values u = 0 and v = 0 and for which
N = (0, 1, 0). It is also clear that the surface must lie in −1 ≤ y ≤ 1
so the outward normal ought to have a positive y-component, which is the
case.
Now we can calculate the flux
Z ZZ
F · dS = F · N dudv
S D
1 1 2π
Z Z
= (v, (1 − v2 ) cos u, (1 − v2 ) sin u) · (1 − v2 )(2v, cos u, sin u)dudv
3 v=−1 u=0
1 1
Z 2π
2π 1 16π
Z Z
= (1 + v2 )(1 − v2 )dudv = (1 − v4 )dv = .
3 v=−1 u=0 3 −1 15
Theorem 5.8. Let r describe the path of a particle moving through a force
field F. The work done is the circulation of F along C:
Z
W= F · dr.
C
∇ f = (∂ x f , ∂y f , ∂z f ) = (3x2 y + yz + 2, x3 + xz + 1, xy + 2z).
F( x, y, z) = (3x2 y + yz + 2, x3 + xz + 1, xy + 2z)
done
Z Z 1 Z 1
F · dr = F(1, 1, t) · (0, 0, 1) dt = (5 + t, 2 + t, 1 + 2t) · (0, 0, 1) dt
A2 B 0 0
= [t + t2 ]10 = 2.
Notice that the sum of these three parts is 6; exactly the same as
we calculated for the direct path from A to B. (Incidentally, note
that the fact that the three results for the individual steps are each
equal to two is pure coincidence; in general we would expect the
work done on the various parts of the path to be different.)
The fact that the results for the two routes taken from A to be B
are the same is not fortuitous for F( x, y, z) is an example of what is
known as a conservative field. As an exercise, design some other path
from A to B and then compute the total work done by F in moving
from A to B. You should get the same answer 6 irrespective of how
simple or complicated your choice of path is (though, of course, if
you choose a very intricate way of moving from A to B you will
have a lot of calculations to do and are more likely to make a mis-
take!). The conclusion we draw is that it does not matter how one
moves between A and B the work done is the same independent
of the path traversed. However, the work done does depend on the
endpoints A and B; if one (or both) of them moves then the work
done is likely to change.
We can encapsulate these ideas in the following definition:
F = ∇ψ
or, in other words, its value depends only on the locations of the
endpoints S and F. Note that in working out this integral we have
not had to say anything about the details of the path C, except
for the very weak assumption that we can parameterise it in an
appropriate way. Therefore the vector field F satisfies our earlier
definition of a conservative field that the circulation depends only
on the choice of endpoints and is independent of the path taken
between the two.
So what do we know now? Well, we can be assured that any
field that can be written as the gradient of a scalar function is con-
servative. What we do not yet know is whether there are any other
conservative fields; i.e. are there conservative fields that cannot be
expressed as the gradient of a suitable scalar function?
So suppose that some vector field F is known to be conservative
and suppose
where C is any path from the origin 0 to the point ( x, y, z); this is
well defined since F is conservative. Notice ψ is a scalar function
(from R3 to R) since the circulation of a field along a curve is a real
number.
Then, for a fixed point ( x, y, z) and a fixed (small) real number
∆x, Z
ψ( x + ∆x, y, z) − ψ( x, y, z) = F · dr
C0
82 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Example 6.5. Find the work done in moving a particle of mass m from
(1, 1, 1) along a straight line C to (2, 2, 2) by the gravity field
( x, y, z)
F( x, y, z) = − GmM ,
( x2 + y2 + z2 )3/2
where G is the gravitational constant and M the mass of a second body
located at the origin.
Solution: First we need to check if the field is conservative. We note that
!
∂ 1 x
p =− 2 2 + z2 )3/2
∂x 2
x +y +z2 2 ( x + y
The same pattern continues for the second and third components of F so
that we have
GmM
F = ∇ψ where ψ( x, y, z) = p .
x2 + y2 + z2
need to do is find the difference in the values of ψ at the ends of the line C
so that the work done is just
( ∂ y f 3 − ∂ z f 2 , ∂ z f 1 − ∂ x f 3 , ∂ x f 2 − ∂ y f 1 ).
We denote it by ∇ × F.
∇ × F = (0, 0, Q x − Py ).
mathematical methods 2 85
∂y f 3 − ∂z f 2 = ∂y (e xy ) − ∂z (z2 ) = xe xy − 2z,
∂z f 1 − ∂ x f 3 = ∂z ( xy) − ∂ x (e xy ) = 0 − ye xy = −ye xy ,
∂ x f 2 − ∂y f 1 = ∂ x (z2 ) − ∂y ( xz) = 0 − 0 = 0.
Note carefully what this statement says (and, perhaps more im-
portantly, what it does not). It tells us that the curl of a conservative
vector field is zero but not the converse; the argument above does
not allow us to conclude that if the curl of a vector field is zero then
the field is conservative. In fact this statement is true for a certain
kind of solids (regions in R2 ), but the proof is relatively long and
complicated so is not part of this unit.
1. If F is conservative, then ∇ × F = 0.
g1 = ∂y (0) − ∂z (1 + 2xe2y ) = 0,
g2 = ∂z (e2y ) − ∂ x (0) = 0 − ye xy = 0
g3 = ∂ x (1 + 2xe2y ) − ∂y (e2y ) = 2e2y − 2e2y = 0.
1. If F is conservative, then Q x = Py in U.
for some apparently arbitrary function g(y, z). However we are not
at liberty to choose g(y, z) as we wish, because we are yet to satisfy
the other two equations
∂ψ ∂ψ
f2 = and f3 = .
∂y ∂z
Let us take the first of these. If we substitute in (♣) then
Z
f 2 = ∂y f 1 dx + ∂y g(y, z).
To find h(z) we need to ensure that the third (and last) equation
f 3 = ∂ψ/∂z holds. Substituting in the form of ψ leads to an equa-
tion for h0 (z) which can be integrated to find h(z) within an arbi-
trary constant. Remember (Remark 6.4) that a potential can only be
tied down within a constant.
The formulae above are somewhat involved (and need not be
remembered). In practice the process of finding potentials is less
daunting as long as it is done in a systematic manner. Let us illus-
trate the method with a couple of examples.
x2 x2
ψ( x, y, z) = + g(y, z) = + ey sin z + h(z). (♦)
2 2
Now we turn to the last equation ∂ψ/∂z = ey cos z. With ψ given by (♦)
it follows that
a constant. Given that this potential does not arise from any particular
situation we are free to choose the value of C and the conventional choice is
to take C = 0. Then h(z) = 0 and (♦) gives the final potential as
x2
ψ( x, y, z) = + ey sin z.
2
Example 6.14. Let F = (ye xy , xe xy + z, y + 2) be a force field. Show that
F is conservative and find the corresponding potential function φ( x, y, z).
Find the work done in moving a particle from A(0, 0, 1) to B(1, 1, 1).
Verify this result by calculating the work done by F in moving the particle
along the straight-line path between A and B.
Solution: First we need to check that ( g1 , g2 , g3 ) = ∇ × F is zero. As
g1 = ∂y (y + 2) − ∂z ( xe xy + z) = 1 − 1 = 0,
g2 = ∂z (ye xy ) − ∂ x (y + 2) = 0
g3 = ∂ x ( xe xy + z) − ∂y (ye xy ) = (1 + xy)e xy − (1 + xy)e xy = 0,
∂φ ∂φ ∂φ
= ye xy = xe xy + z and = y + 2.
∂x ∂y ∂z
note that integrating with respect to z means that our arbitrary function
depends on x and y. If we substitute this form into ∂φ/∂y = xe xy + z
then
∂g ∂g
z+ = xe xy + z ⇒ = xe xy
∂y ∂y
mathematical methods 2 89
φ = (y + 2)z + e xy + h( x ).
∂φ dh dh
= ye xy + = ye xy ⇒ =0
∂x dx dx
so h( x ) is a constant that we choose to be zero. Hence
φ = (y + 2)z + e xy
is a suitable potential.
To find the work done in moving the particle from A(0, 0, 1) to B(1, 1, 1)
R
we can use the result that the circulation C F · dr of a conservative field is
just the difference in the values of the potential at the endpoints of C. Thus
in this case the work done is simply
φ(1, 1, 1) − φ(0, 0, 1) = (3 + e) − (2 + 1) = e.
7.1 Introduction
Periodic function
function for its sketch takes the form of a curve that clearly repeats
itself after an interval of length P. The rather peculiar function on
the right possesses discontinuities at t = (2n + 1)π for n ∈ Z but 0
−3π −2π −π 0 π 2π 3π
92 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
−1
−π −π/2 0 π/2 π
−1
−π −π/2 0 π/2 π
−1
−3π −2π −π 0 π 2π 3π
partial sum N=11 partial sum N=21 partial sum N=101 Figure 7.2: Graph of S N f for different
values of N, on [−π, π ].
1 1 1
0 0 0
−1 −1 −1
−π 0 π −π 0 π −π 0 π
a0
Z π Z π
FS f (t)dt = dt = πa0
−π −π 2
because
Z π Z π
cos(nt) dt = 0 and sin(nt) dt = 0.
−π −π
So we set
1
Z π
a0 = f (t)dt. (7.2a)
π −π
sin( A + B) + sin( A − B)
sin( A) cos( B) = ,
2
or using integration by parts (twice).
This allows us to calculate the other Fourier coefficients. To ob-
tain an we multiply (7.1) by cos(nt) and integrate from −π to π:
Z π
FS f (t) cos(nt)dt = πan
−π
mathematical methods 2 95
so we set
1
Z π
an = f (t) cos(nt)dt. (7.2b)
π −π
Note that only one term on the right-hand-side survives the in-
tegration. To obtain bn we multiply (7.1) by sin(nt) and integrate
from −π to π: Z π
FS f (t) sin(nt)dt = πbn ,
−π
so we set
1
Z π
bn = f (t) sin(nt)dt. (7.2c)
π −π
Again, only one term on the right-hand-side survives. The above
process is called expanding the function f as an infinite sum of
orthogonal1 functions. 1
We can think of functions as vectors,
R π define a dot product f 1 . f 2 =
and
−π f 1 ( t ) f 2 ( t ) dt, so that the functions
The expressions (7.2a, b, c) are called Euler’s formulae. in the Fourier series are mutually
orthogonal (dot product equal to 0).
1
Z π
a0 = f (t)dt Notice we need to assume here that
π −π f is sufficiently regular for all these
integrals to be defined (for instance
1
Z π
it is sufficient for f to be piecewise
an = f (t) cos(nt)dt
π −π continuous).
1
Z π
bn = f (t) sin(nt)dt
π −π
−3π −2π −π 0 π 2π 3π
or, for k = 1, 2, . . .,
2/ π (2k − 1)2
a2k−1 =
a2k = 0.
A similarly complicated calculation yields
1 1
Z π
bn = (π − t) sin(nt)dt = .
π 0 n
Hence the Fourier series of the above function is
∞ ∞
π 2 cos(nt) sin(nt)
FS f (t) =
4
+
π ∑ n 2
+ ∑ n
n=1, n odd n =1
π 2 cos(3t) cos(5t)
= + cos(t) + + +···
4 π 9 25
sin(2t) sin(3t)
+ sin(t) + + +···
2 3
∞ ∞
a0 nπt nπt
FS f (t) = + ∑ an cos + ∑ bn sin
2 n =1
L n =1
L
where
1 L
Z
a0 = f (t)dt,
L −L
Z L
1 nπt
an = f (t) cos dt,
L −L L
Z L
1 nπt
bn = f (t) sin dt.
L −L L
mathematical methods 2 97
Theorem 7.5. Provided that f (t) and f 0 (t) are bounded and piece-
wise continuous on [− L, L], the Fourier series will converge to (be Recall Definition 1.13, of a piecewise
continuous function.
equal to) f (t) except at points of discontinuity, where it will con-
verge to the average of the right- and left-hand limits of f (t) at that
point, i.e.
f (t+ ) + f (t− )
2
where f (t+ ) is the right-hand limit and f (t− ) is the left-hand limit.
Example 7.6. Example 7.3 revisited. We see that f (t) and f 0 (t) are
bounded and piecewise continuous on [−π, π ] (with the only disconti-
nuity point being 0).2 Since for this function we have f (0− ) = 0 and 2
Note that f 0 (t) considered on its full
f (0+ ) = π, then by Theorem 7.5 the Fourier series FS f (t) converges domain also has discontinuity points
in the odd multiples of π.
to the average of these values, ie. π/2 at t = · · · , −2π, 0, 2π, · · · . The
graph of the Fourier series is then
98 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
−3π −2π −π 0 π 2π 3π
Note that it is identical to the graph of f (t) except that it takes the
value π/2 at integer multiples of 2π whereas the function itself is 0 at
these points.
Now φ(t) is defined for all values of t and is naturally a periodic Exponential periodic extension
3
function of period 2L. Hence we are able to apply the theory of
Fourier series to φ(t).
2
Example 7.8. The graph of the Fourier series of the periodic extension of
et , −1 < t ≤ 1 is on the right. 1
−1, 1, 3, 5,· · · .
together with the function in Example 7.2. In particular all the odd
power functions f (t) = t2n+1 are odd.
Even and Odd functions
1
−1
−1 0 1
In words, this tells us that the sum of two even (odd) functions is
itself even (odd). The product of two even or two odd functions is
even while the product of an odd and an even function is odd. The
integral results are particularly important for they facilitate some
great simplifications in the calculation of Fourier series.
Example 7.11. Determine the Fourier series of the even (‘Hats’3 ) function 3
Note that the name ‘hats’ function
derives from the from of its graph
π + t, −π < t ≤ 0 sketched in the Figure.
f (t) = f (t + 2π ) = f (t).
π − t, 0 < t ≤ π
Hats
−3π −2π −π 0 π 2π 3π
2 2(1 − (−1)n ) 0 if n even,
Z π
an = (π − t) cos(nt)dt = =
π 0 n2 π 24 if n odd,
n π
Odd functions must have odd Fourier series and hence an = 0 for
all n leading to a Fourier sine series. Again it is relatively straightfor-
ward to check that the an = 0 because from the definition
1 L 1 L
nπt
Z Z
an = f (t) cos dt = (odd)(even)dt
L −L L L −L
Z L
1
⇒ an = (odd)dt = 0.
L −L
where Z L
2 nπt
bn = f (t) sin dt
L 0 L
−1
−3π −2π −π 0 π 2π 3π
102 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
−1
−2
−2 −1 0 1 2
0 0
−1 −1
−2 −2
−2 −1 0 1 2 −2 −1 0 1 2
f ( t ) = t2 , 0 ≤ t ≤ 1.
Solution:
Even expansion: The even expansion is just g(t) = t2 , −1 ≤ t ≤ 1.
We find the Fourier coefficients (as g(t) is even, bn = 0):
Z 1
2
a0 = 2 t2 dt = ,
0 3
4(−1)n
Z 1
an = 2 t2 cos(nπt) dt = .
0 π 2 n2
Then the Fourier cosine series is
∞
1 4 (−1)n
+
3 π2 ∑ 2
cos(nπt).
n =1 n
then
a20 ∞
1
Z π
[ f (t)]2 dt = + ∑ a2n + bn2 .
π −π 2 n =1
π −π
12 dt +
π 0
(−1)2 dt = ∑ −
nπ
n=1, n odd
∞ ∞
16 1 1 π2
⇒ 1+1 =
π2 ∑ n2
⇒ ∑ n 2
=
8
n=1, n odd n=1, n odd
We can use the result of the previous example to find the value
of ∑∞ 1
n=1 n2 by noting that
∞ ∞ ∞
1 1 1
∑ 2
= ∑ 2
+ ∑ 2
n =1 n n=1, n odd
n n=1, n even n
π 0
2
(π − t) dt =
8
+ ∑ πn2
+∑
n
n=1, n odd n =1
∞ ∞
π2 π2 4 1 1
⇒
3
=
8
+ 2
π ∑ n 4
+∑ 2
n
n=1, n odd n =1
∞
5π 2 4 1 π2
⇒
24
=
π2 ∑ n4
+
6
.
n=1, n odd
∞
1 π4
⇒ ∑ n4
=
96
.
n=1, n odd
then
∞ ∞
(FS f )0 (t) = − ∑ nan sin(nt) + ∑ nbn cos(nt).
n =1 n =1
Moreover, the following theorem tells us what f 0 (t) is, if f is contin-
uous.
The graph of the original function (and of its Fourier series) appeared 1
in Example 7.11. We recognise from Example 7.12 that (FS f )0 (t) is the
Fourier series of the function from Example 7.2, so the graph of (FS f )0 (t)
0
is the one in the margin. Note that it has the value 0 at multiples of π,
the average of the left and right limits, but that f 0 (t) is not defined at
multiples of π. −1
−1
−2
−3π −2π −π 0 π 2π 3π
Since f is an odd function, it has a sine Fourier series and we can show
(exercise) that
∞
(−1)n+1
FS f = ∑ sin(nt).
n =1
n
mathematical methods 2 107
zero as n → ∞ (although you should also remember that the fact that the
nth term → 0 is not in itself sufficient to guarantee the convergence of a
series). For example, if we try to evaluate (♦) at t = 0 we have
1−1+1−1+...
Actual derivative of slopes
and this series clearly does not converge. Moreover, note that t = 0 is 1
not a problem point for f and f 0 (0) = 1/2 so it is not the case that the
differentiated series only fails at points where f (t) is not differentiable or
has some other problem. If we evaluate (♦) at t = π we have −1 − 1 − 0
It turns out that the Integration of Fourier series is more stable than
differentiation in the sense that fewer potential problems tend to
arise.
t2 − π 2 2
Note that the average value of function 4 is − π6 , and that the
integral function is continuous.
−1
−2
−3π −2π −π 0 π 2π 3π
8
Laplace Transforms
defined for all s ∈ R for which the above integral is convergent. We often
write F (s) as L( f ), or, more precisely L( f )(s).
1
L(1) = .
s
Notice that here F (s) is not defined for all real values of s, just for s > 0.
The definition of the ILT now implies that
1
L −1 =1.
s
Substituting u = ts gives
Z ∞ u n du Z ∞
1 n!
L(tn ) = e−u = n +1 un e−u du = n+1 ( for s > 0)
0 s s s 0 s
Z ∞
" #∞
−st at e( a−s)t
F (s) = e e dt =
0 a−s
0
1 1 1
=− + lim e(a−s)t = .
a−s a − s t→∞ s−a
Thus, the integral exists for s > a (note F (s) does not exist for s ≤ a) and
1
L(e at ) = (s > a) .
s−a
Example 8.5. Let f (t) = sin( at) for some a 6= 0, and let F = L( f ).
Notice that for s > 0 we have lim e−st sin( at) = 0 (by the Squeeze Theo-
t→∞
rem); similarly, lim e−st cos( at) = 0. Using this and two integrations by
t→∞
mathematical methods 2 111
parts, we get
Z ∞ Z ∞
1
F (s) = e−st sin( at) dt = − (e−st )0 sin( at) dt
0 s 0
1 a ∞ −st
Z
= − [e−st sin( at)]0∞ + e cos( at) dt
s s 0
Z ∞
a
= 0− 2 (e−st )0 cos( at) dt
s 0
a a2 ∞ −st
Z
= − 2 [e−st cos( at)]0∞ − 2 e sin( at) dt
s s 0
a a2
= 2 − 2 F (s) .
s s
This gives an equation for F (s):
a a2
F (s) = 2
− 2 F (s) .
s s
It is a matter of simple algebra to rearrange to deduce that
a
L(sin( at)) = F (s) = ( s > 0).
s2 + a2
It is then immediately obvious that
a
L −1 = sin( at) .
s + a2
2
for s ≥ a.
1 t 1
cosh(t) = (e + e−t ) and sinh(t) = (et − e−t ).
2 2
112 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Theorem 8.9. Let F (s) and G (s) be functions. If the inverse Laplace
transforms f (t) = L−1 ( F )(t) and g(t) = L−1 ( G )(t) exist, then for
any constants α ∈ R and β ∈ R we have
Example 8.12. For f (t) = e at , we saw in Example 8.4 that the transform
exists for s > a. This is consistent with Theorem 8.11 since f (t) is of
exponential order ≤ a: taking M = 1 and γ = a we see that | f (t)| ≤
Me at for all t ≥ 0.
2 2
Example 8.13. For f (t) = et there are no M and γ for which et ≤
2
Meγt for all t ≥ 0. In very informal terms, et grows more quickly than
2
eγt for any γ. The Laplace transform L(et ) does not exist in this case.
This example proves that not every well-defined function necessarily has a
Laplace transform.
8.1.3 Exercises
Based on the examples above you should be able to tackle the fol-
lowing problems:
Exercise 8.1.2. Use integration by parts to show that for any constants
a > 0 and ω ∈ R we have
ω
(a) L(e− at sin(ωt)) = , for s > 0
( s + a )2 + ω 2
s+a
(b) L(e− at cos(ωt)) = , for s > 0
( s + a )2 + ω 2
(Hint: Write down the definition of the Laplace transform in each case.
A suitable substitution will reduce the integrals to those in Example 8.5
thereby circumventing the need to do pages of laborious calculation.)
Exercise 8.1.3. Use the linearity of the Laplace transform and some of the
above examples to find the Laplace transforms of:
(a) f (t) = cos t − sin t,
(b) f (t) = t2 − 3t + 5,
(c) f (t) = 3e−t + sin(6t).
Exercise 8.1.4. Use the linearity of the inverse Laplace transform and
some of the above examples to find the inverse Laplace transforms of:
2
(a) F (s) = − , s > −16
s + 16
4s
(b) F (s) = 2 , s>3
s −9
3 1
(c) F (s) = + 2 , s > 7.
s−7 s
114 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
2s − 1 A B C
F (s) = = + + .
(s − 1)(s + 1)(s + 3) s−1 s+1 s+3
2s2 − s + 4
Example 8.15. Suppose F (s) = , for s ≥ 0. To find f (t) =
s3 + 4s
− 1
L ( F ), we first use partial fractions:
2s2 − s + 4 A Bs + C
F (s) = = + 2 .
s ( s2 + 4) s s +4
mathematical methods 2 115
1
= 1 + cos(2t) − sin(2t).
2
8.2.1 Exercise
Exercise 8.2.1. Use partial fractions to find the inverse Laplace trans-
forms of:
2s
(a) F (s) = −
(s + 1)(s2 + 1)
1
(b) F (s) = 4 .
s − 16
so that
L( f 00 )(s) = s2 L( f )(s) − s f (0) − f 0 (0) .
Similarly,
Example 8.17. The above can be used to find L(sin( at)) by an alterna-
tive route than that taken in Example 8.5. Let f (t) = sin( at). Then f (t)
is continuous and of exponential order ≤ 0 (since | f (t)| ≤ 1e0t ) and f 0 (t)
is continuous, so we can apply Theorem 8.16. We have f 00 (t) = − a2 f (t)
and so
1
L( g)(s) = L( f )(s) .
s
In other words, if the Laplace transform of f (t) is F (s), then
Z t
F (s) F (s)
L( g)(s) = and L−1 = L−1 ( F )(u) du .
s s 0
8.3.1 Exercises
Exercise 8.3.1. Use the formula for the Laplace transform of a double
derivative to show that
2ωs
(a) L(t sin(ωt)) = 2
( s + ω 2 )2
s2 − ω 2
(b) L(t cos(ωt)) = 2 .
( s + ω 2 )2
Exercise 8.3.2. Use the formula for the Laplace transform of a derivative
to find:
(a) L(te at )
(b) L(tn e at )
Use the
Exercise 8.3.3. formula for
the Laplace transform of an integral
−1 1 −1 1
to find L and L (no partial fractions
s ( s + 3) s2 ( s + 3)
required).
where a, b are constants and r (t) is a given function, such that y(t)
satisfies the initial conditions
y ( 0 ) = K0 , y 0 ( 0 ) = K1 .
where R(s) = L(r )(s). In terms of Y (s) = L(y), this gives the
equation
so we have
R ( s ) + ( s + a ) K0 + K1
Y (s) = .
s2 + as + b
Therefore
−1 −1 R ( s ) + ( s + a ) K0 + K1
y(t) = L (Y ) = L .
s2 + as + b
1
( s 2 − 1 )Y = s + 1 + .
s2
Solving for Y (s) gives
s+1 1 1 1 1
Y (s) = + = + − ,
s2 − 1 s2 ( s2 − 1) s − 1 s2 − 1 s2
Hence
1 −1 1 1 1 1 1
y ( t ) = L − 1 (Y ) = L ( ) − L −1 ( ) + L −1 ( 2 )
4 s−1 4 s+1 2 s +1
1 t 1 −t 1
= e − e + sin t .
4 4 2
8.4.1 Exercises
Exercise 8.4.1. Solve the initial value problems using the Laplace trans-
form :
(a) y0 (t) − 9y(t) = t, y(0) = 5
(b) y00 (t) − 4y0 (t) + 4y(t) = cos t, y(0) = 1, y0 (0) = −1
(c) y00 (t) − 5y0 (t) + 6y(t) = e−t , y(0) = 0, y0 (0) = 2
(d) y(4) (t) − 4y(t) = 0, y(0) = 1, y0 (0) = 0, y00 (0) = −2,
y000 (0) = 0.
We have now seen the general strategy for solving differential equa-
tions using Laplace transforms; we transform the differential prob-
lem to an algebraic one for Y (s) and then, given our knowledge of
inverse Laplace transforms, we attempt to reconstruct the form of
y(t). It is this last step that is potentially the tricky one for there is
always the possibility that Y (s) is of a form we do not recognise.
The situation gets worse. It is relatively straightforward to find
the Laplace transform of a function in as much that given an f (t)
we can, at least theoretically, compute F (s) using the definition of
a Laplace transform but, unfortunately, there is no easy equivalent
definition for going in the reverse direction (ie. given F (s), deduce
f (t)). Thus it is of importance to expand our repertoire of easily
identifiable inverse functions and this is facilitated using two so-
called shift theorems.
Proof. We have
Z ∞ Z ∞
L e at f (t) = e−st e at f (t) dt = e−(s− a)t f (t)dt = F (s − a) ,
0 0
which proves the statement.
2s 2s
G (s) = = .
s2 + 2s + 5 ( s + 1)2 + 4
2( s + 1) − 2
L −1 ( G ) ( t ) = L −1
( s + 1)2 + 4
2s − 2
= e − t L −1
s2 + 4
s 1
= e − t L −1 2 2 − 2
s + 22 s2 + 22
= e−t [2 cos(2t) − sin(2t)] .
0
, t<a,
g(t)[ H (t − a) − H (t − b)] = g(t) , a ≤ t < b ,
0 , t≥b.
mathematical methods 2 123
You can try to prove this result (it involves the definition of the
Laplace transform of H and one change of coordinate).
e− as
L ( H (t − a)) = L ( f (t − a) H (t − a)) = e−as F (s) = , for s > 0.
s
Example 8.30. Find L−1 (e−4s /s3 ).
Solution: We apply the theorem with a = 4 and F (s) = 1/s3 , so that
L−1 (e−4s /s3 ) = f (t − 4) H (t − 4). All we have to do is determine f (t).
From the table we get that f (t) = t2 /2, so that,
( t − 4) 2 0 if t < 4 ,
L−1 e−4s /s3 = H (t − 4) =
2 (t − 4)2 /2 if t ≥ 4 .
This is still not in the form required in order to use Theorem 8.28: in
the second term we have to write 1 − 4t as a function of t − 3. Since
t = (t − 3) + 3, we have 1 − 4t = 1 − 4(t − 3) − 12 = −4(t − 3) − 11.
Thus
g(t) = t − [4(t − 3) + 11] H (t − 3) .
We now apply Theorem 8.28 with f (t) = 4t + 11.
4 11
L([4(t − 3) + 11] H (t − 3)) = L( f (t − 3) H (t − 3))(s) = e−3s F (s) = e−3s ( + ).
s2 s
124 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Thus
1 4 11
L( g)(s) = 2
− e−3s ( 2 + ).
s s s
We now solve a differential equation where the right-hand side
uses Heavside functions.
y00 + y = H (t − 1) − H (t − 2)
e−s e−2s
s2 Y (s) − sy(0) − y0 (0) + Y (s) = −
s s
and solving for Y yields that
1 1
Y (s) = + (e−s − e−2s ) 2
s2 + 1 s ( s + 1)
1 −s −2s 1 s We used partial fractions here.
= + ( e − e ) −
s2 + 1 s s2 + 1
1 −s 1 s −2s 1 s
= + e − − e − .
s2 + 1 s s2 + 1 s s2 + 1
Hence
sin t,
0 ≤ t < 1,
y(t) = sin t + 1 − cos(t − 1), 1≤t<2
sin t − cos(t − 1) + cos(t − 2), 2 ≤ t.
8.5.1 Exercises
Exercise 8.5.1. Find the Laplace transform of the functions
(a) (t3 − 3t + 2)e−2t
(b) e4t (t − cos t)
(c) f (t) = 2t + 1 for 0 ≤ t < 2 and f (t) = 2 − 3t for t > 2.
Consequently,
L −1 F 0 ( s ) = − t f ( t )
( t ≥ 0).
8.7 Convolution
f ∗g = g∗ f commutative
f ∗ ( g1 + g2 ) = f ∗ g1 + f ∗ g2 distributive
( f ∗ g) ∗ h = f ∗ ( g ∗ h) associative
f ∗0 = 0∗ f = 0
This result tells us that if f (t) and g(t) have Laplace transforms
F (s) and G (s) respectively then the function ( f ∗ g)(t) has Laplace
transform F (s) G (s).
1
K (s) =
( s − 1)2 ( s − 3)2
using the Convolution Theorem. Another method would be to use
partial fractions as in Section 8.2.
Solution: Notice that, by Example 8.25,
−1 1 t −1 1
L = te = f (t) , L = te3t = g(t) .
( s − 1)2 ( s − 3)2
Therefore by the Convolution Theorem,
−1 −1 1 1
L (K ) = L · = f (t) ∗ g(t)
( s − 1)2 ( s − 3)2
Z t Z t
= f (u) g(t − u) du = ueu (t − u)e3(t−u) du
0 0
Z t
= e3t (tu − u2 )e−2u du .
0
1 1
Solution: Since F (s) = · and L−1 ( s2 1+1 ) = sin(t), by
( s2 + 1) ( s2 + 1)
the Convolution Theorem,
Z t
f (t) =L−1 ( F ) = sin(t) ∗ sin(t) = sin(u) sin(t − u) du
0
1 t
Z
= [cos(u − (t − u)) − cos(u + (t − u))] du (using the cosine sum formula)
2 0
Z t
1
= [cos(2u − t) − cos t] du
2 0
t
1 1
= sin(2u − t) − u cos t
2 2 0
1 1
= sin t − t cos t.
2 2
Exercise 8.7.1. In the same manner as the previous example, show that
s2
−1 1 1
L = sin t + t cos t .
( s2 + 1)2 2 2
1 1
Y (s) = + (1 − e − s ) 2
s2 + 1 s ( s + 1)
1 1 s e−s 1
= 2 + − 2 − ·
s +1 s s +1 s s2 + 1
Taking the inverse Laplace transform of this and using the Convolution
Theorem for the last term, we get
8.7.1 Exercises
Exercise 8.7.2. Find the inverse Laplace transform of the functions using
the Convolution Theorem:
1
(a) 2
(s + 4)(s2 − 4)
s
(b) 2
(s + a2 )(s2 + b2 )
e−2s
(c) 2 .
s + 16
Exercise 8.7.3. Solve the following initial value problems using Laplace
transforms:
(a) y00 (t) + 4y0 (t) + 13y(t) = f (t), y(0) = y0 (0) = 0, where f (t) = 1
for 0 ≤ t < π and f (t) = 0 for t ≥ π.
(b) y00 (t) + 2y0 (t) + 2y(t) = sin t, y(0) = y0 (0) = 0.
mathematical methods 2 129
1 e− as
1 H (t − a)
s s
1 t n −1
, n ∈ Z+ e− as F (s) f (t − a) H (t − a)
sn ( n − 1) !
1
e− at F (s − a) e at f (t)
s+a
1 t n −1
, n ∈ Z + e− at sF (s) − f (0) f 0 (t)
(s + a)n ( n − 1) !
1 sin(ωt)
s2 F ( s ) − s f (0) − f 0 (0) f 00 (t)
s2 + ω 2 ω
s
cos(ωt) F 0 (s) −t f (t)
s2 + ω 2
1 e−at sin(ωt)
F (n) ( s ) (−t)n f (t)
( s + a )2 + ω 2 ω
Z t
s+a F (s)
e− at cos(ωt) f (u) du
( s + a )2 + ω 2 s 0
1 sin(ωt) − ωt cos(ωt)
F (s) G (s) f ∗ g (t)
( s2 + ω 2 )2 2ω 3
s t sin(ωt)
( s2 + ω 2 )2 2ω
Higher derivatives
L f ( n ) ( t ) = s n F ( s ) − s n −1 f (0 ) − s n −2 f 0 (0 ) − · · · − s f ( n −2) (0 ) − f ( n −1) (0 )
• Hooke’s Law: The elastic force exerted by the free end of the
spring is proportional to the displacement of that end. The factor
k > 0 of proportionality is called the spring constant.
which is equivalent to
where F (s) = L( f )(s). Notice that we can write f (t) in terms of Heavi-
side functions f (t) = F0 [1 − H (t − 1)], so
1 e−s
F (s) = F0 [L(1) − L( H (t − 1))] = F0 − .
s s
Hence
e−s
F (s) 1
X (s) = 2 = F0 − .
s +4 s ( s2 + 4) s ( s2 + 4)
and therefore
1 1 1 s
L −1 = L −1 − L −1 2
s ( s2 + 4) 4 s s +4
1
= (1 − cos 2t).
4
Next, since L−1 (s/(s2 + 4)) = cos 2t, using the (second) Shifting
Theorem, we get
e−s 1 −1 e − s se−s
−1
L = L − 2
s ( s2 + 4) 4 s s +4
1 1
= H (t − 1) − H (t − 1) cos 2(t − 1).
4 4
Hence
e−s
1
x (t) = L−1 ( X (s)) = F0 L−1 − F0 L −1
s ( s2 + 4) s ( s2 + 4)
1 F F
= F0 − cos 2t − 0 H (t − 1) + 0 H (t − 1) cos 2(t − 1) .
4 4 4
132 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
• Kirchoff’s Voltage Law: the algebraic sum of the voltage rises and
drops in any closed loop in a circuit is equal to zero.
d2 q dq 1
L 2
+ R + q = E(t) . (8.4)
dt dt C
On differentiating both sides of this equation, one derives a second
order equation for the current i (t):
d2 i di 1 dE
L 2
+R + i= (t) .
dt dt C dt
mathematical methods 2 133
10e−2s 10e−3s
e(s) = L( E)(s) = 10[L( H (t − 2)) − L( H (t − 3))] = − .
s s
10e−2s 10e−3s
Thus, 106 (s/4 + 1) Q(s) = + , and therefore
s s
4 4
Q(s) = 10−5 e−2s − e−3s
s ( s + 4) s ( s + 4)
Using partial fraction (or a direct observation), one gets
4 1 1
= − .
s ( s + 4) s s+4
This and the (second) Shift Theorem imply
−1 −5 −1 4 −2s −1 4 −3s
q(t) = L ( Q(s)) = 10 L e −L e
s ( s + 4) s ( s + 4)
−2s
e−2s e−3s e−3s
e
= 10−5 L−1 − − L −1 −
s s+4 s s+4
= 10−5 [ H (t − 2) − e−4(t−2) H (t − 2) − H (t − 3) + e−4(t−3) H (t − 3)] .
Output Voltage
Figure 8.3: Output voltage for Example
10 8.41
9
0
0 1 2 3 4 5 6
m1 m2
f1→ f2→
friction). For S1 the restoring force at time t is −k1 x1 (t), while for
S2 the restoring force acting on m1 at time t is k2 [ x2 (t) − x1 (t)]. The
damping force is −c1 x10 (t) for some (damping) constant c1 ≥ 0.
Hence
Given some initial conditions for x1 and x2 (i.e. xi (0) and xi0 (0) for
i = 1, 2), the system can be solved for x1 (t) and x2 (t). One possible
way of doing this is by applying the Laplace transform to each
differential equation, thus obtaining a system of algebraic equations
for X1 (s) = L( x1 (t)) and X2 (s) = L( x2 (t)).
Example 8.42. In the model considered above, assume that the friction is
negligible and the system is at rest initially (i.e. both masses are at their
equilibrium positions). From time t = 0 until time t = 3 seconds a force
of magnitude 2 N acts on the mass m2 in direction to the right. Find the
positions x1 (t) and x2 (t) at any time t ≥ 0, provided m1 = m2 = 1 kg
and all spring constants are the same: k1 = k2 = k3 = k > 0.
If we denote
2 2e−3s
F2 (s) = L( f 2 (t)) = 2L(1) − 2L( H (t − 3)) = − ,
s s
then applying the Laplace transform to each of the differential equations in
the above system leads to the algebraic equations
s2 X (s) + 2kX (s) − kX (s) = 0
1 1 2
(8.6)
s2 X2 (s) + 2kX2 (s) − kX1 (s) = F2 (s) .
Solving the system with respect to X1 (s) and X2 (s) (e.g. by expressing
X2 (s) in terms of X1 (s) from the first equation and then substituting it
into the second), one finds
k s2 + 2k
X1 ( s ) = F2 ( s ) , X 2 ( s ) = F2 (s) .
(s2 + 2k)2 − k2 (s2 + 2k)2 − k2
136 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Thus,
2 2e−3s
1 1 1
X1 ( s ) = − −
2 s2 + k s2 + 3k s s
1 e−3s 1 e−3s
= − − + .
s ( s2 2 2 2
+ k) s(s + k) s(s + 3k) s(s + 3k)
Using partial fractions, one gets
1 1 1 s 1 1 1 s
= − 2 , = − 2 .
s ( s2 + k ) k s s +k s(s2 + 3k ) 3k s s + 3k
Thus,
1 e−3s s e−3s
1 1 s
X1 ( s ) = − 2 − − 2
k s s +k k s s +k
−3s
s e−3s
1 1 s 1 e
− − 2 + − 2 .
3k s s + 3k 3k s s + 3k
Applying the inverse Laplace transform and using the Shift Theorem, one
finds
2 1 √ 1 √
x1 ( t ) = [1 − H (t − 3)] − cos( kt) + H (t − 3) cos( k(t − 3))
3k k k
1 √ 1 √
+ cos( 3kt) − H (t − 3) cos( 3k(t − 3)) .
3k 3k
Similarly, for X2 (s) we have
2 2e−3s
1 1 1
X2 ( s ) = + −
2 s2 + k s2 + 3k s s
1 1
= (1 − e−3s ) + (1 − e−3s )
s ( s2 + k ) s(s2 + 3k )
1 1 s 1 1 s
= − 2 (1 − e−3s ) + − 2 (1 − e−3s ) ,
k s s +k 3k s s + 3k
and therefore
4 1 √ 1 √
x2 ( t ) = [1 − H (t − 3)] − cos( kt) + H (t − 3) cos( k(t − 3))
3k k k
1 √ 1 √
− cos( 3kt) + H (t − 3) cos( 3k(t − 3)) .
3k 3k
Example 8.43. Next, consider an electric circuit containing two closed
loops as in the figure.
Here R1 and R2 are the resistances of the corresponding resistors, q1 (t)
and q2 (t) the charges of the corresponding capacitors, i1 (t), i2 (t) and
i (t) are the currents in the corresponding parts of the circuit, L is the
inductance of the single inductor, and E(t) is the electromotive force.
Using Kirchoff’s current law for the first junction, one sees that i1 −
i2 − i = 0, i.e. i = i1 − i2 . (The same result follows if the current law is
applied at the second junction.)
mathematical methods 2 137
R1 R2
E L
C1 C2
Next, applying Kirchoff’s voltage law to the left loop, it follows that
di 1
R1 i1 ( t ) + L ( t ) + q1 ( t ) = E ( t ) .
dt C1
1
L[q100 (t) − q200 (t)] + R1 q10 (t) + q (t) = E(t) . (8.7)
C1 1
Similarly, it follows from Kirchoff’s voltage law applied to the right loop
that
1
− L[q100 (t) − q200 (t)] + R2 q20 (t) + q2 (t) = 0 . (8.8)
C2
The equations (8.7) and (8.8) form a system of differential equations that
determine q1 (t) and q2 (t) provided some initial conditions for these are
given.
One can deal in a similar way with other electric circuits. Some-
times it is more convenient to consider differential equations in-
volving the currents (i1 and i2 in the case above).
i1→ i2→
↓ i1−i2
E R1=200Ω R2=300Ω
L=0.5H C=5×105F
Using this and Kirchoff’s laws, we obtain the following system of equa-
tions:
1
i10 (t) + 200[i1 (t) − i2 (t)] = 50
2 Z t
1
300i2 (t) + 200[i2 (t) − i1 (t)] +
i2 (u) du = 0 .
50 × 10−6 0
Applying the Laplace transform to each equation, one finds the following
system for I1 (s) = L(i1 (t)) and I2 (s) = L(i2 (t)):
100
sI1 (s) + 400I1 (s) − 400I2 (s) =
s
5I2 (s) − 2I1 (s) + 200 I2 (s) = 0 .
s
Solving, one finds I1 (s) and I2 (s), and then i1 (t) = L−1 ( I1 (s)) and
i2 (t) = L−1 ( I2 (s)). We leave the details to the reader.
8.8.4 Exercises
Exercise 8.8.1. A body of mass m is fixed at the end of a vertical spring
with spring constant k > 0. Initially the mass is at the equilibrium po-
sition. At time t = 0 the upper end of the spring starts to be moved
vertically by means of a rotating eccentric or cam with the periodic dis-
placement y(t) = a sin ωt. Assume that the air resistance is negligible
and that y(t) and the displacement x (t) of the mass at time t are measured
from their initial positions and both x and y are positive for displacements
downward. Find the position x (t) of the mass at any time t ≥ 0.
Exercise 8.8.2. Consider the electric circuit in the figure below. Assume
that the capacitor has zero charge initially and that there is zero initial
current. At time t = 2 seconds the switch is moved from position B to
position A, held there for 1 second, and then switched back to position B.
Find the current i (t) at any time t ≥ 0.
m1=1 m2=1
f1→ f2→
M m
f→
c1
(a) Derive and solve the system of equations for the displacement func-
tions of M and m.
√
(b) Show that if m and k2 are chosen so that ω = k2 /m, then the
mass m cancels the forced vibrations of M. (In this case m is called a
vibration absorber.)
Exercise 8.8.5. Find the currents i1 (t) and i2 (t) for the electric circuit
in the figure, assuming that E(t) = 1 − H (t − 4) sin[2(t − 4)] and the
initial charges and currents are zero.
i1→ i2→
R1=2Ω R2=1Ω
E L=5H
R3=3Ω R4=4Ω
9
Complex Functions — Derivatives
p r
where r = |z| = x2 + y2 is called the modulus of z and θ = arg(z)
its argument. This θ is the angle of the radius vector to the point
( x, y) in the argand diagram and, by convention, is usually chosen θ
x
to be in the interval (−π, π ].
√ √
Example 9.1. i = 1 cis π/2 = eiπ/2 and 1 + i = 2 cis π/4 = 2eiπ/4 .
So z1/m does not take a unique value but actually has up to m dif-
ferent values! Hence simply writing f (z) = z1/m does not immedi-
ately give us a well-defined function although this can be done by
exercising a little more care. We can make z1/m a genuine function
by specifying just one of the roots, say the first, k = 0. This might
seem a little strange but we should recognise this type of behaviour
from a simple discussion of the square root of a positive real num-
ber. Everyone knows that the square root of 4 is 2; but the roots
of the equation x2 = 4 are 2 and −2. Thus if z is a real positive
number then z1/2 has two possible values but we almost always
identify the square root of a positive number uniquely by specify-
ing the square root to be positive as well. This makes the square
root function well-defined and it precisely this sort of restriction
we are imposing when we take one particular root of the complex
142 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
log 0 does not exist). But we have just seen that ez , while defined
for all z ∈ C is not one-to-one. This is a similar issue as when we
wanted to find inverse functions for sin and cos in Chapter 1. We
will use here the same trick as we did then: by restricting the range
of log z we can form a well-defined complex function. The standard
range of log z is taken to be
{ x + iy | x ∈ R, −π < y ≤ π }
that is the real part of log z is unrestricted but the imaginary part
is confined to lie in the interval (−π, π ]. Confining the range of
a complex function to something smaller than the whole of the
complex plane is called choosing a branch of the function.
Let w = x + iy and suppose ew = z for z ∈ C. We have
z = ew = e x cos y + i sin y ,
ranges.
Given this definition, we are now in a position to define ab for
a, b ∈ C, a 6= 0. We have a = elog a from the above definition, so
ab = eb log a . Note that in general this is multi-valued, because we
can choose different branches of the log function.
Thus
cos2 z + sin2 z = 1,
• zb = eb log z , for z 6= 0
• sin z = 1 iz
2i ( e − e−iz )
f : R2 → R2 : ( x, y) 7→ (u, v),
D ( z0 , r ) = { z ∈ C : | z − z0 | < r }.
lim f (z) = a
z → z0
means that we can make the values of f (z) arbitrarily close to a by taking More precisely, for the interested
student: for every e > 0, there is
z to be sufficiently close to z0 , but not equal to z0 . This means that f (z) a δ > 0 such that z ∈ D (z0 , r ),
is close to a when z is close to z0 in all directions. Note that the limit a is z 6= z0 , and |z − z0 | < δ imply that
| f (z) − a| < e.
unique.
The following equations hold (as for real functions of real vari-
ables). Let limz→z0 f (z) = a and limz→z0 g(z) = b, then
1. limz→z0 f (z) ± g(z) = a ± b.
2. limz→z0 f (z) g(z) = ab.
3. limz→z0 f (z)/g(z) = a/b, if b 6= 0.
f ( z ) − f ( z0 )
lim
z → z0 z − z0
exists. This limit is denoted f 0 (z0 ) (or sometimes d f /dz (z0 )) and is a
complex number. When f is complex differentiable at every point in A,
we say f is analytic on A. Note that the phrase “analytic at z0 ” means
analytic in a neighbourhood of z0 .
d exp(z)
= exp(z),
dz
d sin(z) d cos(z)
= cos(z), = − sin(z).
dz dz
Example 9.9. The logarithmic function log z is analytic on any open set
not containing any real negative numbers, and
d log(z)
= 1/z.
dz
Theorem 9.10. Differentiability implies continuity. If f 0 (z0 ) exists then
f is continuous at z0 .
( a f + b g ) 0 ( z ) = a f 0 ( z ) + b g 0 ( z ).
( f g ) 0 ( z ) = f 0 ( z ) g ( z ) + f ( z ) g 0 ( z ).
−1
Similarly f /g = f (z) g(z) , where g(z) 6= 0 ∀z, is analytic on
A, and the quotient rule holds:
−1 −2 0
( f /g)0 (z) = f 0 (z) g(z) − f (z) g(z) g (z)
f 0 (z) g(z) − f (z) g0 (z)
= 2 .
g(z)
146 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
a0 + a1 z + · · · + a n z n
b0 + b1 z + · · · + bm zm
is analytic on C \ {zeroes of the denominator} (there are at most
m zeroes).
d
( h ◦ f )(z) = h0 f (z) f 0 (z).
dz
f ( x + iy) = u( x, y) + iv( x, y)
∂u ∂v ∂u ∂v
= e x cos y = and = −e x sin y = − .
∂x ∂y ∂y ∂x
We now have the two expressions (♥) and (♠) for the value of
f 0 (z0 ) depending on how we allow z → z0 . But from the definition
of complex differentiability, if the function f (z) is differentiable at
z = z0 then the values of f 0 (z0 ) must be the same however we let
z → z0 . Hence the two forms (♥) and (♠) must be identical so
comparing the real and imaginary parts gives us
∂u ∂v ∂v ∂u
( x0 , y0 ) = ( x0 , y0 ) and ( x0 , y0 ) = − ( x0 , y0 ).
∂x ∂y ∂x ∂y
1. Let z0 = x0 + iy0 .
If f 0 (z0 ) exists then f is differentiable in the sense of real variables
at ( x0 , y0 ) and u, v satisfy
∂u ∂v ∂u ∂v
( x0 , y0 ) = ( x0 , y0 ) and ( x0 , y0 ) = − ( x0 , y0 ).
∂x ∂y ∂y ∂x
∂u ∂v ∂v ∂u ∂u ∂u ∂v ∂v
f 0 ( z0 ) = +i = −i = −i = +i .
∂x ∂x ∂y ∂y ∂x ∂y ∂y ∂x
−uv x + vu x = 0. (3)
( u2 + v2 ) u x = 0 and (u2 + v2 )v x = 0.
∂2 u ∂2 u
+ 2 = 0;
∂x2 ∂y
in other words the real part of an analytic function satisfies Laplace’s
equation. On the other hand, we can also differentiate the two
Cauchy–Riemann equations with respect to y and x in turn to
obtain u xy = vyy and vyy = −u xy . Now if we eliminate u xy we
conclude that
∂2 v ∂2 v
+ =0
∂x2 ∂y2
as well.
Arbitrary sums of these functions are also solutions, because the
Laplace equation is linear and homogeneous.
Example 9.14. Show that the real and imaginary parts of z2 satisfy
Laplace’s equation.
Solution: As z2 = ( x + iy)2 = x2 − y2 + 2ixy we have that u( x, y) =
x2 − y2 and v( x, y) = 2xy. Simple differentiation gives that u xx = 2 and
150 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Exercise 9.5.2. Show that the real and imaginary parts of the function
f (z) = z3 satisfy Laplace’s equation.
Theorem 9.15. Both the real and imaginary parts of any analytic
function, considered as functions of x and y, are solutions to Laplace’s
equation:
∂2 u ∂2 u ∂2 v ∂2 v
+ 2 =0 and + 2 = 0.
∂x2 ∂y ∂x 2 ∂y
x
Example 9.16. Show that u( x, y) = x 2 + y2
satisfies Laplace’s equation
for ( x, y) 6= (0, 0). Find the corresponding analytic f (z) having real part
u( x, y).
x
Solution: If u = x 2 + y2
then elementary differentiation gives
∂u y2 − x 2 ∂2 u 2x3 − 6xy2
= 2 ⇒ = 2
∂x ( x + y2 )2 ∂x 2 ( x + y2 )3
and
∂u 2xy ∂2 u −2x3 + 6xy2
=− 2 ⇒ = .
∂y ( x + y2 )2 ∂y2 ( x 2 + y2 )3
As u xx = −uyy it follows that u xx + uyy = 0 so u( x, y) satisfies Laplace’s
equation.
To find the complex function f = u( x, y) + iv( x, y), we attempt to find
v( x, y) so that the Cauchy-Riemann equations are satisfied. This requires
∂v ∂u y2 − x 2
= = 2 (♣)
∂y ∂x ( x + y2 )2
and
∂v ∂u 2xy
=− = 2 (♦).
∂x ∂y ( x + y2 )2
The method to solve this problem is similar to the method for finding
a potential, seen in Chapter 6. Notice that we will determine v up to a
constant, that we can then choose. If we integrate (♦) with respect to x we
have that
y
v( x, y) = − 2 + Q(y)
x + y2
mathematical methods 2 151
∂v y2 − x 2 dQ ∂u y2 − x 2
= 2 + = = .
∂y ( x + y2 )2 dy ∂x ( x 2 + y2 )2
x − iy
f = u + iv =
x 2 + y2
and if z = x + iy then
z̄ z̄ 1
f (z) = = = .
| z |2 zz̄ z
Hence the differentiable function 1/z (for z 6= 0) is analytic with real part
the given u( x, y).
P ( A1 ∪ A2 ∪ . . . ∪ A n ) = P ( A1 ) + P ( A2 ) + . . . + P ( A n ).
Axiom 3 is more subtle than the first two, and is known as the
additivity property of probability. It says we can calculate probabili-
ties of complicated events by adding up the probabilities of smaller
events provided the smaller events are disjoint and together make
up the entire complicated event. When we say disjoint we mean the
events do not intersect. This axiom can be extended to a countable
sequence of disjoint events A1 , A2 , . . . and is needed in general but
not for this unit.
Axioms 1-3 imply other basic properties or theorems that are
true for any probability function. We list four now.
Theorem 10.7. P(∅) = 0.
Proof. S ∪ ∅ = S and thus P(S ∪ ∅) = P(S) = 1. But S ∩ ∅ = ∅ so
that P(S ∪ ∅) = P(S) + P(∅) = 1 + P(∅) by Axioms 3 and 1. Thus
1 + P(∅) = 1 and P(∅) = 0.
P( A ∪ ( A ∩ B)) = P( A) + P( A ∩ B)
= P( A) + P( B) − P( A ∩ B)
P ( A ∪ B ) = P ( A ) + P ( B ) − P ( A ∩ B ).
1. A1 ∪ A2 ∪ . . . ∪ Ak = S
2. Ai ∩ A j = ∅ for all i 6= j
P( B| A j ) P( A j )
P( A j | B) = , j = 1, 2, . . . , k.
∑ik=1 P ( B | Ai ) P ( Ai )
P ( B | A1 ) P ( A1 )
P ( A1 | B ) =
P ( B | A1 ) P ( A1 ) + P ( B | A2 ) P ( A2 )
30/100 · 1/3
=
30/100 · 1/3 + 2/100 · 2/3
= 30/34 ≈ 0.882.
Example 10.18. We throw a red die and a blue die. The outcome of the
experiment is a pair of numbers ( a, b), where a is the number on the red
die and b the number on the blue die. The sample space S is the set
A( x ) = {s : X (s) = x } ⊆ S.
p X ( x ) ≥ 0 and ∑ p X ( x ) = 1.
all x
From now on we shall write P( X = x ) for P( X (s) = x ) and
suppress the functional dependence of X on the elements S un-
less necessary. We also define the cumulative distribution function
(c.d.f.).
FX (t) = P( X ≤ t) = ∑ p X ( x ).
x ≤t
1. 0 ≤ FX (t) ≤ 1
3. limt→∞ FX (t) = 1
4. limt→−∞ FX (t) = 0
160 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
E[ X ] = ∑ xp X ( x ).
range of X
E[ H ( X )] = ∑ H ( x ) p X ( x ).
range of X
1. E[c] = c,
2. E[cX ] = cE[ X ],
3. E[ X + Y ] = E[ X ] + E[Y ].
Var( X ) = E[( X − µ X )2 ] = ∑ ( x − µ X )2 p X ( x ),
range of X
p
and the standard deviation is defined by σX = Var( X ).
Intuitively, σX2 and σX are measures of how spread out the distri-
bution of X is or how much it varies. As a measure of variability
the variance is not so intuitive because it is measured in different
units than the random variable, but it is convenient mathemati-
cally. Therefore, the standard deviation σX is also defined and is the
square root of the variance, providing a measure of variability in
the same units as the random variable. The variance of a random
variable plays a very important role in probability theory and statis-
tical inference so we pause to present some properties of variances
(without proof).
P( X = x and Y = y) = P( X = x ) P(Y = y)
for all x, y.
1. Var( X ) ≥ 0.
3. Var( X ) = E( X 2 ) − µ2X = E[ X 2 ] − E[ X ]2 .
E [Y ] = 1 · p + 0 · q = p
(clear).
1.0
● ● ● ● ● ● ●● ●● ● ● ● ● ●
●
●
●
●
0.20
0.8 ●
●
0.15
0.6
●
p
y
0.10
0.4
● ●
●
0.05
0.2
●
●
0.00
●
0.0
●
● ● ● ● ●
0 5 10 15 20 0 5 10 15 20
x x
From Theorem 10.24, Theorem 10.27 and the fact that the Bernoulli
trials are independent it can be shown that if X ∼ Binom(n, p) then
µ X = np and the variance is σX2 = npq by
" #
n n n
µ X = E[ X ] = E ∑ Yi = ∑ E[Yi ] = ∑ p = np.
i =1 i =1 i =1
" #
n n n
σX2 = Var ∑ Yi = ∑ Var[Yi ] = ∑ pq = npq.
i =1 i =1 i =1
To summarise:
Name p.m.f. ( Mean Variance
p y q 1− y y = 0, 1.
Y ∼ Bern( p) PY (y) = p pq
0 otherwise.
(
(nx) p x qn− x x = 0, 1, 2, . . . , n.
X ∼ Binom(n, p) PX ( x ) = np npq
0 otherwise.
are working the next morning the factory can’t meet customer demand
and incur severe financial losses. How many machines do we expect to
be working the next day? What is the probability that at least 9 out of 10
machines will work? Provide a measure of variability (standard deviation)
for how many machines are still operating the next day.
Solution: Define X to be the number of machines working the next day
out of a total of 10 machines. Then X ∼ Binom(10, 0.95). We would
therefore expect E[ X ] = 10 × 0.95 = 9.5 machines to be operating the
next day. The probability that at least 9 machines out of 10 are working
the next day is
10 10
P( X ≥ 9) = p X (9) + p X (10) = 0.959 0.051 + 0.9510 0.050
9 10
≈ 0.3151 + 0.5987 = 0.9138.
Observe that the ratio is close to 3 for small p and close to 1 for large p
illustrating the general principle of reliability: that redundancy improves
system reliability when components are ‘unreliable’ but there is little
advantage in having redundancy when the components are highly reliable.
by
Z t
FX (t) = P( X ≤ t) = f X ( x )dx.
−∞
Note that Theorem 10.24 and Theorem 10.27 hold also for contin-
uous random variables.
density function (
λe−λx x>0
f T (x) = (10.1)
0 x≤0
Clearly, f T ( x ) ≥ 0 for all x and
Z ∞ Z ∞
f T ( x )dx = λe−λx dx = [−e−λx ]0∞ = 1
−∞ 0
0.8
0.6
1.0
F
f
0.4
0.5
0.2
0.0
0.0
0 1 2 3 4 5 0 1 2 3 4 5
t t
so
2 1 1
Var( T ) = E[ T 2 ] − E[ T ]2 = 2
− 2 = 2
λ λ λ
and
1
q
σT = Var( T ) = .
λ
To summarise:
For T ∼ Exp(λ):
(
λe−λx x>0
f T (x) =
0 x≤0
(
1 − e−λt t>0
FT (t) =
0 t≤0
1 1
E[ T ] = and Var( T ) = 2 .
λ λ
Example 10.38. Past experience tells us the life time in hours of a certain
type of lightbulb produced by a manufacturer follows an Exp(λ) distribu-
tion, where λ = 0.001. How many hours do we expect a light bulb to last?
What is the standard deviation of the life time of a lightbulb? What is the
probability that a light bulb lasts more than 900 hours?
Solution: Let T denote the life time in hours of a light bulb. We know
T ∼ Exp(λ) and λ = 0.001. Therefore, the expected lifetime is 1/λ =
1/0.001 = 1000 hours. Also, the standard deviation is 1/λ = 1000
hours.
Finally, because T ∼ Exp(0.001)
Therefore the probability a randomly selected light bulb lasts longer than
900 hours is
1 2 /2σ2
f X (x) = √ e−( x−µ)
σ2 2π
1.0
0.8
0.6
0.6
0.4
f
0.4
0.2
0.2
0.0
0.0
-4 -2 0 2 4 -4 -2 0 2 4
x x
Theorem 10.40 and Corollary 10.41 say that we can write proba-
bility statements about any normally distributed random variable
in terms of the standard normal c.d.f. Section 10.3.3 at the end of
this chapter contains approximations of the P( Z < z) = FZ (z) for
z > 0 for a standard normal random variable. Some equalities will
be useful to remember: P( Z > z) = 1 − P( Z < z) and if −z < 0,
P ( Z < − z ) = P ( Z > z ) = 1 − P ( Z < z ).
P( Z > zα ) = α.
Exercise 10.2.3. Let Z ∼ N (0, 1). Deduce from the standard normal
table the (1 − α)-quantiles for α = 0.05, α = 0.025 and α = 0.0005.
a−µ
a) = 0.01. We know P( X > a) = P( Z > σ ). We also know from the
table that P( Z > 2.33) ≈ 0.01 so
a − 75
≈ 2.33
0.1
and a ≈ 75.233. Thus if 1% of screws get rejected then the smallest to be
rejected would be approximately 75.233 mm.
10.3.1 Estimation
A very basic concept to statistical inference is that of a random sam-
ple. By sample we mean only a fraction or portion of the whole –
crash testing every car to inspect the effectiveness of airbags is not
a good idea. By random we mean essentially that the portion taken
is determined non-systematically. Heuristically, it is expected that a
haphazardly selected random sample will be representative of the
whole population because there is no selection bias. For the pur-
poses of this course (and often in practice), we will assume every
observation in a random sample is generated by the same probabil-
ity distribution and each observation is made independently of the
others. This motivates the following definition
Definition 10.46. A collection of random variables X1 , X2 , . . . , Xn is
independent and identically distributed (or i.i.d.) if the collection
is independent and if each of the n variables has the same distribution.
The i.i.d. collection X1 , X2 , . . . , Xn is called a random sample from the
common distribution.
Definition 10.47. Any function of the elements of a random sample,
which does not depend on unknown model parameters, is called a statistic.
It is itself a random variable.
172 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
x n = ( x1 + x2 + . . . + xn )/n
2.58 2.58 1.75 0.53 3.29 2.04 3.46 2.92 3.10 2.41
3.89 1.99 0.74 1.59 0.35 0.03 0.52 1.42 0.04 4.02
2.0
creases.
The above outlines some nice properties of the estimator of µ,
1.5
X n , if the random sample is generated from a normal distribution.
What if we wanted to construct an estimator for the mean from
a non-normal random sample? Now the Central Limit Theorem
f(x)
1.0
comes to the rescue. Before we can state the theorem, we need the
following definition.
0.5
Definition 10.50. We say that a sequence Y1 , Y2 , . . . of random variables
converges in distribution to the random variable Y if
0.0
−3 −2 −1 0 1 2 3
lim FYn ( x ) = FY ( x )
n→∞ x
Xn − µ
Zn = √ , n = 1, 2, 3, ...
(σ/ n)
where
n
1
Xn =
n ∑ Xi .
i =1
Loosely, the CLT states that for any random sample as n gets
larger the distribution of the sample average X n , properly nor-
malised, approaches a standard normal distribution. Even more
˙ N (µ X , σX2 /n) where ∼
loosely, the CLT states, for large n, X n ∼ ˙
denotes approximately distributed. The CLT is one of the most
powerful results presented in probability theory. Obviously, know-
ing the distribution of the sample average is important because
it provides a measure of how precise we believe our estimate of
the model parameter to be. However, more powerfully, the CLT
means that for any random variable we can always perform statis-
tical inference on the mean or expected value. Figure 10.4 shows
an approximate p.d.f. of an exponential random variable and the
standardised estimator X n based on random samples of increasing
size from this random variable. The figure shows that although the
original density of each Xi is exponential, the p.d.f. of the standard-
ised estimator approaches the standard normal distribution as n
increases.
mathematical methods 2 175
0.5
approximate p.d.f. of Zn when n = 2 is
shown in the upper right panel, when
n = 20 in the bottom left panel and
1.5
0.4
when n = 2000 in the bottom right
panel. The thick black line shows the
p.d.f. of the standard normal random
0.3
1.0
fZn
fX
0.2
p.d.f.
0.5
0.1
0.0
0.0
0 1 2 3 4 5 6 −2 0 2 4 6 8 10 12
x zn
0.4
0.4
0.3
0.3
0.2
0.2
fZn
fZn
0.1
0.1
0.0
0.0
−2 0 2 4 6 −4 −2 0 2 4
zn zn
P ( L1 < θ < L2 ) = 1 − α
σ σ
P( X n − 1.96 √ ≤ µ ≤ X n + 1.96 √ ) ≈ 0.95.
n n
176 a. bassom, e. cripps, a. devillers, l. jennings, a. niemeyer, t. stemler, l. stoyanov
Notice the confidence interval does not extend below 0.05 and we should
probably doubt the producer’s claim that no more than 5% of their bullets
misfire.
0.05
0.00
xn
• For x n = 1012
• For x n = 1002
0.05
0.00
xn
mathematical methods 2 179
The discussion to this point in the present section has been our
motivation for testing hypotheses; now we formalise some ideas.
By a statistical hypothesis we mean a statement about the value of a
specified parameter in a probability model. For example, based on
the probability model described by
X n ∼ N (µ, σ2 /n)
P-value method
Based on the observed test statistic, z, we determine the P-value,
this being the probability of obtaining a value of the test statistic
at least as extreme as that observed. In determining the P-value we
use one or both tails of the distribution, depending on whether the
alternative hypothesis is one-sided (P( Z > z) or P( Z < z)) or
two-sided (P(| Z | > z)), respectively.
H0 : p = 0.05
H1 : p > 0.05
X −p
qn .
p (1− p )
n
X − 0.05
q n ∼ N (0, 1)
0.05(1−0.05)
60
1000
50
because the sample size is large enough for the CLT to apply.
40
4. Calculate the P-value based on the data
The observed value of the estimator is x n = 70/1000 = 0.07 and
30
n
fX
the observed test statistic is
0.07 − 0.05
20
p ≈ 2.9019.
0.05(1 − 0.05)/1000
10
Hence from tables at the end of this chapter the P-value is
●
0
P( X n > 0.07| p = 0.05) = P( Z > 2.9019) = 1 − P( Z <
0.02 0.03 0.04 0.05 0.06 0.07 0.08
2.9019) ≈ 0.0019, where Z ∼ N (0, 1). See the figure on the right
x/n
for a graphical representation of this result.
The p.d.f. of X n if the misfire rate of
5. State your conclusion the bullets were 0.05. The observed
The P-value ≈ 0.0019 is less than α = 0.05 and we reject H0 in estimate x n = 0.07 is indicated with
a black dot and lies far out in the
favour of H1 . tail of the p.d.f. The probability of
lying further out in the upper tail, the
6. Interpret your conclusion in the relevant context. P-value, is the shaded area.
There is sufficient evidence to reject the manager’s claim that
less than 5% of the bullets misfire.
H0 : µ = 50
H1 : µ 6= 50
The test is two sided because we do not have any extra informa-
tion to specify which way we should test.
X n − 50
q ∼ N (0, 1).
16
25
mathematical methods 2 183
0.5
zα/2 ) = α. From the table zα/2 = z0.025 = 1.96.
0.4
5. State your conclusion
The observed estimate is x n = 51.3 and the observed test
0.3
statistic is
51.3 − 50
fXn
= 1.625.
4/5
0.2
The observed statistic is less than the critical value, so we do not
reject H0 in favour of H1 .
0.1
6. Interpret your conclusion in the relevant context.
0.0
We do not have enough evidence to reject the claim that less ●
xn
z 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
0.0 0.5000 0.5040 0.5080 0.5120 0.5160 0.5199 0.5239 0.5279 0.5319 0.5359
0.1 0.5398 0.5438 0.5478 0.5517 0.5557 0.5596 0.5636 0.5675 0.5714 0.5753
0.2 0.5793 0.5832 0.5871 0.5910 0.5948 0.5987 0.6026 0.6064 0.6103 0.6141
0.3 0.6179 0.6217 0.6255 0.6293 0.6331 0.6368 0.6406 0.6443 0.6480 0.6517
0.4 0.6554 0.6591 0.6628 0.6664 0.6700 0.6736 0.6772 0.6808 0.6844 0.6879
0.5 0.6915 0.6950 0.6985 0.7019 0.7054 0.7088 0.7123 0.7157 0.7190 0.7224
0.6 0.7257 0.7291 0.7324 0.7357 0.7389 0.7422 0.7454 0.7486 0.7517 0.7549
0.7 0.7580 0.7611 0.7642 0.7673 0.7704 0.7734 0.7764 0.7794 0.7823 0.7852
0.8 0.7881 0.7910 0.7939 0.7967 0.7995 0.8023 0.8051 0.8078 0.8106 0.8133
0.9 0.8159 0.8186 0.8212 0.8238 0.8264 0.8289 0.8315 0.8340 0.8365 0.8389
1.0 0.8413 0.8438 0.8461 0.8485 0.8508 0.8531 0.8554 0.8577 0.8599 0.8621
1.1 0.8643 0.8665 0.8686 0.8708 0.8729 0.8749 0.8770 0.8790 0.8810 0.8830
1.2 0.8849 0.8869 0.8888 0.8907 0.8925 0.8944 0.8962 0.8980 0.8997 0.9015
1.3 0.9032 0.9049 0.9066 0.9082 0.9099 0.9115 0.9131 0.9147 0.9162 0.9177
1.4 0.9192 0.9207 0.9222 0.9236 0.9251 0.9265 0.9279 0.9292 0.9306 0.9319
1.5 0.9332 0.9345 0.9357 0.9370 0.9382 0.9394 0.9406 0.9418 0.9429 0.9441
1.6 0.9452 0.9463 0.9474 0.9484 0.9495 0.9505 0.9515 0.9525 0.9535 0.9545
1.7 0.9554 0.9564 0.9573 0.9582 0.9591 0.9599 0.9608 0.9616 0.9625 0.9633
1.8 0.9641 0.9649 0.9656 0.9664 0.9671 0.9678 0.9686 0.9693 0.9699 0.9706
1.9 0.9713 0.9719 0.9726 0.9732 0.9738 0.9744 0.9750 0.9756 0.9761 0.9767
2.0 0.9772 0.9778 0.9783 0.9788 0.9793 0.9798 0.9803 0.9808 0.9812 0.9817
2.1 0.9821 0.9826 0.9830 0.9834 0.9838 0.9842 0.9846 0.9850 0.9854 0.9857
2.2 0.9861 0.9864 0.9868 0.9871 0.9875 0.9878 0.9881 0.9884 0.9887 0.9890
2.3 0.9893 0.9896 0.9898 0.9901 0.9904 0.9906 0.9909 0.9911 0.9913 0.9916
2.4 0.9918 0.9920 0.9922 0.9925 0.9927 0.9929 0.9931 0.9932 0.9934 0.9936
2.5 0.9938 0.9940 0.9941 0.9943 0.9945 0.9946 0.9948 0.9949 0.9951 0.9952
2.6 0.9953 0.9955 0.9956 0.9957 0.9959 0.9960 0.9961 0.9962 0.9963 0.9964
2.7 0.9965 0.9966 0.9967 0.9968 0.9969 0.9970 0.9971 0.9972 0.9973 0.9974
2.8 0.9974 0.9975 0.9976 0.9977 0.9977 0.9978 0.9979 0.9979 0.9980 0.9981
2.9 0.9981 0.9982 0.9982 0.9983 0.9984 0.9984 0.9985 0.9985 0.9986 0.9986
3.0 0.9987 0.9987 0.9987 0.9988 0.9988 0.9989 0.9989 0.9989 0.9990 0.9990
3.1 0.9990 0.9991 0.9991 0.9991 0.9992 0.9992 0.9992 0.9992 0.9993 0.9993
3.2 0.9993 0.9993 0.9994 0.9994 0.9994 0.9994 0.9994 0.9995 0.9995 0.9995
3.3 0.9995 0.9995 0.9995 0.9996 0.9996 0.9996 0.9996 0.9996 0.9996 0.9997
3.4 0.9997 0.9997 0.9997 0.9997 0.9997 0.9997 0.9997 0.9997 0.9997 0.9998
3.5 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998
11
Index