HELM Workbook 13 Integration
HELM Workbook 13 Integration
Integration
13.1 Basic Concepts of Integration 2
Learning outcomes
In this Workbook you will learn about integration and about some of the common techniques
employed to obtain integrals. You will learn that integration is the inverse operation to
differentiation and will also appreciate the distinction between a definite and an indefinite
integral. You will understand how a definite integral is related to the area under a curve.
You will understand how to use the technique of integration by parts to obtain integrals
involving the product of functions. You will also learn how to use partial fractions and
trigonometric identities in integration.
Basic Concepts ◆ ⇣
of Integration ✓13.1 ⌘
Introduction
df
When a function f (x) is known we can di↵erentiate it to obtain its derivative . The reverse process
dx
is to obtain the function f (x) from knowledge of its derivative. This process is called integration.
Applications of integration are numerous and some of these will be explored in subsequent Sections.
First, what is important is to practise basic techniques and learn a variety of methods for integrating
functions.
2 HELM (2008):
Workbook 13: Integration
®
differentiate
x2 2x
integrate
Figure 1
The situation is just a little more complicated because there are lots of functions we can di↵erentiate
to give 2x. Here are some of them: x2 + 4, x2 15, x2 + 0.5
All these functions have the same derivative, 2x, because when we di↵erentiate the constant term we
obtain zero. Consequently, when we reverse the process, we have no idea what the original constant
term might have been. So we include in our answer an unknown constant, c say, called the constant
of integration. We state that the integral of 2x is x2 + c.
d
When we want to di↵erentiate a function, y(x), we use the notation as an instruction to di↵er-
dx
d
entiate, and write y(x) . In a similar way, when we want to integrate a function we use a special
Z dx
notation: y(x) dx.
Z
The symbol for integration, , is known as an integral sign. To integrate 2x we write
!
2x dx = x2 + c
integral
sign
this term is constant of integration
called the
integrand there must always be a
term of the form dx
Note that along with the integral sign there is a term of the form dx, which must always be written,
and which indicates the variable involved, in this case x. We say that 2x is being integrated with
respect to x . The function being integrated is called the integrand. Technically, integrals of this
sort are called indefinite integrals, to distinguish them from definite integrals which are dealt with
subsequently. When you find an indefinite integral your answer should always contain a constant of
integration.
Exercises
1. (a) Write down the derivatives of each of: x3 , x3 + 17, x3 21
Z
(b) Deduce that 3x2 dx = x3 + c.
2. Explain why, when finding an indefinite integral, a constant of integration is always needed.
HELM (2008): 3
Section 13.1: Basic Concepts of Integration
Answers
1. (a) 3x2 , 3x2 , 3x2 (b) Whatever the constant, it is zero when di↵erentiated.
2. Any constant will disappear (i.e. become zero) when di↵erentiated so one must be reintroduced
to reverse the
process.
2. A table of integrals
We could use a table of derivatives to find integrals, but the more common ones are usually found
in a ‘Table of Integrals’ such as that shown below. You could check the entries in this table using
your knowledge of di↵erentiation. Try this for yourself.
constant, k
kx + c
1 2
x 2
x +c
1 3
x2 3
x +c
xn+1
xn + c, n 6= 1
n+1
1
x 1 (or ) ln |x| + c
x
cos x sin x + c
sin x cos x + c
1
cos kx sin kx + c
k
1
sin kx cos kx + c
k
1
tan kx ln | sec kx|+c
k
ex ex + c
e x e x+c
1 kx
ekx e +c
k
When dealing with the trigonometric functions the variable x must always be measured in radians
and not degrees. Note that the fourth entry in the Table, for xn , is valid for any value of n, positive
or negative, whole number or fractional, except n = 1. When n = 1 use the fifth entry in the
Table.
4 HELM (2008):
Workbook 13: Integration
®
Example 1 Z
7
Use Table 1 to find the indefinite integral of x : that is, find x7 dx
Solution
Z
xn+1
From Table 1 note that xn dx = + c. In words, this states that to integrate a power
n+1
of x, increase the power by 1, and then divide the result by the new power. With n = 7 we find
Z
1
x7 dx = x8 + c
8
Example 2 Z
Find the indefinite integral of cos 5x: that is, find cos 5x dx
Solution
Z
sin kx
From Table 1 note that cos kx dx = +c
k
Z
1
With k = 5 we find cos 5x dx = sin 5x + c
5
In Table 1 the independent variable is always given as x. However, with a little imagination you will
be able to use it when other independent variables are involved.
Example
Z 3
Find cos 5t dt
Solution
We integrated cos 5x in the previous example. Now the independent variable is t, so simply use
Table 1 and replace every x with a t. With k = 5 we find
Z
1
cos 5t dt = sin 5t + c
5
It follows immediately that, for example,
Z Z
1 1
cos 5! d! = sin 5! + c, cos 5u du = sin 5u + c and so on.
5 5
HELM (2008): 5
Section 13.1: Basic Concepts of Integration
Example 4 Z
1 1
Find the indefinite integral of : that is, find dx
x x
Solution
This integral deserves special mention. You may be tempted
Z to try to write the integrand as x 1
n+1
x
and use the fourth row of Table 1. However, the formula xn dx = + c is not valid when
n+1
n = 1 as Table 1 makes clear. ZThis is because we can never divide by zero. Look to the fifth
entry of Table 1 and you will see x 1 dx = ln |x| + c.
Example
Z 5 Z
Find 12 dx and 12 dt
Solution
In this Example we are integrating a constant, 12. Using Table 1 we find
Z Z
12 dx = 12x + c Similarly 12 dt = 12t + c.
Task Z
Find t4 dt
Your solution
Answer
Z
1
t4 dt = t5 + c.
5
6 HELM (2008):
Workbook 13: Integration
®
Task Z
1 5
Find dx using the laws of indices to write the integrand as x and then use
x5
Table 1:
Your solution
Answer
1 4 1
x +c= + c.
4 4x4
Task Z
2x
Find e dx using the entry in Table 1 for integrating ekx :
Your solution
Answer
R 2x 1 2x
With k = 2, we have e dx = e + c.
2
Exercises
1. Integrate each of the following functions with respect to x: p
(a) x9 , (b) x1/2 , (c) x 3 , (d) 1/x4 , (e) 4, (f) x, (g) e4x
Z Z Z Z
2
2. Find (a) t dt, (b) 6 dt, (c) sin 3t dt, (d) e7t dt.
Text
Answers
1 10 2 3/2 1 2 1 3
1 (a) x + c, (b) x + c, (c) x + c, (d) x + c, (e) 4x + c,
10 3 2 3
1
(f) same as (b), (g) e4x + c
4
1 3 1 1 7t
2. (a) t + c, (b) 6t + c, (c) cos 3t + c, (d) e +c
3 3 7
HELM (2008): 7
Section 13.1: Basic Concepts of Integration
3. Some rules of integration
To enable us to find integrals of a wider range of functions than those normally given in a table of
integrals we can make use of the following rules.
Key Point 1
Z Z
k f (x) dx = k f (x) dx
Example 6 Z
2
Find the indefinite integral of 11x : that is, find 11x2 dx
Solution
Z Z ✓ 3 ◆
2 2 x 11x3
11x dx = 11 x dx = 11 +c = +K where K is a constant.
3 3
Example 7 Z
Find the indefinite integral of 5 cos x; that is, find 5 cos x dx
Solution
Z Z
5 cos x dx = 5 cos x dx = 5 (sin x + c) = 5 sin x + K where K is a constant.
8 HELM (2008):
Workbook 13: Integration
®
Key Point 2
Z Z Z
[ f (x) + g(x) ] dx = f (x) dx + g(x) dx
Z Z Z
[ f (x) g(x) ] dx = f (x) dx g(x) dx
Example
Z 8
Find (x3 + sin x) dx
Solution
Z Z Z
1
(x + sin x) dx = x dx + sin x dx = x4
3 3
cos x + c
4
Note that only a single constant of integration is needed.
Task Z p
Find (3t4 + t) dt
Your solution
Answer
3 5 2 3/2
t + t +c
5 3
HELM (2008): 9
Section 13.1: Basic Concepts of Integration
Task
The hyperbolic sine and cosine functions, sinh x and cosh x, are defined as follows:
ex e x
ex + e x
sinh x = cosh x =
2 2
Note that they are combinations of the exponential functions ex and e x .
Find the indefinite integrals of sinh x and cosh x.
Your
Z solution Z ✓ ◆
ex e x
sinh x dx = dx =
2
Z Z ✓ ◆
ex + e x
cosh x dx = dx =
2
Answer
Z Z Z
1 x 1 1 1 1 x
sinh x dx = e dx e x dx = ex + e x
+c= e +e x
+ c = cosh x + c.
2 2 2 2 2
Z
Similarly cosh x dx = sinh x + c.
Further rules for finding more complicated integrals are dealt with in subsequent Sections.
10 HELM (2008):
Workbook 13: Integration
®
Engineering Example 1
Electrostatic charge
Introduction
Electrostatic charge is important both where it is wanted, as in the electrostatic precipitator plate
systems used for cleaning gases, and where it is unwanted, such as when charge builds up on moving
belts. This Example is concerned with a charged object with a particular idealised shape - a sphere.
However, similar analytical calculations can be carried out for certain other shapes and numerical
methods can be used for more complicated shapes.
The electric field at all points inside and outside a charged sphere is given by
Qr
E(r) = if r < a (1a)
4⇡"0 a3
Q
E(r) = if r a (1b)
4⇡"0 r2
where "0 is the permittivity of free space, Q is the total charge, a is the radius of the sphere, and r
is the radial distance between the centre of the sphere and a point of observation (see Figure 2).
Charged sphere
a r
O
Spherical surface S
The electric field associated with electrostatic charge has a scalar potential. The electric field defined
by (1a) and (1b) shows only a radial dependence of position. Therefore, the electric scalar potential
V (r) is related to the field E(r) by
dV
E(r) = . (2)
dr
Problem in words
A sphere is charged with a uniform density of charge and no other charge is present outside the
sphere in space. Determine the variation of electric potential with distance from the centre of the
sphere.
Mathematical analysis
Equation (2) yields V (r) as the negative of the indefinite integral of E(r).
HELM (2008): 11
Section 13.1: Basic Concepts of Integration
Z Z
dV = E(r) dr. (3)
Q
V (r) = + c2 if r a (5b)
4⇡"0 r
The integration constant c2 can be determined by assuming that the electric potential is zero at an
infinite distance from the sphere:
Q
lim [V (r)] = 0 ) lim + c2 = 0 ) c2 = 0.
r!1 r!1 4⇡"0 r
The constant c1 can be determined by assuming that the potential is continuous at r = a.
From equation (5a)
Qa2
V (a) = + c1
8⇡"0 a3
From equation (5b)
Q
V (a) =
4⇡"0 a
Hence
Q 2Q 3Q
c1 = + = .
4⇡"0 a 8⇡"0 a 8⇡"0 a
Substituting for c1 in (5), the electric potential is obtained for all space is:
✓ 2 ◆
Q 3a r2
V (r) = if r < a.
4⇡"0 2a3
Q
V (r) = if r a
4⇡"0 r
Interpretation
The potential of the electrostatic field outside a charged sphere varies inversely with distance from
the centre of the sphere. Inside the sphere, the electrostatic potential varies with the square of the
distance from the centre.
An Engineering Exercise in 29.3 derives the corresponding expressions for the variation of the
electrostatic field and an Engineering Exercise in 27.4 calculates the potential energy due to
the charged sphere.
12 HELM (2008):
Workbook 13: Integration
®
Exercises
Z
1. Find (2x ex ) dx
Z
2. Find 3e2x dx
Z
1
3. Find (x + cos 2x) dx
Z 3
4. Find 7x 2 dx
Z
5. Find (x + 3)2 dx, (be careful!)
Answers
1. x2 ex + c
3
2. e2x + c
2
1 1
3. x2 + sin 2x + c
6 6
7
4. +c
x
1
5. x3 + 3x2 + 9x + c
3
HELM (2008): 13
Section 13.1: Basic Concepts of Integration
Introduction
When you were first introduced to integration as the reverse of differentiation, the integrals you dealt
with were indefinite integrals. The result of finding an indefinite integral is usually a function plus a
constant of integration. In this Section we introduce definite integrals, so called because the result
will be a definite answer, usually a number, with no constant of integration. Definite integrals have
many applications, for example in finding areas bounded by curves, and finding volumes of solids.
• understand integration as the reverse of
Prerequisites differentiation
Before starting this Section you should . . . • be able to use a table of integrals
• find simple definite integrals
Learning Outcomes
• handle some integrals involving an infinite
On completion you should be able to . . . limit of integration
14 HELM (2008):
Workbook 13: Integration
®
1. Definite integrals
Z
We saw in the previous Section that f (x) dx = F (x) + c where F (x) is that function which, when
dF
differentiated, gives f (x). That is, = f (x). For example,
dx
Z
cos(3x)
sin(3x) dx = − +c
3
Here, f (x) = sin(3x) and F (x) = − 13 cos(3x) We now consider a definite integral which is simply
an indefinite integral but with numbers written to the upper and lower right of the integral sign. The
quantity
Z b
f (x) dx
a
is called the definite integral of f (x) from a to b. The numbers a and b are known as the lower
limit and upper limit respectively of the integral. We define
Z b
f (x) dx = F (b) − F (a)
a
so that a definite integral is usually a number. The meaning of a definite integral will be developed
in later Sections. For the present we concentrate on the process of evaluating definite integrals.
Example 9 Z 4
2
Find the definite integral of x from 1 to 4; that is, find x2 dx
1
Solution
Z
x2 dx = 13 x3 + c
3
Here f (x) = x2 and F (x) = x3 . Thus, according to our definition
Z 4
43 13
x2 dx = F (4) − F (1) = − = 21
1 3 3
HELM (2008): 15
Section 13.2: Definite Integrals
Writing F (b) − F (a) each time we calculate a definite integral becomes laborious so we replace this
b
difference by the shorthand notation F (x) . Thus
a
b
F (x) ≡ F (b) − F (a)
a
Example 10 Z π/2
π
Find the definite integral of cos x from 0 to ; that is, find cos x dx.
2 0
Solution
R
Since cos x dx = sin x + c then
Z π/2 π/2
cos x dx = sin x
0
π 0
= sin − sin 0 = 1 − 0 = 1
2
Always remember, that if you use a calculator to evaluate any trigonometric functions, you must
work in radian mode.
Task Z 2
2
Find the definite integral of x + 1 from 1 to 2; that is; find (x2 + 1) dx
1
Your solution
Answer
2
1 3
x +x .
3 1
16 HELM (2008):
Workbook 13: Integration
®
Now insert the limits of integration, the upper limit first, and hence evaluat the integral:
Your solution
Answer
8 1 10
+2 − +1 = or 3.333 (3 d.p.).
3 3 3
Task Z 1
Find (x2 + 1) dx.
2
This Task is very similar to the previous Task. Note the limits have been interchanged:
Your solution
Answer
1
1 3 1 8 10
x +x = +1 − +2 =− .
3 2 3 3 3
Note from these two Tasks that interchanging the limits of integration, changes the sign of the
answer.
Key Point 3
If you interchange the limits, you must change the sign:
Z b Z a
f (x) dx = − f (x) dx
a b
HELM (2008): 17
Section 13.2: Definite Integrals
Task
When a spring is fixed at one end and stretched at the free end it exerts a restoring
force that is proportional to the displacement of the free end. The constant
of proportionality k N m−1 is known as the stiffness of the spring. Calculate
the work done in stretching a spring with stiffness k from displacement x1 m to
displacement x2 m (x2 > x1 ) given that the work done (W ) is the product of
force and displacement.
Your solution
Answer
The restoring force varies during the displacement. So the work done during the extension cannot
be determined from a single simple product.
Consider a small element ∆x of the extension beyond an arbitrary displacement x. The element is
sufficiently small that the force during the displacement can be regarded as constant and equal to
the force at displacement x is kx. So the work done ∆W in extending the spring from displacement
x to displacement x + ∆x is approximately kx∆x.
Using the idea of integration as a limit of a sum, in this case as ∆x tends to zero,
Z x2 x
1 2 2 1
W = kx dx = kx = k(x22 − x21 )
x1 2 x1 2
Exercises
Z 1 Z 3 Z 2 Z 1
2 1 x
1. Evaluate (a) x dx, (b) dx (c) e dx (d) (1 + t2 ) dt
0 x2 −1
Z π/3 Z 2π 1
Z 3
2. Find (a) cos 2x dx (b) sin x dx (c) e2t dt
0 0 1
Answers
1 1
1. (a) (b) (c) e2 − e1 = 4.671 (d) 2.667
3 6
√
2 (a) 3/4 = 0.4330 (b) 2 (c) 198.019
18 HELM (2008):
Workbook 13: Integration
®
Engineering Example 2
Introduction
For materials such as mild-steel, the relationship between applied shear stress and shear strain (de-
formation) can be described as follows.
• For small values of the shear strain, the shear stress (τ ) and shear strain (ω) are proportional
to one another, i.e.
1
ω= ×τ (1)
G
(where G is the shear modulus). This is known as elastic behaviour.
• There is a maximum shear stress that the material is capable of supporting. If the shear
strain is increased further, the shear stress remains roughly constant. This is known as plastic
behaviour.
Figure 3 summarises the relationship between shear stress and shear strain; the point (ωY , τY ) is
known as the yield point.
τ
6
τY p p p p p p p
p
shear pp
stress pp
p
p
p
p
p
-ω
ωY
shear strain
Figure 3
Now suppose that one end of a bar of circular cross section is twisted through an angle θ, then the
shear strain on the surface is given by
Rθ
ωS = (2)
L
(where R and L are the radius and length of the bar respectively), while the shear strain, at a distance
r from the central core, is given by
rθ
ω= (3)
L
The torque transmitted by a bar is given by the integral
Z R
T = 2π r2 τ (r) dr (4)
0
HELM (2008): 19
Section 13.2: Definite Integrals
As the shear strain is a function of distance from the central axis of the bar, it may be that the shear
strain on the surface is greater than the critical shear strain ωY . In this scenario the shear stress is
given by
τ
Y
ω ω ≤ ωY
ωY
τ= (5)
τY ω > ωY
i.e. the regions near the central axis exhibit elasticity, but in those regions near the surface the elastic
limit has been exceeded and the metal exhibits plasticity (see Figure 4).
R
elastic
zone re
plastic zone
τ
τY
0 re R
Figure 4
Problem in words
Find an expression for the torque transmitted by a bar as a function of the angle θ through which
one end is turned.
Mathematical statement of problem
Using Equations (3) to (5), find a formula for T in terms of the variable θ.
Mathematical analysis
Substituting (3) into (5)
τY r θ
rθ
≤ ωY
ω L
L
Y
τ =
τY rθ
> ωY
L
L ωY
τY r θ
ω L r≤ = re
θ
Y
=
τ L ωY
r> = re
Y
θ
20 HELM (2008):
Workbook 13: Integration
®
For small values of θ, re ≥ R so that the whole of the bar will be in the elastic region, i.e.
τY r θ
τ=
ωY L
Now (4) becomes
Z R Z R R
2 τY r θ τY θ 3 τY θ r 4 π τY θ 4
T = 2π r dr = 2π r dr = 2π = R (6)
0 ωY L ωY L 0 ωY L 4 0 2 ωY L
i.e. the torque is directly proportional to the twist, θ.
For larger θ, re < R, so that (4) becomes
Z re Z R
2 τY r θ
T = 2π r dr + 2π r2 τY dr
0 ωY L re
Z re Z R
τY θ 3
= 2π r dr + 2π τY r2 dr
ωY L 0 re
4 re 3 R
τ θ r r
= 2π Y + 2π τY
ωY L 4 0 3 re
π τY θ 4 2π
τY R3 − re3
= re +
2 ωY L 3
But re = L ωY /θ, so
4 4
π τY θ L ωY 2π 3 2π L3 ω 3
T = + τ R − τ Y
2 ωY L θ4 3 Y 3 Y θ3
3 3
L ω
2π 3 1 2
= τY R + π τY − τY Y
3 2 3 θ3
2π 3 π L3 ω 3
= τ R − τY Y
(7)
3 Y 6 θ3
Equation (6) will apply when re ≥ R, i.e. (L ωY /θ) ≥ R or θ ≤ (L ωY /R), so that combining (6)
and (7) gives overall
π τY θ 4 L ωY
R θ≤
2 ω L
R
Y
T = (8)
3 3
2π π L ω L ω
τ R 3 − τY
Y Y
θ>
3 Y 6 θ3 R
Interpretation and further comment
At the critical value of θ, i.e. when the outer edge begins to exhibit plasticity, both formulae in (8)
give
π
Tcrit = τY R3
2
Furthermore, the first derivatives are both
dT π τY R 4
=
dθ 2 ωY L
HELM (2008): 21
Section 13.2: Definite Integrals
i.e. the curves join smoothly.
The second derivatives, though, are not equal (zero in one case). In the theoretical limit as θ → ∞
2π
T = τ R3
3 Y
so this is the total torsional torque which can be carried by the bar. (The critical torque above is
three-quarters of this value.) However, clearly θ → ∞ is merely a theoretical limit since the bar
would, in fact, shear at a finite value of θ.
Example 11 Z ∞
Find the definite integral of e −x
from 0 to ∞; that is, find e−x dx.
0
Solution
Z ∞ ∞
−x −x
The integral is found in the normal way: e dx = −e
0 0
There is no difficulty in evaluating the square bracket at the lower limit. We obtain simply −e−0 =
−1. At the upper limit we must examine the behaviour of −e−x as x gets infinitely large. This is
where it is important that you are familiar with the properties of the exponential function. If you
refer to the graph (Figure 5) you will see that as x tends to infinity e−x tends to zero.
Consequently the contribution to the integral from the upper limit is zero. So
! ∞ " #∞
−x −x
e dx = −e
0 0
= (−e ) − (−e−0 )
−∞ e−x
= (0) − (−e−0 )
= 1 x
Figure 5
Z ∞
Thus the value of e−x dx is 1.
0
22 HELM (2008):
Workbook 13: Integration
®
Z ∞ Z b
−x
e dx = lim e−x dx
0 b→∞ 0
Z b b
−x −x
= −e−b − −e−0 = −e−b + 1
Now, e dx = −e
0 0
−b
Then as b tends to infinity −e tends to zero, and the resulting integral has the value 1, as before.
Many integrals having infinite limits cannot be evaluated in a simple way like this, and many cannot
be evaluated at all. Fortunately, most of the integrals you will meet will exhibit the sort of behaviour
seen in the last example.
Exercise
Z ∞ Z ∞ Z ∞ Z ∞
−x −2x −3x 4
Evaluate (a) e dx (b) e dx (c) e dx (d) dt
1 0 2 1 t2
Answer
(a) e−1 ∼ 0.368 (b) 1
2
(c) 13 e−6 = 0.0008 (4 d.p.) (d) 4
HELM (2008): 23
Section 13.2: Definite Integrals
The Area Bounded
by a Curve 13.3
Introduction
One of the important applications of integration is to find the area bounded by a curve. Often such
an area can have a physical significance like the work done by a motor, or the distance travelled by
a vehicle. In this Section we explain how such an area is calculated.
' $
• understand integration as the reverse of
differentiation
24 HELM (2008):
Workbook 13: Integration
®
y area is A(x)
y = f (x)
a x x
Figure 6
A(x) is clearly a function of x since as the upper limit changes so does the area. How does the area
change if we change the upper limit by a very small amount δx? See Figure 7 below.
y
y = f (x)
a x x + δx x
Figure 7
To a good approximation the change in the area is:
A(x + δx) − A(x) ≈ f (x)δx
[This is because the shaded area is approximately a rectangle with base δx and height f (x).] This
approximation gets better and better as δx gets smaller and smaller. Rearranging gives:
A(x + δx) − A(x)
f (x) ≈
δx
Clearly, in the limit as δx → 0 we have
A(x + δx) − A(x)
f (x) = lim
δx→0 δx
But this limit on the right-hand side is the derivative of A(x) with respect to x so
dA(x)
f (x) =
dx
Thus A(x) is an indefinite integral of f (x) and we can therefore write:
Z
A(x) = f (x)dx
Now the area under the curve from a to b is clearly A(b) − A(a). But remembering our shorthand
notation for this difference, introduced in the last Section we have, finally
b Z b
A(b) − A(a) ≡ A(x) = f (x)dx
a a
We conclude that the area under the curve y = f (x) from a to b is given by the definite integral of
f (x) from a to b.
HELM (2008): 25
Section 13.3: The Area Bounded by a Curve
2. The area bounded by a curve lying above the x-axis
Consider the graph of the function y = f (x) shown in Figure 8. Suppose we are interested in
calculating the area underneath the graph and above the x-axis, between the points where x = a
and x = b. When such an area R b lies entirely above the x-axis, as is clearly the case here, this area is
given by the definite integral a f (x) dx.
y area required
y = f (x)
a b x
Figure 8
Key Point 4
Z b
The area under the curve y = f (x), between x = a and x = b is given by f (x) dx
a
when the curve lies entirely above the x-axis between a and b.
Example 12
Calculate the area bounded y = x−1 and the x-axis, between x = 1 and x = 4.
Solution
Below is a graph of y = x−1 . The area required is shaded; it lies entirely above the x-axis.
y 1
y=
x area required
1
x
O 1 2 3 4 5
Figure 9
Z 4 4
1
area = dx = ln |x| = ln 4 − ln 1 = ln 4 = 1.386 (3 d.p.)
1 x 1
26 HELM (2008):
Workbook 13: Integration
®
Task
Find the area bounded by the curve y = sin x and the x-axis between x = 0 and
x = π. (The required area is shown in the figure. Note that it lies entirely above
the x-axis.)
y y =sinx
area required
O π x
Your solution
Answer
Z π π
sin x dx = − cos x = 2.
0 0
Task
Find the area under f (x) = e2x from x = 1 to x = 3 given that the exponential
function e2x is always positive.
Your solution
Answer
Z 3 3
2x 1 2x
area = e dx = e = 198 to 3 significant figures.
1 2 1
HELM (2008): 27
Section 13.3: The Area Bounded by a Curve
Example 13
The figure shows the graphs of y = sin x and y = cos x for 0 ≤ x ≤ 21 π. The two
graphs intersect at the point where x = 14 π. Find the shaded area.
y
y =cos x
π
2 x
Figure 10
Solution
To find the shaded area we could calculate the area under the graph of y = sin x for x between
0 and 14 π, and subtract this from the area under the graph of y = cos x between the same limits.
Alternatively the two processes can be combined into one and we can write
Z π/4
shaded area = (cos x − sin x)dx
0
π/4
= sin x + cos x
0
1 1
= sin + cos 4
π− (sin 0 + cos 0)
4
π
1 1 2 √
= ( √ + √ ) − (0 + 1) = √ − 1 = 2 − 1
2 2 2
Exercises
In each question you should check that the required area lies entirely above the horizontal axis.
1. Find the area under the curve y = 7x2 and above the x-axis between x = 2 and x = 5.
2. Find the area bounded by the curve y = x3 and the x-axis between x = 0 and x = 2.
3. Find the area bounded by the curve y = 3t2 and the t-axis between t = −3 and t = 3.
4. Find the area under y = x−2 between x = 1 and x = 10.
Answer
1. 273, 2. 4, 3. 54, 4. 0.9.
28 HELM (2008):
Workbook 13: Integration
®
y
y =−x2 +1
-2 -1 1 2 3 x
area required
Figure 11
The shaded area is bounded
R2 by the x-axis and the curve, but lies entirely below the x-axis. Let us
2
evaluate the integral 1 (−x + 1)dx.
Z 2 3 2
2 x
(−x + 1)dx = − + x
1 3
3 1 3
2 1
= − +2 − − +1
3 3
7 4
= − +1=−
3 3
The evaluation of the area yields a negative quantity. There is, of course, no such thing as a negative
area. The area is actually 43 , and the negative sign is an indication that the area lies below the x-axis.
(However, in applications of integration such as work/energy or distance travelled in a given direction
negative values can be meaningful.)
If an area contains parts both above and below the horizontal axis, care must be taken when calcu-
lating this area. It is necessary to determine which parts of the graph lie above the horizontal axis
and which lie below. Separate integrals need to be calculated for each ‘piece’ of the graph. This idea
is illustrated in the next Example.
HELM (2008): 29
Section 13.3: The Area Bounded by a Curve
Example 14
Find the total area enclosed by the curve y = x3 −5x2 +4x and the x-axis between
x = 0 and x = 3.
Solution
We need to determine which parts of the graph lie above and which lie below the x-axis. To do this
it is helpful to consider where the graph cuts the x-axis. So we consider the function x3 − 5x2 + 4x
and look for its zeros
x3 − 5x2 + 4x = x(x2 − 5x + 4) = x(x − 1)(x − 4)
So the graph cuts the x-axis when x = 0, x = 1 and x = 4. Also, when x is large and positive,
y is large and positive since the term involving x3 dominates. When x is large and negative, y is
large and negative for the same reason. With this information we can sketch a graph showing the
required area:
1 2 3 4 x
area required
Figure 12
From the graph we see that the required area lies partly above the x-axis (when 0 ≤ x ≤ 1) and
partly below (when 1 ≤ x ≤ 3). So we evaluate the integral in two parts: Firstly:
Z 1 4 1
5x3 4x2
3 2 x 1 5 7
(x − 5x + 4x)dx = − + = − + 2 − (0) =
0 4 3 2 0 4 3 12
This is the part of the required area which lies above the x-axis. Secondly:
3 3
x4 5x3 4x2
Z
3 2
(x − 5x + 4x)dx = − +
1 4 3 2 1
81 135 1 5 22
= − + 18 − − +2 =−
4 3 4 3 3
This represents the part of the required area which lies below the x-axis. The actual area is 22
3
.
Combining the results of the two separate calculations we can find the total area bounded by the
curve:
7 22 95
area = + =
12 3 12
30 HELM (2008):
Workbook 13: Integration
®
Task
(a) Sketch the graph of y = sin 2x for 0 ≤ x ≤ π.
(b) Find the total area bounded by the curve and the x-axis between x = 13 π
and x = 43 π.
(a) Sketch the graph and indicate the required area noting where the graph crosses the x-axis:
Your solution
Answer
y
π 3π
2 4 x
π π
3
y =sin 2x
(b) Perform the integration in two parts to obtain the required area:
Your solution
Answer
Z π/2 Z 3π/4
1 1
sin 2x dx = and sin 2xdx = − .
π/3 4 π/2 2
1 1 3
The required area is + = .
4 2 4
HELM (2008): 31
Section 13.3: The Area Bounded by a Curve
Exercises
1. Find the total area enclosed between the x-axis and the curve y = x3 between x = −1 and
x = 1.
6. Find the area enclosed between y = x(x − 1)(x − 2) and the x axis.
Answers
16 1
1. 0.5 2. 0.4207 3. (a) 3
, (b) 9, (c) 4 4. 0.5 5. e2 − 1 6. 2
32 HELM (2008):
Workbook 13: Integration
®
Introduction
Integration by Parts is a technique for integrating products of functions. In this Section you will learn
to recognise when it is appropriate to use the technique and have the opportunity to practise using
it for finding both definite and indefinite integrals.
' $
• understand what is meant by definite and
indefinite integrals
Prerequisites • be able to use a table of integrals
Before starting this Section you should . . .
• be able to differentiate and integrate a range
of common functions
& %
' $
• decide when it is appropriate to use the
method known as integration by parts
HELM (2008): 33
Section 13.4: Integration by Parts
1. Indefinite integration
The technique known as integration by parts is used to integrate a product of two functions, such
as in these two examples:
Z Z 1
(i) 2x
e sin 3x dx (ii) x3 e−2x dx
0
Note that in the first example, the integrand is the product of the functions e2x and sin 3x, and in
the second example the integrand is the product of the functions x3 and e−2x . Note also that we
can change the order of the terms in the product if we wish and write
Z Z 1
(i) 2x
(sin 3x) e dx (ii) e−2x x3 dx
0
What you must never do is integrate each term in the product separately and then multiply - the
integral of a product is not the product of the separate integrals. However, it is often possible to
find integrals involving products using the method of integration by parts - you can think of this as
a product rule for integrals.
The integration by parts formula states:
Key Point 5
Integration by Parts for Indefinite Integrals
For indefinite integrals, given functions f (x) and g(x):
Z Z Z Z
df
f · g dx = f · g dx − · gdx dx
dx
Study the formula carefully and note the following observations. Firstly, to apply the formula we must
df
be able to differentiate the function f to find , and we must be able to integrate the function, g.
dx
Secondly the formula replaces one integral, the one on the left, with a different integral, that on the
far right. The intention is that the latter, whilst it may look more complicated in the formula above,
is simpler to evaluate. Consider the following Example:
34 HELM (2008):
Workbook 13: Integration
®
Example 15 R
Find the integral of the product of x with sin x; that is, find x sin x dx.
Solution
Compare the required integral with the formula for integration by parts: we choose
f =x and g = sin x
It follows that
Z Z
df
=1 and g dx = sin x dx = − cos x
dx
(When integrating g there is no need to worry about a constant of integration. When you become
confident with the method, you may like to think about why this is the case.)
Applying the formula we obtain
Z Z Z Z
df
x sin x dx = f · g dx − · gdx dx
dx
Z
= x(− cos x) − 1(− cos x) dx
Z
= −x cos x + cos x dx = −x cos x + sin x + c
Task Z
Find (5x + 1) cos 2x dx.
Z
df
Let f = 5x + 1 and g = cos 2x. Now calculate and g dx:
dx
Your solution
Answer Z
df 1
= 5 and cos 2x dx = sin 2x.
dx 2
Substitute these results into the formula for integration by parts and complete the Task:
Your solution
HELM (2008): 35
Section 13.4: Integration by Parts
Answer Z
1 1 1 5
(5x + 1)( sin 2x) − 5( sin 2x)dx = (5x + 1) sin 2x + cos 2x + c
2 2 2 4
Sometimes it is necessary to apply the formula more than once, as the next Example shows.
Example
Z 16
Find 2x2 e−x dx
Solution
Z
df
2
We let f = 2x and g = e . Then −x
= 4x and gdx = −e−x
dx
Using the formula for integration by parts we find
Z Z Z
2x e dx = 2x (−e ) − 4x(−e )dx = −2x e + 4xe−x dx
2 −x 2 −x −x 2 −x
Z Z
4xe dx = 4x(−e ) − 4(−e−x )dx
−x −x
Z
= −4xe + 4e−x dx = −4xe−x − 4e−x
−x
Altogether we have
Z
2x2 e−x dx = −2x2 e−x − 4xe−x − 4e−x + c = −2e−x (x2 + 2x + 2) + c
Exercises
In some questions
Z below it will be necessary
Z to apply integration
Z by parts more than once.
3t
1. Find (a) x sin(2x)dx, (b) te dt, (c) x cos x dx.
Z
2. Find (x + 3) sin x dx.
Z
3. By writing ln x as 1 × ln x find ln x dx.
Z Z Z
−1
4. Find (a) tan x dx, (b) −7x cos 3x dx, (c) 5x2 e3x dx,
Z Z
5. Find (a) x cos kx dx, where k is a constant (b) z 2 cos kz dz, where k is a constant.
Z Z
−st
6. Find (a) te dt where s is a constant, (b) Find t2 e−st dt where s is a constant.
36 HELM (2008):
Workbook 13: Integration
®
Answers
1 1 1 1
1. (a) sin 2x − x cos 2x + c, (b) e3t ( t − ) + c, (c) cos x + x sin x + c
4 2 3 9
2. −(x + 3) cos x + sin x + c.
3. x ln x − x + c.
1 7 7 5
4. (a) x tan−1 x − ln(x2 + 1) + c, (b) − cos 3x − x sin 3x + c, (c) e3x (9x2 − 6x + 2) + c,
2 9 3 27
cos kx x sin kx 2z cos kz z 2 sin kz 2 sin kz
5. (a) + + c, (b) + − + c.
k2 k k2 k k3
−e−st (st + 1) −e−st (s2 t2 + 2st + 2)
6. (a) + c, (b) + c.
s2 s3
2. Definite integration
When dealing with definite integrals the relevant formula is as follows:
Key Point 6
Integration by Parts for Definite Integrals
For definite integrals, given functions f (x) and g(x):
Z b Z b Z b Z
df
f · g dx = f · g dx − · gdx dx
a a a dx
Z b b Z b
dv du
Alternatively, given functions u and v: u dx = uv − v dx
a dx a a dx
Example 17
Z 2
Find xex dx.
0
Solution
Z
df
We let f = x and g = e . Then = 1 and g dx = ex . Using integration by parts we obtain
x
dx
Z 2 2 Z 2 2
x x
xe dx = xe − 1.e dx = 2e − ex = 2e2 −[e2 −1] = e2 +1
x 2
(or 8.389 to 3 d.p.)
0 0 0 0
Sometimes it is necessary to apply the formula more than once as the next Example shows.
HELM (2008): 37
Section 13.4: Integration by Parts
Example 18 Z 2
2 x
Find the definite integral of x e from 0 to 2; that is, find x2 ex dx.
0
Solution
Z
df
2
We let f = x and g = e . Then x
= 2x and g dx = ex . Using integration by parts:
dx
Z 2 2 Z 2 Z 2
2 x 2 x x 2
x e dx = x e − 2xe dx = 4e − 2 xex dx
0 0 0 0
The remainingZ integral must be integrated by parts also but we have just done this in the example
2
above. So x2 ex dx = 4e2 − 2[e2 + 1] = 2e2 − 2 = 12.778 (3 d.p.)
0
Task Z π/4
Find (4 − 3x) sin x dx.
0
Your solution
Answer
Take f = 4 − 3x and g = sin x.
Now complete the integral:
Your solution
Z π/4
(4 − 3x) sin x dx =
0
Answer
Z π/4 π/4 Z π/4
(4 − 3x) sin x dx = (4 − 3x)(− cos x) −3 cos x dx
0 0 0
π/4 π/4
= (4 − 3x)(− cos x) −3 sin x
0 0
= 0.716 to 3 d.p.
38 HELM (2008):
Workbook 13: Integration
®
Exercises
Z 1 Z π/2 Z 1
1. Evaluate the following: (a) x cos 2x dx, (b) x sin 2x dx, (c) te2t dt
Z 2 0 0 −1
2. Find (x + 2) sin x dx
1
Z 1
3. Find (x2 − 3x + 1)ex dx
0
Answers
1. (a) 0.1006, (b) π/4 = 0.7854, (c) 1.9488.
2. 3.3533.
3. −0.5634.
HELM (2008): 39
Section 13.4: Integration by Parts
Integration by
Substitution and Using
Introduction
The first technique described here involves making a substitution to simplify an integral. We let
a new variable equal a complicated part of the function we are trying to integrate. Choosing the
correct substitution often requires experience. This skill develops with practice.
Often the technique of partial fractions can be used to write an algebraic fraction as the sum of simpler
fractions. On occasions this means that we can then integrate a complicated algebraic fraction. We
shall explore this approach in the second half of the section.
' $
• be able to find a number of simple definite
and indefinite integrals
Before starting this Section you should . . . • be familiar with the technique of expressing
an algebraic fraction as the sum of its partial
fractions
&
' %
$
• make simple substitutions in order to find
definite and indefinite integrals
40 HELM (2008):
Workbook 13: Integration
®
1. Making a substitution
The technique described here involves making a substitution in order to simplify an integral. We let
a new variable, u say, equal a more complicated part of the function we are trying to integrate. The
choice of which substitution to make often relies upon experience: don’t worry if at first you cannot
see an appropriate substitution. This skill develops with practice. However, it is not simply a matter
of changing the variable - care must be taken with the differential form dx as we shall see. The
technique is illustrated in the following Example.
Example 19
Z
Find (3x + 5)6 dx.
Solution
First look at the function we are trying to integrate: (3x + 5)6 . It looks quite complicated to
integrate. Suppose we introduce a new variable, u, such that u = 3x + 5. Doing this means that
the function we must integrate becomes u6 . Would you not agree that this looks a much simpler
function to integrate than (3x + 5)6 ? There is a slight complication however. The new function of
u must be integrated with respect to u and not with respect to x. This means that we must take
care of the term dx correctly.
du dx 1
Long Method u = 3x + 5 so = 3, or =
dx du 3
Z Z
6
Let I= (3x + 5) dx = u6 dx (substituting for 3x + 5)
Z
dx
= u6 du (to change from x to u)
du
Z
1 dx
= u6 . du (substituting for )
3 du
u7
Z
1
= u6 dx = + constant
3 21
du 1
Short Method u = 3x + 5 so = 3, so dx = du
dx 3
u7
Z Z Z Z
1 1
Let I = (3x + 5)6 dx = u6 dx = u6 . . du = u6 du = + constant
3 3 21
To finish off we must rewrite this answer in terms of the original variable x and replace u by 3x + 5:
(3x + 5)7
Z
(3x + 5)6 dx = +c
21
HELM (2008): 41
Section 13.5: Integration by Substitution and Using Partial Fractions
In practice the short method is generally used but mathematicians don’t like to separate the ‘dx’
1
from the ‘du’ as in the statement ‘dx = du’ as it is meaningless mathematically (but it works!). In
3
the future we will use the short method, with apologies to the mathematicians!
Task Z
By making the substitution u = sin x find cos x sin2 x dx
du
You are given the substitution u = sin x. Find :
dx
Your solution
Answer
du
= cos x
dx
Now make the substitution, simplify the result, and finally perform the integration:
Your solution
Answer
Z Z
1
cos x sin x dx simplifies to u2 du. The final answer is sin3 x + c.
2
3
Exercise
Use suitable substitutions to find
Z Z
(a) (4x + 1) dx (b) t2 sin(t3 + 1)dt
7
(Hint: you need to simplify sin(t3 + 1))
Answer
(4x + 1)8 cos(t3 + 1)
(a) +c (b) − +c
32 3
42 HELM (2008):
Workbook 13: Integration
®
Example 20 Z 3
Find the definite integral t sin(t2 )dt by making the substitution u = t2 .
2
Solution
du du
Note that if u = t2 then = 2t so that dt = . We find
dt 2t
Z t=3 Z t=3
1 t=3
Z
2 du
t sin(t )dt = t sin u = sin u du
t=2 t=2 2t 2 t=2
An important point to note is that the limits of integration are limits on the variable t, not u. To
emphasise this they have been written explicitly as t = 2 and t = 3. When we integrate with respect
to the variable u, the limits must be written in terms of u. From the substitution u = t2 , note that
when t = 2 then u = 4 and when t = 3 then u = 9 so the integral becomes
9
1 u=9
Z
1 1
sin u du = − cos u = (− cos 9 + cos 4) = 0.129 to 3 d.p.
2 u=4 2 4 2
Exercise
Z 2 Z 1
7 3
Use suitable substitutions to find (a) (2x + 3) dx, (b) 3t2 et dt.
1 0
Answer
(a) u = 2x + 3 is suitable; 3.359 × 105 to 4 sig. figs. (b) 1.718 to 3 d.p.
HELM (2008): 43
Section 13.5: Integration by Substitution and Using Partial Fractions
3. Integrals giving rise to logarithms
Example 21
3x2 + 1
Z
Find dx
x3 + x + 2
Solution
Let us consider what happens when we make the substitution z = x3 + x + 2. Note that
dz
= 3x2 + 1 so that we can write dz = (3x2 + 1)dx
dx
Then
3x2 + 1
Z Z
1
3
dx = dz = ln |z| + c = ln |x3 + x + 2|
x +x+2 z
Note that in the last Example, the numerator of the integrand (3x2 + 1) is the derivative of the
denominator (x3 + x + 2). The result is the logarithm of the denominator. This is a special case of
the following rule:
Key Point 7
f 0 (x)
Z
dx = ln |f (x)| + c
f (x)
Note that it is the modulus of f (x) in the answer.
Task
Write
Z down, purely by inspection,
Z the following Z
integrals:
1 2x 1
(a) dx, (b) 2
dx, (c) dx.
x+1 x +8 x−3
Hint: In each case the numerator of the integrand is the derivative of the denominator.
Your solution
Answer
(a) ln |x + 1| + c, (b) ln |x2 + 8| + c, (c) ln |x − 3| + c
44 HELM (2008):
Workbook 13: Integration
®
Task 4
3t2 + 2t
Z
Evaluate the definite integral dt.
2 t3 + t2 + 1
Your solution
Answer
4
3 2
ln |t + t + 1| = ln 81 − ln 13 = 1.83
2
Sometimes it is necessary to make slight adjustments to the integrand to obtain a form for which
the rule in Key Point 7 is suitable. Consider the next Example.
Example 22
x2
Z
Find the indefinite integral dx.
x3 + 1
Solution
In this Example the derivative of the denominator is 3x2 whereas the numerator is just x2 . We
adjust the numerator as follows:
x2 3x2
Z Z
1
dx = dx and integrate by the rule to get 13 ln |x3 + 1| + c
x3 + 1 3 x3 + 1
Note that the sort of procedure in the last Example is only possible because we can move constant
factors through the integral sign. It would be wrong to try to move terms involving the variable x
in a similar way.
Exercise
Write
Z down the result
Z of finding the following
Z integrals. Z
1 2t 1 2
(a) dx, (b) 2
dt, (c) dx, (d) dx.
x t +1 2x + 5 3x − 2
Answer
1 2
(a) ln |x| + c, (b) ln |t2 + 1| + c, (c) 2
ln |2x + 5| + c, (d) 3
ln |3x − 2| + c.
HELM (2008): 45
Section 13.5: Integration by Substitution and Using Partial Fractions
4. Integration using partial fractions
Sometimes expressions which at first sight look impossible to integrate using the techniques already
met may in fact be integrated by first expressing them as simpler partial fractions, and then using
the techniques described earlier in this Section. Consider the following Task.
Task
23 − x
Express as the sum of its partial fractions.
(x − 5)(x + 4)
23 − x
Z
Hence find dx
(x − 5)(x + 4)
A B
First produce the partial fractions. Write the fraction in the form + and find A, B.
x−5 x+4
Your solution
Answer
A = 2, B = −3
Now integrate each term separately:
Your solution
23 − x
Z Z Z
A B
dx = dx + dx =
(x − 5)(x + 4) x−5 x+4
Answer
2 ln |x − 5| − 3 ln |x + 4| + c
46 HELM (2008):
Workbook 13: Integration
®
Exercises
By expressing the following in partial fractions, evaluate each integral:
Z
1
1. 3
dx
x +x
13x − 4
Z
2. dx
6x2 − x − 2
Z
1
3. dx
(x + 1)(x − 5)
Z
2x
4. dx
(x − 1)2 (x + 1)
Answers
1
1. ln |x| − ln |x2 + 1| + c
2
3 2
2. ln |2x + 1| + ln |3x − 2| + c
2 3
1 1
3. ln |x − 5| − ln |x + 1| + c
6 6
1 1 1
4. − ln |x + 1| + ln |x − 1| − +c
2 2 x−1
HELM (2008): 47
Section 13.5: Integration by Substitution and Using Partial Fractions
Integration of
Trigonometric
Functions 13.6
Introduction
Integrals involving trigonometric functions are commonplace in engineering mathematics. This is
especially true when modelling waves and alternating current circuits. When the root-mean-square
(rms) value of a waveform, or signal is to be calculated, you will often find this results in an integral
of the form
Z
sin2 t dt
In this Section you will learn how such integrals can be evaluated.
' $
• be able to find a number of simple definite
and indefinite integrals
Prerequisites • be able to use a table of integrals
Before starting this Section you should . . .
• be familiar with standard trigonometric
identities
&
%
• use trigonometric identities to write
Learning Outcomes integrands in alternative forms to enable
On completion you should be able to . . . them to be integrated
48 HELM (2008):
Workbook 13: Integration
®
Task
Write
Z down the following
Z integrals: Z Z
(a) sin x dx, (b) cos x dx, (c) sin 2x dx, (d) cos 2x dx
Your solution
(a) (b)
(c) (d)
Answer
1 1
(a) − cos x + c, (b) sin x + c, (c) − cos 2x + c, (d) sin 2x + c.
2 2
The basic rules from which these results can be derived are summarised here:
Key Point 8
Z Z
cos kx sin kx
sin kx dx = − +c cos kx dx = +c
k k
In engineering applications it is often necessary to integrate functions involving powers of the trigono-
metric functions such as
Z Z
2
sin x dx or cos2 ωt dt
Note that these integrals cannot be obtained directly from the formulas in Key Point 8 above.
However, by making use of trigonometric identities, the integrands can be re-written in an alternative
form. It is often not clear which identities are useful and each case needs to be considered individually.
Experience and practice are essential. Work through the following Task.
HELM (2008): 49
Section 13.6: Integration of Trigonometric Functions
Task
1
Use the trigonometric identity sin2 θ ≡ (1 − cos 2θ) to express the integral
Z 2
2
sin x dx in an alternative form and hence evaluate it.
Answer Z
1
The integral can be written (1 − cos 2x)dx.
2
Note that the trigonometric identity is used to convert a power of sin x into a function involving
cos 2x which can be integrated directly using Key Point 8.
(b) Now evaluate the integral:
Your solution
Answer
1
x − 12 sin 2x + c = 21 x − 14 sin 2x + K where K = c/2.
2
Task Z
Use the trigonometric identity sin 2x ≡ 2 sin x cos x to find sin x cos x dx
Answer
1
The integrand can be written as 2
sin 2x
(b) Now evaluate the integral:
Your solution
Answer
Z 2π Z 2π 2π
1 1 1 1 1 1
sin x cos x dx = sin 2x dx = − cos 2x + c = − cos 4π + cos 0 = − + = 0
0 0 2 4 0 4 4 4 4
50 HELM (2008):
Workbook 13: Integration
®
Engineering Example 3
Magnetic flux
Introduction
The magnitude of the magnetic flux density on the axis of a solenoid, as in Figure 13, can be found
by the integral:
Z β2
µ0 nI
B= sin β dβ
β1 2
where µ0 is the permeability of free space (≈ 4π × 10−7 H m−1 ), n is the number of turns and I is
the current.
β2 β1
Problem in words
Predict the magnetic flux in the middle of a long solenoid.
Mathematical statement of the problem
We assume that the solenoid is so long that β1 ≈ 0 and β2 ≈ π so that
Z β2 Z π
µ0 nI µ0 nI
B= sin β dβ ≈ sin β dβ
β1 2 0 2
Mathematical analysis
µ0 nI
The factor can be taken outside the integral i.e.
2
π
µ0 nI π
Z
µ0 nI µ0 nI
B= sin β dβ = − cos β = (− cos π + cos 0)
2 0 2 0 2
µ0 nI
= (−(−1) + 1) = µ0 nI
2
Interpretation
The magnitude of the magnetic flux density at the midpoint of the axis of a long solenoid is predicted
to be approximately µ0 nI i.e. proportional to the number of turns and proportional to the current
flowing in the solenoid.
HELM (2008): 51
Section 13.6: Integration of Trigonometric Functions
2. Orthogonality relations
In general two functions f (x), g(x) are said to be orthogonal to each other over an interval a ≤ x ≤ b
if
Z b
f (x)g(x) dx = 0
a
It follows from the previous Task that sin x and cos x are orthogonal to each other over the interval
0 ≤ x ≤ 2π. This is also true over any interval α ≤ x ≤ α + 2π (e.g. π/2 ≤ x ≤ 5π, or
−π ≤ x ≤ π).
More generally there is a whole set of orthogonality relations involving these trigonometric functions
on intervals of length 2π (i.e. over one period of both sin x and cos x). These relations are useful
in connection with a widely used technique in engineering, known as Fourier analysis where we
represent periodic functions in terms of an infinite series of sines and cosines called a Fourier series.
(This subject is covered in 23.)
We shall demonstrate the orthogonality property
Z 2π
Imn = sin mx sin nx dx = 0
0
1 2π
Z
Imn = [cos(m − n)x − cos(m + n)x] dx
2 0
2π
1 sin(m − n)x sin(m + n)x
= −
2 (m − n) (m + n) 0
= 0
follows by use of a similar identity to that just used. Here again m and n are integers such that
m 6= n.
52 HELM (2008):
Workbook 13: Integration
®
Example 23
1
Use the identity sin A cos B ≡ (sin(A + B) + sin(A − B)) to show that
2
Z 2π
Kmn = sin mx cos nx dx = 0 m and n integers, m 6= n.
0
Solution
1 2π
Z
Kmn = [sin(m + n)x + sin(m − n)x] dx
2 0
2π
1 cos(m + n)x cos(m − n)x
= − −
2 (m + n) (m − n) 0
1 cos(m + n)2π − 1 cos(m − n)2π − 1
= − + =0
2 (m + n) (m − n)
Task
Derive the orthogonality relation
Z 2π
Kmn = sin mx cos nx dx = 0 m and n integers, m = n
0
Hint: You will need to use a different trigonometric identity to that used in Example
23.
Your solution
HELM (2008): 53
Section 13.6: Integration of Trigonometric Functions
AnswerZ
2π
Kmn = sin mx cos mx dx
0
Putting m = n 6= 0, and then using the identity sin 2A ≡ 2 sin A cos A we get
Z 2π
Kmm = sin mx cos mx dx
0
1 2π
Z
= sin 2mx dx
2 0
2π
1 cos 2mx 1 1
= − =− (cos 4mπ − cos 0) = − (1 − 1) = 0
2 2m 0 4m 4m
Z 2π
1
Putting m = n = 0 gives K00 = sin 0 cos 0 dx = 0.
2 0
Note that the particular case m = n = 1 was considered earlier in this Section.
3. Reduction formulae
You have seen earlier in this Workbook how to integrate sin x and sin2 x (which is (sin x)2 ). Appli-
cations sometimes arise which involve integrating higher powers of sin x or cos x. It is possible, as
we now show, to obtain a reduction formula to aid in this Task.
Task Z
Given In = sinn (x) dx write down the integrals represented by I2 , I3 , I10
Your solution
I2 = I3 = I10 =
Answer Z Z Z
I2 = sin2 x dx I3 = 3
sin x dx I10 = sin10 x dx
54 HELM (2008):
Workbook 13: Integration
®
Task
In the notation used earlier in this Workbook for integration by parts
Z (Key Point
df
5, page 31) put f = sinn−1 x and g = sin x and evaluate and g dx.
dx
Your solution
Answer
df
= (n − 1) sinn−2 x cos x (using the chain rule of differentiation),
dx
Z Z
g dx = sin x dx = − cos x
Z
Now use the integration by parts formula on sinn−1 x sin x dx. [Do not attempt to evaluate the
second integral that you obtain.]
Your solution
Answer
Z Z Z Z
n−1 n−1 df
sin x sin x dx = sin (x) g dx −
g dx
dx
Z
= sin (x)(− cos x) + (n − 1) sinn−2 x cos2 x dx
n−1
Z
We now need to evaluate sinn−2 x cos2 xdx. Putting cos2 x = 1 − sin2 x this integral becomes:
Z Z
n−2
sin (x) dx − sinn (x) dx
Z is expressible as In−2 − In so finally, using this and the result from the last Task we have
But this
In = sinn−1 (x) sin x dx = sinn−1 (x)(− cos x) + (n − 1)(In−2 − In )
HELM (2008): 55
Section 13.6: Integration of Trigonometric Functions
Key Point 9
Reduction Formula
Z
Given In = sinn xdx
1 n−1
In = − sinn−1 (x) cos x + In−2
n n
This is our reduction formula for In . It enables us, for example, to evaluate I6 in terms of I4 , then
I4 in terms of I2 and I2 in terms of I0 where
Z Z
0
I0 = sin x dx = 1 dx = x.
Task
Use the reduction formula in Key Point 9 with n = 2 to find I2 .
Your solution
Answer
1 1
I2 = − [sin x cos x] + I0
2 2
1 1 x
= − [ sin 2x] + + c
Z 2 2 2
1 x
i.e. sin2 x dx = − sin 2x + + c
4 2
as obtained earlier by a different technique.
56 HELM (2008):
Workbook 13: Integration
®
Task Z
Use the reduction formula in Key Point 9 to obtain I6 = sin6 x dx.
Your solution
Answer
1 5
Using Key Point 9 with n = 6 gives I6 = − sin5 x cos x + I4 .
6 6
1 3
Then, using Key Point 9 again with n = 4, gives I4 = − sin3 x cos x + I2
4 4
Now substitute for I2 from the previous Task to obtain I4 and hence I6 .
Your solution
Answer
1 3 3
I4 = − sin3 x cos x − sin 2x + x+ constant
4 16 8
1 5 5 5
∴ I6 = − sin5 x cos x − sin3 x cos x − sin 2x + x + constant
6 24 32 16
Definite integrals can also be readily evaluated using the reduction formula in Key Point 9. For
example,
Z π/2 Z π/2
n
In = sin x dx so In−2 = sinn−2 x dx
0 0
We obtain, immediately
π/2
1 n−1 n−1
In = − sin (x) cos x + In−2
n 0 n
π (n − 1)
or, since cos = sin 0 = 0, In = In−2
2 n
This simple easy-to-use formula is well known and is called Wallis’ formula.
HELM (2008): 57
Section 13.6: Integration of Trigonometric Functions
Key Point 10
Reduction Formula - Wallis’ Formula
Z π/2 Z π/2
n
Given In = sin x dx or In = cosn x dx
0 0
(n − 1)
In = In−2
n
Task Z π/2
If In = sinn x dx calculate I1 and then use Wallis’ formula, without further
0
integration, to obtain I3 and I5 .
Your solution
Answer
Z π/2 π/2
I1 = sin x dx = − cos x =1
0 0
58 HELM (2008):
Workbook 13: Integration
®
Task
The total power P of an antenna is given by
Z π
ηL2 I 2 π
P = 2
sin3 θ dθ
0 4λ
where η, λ, IZ are constants as is the length L of antenna. Using the reduction
formula for sinn x dx in Key Point 9, obtain P .
Your solution
Answer
Ignoring the constants for the moment, consider
Z π
I3 = sin3 θ dθ which we will reduce to I1 and evaluate.
0
Z π π
I1 = sin θ dθ = − cos θ =2
0 0
HELM (2008): 59
Section 13.6: Integration of Trigonometric Functions
4. Harder trigonometric integrals
The following seemingly innocent integrals are examples, important in engineering, of trigonometric
integrals that cannot be evaluated as indefinite integrals:
Z Z
2
(a) sin(x ) dx and cos(x2 ) dx These are called Fresnel integrals.
Z
sin x
(b) dx This is called the Sine integral.
x
Definite integrals of this type, which are what normally arise in applications, have to be evaluated
by approximate numerical methods.
Fresnel integrals with limits arise in wave and antenna theory and the Sine integral with limits in
filter theory.
It is useful sometimes to be able to visualize the definite integral. For example consider
Z t
sin x
F (t) = dx t>0
0 x
Z 0
sin x sin x
Clearly, F (0) = dx = 0. Recall the graph of against x, x > 0:
0 x x
sin x
x
t π 2π x
Figure 14
For any positive value of t, F (t) is the shaded area shown (the area interpretation of a definite integral
was covered earlier in this Workbook). As t increases from 0 to π, it follows that F (t) increases from
0 to a maximum value
Z π
sin x
F (π) = dx
0 x
whose value could be determined numerically (it is actually about 1.85). As t further increases from
sin x
π to 2π the value of F (t) will decrease to a local minimum at 2π because the curve is below
x
the x-axis between π and 2π. Note that the area below the curve is considered to be negative in
this application.
Continuing to argue in this way we can obtain the shape of the F (t) graph in Figure 15: (can you
60 HELM (2008):
Workbook 13: Integration
®
F (t)
1.85
π
2
π 2π t
Figure 15
Z ∞
sin x π
The result dx = is clearly illustrated in the graph (you are not expected to know
0 x 2
how this result is obtained). Methods for solving such problems are dealt with in 31.
HELM (2008): 61
Section 13.6: Integration of Trigonometric Functions
Exercises
You will need to refer to a Table of Trigonometric Identities to answer these questions.
Z Z π/2 Z
2 2
1. Find (a) cos xdx (b) cos tdt (c) (cos2 θ + sin2 θ)dθ
0
Z
2. Use the identity sin(A + B) + sin(A − B) ≡ 2 sin A cos B to find sin 3x cos 2xdx
Z
3. Find (1 + tan2 x)dx.
4. The mean square value of a function f (t) over the interval t = a to t = b is defined to be
Z b
1
(f (t))2 dt
b−a a
Find the mean square value of f (t) = sin t over the interval t = 0 to t = 2π.
Z
5. (a) Show that the reduction formula for Jn = cosn x dx is
1 (n − 1)
Jn = cosn−1 (x) sin x + Jn−2
n n
(b) Using the reduction formula in (a) show that
Z
1 4 8
cos5 x dx = cos4 x sin x + cos2 x sin x + sin x
5 15 15
Z π/2
n n−1
(c) Show that if Jn = cos x dx, then Jn = Jn−2 (Wallis’ formula).
0 n
Z π/2
5
(d) Using Wallis’ formula show that cos6 x dx = π.
0 32
Answers
3. tan x + c.
1
4. 2
.
62 HELM (2008):
Workbook 13: Integration