Module 11 - Integration 4 (Self Study)
Module 11 - Integration 4 (Self Study)
MODULE 11 INTEGRATION IV
Module Topics
1. Double integration
2. Polar coordinates
3. Triple integration
You have all integrated functions of a single variable, y = f (x) say. In the real world, however, most physical
quantities depend on more than one variable and therefore it is important in many branches of engineering
to be able to integrate these more general functions. The so-called multiple integrals are not covered in
J. and so you should study the sections below.
1. Recall that in Module 3 we defined the definite integral of a function of one variable, f (x), by the limit
Z b n
X
f (x) dx = lim lim f (xi ) δxi ,
a δxi →0 n→∞
i=1
where a = x0 < x1 < x2 < ... < xn = b, δxi = xi − xi−1 and xi−1 ≤ xi ≤ xi . This integral represents the
area between the curve y = f (x) and the x-axis and between x = a and x = b, as shown in Figure 11.1.
.. y
........
.........
..
...
....
..
...
...
...
...
...
... ...................................
... y = f (x) ......... .....
... ........
... ...
........... ..
...
.
... ...... ...
... .............
. ............ ...
.
.....
... .. ......................
. .
.........
. ...
... .. ...........................
. . . .
.......... . ...
... .... .
......... ...
.... ... ............ ....
.. ... ...........
..
... ... ...
... ... ............ ...
... ... ........... ...
... ... ............ .. .
................................................................................................................................................................................................................................................................................................................................
... .
0 a
... x i b x
...
...
..
Figure 11.1
2. Let us now define a double integral. Consider a region R in the xy-plane bounded by a closed curve,
and subdivide this into a large number (N ) of small regions δS1 , δS2 ,... δSN by means of a grid of lines
parallel to the two axes at distances δx, δy apart respectively. Then, except near the boundary of S,
δSi = δx δy (see Figure 11.2).
–1–
..
........................
.... z ..................................... ........................
....... ... ................................................................................................ ..
.
........ ..
. .
....
.... ... ..
...
... ... ...
... ... ..
... ... .....
....................... ..
... ... ..................... ...
... ........... ......
...
... ........ .............................. ...
... ... ....... ... ...
... ... ..
........
. ....
. .
.... ...
. ....
.
... z = f (x, y) .
...
.
........
.....
.
... ... ... ...
..
.. .
...
...
.
... ... .... ... .. ... . ... .. ...
... .. ..... ... . ... . ..... .. ...
... ..... .
. . ..... ... .... ...
. ...
....... .... ...
... . ... ....
...
... .......... ................................. .. .............. ....
... ... ............... .. ....
... . .. ............ ... ...
... ... .. ..................................... ... .. ............
... .........................................................................................................................................................................
... .... ... . .. . ....
... ... ... ... .. ... ...
... ... ... . ... ...
.... ... ... .
. ...
.. ... .
. .. .
. ..
. ...
... ... . . .
. .
. ...
... ... ... . ... . ...
... .. . .
.
...
y ..................
.
.
.
.
.
.
.
. .
. .
.
.
.
.
....
..
... ........... .
.
.
.
.
. . ... .
...
... ..
.... . . . .
. .
.
.. .
...... .
.
. .
. . .
. ........ .. .
... ...
... ......
... .
.
. . .
.
. .......... . ........ .
...
. ...
... ..... . .
. . . .
... δy ......................................................................................................................................................................................................................................................................................
. .
... .
..... . ..... .............
. ..................................................................................................... δS = δxδy
... .
...
.
... .
............................................................................................................................................................................................................................................................... i
... ... ........ .. ..............................
... ...... . .. .......
... ....... .......... ....... ......... .....................
... .
......... ......... .................................................................... ...........
...
.
.....
......
.
.......
.......
... .
......
......
R
.... ....... ...... .......
....... ....... .......
.. ...... ...... .......
... ............. .. ..... .
... ....
.
............................................................................................................................................................................................................................................................................................................................................................................................
..
0 δx x
Figure 11.2
Let (xi , yi ) be a point inside the ith element of area δSi and let fi = f (xi , yi ), so that fi is the height
of the surface above (xi , yi ). Then fi δSi is approximately equal to the volume of a thin vertical tube of
rectangular cross-section constructed on δSi as base, the upper boundary being where this tube intersects
the surface.
Let f1 , f2 , ..., fN be the values taken by the function f at arbitrary points inside the elements of area δS1 ,
δS2 , ..., δSN respectively. Let
N
X
VN = f1 δS1 + f2 δS2 + ... + fN δSN ≡ fi δSi .
i=1
3. To evaluate a double integral you simply have to perform two successive single integrations. The
difficult problem is getting the limits right . The first step is to draw the region R in the xy-plane over which
the double integral is to be evaluated. This is easy when the region is rectangular, as illustrated in the
following example.
Z y=1 Z x=2
Example 1 Evaluate (1 + x + y) dx dy.
y=0 x=1
The region of integration in this example is rectangular, since all limits are constants and hence independent
of x and y.
–2–
When a double integral is written in the form displayed in Example 1, with dx as the inner infinitesimal
and dy as the outer, then the function f (x, y) = 1 + x + y is first integrated with respect to x and the
resulting quantity then integrated with respect to y.
Although the limits in Example 1 are written in their fullest form, you will often encounter integrals written
Z 1Z 2
(1 + x + y) dx dy, where you can assume from the order of the infinitesimals and integrals which limits
0 1
apply to which variables. The situation is completely clear if you
introduce an extra set of brackets. For
Z 1 Z 2
instance, the integral could be written (1 + x + y) dx dy.
0 1
The region of integration for Example 1 is easy to find, but in general you proceed as follows. Draw on a
diagram the lines defined by the four limits appearing explicitly in the integral, i.e. x = 1, x = 2, y = 0 and
y = 1 in the current example. For the inner integral y is constant and x varies between 1 and 2. Similarly,
for the second integral y has outer limits y = 0 and y = 1. The region of integration is then the shaded
area in Figure 11.3, with the thin strip showing the part of the region along which the first integral is taken.
...
.....
y
..........
..
...
....
.. ... ...
1 .
......................................................................................................................................................
.. ... .... .... .... .... ..
... ....... .... .... .... ...
... ..........................................................................
... ....................................................................
... .. ... .. ... ... ..
.... ...... ...... ...... ...... .......
.
.. .... ...... ...... ...... ...... ...
... ...... ..... ..... ..... ....
... ... ..... ..... ...... ........
.. .. .... ... ... .... .. .
.....................................................................................................................................................................................................
... ... ... .
0 ... 1 ... 2
... x
... ... ...
Figure 11.3
Let us now evaluate the double integral. First, consider the inner integral which, as stated above, is an
integration with respect to x along a strip with y kept fixed. The calculation of this inner integral
proceeds as follows:
Z 2 x=2
x2 22 1 5
(1 + x + y) dx = x + + xy = 2+ + 2y − 1 + + y = + y.
1 2 x=1 2 2 2
Note that in evaluating the above integral y remained constant. The full calculation of the double integral
then gives
Z 1 Z 2 Z 1 1
5 5y y 2 5 1
(1 + x + y) dx dy = +y dy = + = + − (0) = 3.
0 1 0 2 2 2 0 2 2
Integrals of the type appearing in Example 1 can have a number of interpretations. The one mentioned
earlier is the calculation of volumes and for this situation you can think of the above answer, 3, as the
volume of a solid body with its base on the xy-plane (i.e. z = 0), with height z = f (x, y) = 1 + x + y and
its sides being the planes x = 1, x = 2, y = 0 and y = 1.
Retaining the volume interpretation, the volume in Example 1 was calculated above by dividing it into thin
slices parallel to the x-axis (i.e. parallel to the xz-plane). The volume of one of these slices was calculated
through the first integral, over the variable x, and then the second integral is equivalent to adding together
the volumes of all the thin slices to give the total volume of the solid.
It is not difficult to see that the total volume could also have been divided into slices parallel to the y-axis
(or yz-plane). In this case the first integral would give the volume of the thin slice parallel to the y-axis,
with the second integral providing the sum of the volumes of these slices. Hence, a given double integral can
also be evaluated, at least in principle, by changing the order of integration. During this process the
limits will normally need to be changed too.
–3–
When the first integral in Example 1 is with respect to y, the picture for the region of integration becomes
that shown in Figure 11.4.
...
......
y
..........
....
...
...
..
...... 1 ...............................................................................
.... ... ...... ......... ..... ..... ...
... .. .. . . .
... ...... ........... .......... ....... ........
. . .. .. ..
... ........... ............. ...... ...... ...
... . ... ... .. . ... .... ...
... ... ...... ............... ...... ...........
... ........ ....... ... ..... ..... ...
... .. ... .. ...... .... .....
.
... .... ...... ............. ...... ....... .... .
.................................................................................................................................................................................................
.. .
0 ...
....
1 2 x
Figure 11.4
The function being integrated (i.e. the height of the solid if it is a volume being calculated) is unaltered. So
the double integral can also be expressed
Z 2 Z 1
(1 + x + y) dy dx.
1 0
Note the order of dx and dy, and the limits of integration, have been changed. The limits are normally
determined from the figure. For this new double integral x is fixed in the inner integral and from the figure
y varies along the thin strip between 0 and 1. The calculation of the integral now proceeds as follows:
Z 2 Z 1 Z 2 y=1 Z 2
y2 1
(1 + x + y) dy dx = y + xy + dx = 1+x+ − 0 dx
1 0 1 2 y=0 1 2
Z 2 2
3 3x x2 3 1
= + x dx = + = (3 + 2) − + = 3.
1 2 2 2 1 2 2
The final answer, obviously, must be the same as the value calculated earlier, but the details of the interme-
diate steps are usually very different.
Changing the order of integration can be very useful in practice, since it is often easier to perform the
calculation one way than the other. As you will see later, in some cases it is only possible to evaluate the
integral one way.
***Do Exercise B Evaluate Exercise A(i) after changing the order of integration.
4. As seen in Section 3 the evaluation of double integrals is relatively straightforward when the region of
integration, R, is rectangular. Consider now the more usual situation in which R is non-rectangular.
Suppose R is bounded by the two curves y = g1 (x) and y = g2 (x) which intersect where x = x1 and
x = x2 .
–4–
... y
.......
.........
....
...
...
.............
... ............. ........
... ............
... ..
.. .....
.
.............. .
...... .. B
y = g (x)... . .
... 2 ........... .... .... ..... .
... ..
....
.............
. .
. .
. ....
....... ...
.. ...
... .......... .
.... ... ..............
... ...
........... . . ..
..
... .........
...............
... .......... .
. ..
.... . .
.........
R A ...
... ....... ....
.....
...... ...
... ..................................................................
... ... 1 y = g (x) ...
... .
.
. ...
...
... ... ...
... ... ...
...
... ... ...
... ..
. ...
...
... ... ...
. . . .
...................................................................................................................................................................................................................................................................................................................................................................................................
. .
.. .
0 ....
.
1 x 2 x x
Figure 11.5
For a given value of x, the variable y can take all values between g1 (x) and g2 (x). The region R is therefore
that for which )
g1 (x) ≤ y ≤ g2 (x)
(1)
x1 ≤ x ≤ x2
Consider now the double integral ZZ
f (x, y) dx dy ,
R
which can represent, for instance, the volume under the surface z = f (x, y) above the region R.
With the region described by equation (1), the integral can be written with the integration with respect to
y being carried out first, followed by integration with respect to x, as:
Z x=x2 (Z y=g2 (x) )
I= f (x, y) dy dx . (2)
x=x1 y=g1 (x)
In other words, for a fixed x in the range x1 ≤ x ≤ x2 we evaluate the area AP QB in Figure 11.6:
Z y=g2 (x)
AreaAP QB = f (x, y) dy ,
y=g1 (x)
which is clearly a function of x only. Multiply this by dx to give the contribution to the volume from the
thin slice AP QBA0 P 0 Q0 B 0 , and then integrate with respect to x to add the contributions from the various
slices, so giving the total volume.
.. z
.......
..........
... 0
.... Q
..................................................
Q
..
... ... ....... ............................... ..... . ..........
...
...........
...........
..........
......
.... ...
...... ...... ..
...... ...
..........
........... ................. z = f (x, y)
... .......... ...... .... .......... .. .. ..
...................
... ...... . ....
...
.........
......
.
. .
.... .
.
. ..................................
... .... ... ..... ..... .
. ....
...
. ...
... ............
. . . ....... ......... ... .
.
. .....
. ... ...... ...
... ..........
.... .
...... ...... .... .............. ...
... . ..... ........ ......... ... ..............
. ...... ..
....... .
. ....... ....
.. ...
... .. ........
. ..... .... . ......... . . . .
.
. ...
... ..
. ...... ...
...... ........
. ..... ..... ......... ... .
.
. ...
.. ....... ...... .......................
. . .
... . . .
..
. .
... ........................................ .
. .
... .
. .
. ...
... ..
. . . .
.
.
.
.
. ...
P ..... .
... .... .... 0 .
.
. .
. ...
.... ... ... P .... .... .... 0 ...
.. ... ... .
. .
. B
.. .... .... .......
.
. B ...
... ... ... ..
. . . . . ...
... ... ... ..
. .... ............ .... ......... ....
. ...
.... ......... ... .... .
... ... ... ..
. .... . .
... .. . . . ... y = g2 (x)
... .
. .
. .. ... .... ... ...
.
.
. ........
....
...... .... .
...................
.... .... ....
... .... .... ... ..
.
. ........... ..... . .. .. ..
..............
.. .
... .
. . . .... .... ...
. .
.... .........
... ..
. ... . ...... ..
. ... ......
.......... .
. .
. . .
. ....
..
........
... ...... .... ..
... ....... .... .... ..... ............... ....... .............. ......
.... ..... ....... ...... ............... ............
...
y . ................
... ................... ....... ........... ... ...... .. ...... .................. ....
.....................................
....... ....... ... ............ .... ............. ........................................
...
... .......
..
... .. . .
.......
........ ..............
... ........................................
. .. ..............
.......
...
.. ...
.................. y = g1 (x)
... ....... ....... ....... . 0 .......
... ...
........ ........... ...
........
A A ..
........
. . ..........
. ..... . ...... .
...... ...... .
... ....... ....... ....... ....... .......
....... ....... ....... ....... .......
.... ...... ....... ....... .......
.. ....... ...... ...... ....... ......
... ..
........
. ..
.. .
. .
.. .
.. ............ ......
. . . ...... .
... ....... ......
...... ....... ....... .......
...... ...... ......
... ...... ....... ....... ......
...... .......
... ............. ....... ...... ....... ......
....... . ....... ........ ... .. ...... .
....................................................................................................................................................................................................................................................................................................................................................................................
. ... ... ... .
..
.
0 x1 x x + δx x2 x
Figure 11.6
–5–
As in Section 3 the order of integration can be changed, but much more care is now needed to obtain the
correct limits of integration. Consider again the region R in Figure 11.5. This region can equally well be
specified by the inequalities )
h1 (y) ≤ x ≤ h2 (y)
y1 ≤ y ≤ y2
where the equations y = g1 (x) and y = g2 (x) have been solved for x to give x = h2 (y) and x = h1 (y)
respectively, see Figure 11.7. The actual curves are, of course, the same as before. The slices are now taken
parallel to the xz-plane and the picture of the region of integration, the xy-plane, becomes
.. y
......
..........
...
....
y .
2 ....... .... .... .... .... .... .... .... .... .... .... .... .... .... .... .... .... .... .... .... .......................................
... ...
..........
.
...
... ........................................
.... x = h (y) 1 ..................................................
.. ..
............... .
...
.
. .
... ......... ......
........ ......
... R ........ .......
... ........ .........
... ............. ......................
y .... ............
1 ........ .... .... .... .... .... .... ........................................... x = h (y) 2
....
...
...
...
...
...
...
...
...
.. .
..........................................................................................................................................................................................................................................................................................................................
.
0 .... x
..
....
...
. Figure 11.7
in which the integration with respect to x is performed first , keeping y constant and then, after filling in the
limits in x, carrying out the integration with respect to y.
The procedures for calculating a double integral in a non-rectangular region are illustrated in the following
Examples.
Example 2 The volume of a solid body of triangular cross-section with height z = xy 2 and bounded by
the planes z = 0, x = 0, y = 0 and y = 1 − x is given by
Z x=1 Z y=1−x
V = xy 2 dy dx .
x=0 y=0
Sketch the region in which the body intersects the xy-plane and evaluate the integral. Check the result by
reversing the order of integration.
Draw the lines x = 0, x = 1, y = 0 and y = 1 − x as in Figure 11.8, noting that for the purposes of finding
this region the integrand xy 2 is irrelevant.
–6–
...
.... .............
. y
.... .. ....
... ..
.... .. ...
...
..... 1
......
... .....
...
...
... .......
... ..... ....
........
y =1−x ....
...
...
..... ...
.. .
..
. . ...
... ....... ...... ...
... ..... .
.. .... ...
......... ..
... .... ...
.. ...
...
.....
. ..
..
.. ..... ...
... ... ..
x=1 ...
... ......
.
...... ....
. . ...
... ..... ....
. ..
............. ...
...... .
.... .
.... ... ... ..... ...
... ..
... .
.... ... ... ......... ...
... ......
. ... . ........ .....
...
.. .
.
....... .... ....
.
... ...... .... . .
..... ... . ..
.. .... ....
.
....... .
. . .
. .
.. .
..........................................................................................................................................................................
.
. .
0 .... 1 .... ......
x
... ... .....
.. .. ...
Figure 11.8
In order to reverse the order of integration you first have to find the appropriate limits. Simply reversing
everything and writing V as
Z y=1−x Z x=1
xy 2 dx dy
y=0 x=0
....
...
.......
y
.... ........
... ..
.... ..
..... 1
...........
... ......
... ..... .....
.....
x=1−y ....
...... ...
.... .... ....
.... ......
... ....... ....
......... .... .......
.... ... ...
.... ..
..................................................................
.................................................................
... ...... ... . ..
.... ......
.....
. .
..... .... ...
....
....
.. .... ..
.... ......
.... ...... ...
.. . .... ...
... ..... .... .... ..
........ .... . ... .... ....
... . .. .... ... .
...........................................................................................................................................................................
... .... .
0 ... 1 ...
.... x
... ...
Figure 11.9
–7–
and this evaluates as
Z y=1 x=1−y Z 1 Z 1
2 1 2 1 2 2 1 2
V = y x dy = y (1 − y) dy = y (1 − 2y + y 2 ) dy
y=0 2 x=0 0 2 0 2
Z 1 1
1 1 1 3 1 4 1 5 1 1 1 1 1
= (y 2 − 2y 3 + y 4 ) dy = y − y + y = − + − (0) = .
2 0 2 3 2 5 0 2 3 2 5 60
...
.....
y
..........
..
...................................................................................
1 .... ....... ......
. ...
... ........
....
... ..... .. .... .......
..... .... ........... 2
...x=y ...............
... ...... ...
... ..... ........
.
. .
.. ........
...........................
......................
... .......
... .........
...
........
....
..
..
..... .
.................................................................................................................................................................................
.... .
... 0
........
1 x
... .. .
..
.... .......
.. ..... ....
........ .........
... ..... .......
... ..... .
........
.......
. .. .... ............
... ... .........
.... ....
.... ...... ...
....
.... ..............
.. ..... .... .... . ..
−1 ...........................................................................................
...
...
....
Figure 11.10
The inner integral may be evaluated as follows:
Z x=y2 x=y2
1 2 1
(x + y) dx = x + xy = y4 + y3 ,
x=0 2 x=0 2
and hence Z y=1 5 1
1 4 3 y y4 1 1 −1 1 1
V = y +y dy = + = + − + = .
y=−1 2 10 4 −1 10 4 10 4 5
To change the order of integration draw lines parallel to the y-axis. The resulting strip splits into two parts
and the y-integration must now also take place in two parts.
.
....... y
..........
............................................................................
1 ... .....
... .....
.. .
.... .. .. .... .......
.
...... ... .... ........
..
... .... .... .......................
..... .. ... ............
. ...
... ....... .................
......
... .....
..
..
.
√
.....
.....
...
y=
.... .........
. x
........
..........
......
...
.. . .
...................................................................................................................................................................................
.... .
.. 0 ..... 1 x
...........
√
.. ....
... .........
...... ........
y=− x
.. ... ..
.... ...... ..........
... ..... .......
.......
. .... ............
...
.. ....... ...... .........
... ..... .... .. .. .........
... .. .. .... ..............
...... .... .. ..... ..
−1 .........................................................................................
...
...
...
Figure 11.11
–8–
It is apparent from Figure 11.11 that the double integral can also be expressed
Z x=1 (Z y=−√x Z y=1 )
V = (x + y) dy + √
(x + y) dy dx,
x=0 y=−1 y= x
In any given example, the shape of the boundary and the functions in the integrand will determine which
order of integration is the simpler one to use. In some cases the functions in the integrand may mean that
integration one way round is not possible.
Z x=1 Z √
y= x
***Do Exercise C Evaluate the double integral (x2 + y 2 ) dy dx .
x=0 y=x
where ∆ is the inside of the triangle bounded by the lines x = 1, y = 2 and y = x. Express I in two
3 1
different ways, and evaluate each expression in turn to show that I = − ln 2.
4 2
Z 1 Z 1
3
***Do Exercise E Evaluate √
ey dy dx .
0 x
5. When the region of integration is made up of arcs of circles and straight lines through the origin it is
usually more convenient to use polar coordinates.
With a polar representation each point in the plane is identified by its distance from the origin r (by
definition r > 0) and from the angle that this position vector makes with the polar axis (usually coincident
with the x-axis). The conversion formulae are
x = r cos θ, y = r sin θ.
–9–
An element of area in polar coordinates is obtained by considering the region bounded by arcs of radii r
and r + δr centred on O and by radial lines making angles θ and θ + δθ with the x-axis.
.. y
......
..........
...
.... ......
.. .... .. ....
... .... ... .....
... ...... .. .....
... ..... .... .... .... .....
.
.
.. . . .
... .... .. .. .. ... .....
... ... .. . . . . ..
..... .. ... ... .. ... ...
... .... .......... ... ... .... ..............
... .
.... ..
. .
.
... .. .. ....... .
. .
. .
...
.
... ...
... ... .. .........
...
... ....
rδθ
.... .............
......
... ........ .
.......
........
. δr
... ...
....... ......
. .
.
... ..
. ....
..
. ... ..
... ... ... ............
... δθ .... ....
.... ......
... ... .......
... .
..... ............
. .
... .... ..........
... .................. ....
... ................ ...
θ
...................................................................................................................................................................................................................
......
0 x
Figure 11.12
The shaded region is approximately a rectangle if δr and δθ are small enough and its area is therefore
approximately δS = (r δθ) δr. Hence, in the limit as the area becomes smaller, the expression dS = dx dy
used in earlier sections is replaced by dS = r dr dθ, and the double integral of a function g(r, θ) over a
region R bounded by a closed curve lying in the plane can therefore be written
ZZ
g(r, θ) r dr dθ .
R
Example 4 Evaluate
ZZ
y2
I= dx dy ,
R x2
where R is the part of the annulus a2 ≤ x2 + y 2 ≤ b2 lying in the first quadrant and below the line y = x.
The curves x2 + y 2 = a2 and x2 + y 2 = b2 are arcs of circles of radii a and b respectively centred at the
origin. Hence the region of integration is that shown in Figure 11.13.
...
......
y
.........
....
.. .......
... ... .......
...
.... .... ...
...
.. ..... ...
..
... .... ... ...
.... ..... ...... ...
...
.. ..... .... ...
.. ... ...
... ..... . ... ...
... ..... .... ...
... ..... .... ...
... ... . ...
... y=x
... ......... ...
. R
... ...
...
. ...
... ...... ....
...
...
... ....
. .
... ..
... ..... ... ..
... .
.... .
...
...
..
... ......... .. ..
... ..... ... . ..
.. ... .. ...
... ... ....
π/4 ... ...... ..
..
.
..
...
.......... .. .. .. .
............................................................................................................................................................................................................
.
0 a b x
Figure 11.13
The region R is specified in polar coordinates by 0 ≤ θ ≤ π/4 and a ≤ r ≤ b. The integrand becomes
y2 r2 sin2 θ
= = tan2 θ,
x2 r2 cos2 θ
– 10 –
and the area element is dS = dx dy = r dr dθ. Therefore, the value of the integral is found by
Z θ=π/4 Z r=b r=b Z θ=π/4
2 2 r2
I= tan θ rdr dθ = tan θ dθ
θ=0 r=a θ=0 2 r=a
Z θ=π/4 2 2 Z θ=π/4
b − a2 b − a2
= tan2 θ dθ = (sec2 θ − 1) dθ, using 1 + tan2 θ = sec2 θ,
θ=0 2 2 θ=0
2 2
b − a2 π/4 b − a2 h π i (4 − π)(b2 − a2 )
= [tan θ − θ]0 = 1− − (0) = .
2 2 4 8
***Do Exercise F Calculate the area of the cardioid r = a(1 − cos θ) 0 ≤ θ ≤ 2π, using the formula
ZZ
area = r dr dθ .
R
6. There are no conceptual differences between the evaluation of integrals in two-dimensions and integrals in
three or more dimensions. The integration method involves writing the multiple integral as a succession of
single integrals whose integration ranges generally depend on the other integration variables. Work through
the calculation of a triple integral in the staightforward Example below.
ZZZ
Example 5 Evaluate (xy 2 + z 3 ) dx dy dz , where B is the rectangular box 0 ≤ x ≤ a, 0 ≤ y ≤ b,
B
0 ≤ z ≤ c.
z ..
........
.........
...
...
c .
.
...
.................................................................................................................
...... .... .... ..
..
... ... ..... ...
.
... ..... ..
.
..... .... ...... ...
..........
.
. .... ..... ...
..... . .... .
... .
. . . .. ...
...
..... .
.
. . .... ... ...
..
...... .
.
. ........ ...
........ .
.
. .......
. ...
.
..
.... .
.
. ..
....
. ...
.
.... . .
....
...
. .
. .... ...
.
..
. ... .
.
. . ..
..... ...
..............................................................................................................
.. .. .. ...
... .. .
.
.... ...
. .
.
....
.... .
...
...
. .
.... . . . . .... . .. . . . .
.
......
.. .
......
.... ......... .. . . . .. .. . . .... . ..
... .
...........
..
.......... .... .
.
.
. .
..
.
... ....
. .
... .... .
. .... ....
... ..
..... ... .
. .
.
. ................ ..
...... ...
... ......
. . .......
. .
. .............. ...
.... b .... .......
.... . .
... .
. .
...... . .. .
... ...... ...............................................................................................................................................................................................................
.... ..... . .. .............. . y
.. ..
.......
. ...
....... .... .
. ..
.. .
..
. .. . ...... .........
. .
... ...
..... .....
.....
. .. ............ .......... ......
..... ...... ........... . ......
.... ..... ............. .......... ......
.. .......... ...... ........... .. . .....
........... ...
....... .
.................. .... .......... . ........
.. .. .. .. ..
................................................................................................................................... .....
. ... .....
... ..... ......
... .......... ... ..........
... ....... ... .......
a . .......
..
.................................................................................................................
...
.............
........
x
Figure 11.14
2. Find the inequalities satisfied by a second variable, for example y, compatible with the chosen value
of z. Since the slice obtained in the previous step is rectangular the limiting values of y are independent of
the variable z: 0 ≤ y ≤ b.
3. Find the inequalities satisfied by the final variable, compatible with the values of y and z. Since the
slice is rectangular the limiting values of x are independent of both y and z: 0 ≤ x ≤ a.
– 11 –
The calculation proceeds as follows (clearly the order of integration could be changed)
ZZZ Z c Z b Z a Z c x=a
Z b
2 3 2 3 x2 y 2 3
I= (xy + z ) dx dy dz = (xy + z ) dx dy dz = + xz dy dz
B 0 0 0 0 0 2 x=0
Z cZ b 2 2 Z c 2 3 y=b Z c 2 3
a y 3 a y 3 a b 3
= + az dy dz = + ayz dz = + abz dz
0 0 2 0 6 y=0 0 6
2 3 z=c
a b z abz 4 a2 b3 c abc4
= + = + .
6 4 z=0 6 4
Note that the region of integration in the above Example is a volume, but it is not a volume which is being
calculated. One interpretation of the answer is the mass of a rectangular solid box when the density of the
solid is xy 2 + z 3 .
Z 3 Z 2 Z 1
***Do Exercise G Evaluate x dx dy dz.
0 0 0
7. This module has concentrated on the techniques of evaluating double and triple integrals. Before finishing
we list some of the many geometrical and physical applications of these integrals, a few of which have been
discussed above.
(b) The volume V beneath the surface z = f (x, y) (> 0) and above a region R in the xy-plane is
ZZ
V = f (x, y) dx dy .
R
(c) Let f (x, y) be the density (= mass per unit area) of a thin lamina of area R in the xy-plane. Then the
total mass M of the lamina is ZZ
M= f (x, y) dx dy ,
R
the centre of mass of the lamina has coordinates (x, y), where
ZZ ZZ
1 1
x= x f (x, y) dx dy, y= y f (x, y) dx dy,
M R M R
the moments of inertia Ix and Iy of the lamina about the x- and y-axes respectively are
ZZ ZZ
2
Ix = y f (x, y) dx dy, Iy = x2 f (x, y) dx dy,
R R
and the polar moment of inertia I0 of the lamina about the origin is
ZZ
I0 = Ix + Iy = (x2 + y 2 )f (x, y) dx dy.
R
– 12 –
(e) Let f (x, y, z) be the density (= mass per unit volume) of a solid body of volume V in three dimensional
space. Then the total mass M of the solid is
ZZZ
M= f (x, y, z) dx dy dz,
V
and the centre of mass of the solid has coordinates (x, y, z), where
ZZZ ZZZ ZZZ
1 1 1
x= x f (x, y, z) dx dy dz, y= y f (x, y, z) dx dy dz, z= z f (x, y, z) dx dy dz.
M V M V M V
Corresponding expressions for the moments of inertia of the solid could also be written down.
– 13 –
Specimen Test 11
...
......
y
..........
..
...
....
..
......
1 ....
... ..........
... 2.............................
x=y
..
.... ..
..
........ ......
.
.
.. ..... .... ...
..... .......
... ........
... ..... .............
... ..... .
.
.
.
... ..... ........ ......
y=x ..
.
... ..... ......
... ........ .....
...... ........
...... .....
.. . .
.................................................................................................................................................................
.
0 .... 1 x
Z x=1 Z y=1−x2
3. A double integral is defined by I= xy dy dx .
x=0 y=0
4. (i) State the expression for an element of area in terms of plane polar coordinates (r, θ).
(ii) Evaluate the double integral
Z θ=π Z r=a
sin θ r dr dθ.
θ=0 r=0
– 14 –