General Mathematics II: Dr. Mohamad Alghamdi
General Mathematics II: Dr. Mohamad Alghamdi
Contents
List of Abbreviations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1 CONIC SECTIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1 Parabola . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.1 Vertical Parabolas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.2 Horizontal Parabolas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.2 Ellipse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.2.1 Ellipses with the Center at the Origin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.2.2 Ellipses with the Center Not at the Origin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3 Hyperbola . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.3.1 Hyperbola with the Center at the Origin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.3.2 Hyperbola with the Center Not at the Origin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
4 INTEGRATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
4.1 Antiderivatives and Indefinite Integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
4.1.1 Antiderivatives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2 Contents
5 APPLICATIONS OF INTEGRATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
5.1 Areas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
5.1.1 Region Bounded by a Curve and x-axis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
5.1.2 Region Bounded by a Curve and y-axis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
5.1.3 Region Bounded by Two Curves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
Appendix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
Appendix (1): Integration Rules and Integrals Table . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
List of Abbreviations
|x| absolute value the absolute value of x is always either positive or zero
a∼
=b approximate a is approximately equal to b
π constant π ≈ 3.141592654
Chapter 1
CONIC SECTIONS
1.1 Parabola
Definition 1.1 A parabola is a set of all points in a plane that are equidistant from a fixed point F (called the focus) and a
fixed line D (called the directrix) in the same plane.
Let the focus lies along the x-axis at F = (a, 0) and let the directrix
be the line x = −a. From Definition 1.1, we have |PF| = |PD|. Then
from the distance formula, we obtain
q q
(x − a)2 + (y − 0)2 = (x + a)2 + (y − y)2
q
⇒ (x − a)2 + y2 = (x + a)
⇒ (x − a)2 + y2 = (x + a)2
⇒ y2 = (x + a)2 − (x − a)2
⇒ y2 = 4ax .
The result is the equation of a parabola with vertex at the origin, that opens to the right. Similarly, we can extract the other equations of
the parabola. In each case, a > 0 which represents the distance from the vertex to the focus. The axis of symmetry of the parabola is a
line that passes through the vertex and is perpendicular to the directrix.
When a parabola opens right or left, it has a vertical axis of symmetry. In this case, the parabola is called a vertical parabola. We study
the special and general cases of the vertical parabolas. In the special case, we assume that the vertex of the parabola is at the origin. In
the general case, we assume that the vertex is at V (h, k).
The equation of the vertical parabola with the vertex at the origin is x2 = ±4ay, where a > 0.
8 CONIC SECTIONS
The equation of the vertical parabola with the vertex at V (h, k) is (x − h)2 = ±4a(y − k), where a > 0. The previous form is the general
formula of the vertical parabolas.
Example 1.1 Find the focus and the directrix of the parabola x2 = 4y, and sketch its graph.
Solution:
The equation x2 = 4y takes the form x2 = 4ay with a = 1.
Example 1.2 Find the focus and the directrix of the parabola (x + 1)2 = −4(y − 1), and sketch its graph.
(x − h)2 = −4a(y − k) .
• The directrix of the parabola is y = 2. Figure 1.7: The graph of the parabola (x + 1)2 = −4(y − 1).
Example 1.3 Find the equation of the parabola with vertex (2, 1) and focus F(2, 3). Then, sketch the graph.
10 CONIC SECTIONS
Solution:
Since the vertex and focus are in the same line x = 2, then the axis
of symmetry of the parabola is parallel to the y-axis. Also, from the
y-coordinate of the vertex and focus, the parabola opens upwards.
Thus, the equation of the parabola takes the form
(x − h)2 = 4a(y − k) .
F(h, k + a) = (2, 3) ⇒ a = 2
By substituting the values of a, h and k, the equation of the parabola
becomes (x − 2)2 = 8(y − 1).
Figure 1.8: The graph of the parabola (x − 2)2 = 8(y − 1).
When a parabola opens upwards or downwards, it has a horizontal axis of symmetry. In this case, the parabola is called a horizontal
parabola. We consider the two cases: the vertex at the origin and the vertex at V (h, k).
(A) Horizontal Parabolas with the Vertex at the Origin. The equation of the horizontal parabola with the vertex at the origin is
y2 = ±4ax, where a > 0.
(B) Horizontal Parabolas with the Vertex at V (h, k). The general equation of the horizontal parabola with the vertex at V (h, k) is
(y − k)2 = ±4a(x − h), where a > 0.
Parabola 11
Example 1.4 Find the focus and the directrix of the parabola y2 = −8x, and sketch its graph.
Solution:
Example 1.5 Find the equation of the parabola with focus (3, 0) and directrix x = −3. Then, sketch the graph.
Solution:
Since the focus of the parabola is F(3, 0) = F(a, 0), the axis of
symmetry of the parabola is x-axis and it opens to the right.
Example 1.6 Find the focus and the directrix of the parabola 2y2 − 4y + 8x + 10 = 0, and sketch its graph.
Solution:
Since the quadrature is on the y−term, then the parabola takes the form (y − k)2 = ±4a(x − h).
1.2 Ellipse
Definition 1.2 An ellipse is a set of all points in a plane such that the sum of the distances from each point to two fixed points
(called foci) is constant.
Ellipse 13
Let c2 be a circle with midpoint F2 and radius 2a. From Figure 1.17, the distance of the point P to the circle c2 equals the distance to the
focus F1 . Therefore, if the point P = W1 (0, b), then |PF1 | = |Pc2 | = a. From Pythagoras’ theorem, we have a2 = b2 + c2 .
Figure 1.17
Figure 1.18
2
x2
The equation of the ellipse is a2
+ by2 = 1.
14 CONIC SECTIONS
• The foci of the ellipse are F1 (c, 0), F2 (−c, 0), where
p
c = a2 − b2 .
• The minor axis endpoints (co-vertices) are W1 (0, b), W2 (0, −b).
• The foci of the ellipse are F1 (0, c), F2 (0, −c), where
p
c = b2 − a2 .
• The minor axis endpoints (co-vertices) are W1 (a, 0), W2 (−a, 0).
• The minor axis of the ellipse is y-axis with length 2a. Figure 1.20: The graph of the ellipse x2
+ y2
= 1, where
a2 b2
b > a.
Example 1.7 Identify the features of the ellipse and sketch its graph.
(1) 9x2 + 25y2 = 225 (2) 16x2 + 9y2 = 144
Solution:
2 2
x2 x2
1. By dividing both sides by 225, we have 25 + y9 = 1. The result takes the form a2
+ by2 = 1, where a = 5 and b = 3. Since a > b,
√ √
then c = 25 − 9 = 16 = 4.
• The minor axis endpoints (co-vertices) are W1 (0, 3), W2 (0, −3).
• The minor axis endpoints (co-vertices) are W1 (3, 0), W2 (−3, 0).
Example 1.8 Find an equation of an ellipse if the center is at the origin and
(1) Major axis on x-axis (2) Major axis on y-axis
Major axis length = 14 Minor axis length = 14 √
Minor axis length = 10 Distance of foci from center = 10 2
Solution:
2 2
(1) Since the major axis is on x-axis, then the equation of the ellipse takes the form ax2 + by2 = 1, where a > b. Also, the major axis
length is 2a = 14 and this implies a = 7. The minor axis length is 2b = 10, so b = 5. From this, the equation of the ellipse
x2 y2
becomes 49 + 25 = 1.
2 2
(2) Since the major axis is on y-axis, then the equation of the ellipse takes the form ax2 + by2 = 1, where b > a. From the minor axis
√
length 2a = 14, we have a = 7. Also, the distance of foci from the center is c = 10 2. Since c2 = b2 − a2 , then b2 = 249. By
x2 y2
substituting the values of a and b, we have 49 + 249 = 1.
(x−h)2 (y−k)2
The equation of an ellipse of the form is a2
+ b2
= 1.
• The minor axis of the ellipse is y-axis with length 2b. (x−h)2 (y−k)2
Figure 1.23: The graph of the ellipse a2
+ b2
= 1,
where a > b.
16 CONIC SECTIONS
• The foci of the ellipse are F1 (h, k + c), F2 (h, k − c), where
p
c = b2 − a2 .
• The minor axis of the ellipse is x-axis with length 2a. (x−h)2 (y−k)2
Figure 1.24: The graph of the ellipse a2
+ b2
= 1,
where b > a.
Example 1.9 Find the equation of the ellipse with foci at (−3, 1), (5, 1) and one of its vertice is (7, 1), then sketch its graph.
Solution:
(x−h)2 (y−k)2
Since the y−term in the foci is constant, the equation of the ellipse is of the form a2
+ b2
= 1 where a > b.
(x − 1)2 (y − 1)2
+ =1.
36 20
Example 1.10 Find the equation of the ellipse with foci at (2, 5), (2, −3) and the length of its minor axis equals 6, then and sketch
its graph.
Ellipse 17
Solution:
(x−h)2 (y−k)2
Since the x−term in the foci is constant, the equation of the ellipse is of the form a2
+ b2
= 1 where b > a.
Also, the length of its minor axis equals 2a = 6, hence a = 3. From the formula c2 = b2 − a2 , we have b2 = 16 + 9 = 25, so b = 5. Thus,
the equation of the ellipse is
(x − 2)2 (y − 1)2
+ =1.
9 25
Example 1.11 Identify the features of the ellipse 4x2 + 2y2 − 8x − 8y − 20 = 0, then sketch its graph.
Solution:
4x2 + 2y2 − 8x − 8y − 20 = 0
2x2 + y2 − 4x − 4y − 10 = 0
2x2 − 4x + y2 − 4y = 10 isolate x any y terms
2 2
2(x − 2x + 1) + (y − 4x + 4) = 10 + 2 + 4 completing square: (a ± b)2 = a2 ± 2ab + b2
2(x − 1)2 + (y − 2)2 = 16
(x − 1)2 (y − 2)2
+ =1 divide by 16 .
8 16
where
√ √ √
h = 1, k = 2, a = 2 2, and b = 4, then c = 16 − 8 = 2 2 .
18 CONIC SECTIONS
1.3 Hyperbola
Definition 1.3 A hyperbola is the set of all points in a plane such that the absolute value of the difference of the distances of
each point from two fixed points (called foci) is constant.
• The point midway between the foci is called the center. The
line containing the foci is the transverse axis.
Let c2 be the circle with midpoint F2 and radius 2a. The distance of a point P of the right branch to the circle c2 equals the distance to
the focus F1 : |PF1 | = |Pc2 |.
Hyperbola 19
Figure 1.29
• The foci of the ellipse are F1 (c, 0), F2 (−c, 0), where
p
c = a2 + b2 .
y2 2
(B) The equation of the hyperbola is a2
− bx2 = 1.
• The foci of the ellipse are F1 (0, c), F2 (0, −c), where
p
c = a2 + b2 .
Example 1.12 Identify the features of the hyperbola and sketch its graph.
Solution:
2 2 2
x2
x
(1) By dividing both sides by 64, we have 16 − y4 = 1. The result takes the form a2
− by2 = 1.
√ √
Since a = 4 and b = 2, then c = 16 + 4 = 2 5.
2 2 y2 2
(2) Divide both sides by 36 to have y9 − x4 = 1. The result takes the form b2
− ax2 = 1.
√ √
Since a = 2 and b = 3, then c = 4 + 9 = 13.
y2 2
Figure 1.33: The graph of the hyperbola 9 − x4 = 1.
Example 1.13 Find an equation of the hyperbola if its vertices are V1 (3, 0) and V2 (−3, 0), and one of its foci (4, 0), then sketch its
graph.
Solution:
2
x2
Since the y−term in the vertices is constant, the equation of the hyperbola takes the form a2
− by2 = 1. Also, V1 (a, 0) = V1 (3, 0) implies
a = 3 and F1 (c, 0) = F1 (4, 0) implies c = 4.
Hyperbola 21
√ √
From the formula c2 = a2 + b2 , we have b = 16 − 9 = 7.
2 2
Thus, the equation of the hyperbola is x9 − y7 = 1.
The hyperbola has the following properties:
• The center of the ellipse is P(0, 0)
(x−h)2 (y−k)2
(A) The equation of the hyperbola is a2
− b2
= 1.
(y−k)2 (x−h)2
(B) The equation of the hyperbola is b2
− a2
= 1.
(y−k)2 2
Figure 1.36: The graph of the hyperbola b2
− (x−h)
a2
= 1.
22 CONIC SECTIONS
Example 1.14 Find the equation of the hyperbola with foci at (−2, 2), (6, 2) and one of its vertices is (5, 2), then sketch its graph.
Solution:
Since the y−term in the foci is constant, then the equation of the
(x−h)2 (y−k)2
hyperbola takes the form a2 + b2 = 1. Illustration:
From the given foci, we have
F1 (h + c, k) = (6, 2) ⇒ h + c = 6, k = 2 h+c = 6 → 1
F2 (h − c, k) = (−2, 2) ⇒ h − c = −2, k = 2 h − c = −2 → 2
By doing some calculation, we obtain h = 2 and c = 4.
Also, from the given vertex V1 (h + a, k) = (5, 2), we have h + a = 5. By substituting the value of h, we obtain a = 3. From the formula
√ (x−2)2 (y−2)2
c2 = a2 + b2 , we find b2 = 16 − 9 = 7 and this implies b = 7. Thus, the equation of the hyperbola is 9 + 7 = 1.
(x−2)2 2
Figure 1.37: The graph of the hyperbola 9 + (y−2)
7 = 1.
Example 1.15 Find the equation of the hyperbola with foci at (−1, −6) , (−1, 4) and the length of its transverse axis is 8, and sketch
its graph.
Solution:
(y−k)2 (x−h)2
Since the x−term in the foci is constant, the equation of the hyperbola takes the form b2
− a2
= 1.
Also, the length of the transverse axis is 2b = 8 and this implies b = 4. From the formula c2 = a2 + b2 , we have a2 = 25 − 16 = 9, so
a = 3.
(y+1)2 2
Figure 1.38: The graph of the hyperbola 16 + (x+1)
9 = 1.
Example 1.16 Identify the features of the hyperbola 2y2 − 4x2 − 4y − 8x − 34 = 0. Then, sketch its graph.
Solution:
2y2 − 4x2 − 4y − 8x − 34 = 0,
2y2 − 4y − 4x2 − 8x = 34
2(y2 − 2y) − 4(x2 − 2x) = 34 Rearranging x-terms and y-terms
2(y2 − 2y + 1) − 4(x2 − 2x + 1) = 34 + 2 + 4 completing the square
2 2
2(y − 1) − 4(x + 1) = 40
(y − 1)2 (x + 1)2
− =1 dividing both sides by 40
20 10
.
(y−k)2 (x−h)2 √ √
From the standard form b2
− a2
= 1, we have h = −1, k = 1, a = 10, and b = 2 5. Thus, from the formula c2 = a2 + b2 , we
√
have c = 30.
Exercises
1 - 13 Write an equation of the parabola with the given elements, then sketch the graph.
2 Vertex at (2, 6), passes through (−1, 4) and opens to the left.
5 Vertex at (0, −2), passes through (−2, 0) and opens to the left.
6 Vertex at (5, 3), passes through (3, −1) and opens downwards.
8 Vertex at (4, −3), passes through (10, −6) and opens downwards.
9 Vertex at (−3, −5), passes through (−5, −6) and opens downwards.
12 Vertex at (−8, 2), passes through (2, −3) and opens downwards.
14 - 21 Write an equation of the ellipse with the given elements, then sketch the graph.
14 Center at the origin and major axis on x-axis and its length equals 8 and minor axis length equals 6.
15 One of its vertices (3, 0), one of its foci at (2, 0).
√
16 Center at (2, 2), one of its vertices (4, 2), and one of its foci (2 + 3, 2).
√
17 Center at (1, −1), one of its vertices (4, −1), and one of its foci (1 + 5, −1).
18 Center at (−2, 3), major axis is parallel to y-axis, and its length equals 8 and minor axis length equals 4.
19 Vertices (2, 3) and (2, −2), and minor axis is parallel to x-axis and its length equals 2.
20 Vertices (−1, −1) and (−1, 9), and minor axis is parallel to x-axis with length 8.
21 Foci (10, −2), (4, −2) and one of its vertices (12, −2).
22 - 31 Write an equation of the hyperbola with the given elements, then sketch the graph.
√
22 Vertices (0, −2) and (0, 2) and one of its foci (0, 13).
23 Vertices (0, −6) and (0, 6), and one of its foci (0, −8).
24 Vertices (1, 1) and (11, 1), and one of its foci (12, 1).
Exercises 25
31 Foci (4, −2) and (10, −2), and one of its vertices (8, −2).
32 - 75 Determine the elements of the conic section and sketch its graph.
33 y = −(x − 2)2 − 2 55 y = x2 + 6x + 5
58 x − 5 = (y − 3)2
36 y = −4(x − 1)2 + 1
59 x = −2y2 − 4y − 5
37 y = 4x2 + 24x + 25
60 x2 + 5y2 + 6x − 40y + 84 = 0
38 y = 4(x − 5)2 − 7
61 x2 + 2y + 2x = 2
39 y = −5(x + 4)2 + 9
2
x2
40 y = x2 − 8x + 7 62 25 − y9 = 1
2
x2
41 4x2 + 10y2 = 100 63 16 − y9 = 1
42 x2 + 9y2 = 36 64 y2 x
− 25
2
=1
49
2
x2 y
43 100 + 49 =1 x2 y 2
65 4 − 49 =1
2
x2
44 5 + y7 = 1 x2 y 2
66 25 − 81 =1
2
x2 y
45 49 + 36 =1 y2 x 2
67 64 − 25 =1
(x+3)2 (y−2)2
46 + =1 x2 y 2
16 9 68 36 − 20 =1
2
x2 y
47 36 + 81 =1 (x+3)2 (y−2)2
69 16 − 9 =1
x2 y2
48 15 + 30 =1 (x−2)2 y 2
70 25 − 16 =1
2
x2 y
49 55 + 27 =1 (y−5)2 (x−6)2
71 64 − 25 =1
x2 y2
50 64 + 10 =1 (x+4)2 (y−5)2
72 81 − 55 =1
51 25(x − 3)2 + 10(y + 2)2 = 100
73 −4x2 + 10y2 = 100
2 2
x (y−2)
52 9 + 25 =1 74 10y2 + 49x2 = 490
53 49x2 + 4y2 = 196 75 y2 − 5x2 + 6y − 40x − 76 = 0
26 Matrices and Determinants
Chapter 2
2.1 Matrices
Definition 2.1 A matrix A of order m × n is a set of numbers or expressions arranged in a rectangular array of m rows and n
columns.
···
a11 a12 a13 a1,n−1 a1n
a21 a22 a23 ··· a2,n−1 a2n
a31 a32 a33 ··· a3,n−1 a3n
A= .
.. .. .. .. .. ..
. . . . . .
a
m−1,1 am−1,2 am−1,3 ··· am−1,n−1 am−1,n
am1 am2 am3 ··· am,n−1 amn
Notes :
(1) The horizontal arrays of a matrix are called its rows and the vertical arrays are called its columns.
(2) ai j represents the element of the matrix A that lies in row i and column j.
Example 2.1 Find the order of each matrix, then find the given elements.
2 −4
(1) A = , a11 and a22
1 0
1 3 5
(2) B = , a12 , a21 and a23
2 1 0
Solution:
Definition 2.2 Two matrices A = [ai j ] and B = [bi j ] having the same order m × n are equal if ai j = bi j for each i = 1, 2, ..., m
and j = 1, 2, ..., n.
Solution:
Since the matrices A = B, then from Definition 2.2, we have 4x − 1 = 11. This implies x = 3.
1. Row vector. A row vector of order n is a matrix of order 1 × n written as A = [a1 a2 ... an ]. For example, A = [2 7 0 − 1 9] is
a row vector of order 5.
a1
a2 1
2. Column vector. A column vector of order n is a matrix of order n × 1 written as A = . . For example, A = 7 is a column
..
3
an
vector of order 3.
3. Null matrix. The matrix A = [ai j ]m×n is called a null matrix if ai j = 0 for all i and j i.e.
···
0 0 0
0 0 ··· 0
A = .
.. .. ..
.. . . .
0 0 ··· 0
0 0 0
For example, is a null matrix of order 2 × 3.
0 0 0
4. Square matrix. If the number of rows equals the number of columns (m = n), then the matrix is called a square matrix of order
n. In a square matrix A = [ai j ], the set of elements of the form aii is called the diagonal of the matrix. For example, the diagonal
2 −7 3
of the following square matrix is highlighted in red 1 0 9. = The main diagonal of A
−1 6 8
5. Upper triangular matrix. The square matrix A = [ai j ] of order n is called an upper triangular matrix if ai j = 0 for all i > j:
···
a11 a12 a13 a1n
0
a22 a23 ··· a2n
A= 0
0 a33 ··· a3n
. .. .. ..
..
. . .
0 0 0 ··· ann
2 3 1
For example, 0 −1 4 is an upper triangular matrix of order 3.
0 0 5
6. Lower triangular matrix. The square matrix A = [ai j ] of order n is called a lower triangular matrix if ai j = 0 for all i < j:
···
a11 0 0 0
a21
a22 0 ··· 0
A = a31
a32 a33 ··· 0 .
. .. .. ..
..
. . .
an1 an2 an3 ··· ann
28 Matrices and Determinants
4 0 0
For example, 1 −1 0 is a lower triangular matrix of order 3.
2 3 5
7. Diagonal matrix. The square matrix A = [ai j ] of order n is called a diagonal matrix if ai j = 0 for all i j:
···
a11 0 0 0
0
a22 0 ··· 0
A= 0
0 a33 ··· 0
. .. .. ..
..
. . .
0 0 0 ··· ann
1 0 0 ··· 0
0 1 0 · · · 0
In = 0 0 1 · · · 0
. . . .
.. .. .. ..
0 0 0 ··· 1
1 0 0
For example, 0 1 0 is an identity matrix of order 3.
0 0 1
Definition 2.3 Let A = [ai j ] and B = [bi j ] be two matrices of order m × n. Then,
1. A + B = C with ci j = ai j + bi j .
2. A − B = C with ci j = ai j − bi j .
From Definition 2.3, if A = [ai j ] and B = [bi j ] are two matrices of order m × n, then
···
a11 + b11 a12 + b12 a1n + b1n
a21 + b21 a22 + b22 ··· a2n + b2n
A+B = .
.. .. .. ..
. . . .
am1 + bm1 am2 + bm2 ··· amn + bmn
Also,
a11 − b11 a12 − b12 ··· a1n − b1n
a21 − b21 a22 − b22 ··· a2n − b2n
A−B = .
.. .. .. ..
. . . .
am1 − bm1 am2 − bm2 ··· amn − bmn
Matrices 29
1 3 2 5 0 8
Example 2.3 If A = 5 −4 6 and B = 1 4 −1, find A + B and A − B.
0 9 2 10 11 −2
Solution:
1+5 3+0 2+8 6 3 10
A+B = 5 + 1 −4 + 4 6 + (−1) = 6 0 5 .
0 + 10 9 + 11 2 + (−2) 10 20 0
1−5 3−0 2−8 −4 3 −6
A-B = 5 − 1 −4 − 4 6 − (−1) = 4 −8 7 .
0 − 10 9 − 11 2 − (−2) −10 −2 4
(2) Multiplying a matrix by a scalar:
Definition 2.4 Let A = [ai j ] be a matrix of order m × n. Then, for any k ∈ R, kA is a matrix C = [ci j ] with ci j = kai j .
···
ka11 ka12 ka1n
ka21 ka22 ··· ka2n
kA = . .. .
.. ..
.. . . .
kam1 kam2 ··· kamn
1 3 2
Example 2.4 If A = , find 3A.
0 9 2
Solution:
3×1 3×3 3×2 3 9 6
3A = =
3×0 3×9 3×2 0 27 6
1 6 2 3
Example 2.5 If A = and B = , find −2A + 3B.
−2 4 0 8
Solution:
1 6 2 3 −2 −12 6 9 4 −3
−2A + 3B = −2 +3 = + = .
−2 4 0 8 4 −8 0 24 4 16
4. (k + `)A = kA + `A.
5. k(`A) = (k`)A.
30 Matrices and Determinants
Definition 2.6 Let A = [ai j ] be a matrix of order m × n. There exists a matrix B such that A + B = 0. This matrix B is called
the additive inverse of A and it is denoted by −A = (−1)A.
Definition 2.7 Let A = [ai j ] be a matrix of order m × n and B = [bi j ] be a matrix of order n × p. The multiplication AB is a
matrix C = [ci j ] of order m × p, where
n
ci j = ∑ aik bk j = ai1 b1 j + ai2 b2 j + ... + ain bn j .
k=1
Note: the multiplication AB is defined if and only if the number of columns of A equals the number of rows of B; otherwise, we say the
multiplication is undefined.
1 6 2 3
Example 2.6 If A = and B = , find AB.
−2 4 0 8
Solution:
1 6 2 3 1 × 2 + 6 × 10 1×3+6×8 2 51
AB = = = .
−2 4 0 8 −2 × 2 + 4 × 0 −2 × 3 + 4 × 8 −4 26
2 3 0
1 6 2
Example 2.7 If A = and B = −1 4 2, find AB.
−2 4 1
0 1 7
Solution:
2 3 0
1 6 2 −4 29 26
AB = −1 4 2 = .
−2 4 1 −8 11 15
0 1 7
A special case of multiplication of matrices is multiplying a row vector by a column vector. Let A = [a1 a2 ... an ] be a row vector of
b1
b2
order n and B = . be a column vector of order n. Then the multiplication AB is a matrix C = [c] of order 1, where
..
bn
n
c= ∑ ak bk = a1 b1 + a2 b2 + ... + an bn .
k=1
8
Example 2.8 If A = 2 1 4 and B = −3, find AB
5
Solution:
Matrices 31
AB = [2 × 8 + 1 × (−3) + 4 × 5] = [33].
1
3
Example 2.9 If A = 2 1 4 −6 and B =
−4, find AB
Solution:
AB = [2 × 1 + 1 × 3 + 4 × (−4) + (−6) × 7] = [−53].
2 2 3
Example 2.10 If A = 1 3 5 , B = −1 and C = −1 4 , Compute (if possible) 1. AB 2. BC.
0 0 1
Solution:
2
(1) AB = A = 1 3 5 −1 = [−1].
0
(2) The multiplication BC is not possible since the matrix B of order 3 × 1 and the matrix C of order 3 × 2.
Theorem 2.8 Let A be a matrix of order m × n, B be a matrix of order n × p and C be a matrix of order p × q. Then,
1. The multiplication of matrices is not commutative: AB 6= BA.
Theorem 2.9 Let A and B be any two matrices of order m × n. The multiplication of matrices is distributive:
1. (A + B)C = AC + BC, where C is a matrix of order n × p.
4 3 9 2 3
Example 2.11 If A = and B = Compute (if possible) 1. AB 2. BA.
−1 2 0 0 1
Solution:
(1) The multiplication BC is not possible since the matrix B of order 3 × 1 and the matrix C of order 3 × 2.
2 3 4 3 9 5 12 13
(2) BA = = .
0 1 −1 2 0 −1 2 0
3 4 7 5
Example 2.12 If A = and B = , compute (if possible) 1. AB 2. BA.
−1 2 1 0
Solution:
32 Matrices and Determinants
3 4 7 5 25 15
(1) AB = = .
−1 2 1 0 −5 −5
7 5 3 4 16 38
(2) BA = = .
1 0 −1 2 3 4
From this example, we find AB 6= BA and this means the multiplication of matrices is not commutative.
4 3
Example 2.13 If A = , find AI2
−1 2
4 3 1 0 4 3
Solution: AIn = = .
−1 2 0 1 −1 2
Definition 2.10 Let A = [ai j ] be a matrix of order m × n. Then, the transpose of A is At = [a ji ]n×m .
2. (A + B)t = At + Bt .
3. (kA)t = kAt .
4. (AB)t = Bt At .
Remark 2.12
1. The transpose of a row vector is a column vector and vice-versa.
2. The transpose of a lower triangular matrix is an upper triangular matrix and vice-versa.
5 3
3 −1 0 2 1 −2
Example 2.14 If A = ,B= and C = 1 4 , Compute
2 5 1 −1 0 1
−1 2
(1) (At )t (3) (3A)t
Solution:
3 2
3 −1 0
(1) (At ) = −1 5, so (At )t = .
2 5 1
0 1
5 1
3 −1 0 2 1 −2 5 0 −2
(2) (A + B) = + = . From this, (A + B)t = 0 5.
2 5 1 −1 0 1 1 5 2
−2 2
3 2 9 6
(3) (3A)t = 3At = 3 −1 5 = −3 15.
0 1 0 3
Matrices 33
3 2
5 1 −1 14 3
(4) (AC)t = Ct At = −1 5 = .
3 4 2 5 50
0 1
Definition 2.13 Let A = [ai j ] be a square matrix of order n. Then, the determinant of A can be defined as follows:
a : A = [a]
det(A) =
∑nj=1 (−1)i+ j ai j Ai j (i = 1, ..., n) : otherwise,
Example
2.15
Find the determinant of the matrix.
1 5 4 −1
(1) A = (2) B =
3 7 2 9
Solution:
1 5
1. det(A) = = 1 × 7 − 3 × 5 = 7 − 15 = −8.
3 7
4 −1
2. det(B) = = 4 × 9 − (−1) × 2 = 36 + 2 = 38.
2 9
Before starting evaluating the determinant of an n × n matrix, we first need to define the minor and cofactor of that matrix. The minor
Mi j is the determinant of the matrix obtained by eliminating the ith row and the jth column of A.
1 3 1
Example 2.16 If A = −2 −1 2, find the minors M11 , M12 and M13 .
2 4 5
Solution:
−1 2
M11 = = −1 × 5 − 2 × 4 = −13.
4 5
−2 2
M12 = = −2 × 5 − 2 × 2 = −14.
2 5
−2 −1
M13 = = −2 × 4 − (−1) × 2 = −6.
2 4
34 Matrices and Determinants
Ci j = (−1)i+ j Mi j .
Example 2.17 In the previous example, calculate the corresponding cofactors of the minors M11 , M12 and M13 .
Solution:
• Multiply each of the elements ai j of the row (or column) by its corresponding cofactor Ci j .
To calculate the determinant, choose the first row of A and multiply each of its elements by the corresponding cofactor:
Example
2.18 Findthe determinant of the matrix.
1 6 3 4 1 5
(1) A = 5 −1 4 (2) B = 2 1 −2
−2 9 7 1 8 7
Matrices 35
Solution:
1 6 3
(1) det(A) = 5 −1 4 = 1C11 + 6C12 + 3C13
−2 9 7
= 1(−1)(1+1) M11 + 6(−1)(1+2) M12 + 3(−1)(1+3) M13
−1 4 5 4 5 −1
=1 −6 +3
9 7 −2 7 −2 9
= 1(−7 − 36) − 6(35 + 8) + 3(45 − 2) = −42 − 258 + 129 = −171 .
4 1 5
(2) det(B) = 2 1 −2 = 4C11 + 1C12 + 5C13
1 8 7
= 4(−1)(1+1) M11 + 1(−1)(1+2) M12 + 5(−1)(1+3) M13
1 −2 2 −2 2 1
=4 −1 +5
8 7 1 7 1 8
= 4(7 + 16) − 1(14 + 2) + 5(16 − 1) = 92 − 16 + 75 = 151 .
Sarrus rule
Let A be a square matrix of order 4 as follows:
a11 a12 a13 a14
a21 a22 a23 a24
A=
.
a31 a32 a33 a34
a41 a42 a43 a44
Then
where
a22 a23 a24 a21 a23 a24
M11 = a32 a33 a34 , M12 = a31 a33 a34 ,
a42 a43 a44 a41 a43 a44
a21 a22 a24 a21 a22 a23
M13 = a31 a32 a34 , M14 = a31 a32 a33 .
a41 a42 a44 a41 a42 a43
Example
2.19 Find the determinant
of the matrix.
1 6 3 2 4 1 5 2
5 −1 4 1 2 1 0 1
(1) A =
−2
(2) B =
9 7 3 0 1 3 9
7 1 3 −6 1 7 4 6
36 Matrices and Determinants
Solution:
−1 4 1
7 3 9 3 9 7
(1) det(A1 ) = 9 7 3 = −1 −4 +1
3 −6 1 −6 1 3
1 3 −6
= −1(−42 − 9) − 4(−54 − 3) + 1(27 − 7) = 51 + 228 + 20 = 299 .
5 4 1
7 3 −2 3 −2 7
det(A2 ) = −2 7 3 =5 −4 +1
3 −6 7 −6 7 3
7 3 −6
= 5(−42 − 9) − 4(12 − 21) + 1(−6 − 49) = −255 + 36 − 55 = −274 .
5 −1 1
9 3 −2 3 −2 9
det(A3 ) = −2 9 3 =5 +1 +1
1 −6 7 −6 7 1
7 1 −6
= 5(−54 − 3) + 1(12 − 21) + 1(−2 − 63) = −285 − 9 − 65 = −359 .
5 −1 4
9 7 −2 7 −2 9
det(A4 ) = −2 9 7 =5 +1 +4
1 3 7 3 7 1
7 1 3
= 5(27 − 7) + 1(−6 − 49) + 4(−2 − 63) = 100 − 55 − 260 = −215 .
Thus,
det(A) = 1 det(A1 ) − 6 det(A2 ) + 3 det(A3 ) − 2 det(A4 ) = 1(299) − 6(−274) + 3(−359) − 2(−215) = 1296 .
1 0 1
3 9 1 9 1 3
(2) det(B1 ) = 1 3 9 =1 −0 +1
4 6 7 6 7 4
7 4 6
= 1(18 − 36) − 0 + 1(4 − 21) = −18 − 17 = −35 .
2 0 1
3 9 0 9 0 3
det(B2 ) = 0 3 9 =2 −0 +1
4 6 1 6 1 4
1 4 6
= 2(18 − 36) − 0 + 1(0 − 3) = −36 − 3 = −39 .
2 1 1
1 9 0 9 0 1
det(B3 ) = 0 1 9 =2 −1 +1
7 6 1 6 1 7
1 7 6
= 2(6 − 63) − 1(0 − 9) + 1(0 − 1) = −114 + 9 − 1 = −106 .
2 1 0
1 3 0 3 0 1
det(B4 ) = 0 1 3 =2 −1 +0
7 4 1 4 1 7
1 7 4
= 2(4 − 21) − 1(0 − 3) + 0 = −35 + 3 = −32 .
Thus,
det(B) = 4 det(B1 ) − 1 det(B2 ) + 5 det(B3 ) − 2 det(B4 ) = 4(−35) − 1(−39) + 5(−106) − 2(−32) = −773 .
Matrices 37
Theorem 2.14
1. If A is a square matrix having a zero row (or a zero column), then det(A) = 0.
2. If A is a square matrix having two equal rows (or two equal columns), then det(A) = 0.
3. If A is a square matrix having a row which is a multiple of another row (or a column which is a multiple of another
column), then det(A) = 0.
4. If A is a diagonal matrix or an upper triangular matrix or a lower triangular matrix, then det(A) is the product of the
elements of the main diagonal.
5. The determinant of the null matrix is 0 and the determinant of the identity matrix is 1.
7. If B is obtained from A by interchanging two rows (or two columns), then det(B) = −det(A).
8. If B is obtained from A by multiplying a row by a non-zero constant and adding the result to another row (or multiplying
a column by a non-zero constant and adding the result to another column), then det(B) = det(A).
Example
2.20 Find
the determinant the matrix.
1 2 −2 1 2 −2
(1) A = 0 0 0 (3) C = 4 7 5
3 4 7 3 6 −6
1 2 1 3 1 4
(2) B = 6 5 6 (4) D = 0 −1 2
3 4 3 0 0 5
Solution:
(1) The matrix A contains a zero row, so from (item 1) in Theorem 2.14, we have det(A) = 0.
(2) The matrix A contains two equal columns, so from (item 2) in Theorem 2.14, we have det(A) = 0.
(3) The third row in matrix A is a multiple of the first row by 2, so from (item 3) in Theorem 2.14, we have det(A) = 0.
(4) The matrix A is an upper triangular matrix, so from (item 4) in Theorem 2.14, we have det(A) = −15.
Example
2.21 Findthe determinant the matrix.
1 3 1 4 2 −1
(1) A = 4 2 −1 (3) C = 1 3 1
0 −3 2 0 −3 2
2 6 2 1 3 1
(2) B = 4 2 −1 (4) D = 6 8 2
0 −3 2 0 −3 2
Solution:
Theorem 2.15
1. Let A and B be square matrices of order n. Then det(AB) = det(A)det(B).
2. Let A be a square matrix. Then det(A) = det(At ).
1 3 −2 0
Example 2.22 If A = and B = , find (1) det(AB) (2) det(At ).
4 2 3 1
Solution:
1 3 −2 0
det(A) = = −10 and det(B) = = −2.
4 2 3 1
Exercises
1 3 2 5 0 −2 0
1 - 10 If A = 5 −4 6 , B = 1 4 and C = 0 7, compute the following (if possible):
0 9 2 10 11 5 3
1 B +C 6 BA
2 2B + 3C 7 At
3 C−B 8 (3A)t
4 A −C 9 det(A)
5 AB 10 det(2A)
4 −1 −2 1
1 2 1
11 - 20 If A = 1 5 ,B= 3 6 and C = , compute the following (if possible):
9 5 3
2 7 1 4
11 A + B 16 BC
12 5A 17 AB
13 −3A + 2B 18 Bt
14 A −C 19 (2A)t
15 AC 20 (At )t
21 7det(A) 23 det(At )
22 det(AB) 24 det((AB)t )
25 - 32 Find the determinant.
1 −2 2 3
25 29
2 7 1 4
3 1 3 2 −2 2
26 7 2 9 30 −3 10 1
1 7 4 5 1 1
4 −2 3 1 −2 3
27 3 7 2 31 4 0 1
6 9 5 2 7 0
4 −2 3 1 1 5 3 6
3 7 2 2 1 0 1 2
28 32
6 9 5 −1 2 7 0 1
1 0 5 1 2 1 0 1
40 Chapter 3. SYSTEMS OF LINEAR EQUATIONS
Chapter 3
Definition 3.1 A linear system of m equations in n unknowns x1 , x2 , ..., xn is a set of equations of the form
···
a11 a12 a1n x1 b1
a21 a22 ··· a2n x2 b2
A= . .. , X = .. , and B = .. .
.. ..
.. . . . . .
an1 an2 ··· ann xn bn
X is called the column vector of the variables (or column vector of the unknowns)
B is called the column vector of constants (or column vector of the resultants)
A special case of the linear system of equations is a system of two different variables x1 and x2 :
a11 x1 + a12 x2 = b1
a21 x1 + a22 x2 = b2
a11 a12 x b
The above system of linear equations can be written as AX = B where A = , X = 1 , and B = 1 .
a21 a22 x2 b2
Solution of Linear Equations Systems 41
A solution of the linear system Ax = B is a column vector Y with entries y1 , y2 , ..., yn such that the linear system (3.1) is satisfied if we
replace yi with xi i.e., AY = B holds where Y t = [y1 , y2 , ..., yn ]. Note that for the linear system of equations AX = 0, the column vector
X t = [0, 0, ..., 0] is always solution and it is called the trivial solution.
In this chapter, we present three methods to solve the system of linear equations (3.1): Cramer’s method, Gauss elimination method,
and Gauss-Jordan method.
Theorem 3.2 Let AX = B be a linear system with n equations in n variables. The system has a solution if det(A) 6= 0.
Theorem 3.3 Let AX = B be a linear system with n equations in n variables. If det(A) 6= 0, then the unique solution to this
system is
det(Ai )
xi = for everyi = 1, 2, ..., n ,
det(A)
where Ai is the matrix formed by replacing the ith column of A by the column vector of constants B.
The matrix A1 is formed by replacing the first column of A by the column vector of constants B:
b1 a12 · · · a1n
b2 a22 · · · a2n
A1 = . .. .
.. ..
.. . . .
bn an2 · · · ann
The matrix A2 is formed by replacing the second column of A by the column vector of constants B:
a11 b1 · · · a1n
a21 b2 · · · a2n
A2 = . .. .
.. ..
.. . . .
an1 bn · · · ann
By continuing doing so, the matrix An is formed by replacing the last column of A by the column vector of constants B:
a11 a12 · · · b1
a21 a22 · · · b2
An = . .. .
.. ..
.. . . .
an1 an2 · · · bn
(2) 2x + y + z = 3 (4) x + y + z = 12
4x + y − z = −2 x−y = 2
2x − 2y + z = 6 x−z = 4
42 Chapter 3. SYSTEMS OF LINEAR EQUATIONS
Solution:
2 3
(1) A = ⇒ det(A) = 5.
−1 1
7 3
A1 = ⇒ det(A1 ) = −5.
4 1
2 7
A2 = ⇒ det(A2 ) = 15.
−1 4
Hence,
det(A1 ) −5 det(A2 ) 15
x1 = = = −1 and x2 = = =3
det(A) 5 det(A) 5
−1
The column vector of variables is X = .
3
2 1 1
(2) A = 4 1 −1 ⇒ det(A) = −18.
2 −2 1
3 1 1
A1 = −2 1 −1 ⇒ det(A1 ) = −9.
6 −2 1
2 3 1
A2 = 4 −2 −1 ⇒ det(A2 ) = 18.
2 6 1
2 1 3
A3 = 4 1 −2 ⇒ det(A3 ) = −54.
2 −2 6
Hence,
det(A1 ) −9 1 det(A2 ) 18 det(A3 ) −54
x= = = , y= = = −1 and z = = =3.
det(A) −18 2 det(A) −18 det(A) −18
1
2
The column vector of variables is X = −1.
3
1 2
(3) A = ⇒ det(A) = −3.
2 1
1 2
A1 = ⇒ det(A1 ) = 3.
−1 1
1 1
A2 = ⇒ det(A2 ) = −3.
2 −1
From this,
det(A1 ) 3 det(A2 ) −3
x1 = = = −1 and x2 = = =1.
det(A) −3 det(A) −3
−1
The column vector of variables is X = .
1
1 1 1
(4) A = 1 −1 0 ⇒ det(A) = 3.
1 0 −1
12 1 1
A1 = 2 −1 0 ⇒ det(A1 ) = 18.
4 0 −1
1 12 1
A2 = 1 2 0 ⇒ det(A2 ) = 12.
1 4 −1
1 1 12
A3 = 1 −1 2 ⇒ det(A3 ) = 6.
1 0 4
Therefore,
det(A1 ) 18 det(A2 ) 12 det(A3 ) 6
x= = = 6, y = = = 4 and z = = =2.
det(A) 3 det(A) 3 det(A) 3
Gauss Elimination Method 43
6
The column vector of variables is X = 4.
2
Definition 3.5 Two linear systems are said to be equivalent if one can be obtained from the other by a finite number of
elementary operations.
The elementary operations help us in getting a linear system from the system (3.1), which is easily solvable.
Theorem 3.6 Let CX = D be the linear system obtained from the linear system AX = B by a finite number of elementary
operations. Then, the linear systems AX = B and CX = D have the same set of solutions.
Elementary Operations
Multiply ist equation by λ and add the result Multiply ith row (Ri ) by λ and add the result
λRi + R j
to jth equation to jth row (R j ): −−−−−→
Replace ith equation by jth equation Replace ith row (Ri ) by jth row (R j ):
Ri ↔ R j
Table 3.1
x + y = 11 → 1
2x + y = 25 → 2
44 Chapter 3. SYSTEMS OF LINEAR EQUATIONS
Multiply the 1st row by −2 and add the result to 2nd row
−2x − 2y = −22
2x + y = 25
Definition 3.7 Two matrices are said to be row-equivalent if one can be obtained from the other by a finite number of
elementary row operations.
Example 3.2 The three matrices given below are row equivalent.
2 −3 1 −1 5 2 −2 10 4
−1 R 1 ↔R2 2R 1
5 2 −−−−−→ 2 −3 1 −−→ 2 −3 1 .
1 −2 −7 1 −2 −7 1 −2 −7
Definition 3.8 Gaussian elimination is a method of solving a linear system AX = B by constructing the augmented matrix
[A|B] and transforming the matrix A to an upper triangular matrix [C|D].
The Method:
2. Use the elementary row operations on the augmented matrix to transform the matrix A to an upper triangular matrix with a leading
coefficient of each row equals 1:
1 ···
c12 c13 c14 c1n d1
0 1 c33 c24 ··· c2n d2
.. .. .. .. .. ..
.. .
. . . . . . .
0 0 0 ··· 1 c(n−1)n dn−1
0 0 0 ··· 1 1 dn
3. From the last augmented matrix, we have xn = dn and the rest of the unknowns can be calculated by backward substitution.
Solution: For each system, construct the augmented matrix [A|B]. Then, use elementary row operations on the augmented matrix to
transform the matrix A to an upper triangular matrix with leading coefficient of each row equals 1.
1
3 1 9 R ↔R2 1 2 8 -3 R1 + R2 1 2 8 − R1 1 2 8
(1) [A|B] = −−1−−−→ −−−− −−→ −−5−→ .
1 2 8 3 1 9 0 −5 −13 0 1 13 5
Thus, x2 = 13 and x1 + 2x2 = 8. By substituting the value of x2 , we obtain x1 = 14
5 . Therefore, the column vector of variables is
14 5
X = 13 5 .
5
1 −2 1 4 1 −2 1 4 1 −2 1 4 1
R2
1R + R R2 ↔ R3
(2) [A|B] = −1 2 1 1
−2 −−−−−−→ 0 2
0 2 2 −−−−−→ 0 5 −5 −20 −5−→
−4R1 + R3 1
4 −3 −1 −4 0 5 −5 −20 0 0 2 2 2 R3
1 −2 1 4
0 1 −1 −4 .
0 0 1 1
y − z = −4 and x − 2y + z = 4. By doing some substitution, we obtain y = −3 and x = −3. The column vector of variables
z = 1,
Hence,
−3
is X = −3 .
1
1 1 1 2 1 1 1 2 1 1 1 2
−1R 1 + R 2 −1R2 + R3
(3) [A|B] = 1 −1 2 0 −−−− −−→ 0 −2 1 −2 −−−−−−→ 0 −2 1 −2
−2R1 + R3
2 0 1 2 0 −2 −1 −2 0 0 −2 0
1 1 1 2
− 1 R2
−−2−→ 0 1 − 12 1 .
− 12 R3
0 0 1 0
Thus, z = 0, y − 12 z = 1 and x + y + z = 2. By substituting the value of z and then y, we have y = 1 and x = 1. The column vector of
1
variables is X = 1 .
0
1 2 3 14 1 2 3 14 − 12 R3
1 2 3 14
−2R 1 + R 2
(4) [A|B] = 2 1 2 10 −−−−−−→ 0 −3 −4 −18 −−−→ 0 −3 −4 −18
−3R1 + R3
3 4 −3 2 0 −2 −12 −40 0 1 6 20
1 2 3 14 1 2 3 14 1
R
1 2 3 14
R ↔ R3 3R2 + R3 3
−−2−−−→ 0 1 6 20 −−−−−→ 0 1 6 20 −− 14
−→ 0 1 6 20 .
0 −3 −4 −18 0 0 14 42 0 0 1 3
Definition 3.9 Gauss-Jordan elimination is a method of solving a linear system AX = B by constructing the augmented matrix
[A|B] and transforming the matrix A to an identity matrix [In |D].
Exercises 46
The Method:
2. Use the elementary row operations on the augmented matrix [A|B] to transform the matrix A to the identity matrix In .
···
1 0 0 0 0 d1
0 1 0 0 ··· 0 d2
.. .. .. .. ..
.
. . . . .
0 0 0 ··· 1 0 dn−1
0 0 0 0 ··· 1 dn
(1) x + y = 2 (3) x − 2y + 2z = 5
2x + y = 1 5x + 3y + 6z = 57
x + 2y + 2z = 21
(2) x + y + z = 2
x − y + 2z = 0
2x + z = 2
Solution: For each linear system, construct the augmented matrix [A|B]. Then, use the elementary row operations on the augmented
matrix to transform the matrix A to the identity matrix In .
1 1 2 −2R + R
1 2 1 1 2 1R + R
2 1 0
1 −1 −1R2 1 0 −1
(1) [A|B] = −−−− −−→ −−− −−→ −−−→ .
2 1 1 0 −1 −3 0 −1 −3 0 1 3
−1
Thus, x = −1 and y = 3. The column vector of variables is X = .
3
1 1 1 2 1 1 1 2 1 1 1 2 − 1 R2
−1R 1 + R 2 −1R2 + R3
(2) [A|B] = 1 −1 2 0 −−−−−−→ 0 −2 1 −2 −−−−−−→ 0 −2 1 −2 −−2−→
−2R1 + R3 − 12 R3
2 0 1 2 0 −2 −1 −2 0 0 −2 0
1 1 1 2 1 1 1 1 2 1 0 1 1 1 0 0 1
R 3 + R 2 −1R + R −1R + R
0 1 − 1 1 −2−−−−→ 0 1 0 1 −−−− 2
−−→ 1
0 1 0 1 −−−− 3
−−→ 1
0 1 0 1 .
2
0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0
1
Hence, x = 1, y = 1 and z = 0. The column vector of variables is X = 1 .
0
1 −2 2 5 1 −2 2 5 1
R
1 −2 2 5
−5R1 + R2 3 R ↔ R3
(3) [A|B] = 5 3 6 57 −−−− −−→ 0 13 −4 32 −4−→ 0 13 −4 32 −−2−−−→
−1R1 + R3
1 2 2 21 0 4 0 16 0 1 0 4
1 −2 2 5 1 −2 2 5 1
− 4 R3
1 −2 2 5 1 0 2 13
−13R2 + R3 2R 2 + R 1 −2R3 + R1
0 1 0 4 −−−−−−−→ 0 1 0 4 −−−→ 0 1 0 4 −−−−−→ 0 1 0 4 −−−− −−→
0 13 −4 32 0 0 −4 −20 0 0 1 5 0 0 1 5
1 0 0 3
0 1 0 4 .
0 0 1 5
3
Hence, x = 3, y = 4 and z = 5. The column vector of variables is X = 4 .
5
Exercises 47
Exercises
1 x + y + z = 18 5 x+y−z = 1
x−y+z = 6 x − y + 2z = 2
x+y−z = 4 4x + y − z = 2
2 x1 + x2 − 2x3 = 1 6 x1 + 3x2 − x3 = 2
2x1 − 3x2 + x3 = −8 x1 − x2 + 5x3 = 3
3x1 + x2 + 4x3 = 7 3x1 + x2 − x3 = 1
3 x + 2y − z = −1 7 x + y − z = −3
x − 3y + 2z = 3 x − 6y + 5z = 1
x+y−z = 2 x + 4y + z = 1
4 2x1 + 2x2 − x3 = 1 8 2x − 4y + 3z = 10
2x1 − x2 + 2x3 = 1 3x + y − 2z = 6
x1 + x2 − 4x3 = 5 x + 3y − z = 20
9 x + y + z = 18 13 x1 + 4x2 − x3 = 1
x−y+z = 6 x1 − 2x2 + 7x3 = 9
x+y−z = 4 x1 + 2x2 + x3 = 3
10 x1 + 3x2 − x3 = 2 14 x + 2y + z = 11
x1 − x2 + 5x3 = 3 x + y − z = 20
3x1 + x2 − x3 = 1 x−y+z = 5
11 x − y − z = 1 15 x1 + 2x2 + x3 = 15
x − 6y + 5z = 4 x1 + 3x2 + x3 = 5
2x + y + z = 6 − 3x1 − x2 + 2x3 = 1
12 x1 − x2 − x3 = 4 16 x + y + z = 12
− 2x1 − x2 + x3 = 2 x−y = 2
x1 + x2 + 3x3 = 7 x−z = 4
17 x − 3y + z = 21 21 2x + y + 3z = 5
4x + 2y + z = 14 − 5x − 3y + z = 13
3x + 3y + z = 7 x + y + 2z = 7
18 x1 + 3x2 − x3 = 2 22 x1 + x2 + x3 = 1
x1 − x2 + 5x3 = 3 − x1 + x2 + x3 = 3
3x1 + x2 − x3 = 1 2x1 + x2 − 3x3 = 5
19 2x + 2y + 6z = 8 23 x1 + 6x2 + 3x3 = 4
x + 2y − z = 1 2x1 + x2 + 2x3 = 1
x + y − 3z = 1 3x1 + x2 + x3 = 5
20 x1 + x2 − 2x3 = 1 24 2x − 4y + 3z = 10
2x1 − 3x2 + x3 = −8 3x + y − 2z = 6
3x1 + x2 + 4x3 = 7 x + 3y − z = 20
49
Chapter 4
INTEGRATION
We begin with a definition of antiderivatives and indefinite integrals. Then, we provide basic integration rules.
4.1.1 Antiderivatives
Example 4.1
(1) Consider the functions F(x) = x3 + 4x2 − x + 5 and f (x) = 3x2 + 8x − 1.
0
Since F (x) = 3x2 + 8x − 1 = f (x), then the function F(x) is an antiderivative of f (x).
(2) Consider the functions G(x) = tan x + x2 − 1 and g(x) = sec2 x + 2x.
0
Since G (x) = sec2 x + 2x = g(x), then the function G(x) is an antiderivative of g(x).
Now, assume that F(x) is an antiderivative function of a function f (x), then every function F(x) + c is also antiderivative of f (x), where
c is a constant. For different values of the constant c, we have different antiderivatives, but they are very similar geometrically. The
upcoming theorem states that any antiderivative G(x), which is different from F(x) can be expressed as F(x) + c where c is an arbitrary
constant. In particular, if F(x) and G(x) are antiderivative functions of f (x), then
G(x) = F(x) + c .
F(x) = x3 + 2,
G(x) = x3 − 21 , F(x) − G(x) = x3 + 2 − (x3 − 12 ) = 52 = c
√ √
√ F(x) − H(x) = x3 + 2 − (x3 + 3 2) = 2 − 3 2 = c
H(x) = x3 + 3 2
50 INTEGRATION
are antiderivative functions of the function f (x). The difference between any two antiderivative functions is a constant c. Therefore,
F(x) = x3 + c is a general form of the antiderivatives of the function f (x) = 2x.
From Theorem 4.2, if F(x) + c is an antiderivative function of f (x), then there exist no antiderivative R
functions in different forms for
the function f (x). This leads us to define the indefinite integral. We introduce a symbol, namely, f (x) dx which will represent the
antiderivative the function f (x) and it is read as the indefinite integral of the function f with respect to x.
Definition 4.3 Let f be a continuous function on an interval I. The indefinite integral of f is the general antiderivative of f on
I: Z
f (x) dx = F(x) + c.
Z
The function f is called the integrand, the symbol is the integral sign, x is called the variable of the integration and c is the
constant of the integration.
Now, by using the previous definition, the general antiderivatives in Example 4.1 are
Z
(1) (3x2 + 8x − 1) dx = x3 + 4x2 − x +c.
| {z }
= F(x)
Z
(2) (sec2 x + 2x) dx = tan x + x2 +c.
| {z }
= F(x)
d
Z Z Z Z
1. f (x) dx = f (x). 3. f (x) ± g(x) dx = f (x) dx ± g(x) dx.
dx
Z Z
d
Z
2. (F(x)) dx = F(x) + c. 4. k f (x) dx = k f (x) dx, where k is a constant.
dx
Notes
The properties 3. and 4. can be generalized to a finite number of functions f1 , f2 , ..., fn and real numbers, k1 , k2 , ..., kn :
Z Z Z Z
k1 f1 (x) + k2 f2 (x) + ... + kn fn (x) dx = k1 f1 (x) dx + k2 f2 (x) dx + ... + kn fn (x) dx .
For property 3., we shall write only one constant of integration in the final answer.
In this section, we are given the derivative of certain functions and asked to find their antiderivatives. We will notice that the antiderivative
functions are obtained directly from the corresponding formula for differentiation. That is, by knowing the derivation formulas, we can
write the corresponding formulas for integration.
Properties of Indefinite Integrals 51
Rule 1: Power of x.
d xn+1 xn+1
Z
( ) = xn , so xn dx = + c for n 6= −1 . Note
dx n + 1 n+1
Note that Rule 1 cannot be applied
In words, to integrate the function xn , we add 1 to the power (i.e., for n = −1.
n + 1) and divide the function by n + 1.
For this value, the formula gives
For n = 0, we have a special case
x0
Z
Z x−1 dx = =∞.
1 dx = x + c. 0
d xn+1
Recall, if c ∈ R, then dx n+1 + c = xn . Thus, antiderivatives are
not unique.
1
Z
(3) x2 − dx
x3
Solution:
x2
Z Z Z
(1) (x + 1) dx = x dx + 1 dx = +x+c .
2
4x4 2 3 2
Z Z Z Z
(2) (4x3 + 2x2 + 1) dx = 4x3 dx + 2x2 dx + 1 dx = + x + x + c = x4 + x3 + x + c .
4 3 3
1 x3 x−2
Z Z Z
(3) x2 − dx = x2 dx − x−3 dx = + +c .
x3 3 2
Z Z
d
• dx (sin x) = cos x ⇒ cos dx = sin x + c • d
dx (cot x) = − csc x ⇒ − csc2 x dx = cot x + c
Z Z
d d
• dx (cos x) = − sin x ⇒ − sin x dx = cos x + c • dx (sec x) = sec x tan x ⇒ sec x tan x dx = sec x + c
Z Z
• d
dx (tan x) = sec2 x ⇒ sec2 x dx = tan x + c • d
dx (csc x) = − csc x cot x ⇒ − csc x cot x dx = csc x + c
1
Z
(2) + cos x dx
sec x
tan x
Z
(3) − x2 dx
cos x
Solution:
Z Z Z
(1) (cos x − csc x cot x) dx = cos x dx − csc x cot x dx = sin x + csc x + c .
52 INTEGRATION
1
Z Z Z
(2) − sin x dx = cos x dx − sin x dx = sin x + cos x + c .
sec x
tan x x3
Z Z Z
(3) − x2 dx = tan x sec x dx − x2 dx = sec x − + c .
cos x 3
ln : (0, ∞) −→ R ,
If u = g(x) is a differentiable function, then
0
Z 0
u e : R −→ (0, ∞) ,
d
ln |u| = uu =⇒
• dx dx = ln |u| + c
u y = ex ⇔ ln y = x
3 y
As special case:
1
Z
d
1 2
dx ln |x| = x =⇒ dx = ln |x| + c
x y = ex
1 y = ln x
Z
d
eu = eu .u0 =⇒ eu u0 dx = eu + c
• dx x
1 2 3 4
As special case: Z −1
d x x ex dx = ex + c
dx e = e =⇒
−2
Solution:
3 1
Z Z
(1) dx = 3 dx = 3 ln |x| + c.
x x
1 x4
Z Z Z Z
x3 + x−1 + ex dx = x3 dx + dx + ex dx = + ln |x| + ex + c.
(2)
x 4
1 1 1 1
Z Z Z
(3) −x
+ 2 dx = ex dx + x−2 dx = ex − x−1 + c.
3e x 3 3
1
Z
(4) dx = ln |x + 2| + c.
x+2
3x2 − 1 1 2(3x2 − 1)
Z Z
(5) dx = dx = ln |2x3 − x2 + 1| + c.
2x − x2 + 1
3 2 2x3 − x2 + 1
Properties of Indefinite Integrals 53
sin x − sin x
Z Z
(6) dx = − dx = ln | cos x| + c.
cos x cos x
Z
2 2
(7) 2x e(x +1)
dx = e(x +1)
+c
Z Z
(8) sin x ecos x dx = − − sin x ecos x dx = ecos x + c.
√
e x √ 1 √
Z Z
√ dx = 2 x √ dx = 2e x + c.
(9) e
x 2 x
etan x 1
Z Z Z
(10) dx = etan x
dx = etan x
sec2 x dx = etan x + c.
cos2 x cos2 x
d
sec−1 √1
d
cos−1 1 •
• = −√ dx =
dx 1−x2 x x2 −1
d
csc−1 √1
d
tan−1 1 • =−
• dx = 1+x2 dx x x2 −1
In general, we have
1 x
Z
• √ dx = sin−1 +c
2
a −x 2 a
1 1 x
Z
• dx = tan−1 +c
a2 + x2 a a
1 1 x
Z
• √ dx = sec−1 +c
2
x x −a2 a a
6
Z
(1) dx
4 + x2
1 1
Z
(2) √ + √ dx
1−x 2 x
1
Z
(3) 5x + √ dx
2
x x −5
Solution:
6 1 x
Z Z
(1) 2
dx = 6 2
dx = 3 tan−1 ( ) + c.
4+x 4+x 2
Z
1 1 Z
1
Z
1 √
(2) √ + √ dx = √ dx + x− 2 dx = sin−1 x + 2 x + c.
1−x 2 x 1−x 2
1 1 5x2 1 x
Z Z Z
(3) 5x + √ dx = 5x dx + √ dx = + √ sec−1 ( √ ) + c.
x x2 − 5 x x2 − 5 2 5 5
Notes
54 INTEGRATION
When evaluating an integral, we can always verify our answer by deriving the result of the integral.
In the previous examples, we use x as a variable of the integral. However, for this role, we can use any variable such as y, z, t,
etc . That is, instead of f (x) dx, we can integrate f (y) dy or f (t) dt.
Solution:
y3 y2
Z
(1) (y2 + y + 1) dy = + +y+c .
3 2
Z
(2) (cos t + sec2 t) dt = sin t + tan t + c .
Summation (or sigma notation) is a simple form used to give a concise expression for a sum of values.
n
Definition 4.5 Let {a1 , a2 , ..., an } be a set of numbers. The symbol ∑ ak represents their sum:
k=1
n
∑ ak = a1 + a2 + ... + an .
k=1
Solution:
3
(1) ∑ i2 = 12 + 22 + 32 = 1 + 4 + 9 = 14.
i=1
5
(2) ∑ ( j + 1) = (1 + 1) + (2 + 1) + (3 + 1) + (4 + 1) + (5 + 1) = 2 + 3 + 4 + 5 + 6 = 20.
j=1
3
(3) ∑ (k + 1)k2 = (1 + 1)(1)2 + (2 + 1)(2)2 + (3 + 1)(3)2 = 2 + 12 + 36 = 50.
k=1
Theorem 4.6 Let {a1 , a2 , ..., an } and {b1 , b2 , ..., bn } be sets of real numbers. If n is any positive integer, then
n
1. ∑ c = c + c + ... + c = nc for any c ∈ R.
k=1 | {z }
n-times
n n n
2. ∑ (ak ± bk ) = ∑ ak ± ∑ bk .
k=1 k=1 k=1
n n
3. ∑ c ak = c ∑ ak for any c ∈ R.
k=1 k=1
Definite Integrals 55
Solution:
10
(1) ∑ 6 = (10)(6) = 60.
k=1
3 3 3
(2) ∑ (k2 + 2k) = ∑ k2 + 2 ∑ k = (12 + 22 + 32 ) + 2(1 + 2 + 3) = 14 + 12 = 26.
k=1 k=1 k=1
4 4
(3) ∑ 3(k + 1) = 3 ∑ (k + 1) = 3(2 + 3 + 4 + 5) = 42.
k=1 k=1
Riemann sum is a mathematical model where one of its applications is to approximate the areas of the regions bounded by the graphs of
the functions. In the following, we will provide some basic definitions that we need to define the definite integral.
Definition 4.7 A set P = {x0 , x1 , x2 , ..., xn } is called a partition of a closed interval [a, b] if for any positive integer n,
Notes
The division of the interval [a, b] by a partition P generates n subintervals: [x0 , x1 ], [x1 , x2 ], [x2 , x3 ], ..., [xn−1 , xn ].
The largest length among ∆x1 , ∆x2 , ∆x3 , ..., ∆xn is called the norm of the partition P:
The partition P of the interval [a, b] is regular if ∆x0 = ∆x1 = ∆x2 = ... = ∆xn = ∆x.
x1 = x0 + ∆x ,
x2 = x1 + ∆x = (x0 + ∆x) + ∆x = x0 + 2∆x ,
x3 = x2 + ∆x = (x0 + 2∆x) + ∆x = x0 + 3∆x.
Riemann Sum
Definition 4.8 Let f be a bounded and defined function on a closed bounded interval [a, b] and let P = {x0 , x1 , ..., xn } be
a partition of [a, b]. Let ω = (ω1 , ω2 , ..., ωn ) be a mark on the partition P where ωk ∈ [xk−1 , xk ], k = 1, 2, 3, ..., n. Then the
Riemann sum of f with respect to the partition P and the mark ω is
n
R( f , P, ω) = ∑ f (ωk )∆xk .
k=1
As shown in Figure 4.3, the amount f (ω1 )∆x1 is the area of the rectangle A1 , f (ω2 )∆x2 is the area of the rectangle A2 and so on. The
sum of these areas approximates the area of the whole region under the graph of the function f (x) from x = a to x = b. This indicates
that if the function f is bounded and non-negative on a closed bounded interval [a, b] and P = {x0 = a, x1 , . . . , xn = b} is a partition of
that interval where ω = (ω1 , ω2 , ..., ωn ) is a mark on the partition P, then the Riemann sum estimates the area of the region under the
function f (x) from x = a to x = b. As the number of the subintervals increases n → ∞ (i.e., ||P|| → 0), the estimation becomes better.
Therefore,
n
A = lim R( f , P, ω) = lim ∑ f (ωk )∆xk (4.1)
kPk→0 kPk→0 k=1
The following definition shows that the definite integral of a defined function f on a closed bounded interval [a, b] is a Riemann sum
when k P k→ 0.
Definite Integrals 57
Definite Integrals
Definition 4.9 For any function f bounded and defined on a closed bounded interval [a, b], the definite integral of f from a to
b is Z b
a
f (x) dx = lim
n→∞ ∑ f (ωk )∆xk , (k P k→ 0)
k
if the limit exists. The numbers a and b are called the limits of the integration.
Notes
The limit in the previous definition is over all points in the partition P = {([xk−1 , xk ], ωk )}1≤k≤n . When the limit exists, we say
that f is Riemann integrable (or integrable) on [a, b].
Any continuous function is Riemann integrable on a closed bounded interval [a, b].
Theorem
Z 4.10 b
1. c dx = c(b − a).
a
Z a
2. f (x) dx = 0 if f (a) exists.
a
Solution:
Z 2
(1) 5 dx = 5(2 − 0) = 10.
Z03
(2) (x2 − 1) dx = 0.
3
Theorem 4.11
1. Let f and g be integrable functions on [a, b], then the functions f + g and f − g are integrable on [a, b] and
Z b Z b Z b
f (x) ± g(x) dx = f (x) ± g(x) dx.
a a a
2. Let f be integrable function on [a, b] and k ∈ R, then the function k f is integrable on [a, b] and
Z b Z b
k f (x) dx = k f (x) dx.
a a
Z b Z b Z b
g(x)
Example 4.11 If f (x) dx = 5 and g(x) dx = 9, find the value of the integral 4 f (x) − dx.
a a a 3
Solution:
58 INTEGRATION
Z b Z b Z b
g(x) 1
4 f (x) − dx = 4 f (x) dx − g(x) dx (Using Theorem 4.11 )
a 3 a 3 a
1
= 4(5) − (9) = 17 .
3
Theorem 4.12
1. If f and g are integrable on [a, b] and f (x) ≥ g(x) for all x ∈ [a, b], then
Z b Z b
f (x) dx ≥ g(x) dx.
a a
Z 2 Z 2
Example 4.12 Prove that (x3 + x2 + 2) dx ≥ (x2 + 1) dx without evaluating the integrals.
0 0
Solution: Let f (x) = x3 + x2 + 2 and g(x) = x2 + 1. We can find that f (x) − g(x) = x3 + 1 > 0 for all x ∈ [0, 2]. This implies that
f (x) > g(x) and from Theorem 4.12, we have
Z 2 Z 2
(x3 + x2 + 2) dx ≥ (x2 + 1) dx .
0 0
Theorem 4.13 If f is integrable on the intervals [a, c] and [c, b], then f is integrable on [a, b] and
Z b Z c Z b
f (x) dx = f (x) dx + f (x) dx.
a a c
Theorem 4.15 Suppose that f is continuous on the closed interval [a, b].
Z x
1. If F(x) = f (t) dt for every x ∈ [a, b], then F(x) is an antiderivative of f on [a, b].
a Z b
2. If F(x) is any antiderivative of f on [a, b], then f (x) dx = F(b) − F(a).
a
Z b
From the previous theorem, the definite integral f (x) dx is evaluated by two steps:
a
Techniques of Integration 59
Z
Step 1: Find the indefinite integral f (x) dx = F(x) and no need to put a constant c.
Step 2: Evaluate the antiderivative F at upper and lower limits by substituting x = b and x = a then calculate F(b) − F(a).
Solution:
Z 2 h i2
(2x + 1) dx = x2 + x = 4 + 2 − (−1)2 + (−1) = 6 − 0 = 6.
(1)
−1 −1
Z 3 h x3
27 i3
(2) (x2 + 1) dx = + 3) − 0 = 12.
+x =(
0 33 0
√
Z 2
1 h −2 i2 −2 −2 + 2 2 √
(3) √ dx = √ = √ − (−2) = √ = − 2 + 2.
1 x 3 x 1 2 2
Z π iπ
2
h
2 π π π
(4) (sin x + 1) dx = − cos x + x = (− cos + ) − (− cos 0 + 0) = + 1.
0 0 2 2 2
Z π iπ
h π π
(5) π (sec2 x − 4) dx = tan x − 4x π = tan π − 4π − tan − 4 = −4π − (1 − π) = −3π − 1.
4 4 4 4
x2 i π3 π ( π )2 π2 π2
π
0
Z h
3
(6) (sec x tan x + x) dx = sec x + = sec + 3 − sec 0 + = 2+ −1 = 1+ .
0 2 0 3 2 2 18 18
In the previous section, we demonstrated the importance of basic integral rules and the integral properties in evaluating several integrals.
Those rules are elementary for simple functions, but in this section we will study new techniques that will enable us to evaluate more
complex functions.
The integration by substitution (known as u-substitution) is one of the important methods for evaluating the integrals. The method is
based on changing the variable of the integral to obtain a simple integral.
Theorem 4.16 Let g be a differentiable function on an interval I where the derivative is continuous. Let f be continuous on
the interval J contains the range of the function g. If F is an antiderivative of the function f on J, then
Z
0
f (g(x))g (x) dx = F(g(x)) + c, ∀x ∈ I.
0
In the following examples, we will check if the integrand has the form f (g(x))g (x) so that we can use Theorem 4.16.
Z
Example 4.14 Evaluate the integral 3x2 (x3 + 1)5 dx.
60 INTEGRATION
5 Note
theorem by assuming f (x) = x
Solution: We can use the previous
and g(x) = x3 + 1, then f g(x) = (x3 + 1)5 . Z
f (x) dx = F(x) + c
0
Since g (x) = 3x2 , then
x6
Z
⇒ x5 dx = +c
(x3 + 1)6 6
Z
3x2 (x3 + 1)5 dx = + c.
6
Note that in the previous example, we can end with the same solution by using the following steps of the substitution method.
Steps of the integration by substitution:
Step 1: Choose a new variable u. Observe the integrand f (x) and choose an inside function u as a function of x. Then check if f (x) can
be decomposed into
f (x) = (function of u).constant multiple of u0 (x)
Step 2: Determine the value of du.
Step 3: Make the substitution i.e., eliminate all occurrences of x in the integral by making the entire integral in terms of u.
Step 4: Evaluate the new integral.
Step 5: Return the evaluation to the initial variable x.
In Example 4.14, let u = x3 + 1, this implies du = 3x2 dx. Now apply the substitution by substituting all x-terms into u-terms:
u6
Z Z
(x3 + 1)5 |{z}
3x2 dx = u5 du = +c .
| {z } 6
u du
Z p
Example 4.15 Evaluate the integral (3x − 1) 3x2 − 2x + 1 dx.
√ √
Solution: Assume f (x) = x and g(x) = 3x2 − 2x + 1, then f g(x) = 3x2 − 2x + 1. Since g0 (x) = 6x − 2 = 2(3x − 1), then from
Theorem 4.16,
1
Z p Z p
(3x − 1) 3x2 − 2x + 1 dx = 2(3x − 1) 3x2 − 2x + 1 dx
2
3
1 2 3 (3x2 − 2x + 1) 2
= (3x2 − 2x + 1) 2 + c = +c .
2 3 3
By using the steps of the substitution method, let u = 3x2 − 2x + 1 and this implies du = (6x − 2)dx = 2(3x − 1)dx. By substitution, we
have 3 3
√ du 1 Z 1 u2 (3x2 − 2x + 1) 2
Z
u = u 2 du = +c = +c .
2 2 3 3
√
cos x
Z
Example 4.16 Evaluate the integral √ dx.
x
√ √ 0 √
Solution: Assume f (x) = cos x and g(x) = x, then f g(x) = cos x. Since g (x) = 1/(2 x), then from Theorem 4.16, we have
√ √
Z
cos x
Z
cos x √
√ dx = 2 √ dx = 2 sin x + c .
x 2 x
√ 1
By using the steps of the substitution method, let u = x and this implies du = √ dx. By substitution, we obtain
2 x
Z
√
2 cos u du = 2 sin u + c = 2 sin x + c.
Techniques of Integration 61
Z
Example 4.17 Evaluate the integral tan x dx.
sin x
Solution: Write tan x = and assume u = cos x. This implies du = − sin x dx, then −du = sin x dx.
cos x
Hence,
sin x 1 Rules
Z Z Z
tan x dx = dx = . sin x dx
cos x cos x Z
• tan x dx = − ln | cos x| + c
1
Z
=− du = − ln |u| + c (Using Rule 3 on page 52 ) . = ln | sec x| + c
u Z
• cot x dx = ln | sin x| + c
By returning the evaluation to the initial variable x, we have
= − ln | csc x| + c
Z
tan x dx = − ln | cos x| + c .
x2 − 1
Z
Example 4.18 Evaluate the integral dx.
(x3 − 3x + 1)6
Solution:
Let u = sin x, then du = cos x dx. By substitution, we
have Note
Any power of a trigonometric
u4
Z
u3 du = +c function can be integrated by Rule
4 1 on page 51 when accompanied by
sin4 x its differential.
= + c.
4
sin−1 x
Z
Example 4.20 Evaluate the integral √ dx.
1 − x2
Z
Corollary 4.17 If f (x) dx = F(x) + c, then for any a 6= 0,
1
Z
f (ax ± b) dx = F(ax ± b) + c.
a
Z
(2) sec2 (5x + 4) dx
Notes
Z b
When using the substitution method to evaluate the definite integral f (x) dx, we have two options:
a
Option 1: Change the limits of integration to the new variable. In this case, we evaluate the integral without returning to the
original variable.
Option 2: Leave the limits in terms of the original variable. In this case, we evaluate the integral and return the result to the
original variable. After that, we substitute x = b and x = a into the antiderivative.
Z 1 p
Example 4.22 Evaluate the integral 2x x2 + 1 dx.
0
Solution:
Option 1: Let u = x2 + 1, this implies du = 2x dx. Change the limits u(0) = 1 and u(1) = 2. By substitution, we have
Z 2
2 3 i2 2 3 3
2 √
u1/2 du = u 2 = 22 −12 = 2 2−1 .
1 3 1 3 3
2 3 2
Z
3
Option 2: Let u = x2 + 1, then du = 2x dx. By substitution, we have u1/2 du = u 2 = (x2 + 1) 2 + c. Thus,
3 3
Z 1 p
2 h 3
i1 2 √
2x x2 + 1 dx = (x2 + 1) 2 = 2 2−1 .
0 3 0 3
π
sin x
Z
2
Example 4.23 Evaluate the integral dx.
0 1 + cos2 x
Solution: Let u = cos2 x, this implies du = − sin x dx. Change the limits u(0) = 1 and u( π2 ) = 0. By substitution, we have
Z 0 Z 1 i1
1 1 h π π
− du = du = tan−1 u = tan−1 (1) − tan−1 (0) = ( − 0) = .
1 1 + u2 0 1+u 2 0 4 4
In this section we will learn another important technique of integration, called integration by parts. This technique depends on the
product rule in differentiation, so it can be thought as the product formula for integration. In practice, the integrand is divided into two
parts u and dv, then we find du by deriving u and v by integrating dv. This method transfers a product integral (the original integrand)
into another product integral that can be evaluated.
Theorem 4.18 If u = f (x) and v = g(x) such that f 0 and g0 are continuous, then
Z Z
u dv = uv − v du.
Techniques of Integration 63
Z Z
Theorem 4.18 shows that the integration by parts transfers the integral u dv into an easier integral v du. The question here is, what
we choose as u and what we choose as dv = v0 dx.
(1) Choose u a portion of the integrand that can be easily differentiated. We can choose u to be the function that comes first in this list:
(a) Inverse trigonometric function.
(b) Logarithmic function.
(c) Algebraic function.
(d) Exponential function.
(e) Trigonometric function.
(2) Choose dv the most complicated portion of the integrand that can be easily integrated.
This guideline is useful, but is not enough to solve all product integrals. Sometimes we need to analyze the integrand and examine the
best way of using integration by parts.
Z
Example 4.24 Evaluate the integral x ex dx.
Solution: Note
The integrand xex is a product of two functions x and ex . Now we • We choose u = x because it can be
need to identify one function as u and the other one as dv such that differentiated to a constant. Thus the
R
the new product integral v du is easier than the original integrand. new product integral will not involve
a product anymore.
Z • Try to choose
Let I = x ex dx and choose u = x, and dv = ex dx. Then, u = ex and dv = x dx
u = x ⇒ du = dx , We will obtain
Z ttttttttttttttttttttttttttttttttttttttttttttttttttt
2
x x x2 x
Z
dv = ex dx ⇒ v = ex dx = ex . I= e − e dx.
2 2
From Theorem 4.18, we have x2 x
Z
However, the integral e dx is
Z 2
I = x ex − ex dx = x ex − ex + c. more
Z difficult than the original one
xex dx.
Z
Example 4.25 Evaluate the integral x cos x dx.
Z
Solution: In the same manner as in the preceding example, let I = x cos x dx. Set u = x and dv = cos x dx. Hence,
u = x ⇒ du = dx ,
Z
dv = cos x dx ⇒ v = cos x dx = sin x. Try to choose
u = cos x and dv = x dx
From Theorem 4.18, we have
Z Do you have the same result?
I = x sin x − sin x dx
= x sin x + cos x + c .
Z
Example 4.26 Evaluate the integral ln x dx.
64 INTEGRATION
Z
Solution: Let I = ln x dx and choose u = ln x, and dv = dx. Then,
Remember
1 If u = g(x) is differentiable, then
u = ln x ⇒ du = dx ,
x
d u0
Z
dv = dx ⇒ v = 1 dx = x. (ln |u|) =
dx u
Z
Example 4.27 Evaluate the integral x3 ln x dx.
Solution:
Z
Let I = x3 ln x dx and choose u = ln x, and dv = x3 dx. Then,
Rule
1
u = ln x ⇒ du = dx , In the same manner
Z as in Example
x
4.28, to evaluate xn ln x dx, let
x4
Z
dv = x3 dx ⇒ v = x3 dx = .
4 u = ln x ⇒ du = 1x dx
From Theorem 4.18, we have xn+1
dv = xn dx ⇒ v = xn dx =
R
n+1
x4 x4 1
Z
I= ln x − dx Hence,
4 4 x
x4 1 xn+1 xn+1
Z Z
= ln x − x3 dx xn ln x dx = ln x+ +c
4 4 n+1 (n + 1)2
x4 x4
= ln x − +c .
4 16
Notes
Remember when we consider the integration by parts, we want to obtain an easier integral. As we saw in Example 4.28, if we
Z 2
x x
choose u = ex and dv = x dx, we have e dx which is more difficult than the original integral.
2
When considering the integration by parts, we have to choose dv a function that can be integrated (see Example 4.28 ).
Z
Example 4.28 Evaluate the integral sin x ln(cos x) dx.
Solution:
Z
Let I = sin x ln(cos x) dx and choose u = ln(cos x) for cos x > 0, and dv = sin x dx. Then,
− sin x
u = ln(cos x) ⇒ du = dx = − tan x dx,
Z
cos x
dv = sin x dx ⇒ v = sin x dx = − cos x.
Hence,
Techniques of Integration 65
Z 1
Example 4.29 Evaluate the integral tan−1 x dx.
0
Solution:
Z
Let I = tan−1 x dx and choose u = tan−1 x, and dv = dx. Hence, Rule
Any inverse trigonometric function
1 such that the differential is not
u = tan−1 x ⇒ du = dx ,
x2 + 1 accompanied can be integrated in
Z
the same manner of Example 4.29.
dv = dx ⇒ v = 1 dx = x.
Therefore, Z 1 h 1 i1 1 1 π √
tan−1 x dx = x tan−1 x − ln(x2 + 1) = (tan−1 (1) − ln 2) − (0 − ln 1) = − ln 2.
0 2 0 2 2 4
Note that sometimes we need to use the integration by parts twice as in the following examples.
Z
Example 4.30 Evaluate the integral x2 ex dx.
Solution:
Z
Let I = x2 ex dx and choose u = x2 , and dv = ex dx. Then, Note
In successive application of the
u = x2 ⇒ du = 2x dx ,
Z integration by parts, do not switch
dv = ex dx ⇒ v = ex dx = ex . choices for u and dv. In Example
4.30, we choose dv = ex in integrals
Z I and J.
This implies I = x2 ex − 2 xex dx.
Z Z
We use the integration by parts again for the integral xex dx where we assume J = xex dx.
Therefore, Z
J = xex − ex dx = xex − ex + c .
By substituting the result of J into I, we have
I = x2 ex − 2(xex − ex ) + c
= ex (x2 − 2x + 2) + c .
66 INTEGRATION
Z
Example 4.31 Evaluate the integral ex cos x dx.
Z
Solution: Let I = ex cos x dx and choose u = ex , and dv = cos x dx. Note
Then, In Example 4.31,
u = ex ⇒ du = ex dx , try to choose
Z
dv = cos x dx ⇒ v = cos x dx = sin x. u = cos x and dv = ex
Z
Do you have the same result?
Hence, I = ex sin x − ex sin x dx.
Z Z
The integral ex sin x dx cannot be evaluated. Therefore, we use the integration by parts again where we assume J = ex sin x dx.
u = ex ⇒ du = ex dx ,
Z
dv = sin x dx ⇒ v = sin x dx = − cos x.
Hence, Z
J = −ex cos x + ex cos x dx .
I = ex sin x − J
Z
= ex sin x + ex cos x − ex cos x dx
⇒ I = ex sin x + ex cos x − I.
This implies
1
2I = ex sin x + ex cos x ⇒ I = (ex sin x + ex cos x)
2
ex
Z
⇒ ex cos x dx = (sin x + cos x) + c .
2
f (x)
A rational function is a quotient of two polynomials of the form q(x) = g(x)
. A Polynomial f (x) is a linear sum of powers of x, for
example f (x) = 5x3 + x2 + x + 1
or g(x) = x4 − x.
The degree of a polynomial is the highest power occurring in the polynomial, for
example the degree of f (x) is 3 and the degree of g(x) is 4.
ä Step 1: If the degree of f (x) is equal or greater than the degree of g(x), we do polynomial long-division; otherwise we move to step 2.
By doing the long-division, we reduce the fraction to a
mixed quantity.
h(x)
f (x) r(x) g(x) f(x)
q(x) = = h(x) + ,
g(x) g(x) - ...
...
where h(x) is the quotient and r(x) is the remainder. r(x)
The degree of the numerator of the new fraction will
be less than the degree of the denominator.
ä Step 2: Factor the denominator g(x) into irreducible polynomials where the factors are either linear or irreducible quadratic polynomials.
Techniques of Integration 67
f (x) r(x)
ä Step 3: Find the partial fraction decomposition. This step depends on the result of step 2 where the fraction g(x) or g(x) can be written
as a sum of partial fractions:
q(x) = P1 (x) + P2 (x) + P3 (x) + ... + Pn (x) ,
Ak Ak x + Bk
where Pk (x) = , n ∈ N or Pk (x) = if b2 − 4ac < 0. The constants Ak and Bk are real numbers and computed
(ax + b)n (ax2 + bx + c)n
later. Note that the denominators of the fractions Pk (x) are the factors of the original denominator obtained in step 2.
f (x) A1 A2 A3 An
= + + + ... +
g(x) a1 x + b1 a2 x + b2 a3 x + b3 an x + bn
Case 2: The denominator g(x) has repeated linear factors of the form (ai x + bi )k where k > 1. Then,
f (x) A1 A2 A3 An
= + + + ... + .
g(x) (ai x + bi ) (ai x + bi )2 (ai x + bi )3 (ai x + bi )k
Case 3: The denominator g(x) has factors which are irreducible quadratics of the form ai x2 + bi x + ci where b2i − 4ai ci < 0. In this
Ai x + Bi
case, we include terms of the form .
ai x2 + bi x + ci
Solution:
Illustration:
Step 1: Do the polynomial long-division.
Since the degree of the denominator g(x) is less than the degree of
the numerator f (x), we do the polynomial long-division given on 2x −10
x2 + 3x + 2 2x3 −4x2
−15x +5
the right side. Hence, we have
−(2x3 +6x2 +4x)
11x + 25 −10x2 −19x +5
q(x) = (2x − 10) + . −(−10x2 −30x −20)
x2 + 3x + 2 11x +25
Illustration:
Multiply equation 1 by −2 and
x+1
Z
Example 4.33 Evaluate the integral dx.
x2 − 2x − 8
Solution:
Step 1: This step can be skipped since the degree of f (x) = x + 1 is less than the degree of g(x) = x2 − 2x − 8.
Step 2: Factor the denominator g(x) into irreducible polynomials
g(x) = x2 − 2x − 8 = (x + 2)(x − 4) .
Here we have case 2 in factoring the denominator g(x).
Step 3: Find the partial fraction decomposition
x+1 A B Ax − 4A + Bx + 2B
= + = .
x2 − 2x − 8 x + 2 x − 4 (x + 2)(x − 4)
We need to find the constants A and B by equating the coefficients of like powers of x in the two sides of the equation:
x + 1 = (A + B)x + (−4A + 2B)
constants: − 4A + 2B = 1 → 2 4A + 4B = 4
1
−4A + 2B = 1
By doing some calculation, we obtain A = 6 and B = 56 . Thus,
−−−−−−−−−
x+1 1/6 5/6 6B = 5
= + .
x2 − 2x − 8 x + 2 x − 4
2x2 − 25x − 33
Z
Example 4.34 Evaluate the integral dx.
(x + 1)2 (x − 5)
Solution:
Steps 1 and 2 can be skipped in this example. According to the cases of factoring the denominator, we have cases 1 and 2.
2x2 − 25x − 33 5 1 −3
Z Z Z Z
dx = dx + dx + dx
(x + 1)2 (x − 5) x+1 (x + 1)2 x−5
Z
= 5 ln | x + 1 | + (x + 1)−2 dx − 3 ln | x − 5 |
1
= 5 ln | x + 1 | − − 3 ln | x − 5 | +c.
(x + 1)
x+1
Z
Example 4.35 Evaluate the integral dx.
x(x2 + 1)
Solution:
Steps 1 and 2 can be skipped in this example. Here we have cases 1 and 3 of factoring the denominator.
coefficients of x2 : A+B = 0 → 1
coefficients of x: C=1→ 2
constants: A=1→ 3
We have A = 1, B = −1 and C = 1.
70 INTEGRATION
Exercises
1
Z
Z
√ 46 dx
39 x ln x dx x2 + 2x − 3
x2
Z Z 7
40 2
x sec x dx 47 dx
3 x2 − x − 2
3x2 − 10
Z Z
41 x e−4x dx 48 dx
x2 − 4x + 4
Z
x2 − 9
Z
42 (ln x)2 dx 49 dx
x−1
1
Z
2x4 − 3x3 − 10x2 + 2x + 11
Z
43 dx
(x − 3)(x − 1)2 50 dx
x3 − x2 − 5x − 3
1
Z
1
Z
44 dx 51 dx
x2 + 6x + 8 1 + ex
x
Z
45 dx 2x3 − 18x2 + 29x − 4
Z
x2 − x − 2 52 dx
(x + 1)(x − 2)3
Z b Z c Z b
53 - 58 If f (x) dx = 2, f (x) dx = 2 and g(x) dx = 3 where c ∈ (a, b), evaluate the integral.
a b a
Z a Z a
53 f (x) dx 56 5 f (x) − 3g(x) dx
b b
Z c Z b
54 f (x) dx 1
57 f (x) + 7g(x) dx
a a 3
Z b Z a
55 2 f (x) + g(x) dx 58 4 f (x) + g(x) dx
a a
59 - 64 Use the properties of the definite integrals to prove the inequality without evaluating the integrals.
Z 1 Z 1 Z 3
59 x dx ≥ x2 dx 62 (x2 − 3x + 4) dx ≥ 0
0 0 0
Z 3 Z 3 Z 2√ Z 2√
x
60 dx ≥ x dx 63 5 − x dx ≥ x + 1 dx
0 x3 + 2 0 1 1
Z 4 Z 4 Z 2 p
61 (2x + 2) dx ≥ (3x + 1) dx 64 2 < 1 + x2 dx
1 1 −1
Exercises 73
sin x
Z
65 The value of the integral √ dx is equal to
2 + cos x
√ √
(a) −2
√ 2 + cos x + c (c) −√ 2 + cos x + c
(b) 2 + cos x + c (d) 2 2 + cos x + c
sin (tan x)
Z
66 The value of the integral dx is equal to
cos2 x
(a) cos (tan x) + c (c) − cos (tan x) + c
(b) sin (tan x) + c (d) − sin (tan x) + c
Z p
67 The integral x x2 + 1 dx is equal to
√ 3
(a) 21 x2 x2 + 1 + c (c) − 32 (x2 + 1) 2 + c
3 3
(b) 23 (x2 + 1) + c
2 (d) 13 (x2 + 1) + c
2
x
Z
68 The integral dx is equal to
cos2 x2
(a) 21 tan x2 + c (c) 12 tan x + c
(b) tan x2 + c (d) − cos1 x2 + c
sec2 x
Z
69 The value of the integral dx is equal to
cot2 x
1+cos2 x cot4 x
(a) 3 cos3 x
+c (c) 4 +c
1−3 cos2 x tan3 x
(b) 3 cos3 x
+c (d) 3 +c
cos x
Z
70 The value of the integral √ dx
4 + sin x
√ √
(a) 12 sin x + 4 + c √ x+4+c
(c) 2 sin
√
(b) sin x + 4 + c (d) −2 sin x + 4 + c
Z 1
71 The value of the integral 2 | x |3 dx
−1
(a) 2 (b) 1 (c) 0 (d) −1
74 APPLICATIONS OF INTEGRATION
Chapter 5
APPLICATIONS OF INTEGRATION
5.1 Areas
As shown in Chapter 4, if the function f is bounded and non-negative on a closed bounded interval [a, b] and P = {x0 = a, x1 , . . . , xn = b}
is a partition of that interval where ω = (ω1 , ω2 , ..., ωn ) is a mark on the partition P, then the Riemann sum estimates the area of the
region under the function f (x) from x = a to x = b:
n Z b
A = lim ∑ f (ωk )∆xk = f (x) dx .
kPk→0 k=1 a
In this section we find the area of ??the regions for the following cases:
Figure 5.1
Consider the region between the graph of the function y = f (x), the x-axis, and the ordinates x = a and x = b as shown in Figure 5.2.
Now we want to find the area of the shaded region.
Areas 75
Remember
See the explanation on page 57.
Figure 5.2: The region R bounded by the graph of y = f (x), and the
ordinates x = a and x = b.
The sum of the rectangles areas approximates the area of the whole
region under the graph of the function f from x = a to x = b, where
as the number of the subintervals increases n → ∞ (||P|| → 0), the
estimation becomes better.
Z b
A= f (x) dx
a
Figure 5.3: The region R bounded by the graph of y = f (x), and the
ordinates x = a and x = b.
Example 5.1 Express the area of the shaded region as a definite integral then find the area.
76 APPLICATIONS OF INTEGRATION
(1) (2)
y y
f (x) = x2
f (x) = 2x + 1
R R
x
1 3 x 1 4
Figure 5.4
Solution:
Z 3 h i3 h i
(1) A = (2x + 1) dx = x2 + x = (32 + 3) − (12 + 1) = 12 − 2 = 10.
1 1
Z 4 h i4 h i
(2) A = x2 dx = 1
3 x3 = 31 64 − 1 = 63
3 = 21.
1 1
√
Example 5.2 Sketch the region bounded by the graph of y = x from x = 0 to x = 3, then find its area.
y
Solution:
The region bounded by the graph of the function y =
√
x in the interval [0, 3] is shown in the figure.
√
y= x
The area of the region is
Z 3
√
A= x dx
0
2h i3
= x3/2
3 0
√
= 2 3. x
3
Figure 5.5
Consider the region between the graph of the function x = f (y), the y-axis, and the ordinates y = c and y = d as shown in Figure 5.6.
Now we want to find the area of the shaded region.
Areas 77
Figure 5.6: The region R bounded by the graph of x = f (y), and the ordinates y = c and y = d.
Divide the interval [c, d] into n subintervals and choose y∗i in the ith subinterval. As shown in Figure ??, the area of the rectangle A1 is
f (y∗1 )∆y1 , the area of the rectangle A2 is f (y∗2 )∆y2 and so on.
The sum of the areas of the rectangles approximates the area of the whole region under the graph of the function x = f (y) from y = c to
y = d where as the number of the subintervals increases n → ∞ (||P|| → 0), the estimation becomes better.
Z d
A= f (y) dy
c
Example 5.3 Sketch the region bounded by the graph of x = y + 1 and y-axis over the interval [−1, 0], then find its area.
y
Solution: Figure 5.7 shows the region bounded by the
function x = y + 1 and y-axis over the interval [−1, 0].
Figure 5.7
78 APPLICATIONS OF INTEGRATION
√
Example 5.4 Sketch the region bounded by the graph of x = y from y = 0 to y = 1, then find its area.
y
√
Solution: The region bounded by the function x = y
in the interval [0, 1] is shown in Figure 5.8.
Figure 5.8
If f (x) and g(x) are continuous functions such that f (x) ≥ g(x) ∀x ∈ [a, b], then the area A of the region R bounded by the graphs of
f (x) (the upper boundary of R) and g(x) (the lower boundary of R) from x = a to x = b is subtracting the area of the region under g(x)
from the area of the region under f (x). This can be stated as follows:
Z b
A= f (x) − g(x) dx
a
If f (y) and g(y) are continuous functions such that f (y) ≥ g(y) ∀y ∈ [c, d], then the area A of the region R bounded by the graphs of
f (y) (the right boundary of R) and g(y) (the left boundary of R) from y = c to y = d is subtracting the area of the region bounded by g(y)
from the area of the region bounded by f (y). This can be stated as follows:
Z d
A= f (y) − g(y) dy
c
Areas 79
y y
f (x)
R
R
g(y)
g(x) f (y)
c x
x
a b
Figure 5.9: The area of the region bounded by the graphs of the tow functions f and g.
Example 5.5 Express the area of the shaded region as a definite integral, then find the area.
(2)
(1)
y
y
g(x) = 2x − 1
f (x) = x2 + 2
f (x) = x
g(x) = x + 4
x
x
1 2 3
1 2 3 4
Figure 5.10
Solution:
(1) The area of the region bounded by the two curves f (x) and g(x) is
Z 3 Z 3 h x2 i3 h 9 i 3
A= (2x − 1) − x dx = (x − 1) dx = −x = 9 − 3) − (2 − 2) = .
2 2 2 2 2 2
80 APPLICATIONS OF INTEGRATION
Upper graph: y = x + 4
g(x) = x + 4
Lower graph: y = x2 + 2 f (x) = x2 + 2
Z 2
(x + 4) − (x2 + 2) dx
A1 =
1
Z 2
2 + x − x2 dx
=
1
x
h x2 x3 i2 13
= 2x + − = . 1 2
2 3 1 6
Figure 5.11
y
Region (2) is in the interval [2, 4].
Upper graph: y = x2 + 2
Lower graph: y = x + 4
f (x) = x2 + 2
Z 4
(x2 + 2) − (x + 4) dx
A2 =
2
Z 4
x2 + x − 2 dx
=
2 g(x) = x + 4
h x3 x2 i4 26
= − − 2x = .
3 2 2 3
x
13 26 65
The total area is A = A1 + A2 = + = . 2 4
6 3 6
Figure 5.12
Example 5.6 Sketch the region bounded by the graphs of y = x2 and y = x + 6 over the interval [−2, 3], then find its area.
Solution: The region bounded by the two functions is shown in Figure 5.13.
Areas 81
Z 3
A= (x + 6 − x2 ) dx
−2
h x2 x3 i3 y = x+6
= + 6x −
2 3 −2
h 27 22 i
= +
2 3 y = x2
125
= .
6
x
-2 3
Figure 5.13
Example 5.7 Sketch the region bounded by the graphs of y = x3 and y = x in the interval [−1, 1], then find its area.
Solution:
The figure on the right shows the region bounded by the two functions. The region is divided into two regions as follows:
Region (1) is in the interval [−1, 0]
Upper graph: y = x3 y
Lower graph: y = x
Z 0 h x4 x2 i0 h 1 1 i
A1 = (x3 − x) dx = − = 0−( − )
−1 4 2 −1 4 2
1
= .
4 y = x3
x
Region (2) is in the interval [0, 1] -1 1
Upper graph: y = x
Lower graph: y = x3 y=x
Z 1 h x2 x4 i1 h 1 1 i
A2 = (x − x3 ) dx = − = ( − )−0
0 2 4 0 2 4
1
= .
4
1 1 1
The total area is A = A1 + A2 = + = .
4 4 2
Figure 5.14
Example 5.8 Sketch the region bounded by the graphs of y = x2 and x = y2 over [0, 1], then find its area.
Solution:
The region bounded by the graphs of y = x2 and x = y2 over [0, 1] is displayed in Figure 5.15. We write the two functions in terms of x,
so the upper graph:
√
x = y2 ⇒ y = x.
82 APPLICATIONS OF INTEGRATION
Z 1
√ x = y2
A= ( x − x2 ) dx y = x2
0
h 2 3 x3 i1
= x2 − x
3 3 0 1
h2 1i
= −
3 3
1
= .
3
Figure 5.15
π
Example 5.9 Sketch the region determined by the graphs of y = sin x, y = cos x and y-axis over the interval [0, ]. Then find its
4
area.
Solution:
The figure on the right shows the region bounded by the two functions. Over the interval [0, π4 ], the two curves intersect at π4 .
Figure 5.16
y
Example 5.10 Sketch the region bounded by the graphs of x = 2y and x = 2 + 3, then find its area.
Solution:
Let f (y) = 2y + 3 and g(y) = 2y. To sketch the region bounded by the two functions, we find out whether the two functions are intersected.
Solids of Revolution 83
f (y) = g(y)
y
⇒ + 3 = 2y
2
⇒ y + 6 = 4y
⇒ y = 2.
Substitute y = 2 in f (y) or g(y), we have x = 4. Thus, the two curves intersect at (4, 2).
y
4
Based on the region given in Figure 5.17, the area is
Z 2 3
y
A= ( + 3 − 2y) dy
0 2
Z 2
3 2
= (− y + 3) dy x = 2y
0 2
h 3 i2
= − y2 + 3y
4 0 1
x = y/2 + 3
= −3 + 6 = 3.
x
1 2 3 4 5
Figure 5.17
Definition 5.1 If R is a plane region, the solid of revolution S is a solid generated from revolving R about a line in the same
plane where the line is called the axis of revolution.
Example 5.11 Let y = f (x) ≥ 0 be a continuous function for every x ∈ [a, b]. Let R be a region bounded by the graph of f and the
x-axis from x = a to x = b. The region revolution about x-axis generates a solid given in Figure 5.18 (right).
Figure 5.18: Revolution of the region R about x-axis. The figure on the left shows the region under the continuous function y = f (x) ≥ 0 over the
interval [a, b]. The figure on the right shows the solid S generated by revolving the region about x-axis.
84 APPLICATIONS OF INTEGRATION
Example 5.12 Let y = f (x) be a constant function from x = a to x = b, as in Figure 5.19. The region R is a rectangle and by
revolving it about x-axis, we obtain a circular cylinder.
Figure 5.19: Revolution of the rectangular region R about x-axis. The figure on the left shows the region under the constant function f (x) = c over the
interval [a, b]. The figure on the right shows the circular cylinder generated by revolving the region about x-axis.
Example 5.13 Consider a region R bounded by the graph of x = f (y) from y = c to y = d. Revolution of R about y-axis generates a
solid given in Figure 5.20.
Figure 5.20: Revolution of the region R about y-axis. The figure on the left displays the region under the function x = f (y) over the interval [c, d]. The
figure on the right displays the solid S generated by revolving the region about y-axis.
One of the interesting applications of the definite integral is computing the solids revolution. In this section, we study three methods to
compute the volumes of the revolution solids known as disk method, washer method and cylindrical shells method.
Let y = f (x) ≥ 0 be a continuous function for every x ∈ [a, b] and let R be a region bounded by the graph of f and x-axis form x = a to
x = b. Let S be a solid generated by revolving R about x-axis. Assume that P is a partition of [a, b] and ω = (ω1 , ω2 , ..., ωn ) is a mark
where ωk ∈ [xk−1 , xk ]. From each subinterval [xk−1 , xk ], we form a vertical rectangle with high f (ωk ) and width ∆xk .
Volumes of Revolution Solids 85
V = πr2 h .
Figure 5.21
From Figure 5.22, revolution of the rectangular k generates a circular disk with radius r = f (ωk ) and high h = ∆xk . Thus the volume of
the generated circular disk is
The sum of the volumes of the circular disks approximates the volume of the revolution solid:
n n 2 Z bh i2
V= ∑ Vk = kPk→0
lim ∑ π f (ωk ) ∆xk = π f (x) dx.
k=1 k=1 a
Figure 5.22: The volume of the revolution solid about x-axis by the disk method. The figure on the left shows the region R bounded by the function
y = f (x) ≥ 0 on the interval [a, b] and the figure on the right shows the solid S generated by revolving R about x-axis.
Similarly, we can find the volume of the revolution solid generated by revolving a region R about y-axis. Let f be continuous on [c, d]
and R be a region bounded by the graph of f and y-axis from y = c to y = d. Let S be a solid generated by revolving R about y-axis.
Assume that P is a partition of [c, d] and ω = (ω1 , ω2 , ..., ωn ) is a mark where ωk ∈ [yk−1 , yk ]. From each subinterval [yk−1 , yk ], we form
a horizontal rectangle, its high and width are f (ωk ) and ∆yk , respectively.
The revolution of each horizontal rectangle about y-axis generates a circular disk as shown in Figure 5.23 with radius r = f (ωk ) and high
h = ∆yk . Therefore, the volume of each circular disk is
Figure 5.23: The volume of the revolution solid about y-axis by the disk method. The figure on the left shows the region R bounded by the function f
on the interval [c, d] and the figure on the right shows the solid S generated by revolving R about y-axis.
The sum of the volumes of the circular disks approximates the volume of the revolution solid given in Figure 5.23 (right):
n n
V= lim ∑ π( f (ωk ))2 ∆yk
∑ Vk = kPk→0
k=1 k=1
Z dh i2
=π f (y) dy.
c
Theorem 5.2
1. If R is a region bounded by the graph of f on the interval [a, b], the volume of the revolution solid generated by revolving
R about x-axis is Z h bi 2
V =π f (x) dx.
a
2. If R is a region bounded by the graph of f on the interval [c, d], the volume of the revolution solid generated by revolving
R about y-axis is
Z dh i2
V =π f (y) dy.
c
Example 5.14 Sketch the region R bounded by the graph of y = x + 1 on the interval [0, 2]. Then, find the volume of the solid
generated by revolving R about x-axis.
Solution: First, we sketch the graph of the function y = x + 1 and determine the region R in the interval [0, 2]. Then, we sketch the solid
generated by revolving R about the x-axis.
Volumes of Revolution Solids 87
Figure 5.24
From the figure, we have a vertical disk with radius y = x + 1 and thickness dx. Thus, the volume of the solid S is as follows:
Z 2 i2 π
πh 26π
V =π (x + 1)2 dx = (x + 1)3 = (27 − 1) = .
0 3 0 3 3
√
Example 5.15 Sketch the region R bounded by the graph of y = x from x = 0 to x = 4. Then, find the volume of the solid generated
by revolving R about x-axis.
Solution: Figure 5.25 shows the region R and the solid S generated by revolving the region about the x-axis.
Figure 5.25
√
Since the revolution is about the x-axis, we have a vertical disk with radius y = x and thickness dx.
Thus, the volume of the solid S is as follows:
Z 4
√ Z 4
π h 2 i4 π h i
V =π ( x)2 dx = π x dx = x = 16 − 0 = 8π.
0 0 2 0 2
2
Example 5.16 Sketch the region R bounded by the graph of the function y = x and x-axis from x = −2 to x = 2. Then, find the
volume of the solid generated by revolving R about x-axis.
Solution: The figure on the left shows the region R bounded by the graph of y = x2 in the interval [−2, 2]. The figure to the right shows
the solid S generated by revolving the region about the x-axis.
88 APPLICATIONS OF INTEGRATION
Figure 5.26
From the figure, we have a vertical disk with radius y = x2 and thickness dx. Thus, the volume of the solid S is as follows:
Z 2
π h 5 i2 64π
V =π (x2 )2 dx = x = .
−2 5 −2 5
Example 5.17 Sketch the region R bounded by the graph of the equations x = y + 1, y = 1, y = 3. Then, find the volume of the solid
generated by revolving R about y-axis.
Solution: The figure shows the region R and the solid S generated by revolving the region about the y-axis.
Figure 5.27
From the figure, we have a vertical disk with radius x = y + 1 and thickness dy. Thus, the volume of the solid S is as follows:
Z 3 i3 56π
πh
V =π (y + 1)2 dy = (y + 1)3 = .
1 3 1 3
2
Example 5.18 Sketch the region R bounded by the graph of the equation x = y on the interval [0, 2]. Then, find the volume of the
solid generated by revolving R about y-axis.
Solution: The figure on the left shows the region R bounded by the graph of x = y2 in the interval [0, 2]. The figure to the right shows the
solid S generated by revolving the region about the y-axis.
Volumes of Revolution Solids 89
Figure 5.28
Since the revolution of R is about the y-axis, we have a horizontal disk with radius x = y2 and thickness dy. Thus, the volume of the solid
S is as follows: Z 2
π h 5 i2 32π
V = π (y2 )2 dy = y = .
0 5 0 5
√
Example 5.19 Sketch the region R bounded by the graph of the equation x = y − 1 from y = 1 to y = 3. Then, find the volume of
the solid generated by revolving R about y-axis.
Solution: The figure shows the region R and the solid S generated by revolving the region about the y-axis.
Figure 5.29
√
Since the revolution of R is about the y-axis, we have a horizontal disk with radius x = y − 1 and thickness dy. Thus, the volume of the
solid S is as follows:
Z 3 p h y2 i3
V = π ( y − 1)2 dy = π − y = 2π.
1 2 1
Z bh 2 2 i
V =π f (x) − g(x) dx.
a
Figure 5.30: The volume of the solid generated by revolving the region about the x-axis using the washer method.
Similarly, let R be a region bounded by the graphs of f (y) and g(y) such that f (y) ≥ g(y) for all y ∈ [c, d] as shown in Figure 5.31. The
volume of the solid S generated by revolving R about the y-axis is
Z dh 2 2 i
V =π f (y) − g(y) dy.
c
Figure 5.31: The volume of the solid generated by revolving the region about the y-axis using the washer method.
Theorem 5.3
1. If R is a region bounded by the graphs of f and g on the interval [a, b] such that f ≥ g, the volume of the revolution solid
generated by revolving R about x-axis is
Z bh 2 2 i
V =π f (x) − g(x) dx.
a
2. If R is a region bounded by the graphs of f and g on the interval [c, d] such that f ≥ g, the volume of the revolution solid
generated by revolving R about y-axis is
Z dh 2 2 i
V =π f (y) − g(y) dy.
c
Example 5.20 Let R be a region bounded by the graphs of the functions y = x2 and y = 2x. Evaluate the volume of the solid
generated by revolving R about x-axis.
Solution:
Let f (x) = x2 and g(x) = 2x. First, we check whether the graphs of the two functions are intersecting or not.
Note
f (x) = g(x) ⇒ x2 = 2x ⇒ x2 − 2x = 0 The graphs of two functions f (x)
⇒ x(x − 2) = 0 and g(x) intersect at x = x0 if
f (x0 ) = g(x0 ).
⇒ x = 0 or x = 2.
Substituting x = 0 into f (x) or g(x) gives y = 0. Similarly, substitute x = 2 into any of the two functions gives y = 2. Thus, the two
curves intersect in two points (0, 0) and (2, 4).
Figure 5.32
The figure shows the region R and the solid generated by revolving the region about the x-axis. A vertical rectangle generates a washer
where
the outer radius: y1 = 2x,
the inner radius: y2 = x2 and
the thickness: dx. h i
The volume of the washer is dV = π (2x)2 − (x2 )2 dx.
92 APPLICATIONS OF INTEGRATION
2
Example 5.21 Consider the same region as in Example 5.20 enclosed by the graphs of y = x and y = 2x. Revolve the region about
y-axis instead and find the volume of the generated solid.
Solution: The figure shows the region R and the solid generated by revolving the region about the y-axis.
Figure 5.33
Since the revolution is about the y-axis, we need to rewrite the equations in term of y i.e., x = f (y) and x = g(y).
√ y
y = x2 ⇒ x = y = f (y) and x = = g(y) .
2
The two horizontal rectangles generate a washer where
√
the outer radius: x1 = y,
the inner radius: x2 = 2y and
the thickness: dy. h√ i
The volume of the washer is dV = π ( y)2 − ( 2y )2 dy.
√
Example 5.22 Consider a region R bounded by the graphs of the functions y = x, y = 6 − x and x-axis. Revolve the region about
y-axis and find the volume of the generated solid.
Solution: Since the revolution is about the y-axis, we need to rewrite the functions in terms of y i.e., x = f (y) and x = g(y).
√
y = x ⇒ x = y2 = f (y) and y = 6 − x ⇒ x = 6 − y = g(y).
Volumes of Revolution Solids 93
√
Since y = x, we ignore the value y = −3. By substituting y = 2 into the two functions, we have x = 4. Thus, the two curves intersect in
one point (4, 2). The solid S generated by revolving the region R about y-axis is shown in Figure 5.34.
Also, the revolution is about the y-axis, so we have a horizontal rectangle that generates a washer where
the outer radius: x1 = 6 − y,
the inner radius: x2 = y2 and
the thickness: dy.
The volume of the washer is dV = π (6 − y)2 − (y2 )2 dy.
Z 2 h (6 − y)3 y5 i2
(6 − y)2 − (y2 )2 dy = π −
V =π −
0 3 5 0
h 64 32 216 i
=π − − − − −0
3 5 3
664
= π.
15
Figure 5.34
√
Example 5.23 Consider the same region as in Example 5.22 enclosed by the graphs of y = x, y = 6 − x and x-axis. Revolve the
region about x-axis instead and find the volume of the generated solid.
Solution: From the figure, we find that the solid is made up of two separate regions and each requires its own integral. Hence, we use the
disk method to evaluate the volume of the solid generated by revolving each region.
From revolution of the first region R1 about the x-axis, we have a vertical disk with radius y = x2 and thickness dx. Thus, the volume of
the solid S1 is as follows:
Z 4
√ Z 4
V1 = π ( x)2 dx = π x dx
0 0
π h 2 i4 32
= x = π.
2 0 3
From revolution of the first region R2 about the x-axis, we have a vertical disk with radius y = 6 − x and thickness dx. Thus, the volume
of the solid S2 is as follows:
94 APPLICATIONS OF INTEGRATION
Note
Z 6 Z 6 Use the substitution method to do
V2 = π (6 − x)2 dx = π (6 − x)2 dx the second integral with
4 4
πh i6 32 u = 6 − x and − du = dx
= − (6 − x)3 = π.
3 4 3
V = V1 +V2
8 32
= 8π + π = π.
3 3
Figure 5.35
Figure 5.36
Volumes of Revolution Solids 95
V = V2 −V1
Note
= πr22 h − πr12 h
= π(r22 − r12 )h V = V2 − V1
|{z}
= π(r2 + r1 )(r2 − r1 )h
|{z} |{z}
volume of the cylindrical shell volume of the outer cylinder volume of the inner cylinder
r2 + r1
= 2π( )h(r2 − r1 )
2
= 2πrh∆r.
Now, let R be a region R on the interval [a, b] and S be a solid generated by revolving the region about y-axis (Figure 5.37). Let P be a
partition of the interval [a, b] and let ω = (ω1 , ω2 , ..., ωn ) be a mark on P where ωk is the midpoint of [xk−1 , xk ].
The revolution of the rectangle about the y-axis generates a cylindrical shell where
the high = f (ωk ),
the average radius = ωk ,
the thickness = ∆xk .
A B
Figure 5.37: The volume of the revolution solid about the y-axis by the cylindrical shells method.
Hence, the volume of the cylindrical shell is Vk = 2πωk f (ωk )∆xk . To evaluate the volume of the whole solid, we sum the volumes of all
cylindrical shells. This implies
n n
V= ∑ Vk = 2π ∑ ωk f (ωk )∆xk .
k=1 k=1
This implies
Z b
V = 2π x f (x) dx.
a
Similarly, we can find that if the revolution of the region is about x-axis, the volume of the revolution solid is
Z d
V = 2π y f (y) dy.
c
96 APPLICATIONS OF INTEGRATION
Theorem 5.4
1. If R is a region bounded by the graph of f on the interval [a, b], the volume of the revolution solid generated by revolving
R about y-axis is
Z b
V = 2π x f (x) dx.
a
2. If R is a region bounded by the graph of f on the interval [a, b], the volume of the revolution solid generated by revolving
R about x-axis is Z d
V = 2π y f (y) dy.
c
Note: The importance of the cylindrical shells method appears when solving equations for one variable in terms of another (i.e., solving
x in terms of y). For example, let S be a solid generated by revolving the region bounded by y = 2x2 − x3 and y = 0 about y-axis. By the
washer method, we have to solve the cubic equation for x in terms of y, but this is not simple.
2
Example 5.24 Sketch the region R bounded by the graph of y = 2x − x and x-axis. Then, by the cylindrical shells method, find the
volume of the solid generated by revolving R about y-axis.
Solution:
The figure shows the region R and the solid S generated by revolving the region about the y-axis.
Figure 5.38
Since the revolution is about the y-axis, the rectangle is vertical and by revolving it, we obtain a cylindrical shell where
the high: y = 2x − x2 ,
the average radius: x,
the thickness: dx.
√
Example 5.25 Sketch the region R bounded by the graphs of the equations x = y and y = 4, and y-axis. Then, find the volume of
the solid generated by revolving R about x-axis.
Solution:
Figure 5.39
Since the revolution is about the x-axis, the rectangle is horizontal and by revolving it, we have a cylindrical shell where
√
the high: x = y,
the average radius: y,
the thickness: dy.
√
The volume of the cylindrical shell is dV = 2π y y dy.
Exercises
1 - 36 Sketch the region bounded by the graphs of the given equations, then find its area.
√
2 y = x3 and x-axis over [0, 2] 20 x = y and 2x = y from y = 0 to y = 4
√
3 y = x2 , y = 4 21 y = x + 1, y = x + 1
√
4 x = y3 and y-axis from y = 0 to y = 2 22 y = x2 , y = x
9 y = x, x = 2 − y, y = 0 27 y = e−x from x = −1 to x = 2
√
13 y = x2 , y = x 31 y = cos 2x and x-axis over [0, π4 ]
14 x = y2 , y = x + 1, y = 1, y = 2 32 y = sin x, x = π4 , x = π
2
−π
15 y − 1 = 3x and y − 2 = x from x = 0 to x = 1 33 y = sec2 x, y = 0, x = 4 , x= π
4
√
16 y = x + 1 and y = x − 1 over [1, 3] 34 y = tan x and x-axis from x = 0 to x = π
4
17 y = x3 , y = x2 35 y = x2 − 1, x = 1, x = 2
37 - 54 Sketch the region R bounded by the graphs of the given equations and find the volume of the solid generated by revolving R
about x-axis.
√
46 x = y, x = y
37 y = 2x and x-axis over [0, 1]
47 y = x3 , y = x2 , x = 1, x = 3/2
38 y = x, x + y = 4 and y-axis
48 y = 4x − x2 and x-axis
39 y = x2 , y = 4 − x2
1 49 y = ex over [0, 2]
40 y = x and x-axis over [1, 3]
√ 50 y = x2 , y = 9
41 y = x2 , y = x
51 y = x2 , y = x
42 y = x2 , y = 1 − x2
52 y = ln x over [1, 4]
43 y = x2 , y = x3
π
53 y = sin x and y = cos x from x = 0 to x = 4 (use
44 y = 1 + x3 , x = 1, x = 2, y = 0
cos2 x = 1+cos
2
2x
, sin2 x = 1−cos
2
2x
)
45 x = y + 1, x = 2y − 3, x = 1, x = 3 π
54 y = sin x from x = 0 to x = 2
55 - 72 Sketch the region R bounded by the graphs of the given equations and find the volume of the solid generated by revolving R
about y-axis.
56 x = y2 , x = 2y 65 y = x2 − 1, y = 0, x = 1, x = 2
58 y = x2 , y = 0 and x = 2 67 (x − 2)2 + y = 1, y = 0
59 x = y2 , y = x − 2 68 y = 1 − x2 , y = 1 − x
60 y = cos x, x = 0, x = π
2 69 y = x2 + 1 and x-axis over [0, 1]
π
61 y = cos x, y = sin x, x = 0, x = 4 70 y = 6 − 3x and y-axis
62 y2 = 1 − x, x = 0 71 y = 1 − x2 , x = 0, x = 1
73 The area of the region bounded by the graphs of the functions y = x2 and y = 2 − x2 is equal to
(a) 2 (b) 4 (c) 38 (d) 83
74 The area of the region bounded by the graphs of the functions x = −y2 and x = −1 is equal to
(a) 34 (b) 19 (c) 16 (d) 83
75 The area of the region bounded by the graphs of the functions y = x and y = −x and y = 1 is equal to
(a) 1 (b) 0 (c) 2 (d) 21
76 The area of the region bounded by the graphs of the functions y = 2x and y = x and 0 ≤ x ≤ 1 is equal to
(a) 21 (b) 14 (c) 2 (d) 13
π
77 The area of the region bounded by the graphs of the functions y = cos x, y = sin x, x = 0 and x = is equal to
√ √ √ 4
(a) 2 − 1 (b) 0 (c) 2 + 1 (d) 1 − 2
78 The area of the region bounded by the graphs of the functions x = y2 and x = 2 − y2 is equal to
(a) 31 (b) 8 (c) 1 (d) 83
101
Chapter 6
PARTIAL DERIVATIVES
Definition 6.1
1. A function of two variables is a rule that assigns an ordered pair (x1 , x2 ) to a real number w:
f : R2 −→ R
(x1 , x2 ) −→ w .
2. A function of three variables is a rule that assigns an ordered triple (x1 , x2 , x3 ) to a real number w:
f : R3 −→ R
(x1 , x2 , x3 ) −→ w .
Example 6.1
(1) f (x, y) = x2 + y2 is a function of two variables x and y. The function f (x, y) takes (x, y) ∈ R2 to ω ∈ R. For example, f (1, 2) =
12 + 22 = 5 i.e., the function f takes (1, 2) ∈ R2 to 5 ∈ R.
(2) f (x, y, z) = x2 + y2 + z is a function of three variables x, y and z. The function f (x, y, z) takes (x, y, z) ∈ R3 to ω ∈ R. For example,
the function f takes (1, 2, −1) ∈ R3 to 4 ∈ R.
Definition 6.2 A function of n variables is a rule that assigns an ordered n-tuple (x1 , x2 , ..., xn ) to a real number w:
f : Rn −→ R
(x1 , x2 , ..., xn ) −→ w .
Example 6.2 f (x, y, z, u, v) = x2 + y2 − 7zu + v2 is a function of five variables. The function f (x, y, z, u, v) takes (x, y, z, u, v) ∈ R5 to
ω ∈ R. For example, the function f takes (1, 0, 1, 1, 2) ∈ R5 to −2 ∈ R.
For one variable y = f (x), the derivative dy/dx gives the change rate of y with respect to x. A similar thing occurs with functions of
more than one variable. For example, for a function of two variables ω = f (x, y), the independent variables are x and y while ω is the
dependent variable i.e. as x and y vary the value of ω traces out a surface.
102 PARTIAL DERIVATIVES
2x
Example 6.4 If f (x, y) = y + sin(xy), calculate (1) fx and (2) fy .
Solution:
If u = g(x) is differentiable, then
2
(1) fx = y + y cos(xy) .
d
(sin(u)) = cos(u) u0
(2) fy = − 2x
y2
+ x cos(xy) . dx
∂2 f
• ∂y2
means the second derivative with respect to y holding x constant.
∂2 f
• ∂x∂y
means differentiate first with respect to y and then with respect to x.
∂2 f ∂ ∂f ∂
2. ∂y2
= ( )
∂y ∂y
= f
∂y y
= fyy .
∂2 f ∂ ∂f ∂
3. ∂x∂y
= ( )
∂x ∂y
= f
∂x y
= fyx .
∂2 f ∂ ∂f ∂
4. ∂y∂x
= ( )
∂y ∂x
= f
∂y x
= fxy .
The second partial derivatives of functions of three variables are defined in the same manner given in the previous definition.
Theorem 6.6
(1) Let f (x, y) be a function of two variables. If the second partial derivatives fxy and fyx exist and are continuous, then
fxy = fyx .
(2) Let f (x, y, z) be a function of three variables. If the partial derivatives fxy , fyx , fxz , and fzx exist and are continuous, then
fxy = fyx , fxz = fzx and fyz = fzy .
(3) fxy = 2y, fyz = 3y2 and fzx = 3z2 . At (0, −1, 1), we have fxy = −2, fyz = 3 and fzx = 3.
A chain rule for ordinary derivatives is to differentiate a function of a function (composite functions). If f (x) and g(x) are two
functions, then the composite function of the two functions is ( f ◦ g)(x) = f g(x) . For example, if f (x) = cos x and g(x) = x2 , then
( f ◦ g) = f g(x) = cos x2 . To differentiate such function, we apply the chain rule given in the following definition.
104 PARTIAL DERIVATIVES
Definition 6.7 If g is a differentiable function at x and f is differentiable at g(x), then the composite function F = f ◦ g
defined by F(x) = f g(x) is differentiable at x as follows:
dF d f dg(x)
= .
dx dg(x) dx
dy
Example 6.8 If y = cos x2 , calculate dx .
Solution:
Let f (x) = cos x and g(x) = x2 , then ( f ◦ g)(x) = f g(x) = cos x2 .
It follows that
df dg(x)
= − sin g(x) and = 2x .
dg(x) dx
By applying the chain rule, we have
dy d f dg(x)
=
dx dg(x) dx
= − sin g(x) (2x) = −2x sin x2 .
In the following, we expanded the chain rule for composite functions of two or three functions. Thus, we need to use the chain rule more
than once.
1. If w = f (x, y), x = g(t), and y = h(t) such that f , g and h are differentiable, then
df dw ∂w dx ∂w dy
= = + .
dt dt ∂x dt ∂y dt
2. If w = f (x, y), x = g(t, s), and y = h(t, s) such that f , g and h are differentiable, then
∂f ∂w ∂w ∂x ∂w ∂y
= = + .
∂t ∂t ∂x ∂t ∂y ∂t
∂f ∂w ∂w ∂x ∂w ∂y
= = + .
∂s ∂s ∂x ∂s ∂y ∂s
3. If w = f (x, y, z), x = g(t, s), y = h(t, s), and z = k(t, s) such that f , g, h and k are differentiable, then
∂f ∂w ∂w ∂x ∂w ∂y ∂w ∂z
= = + + .
∂t ∂t ∂x ∂t ∂y ∂t ∂z ∂t
∂f ∂w ∂w ∂x ∂w ∂y ∂w ∂z
= = + + .
∂s ∂s ∂x ∂s ∂y ∂s ∂z ∂s
Note that the previous result can be proven by repeating the chain rule.
∂f ∂f
Example 6.9 If f (x, y) = xy + y2 , x = s2 t, and y = s + t, calculate (1) ∂t
(2) ∂s
.
Solution:
(1) ∂∂tf = ∂ f ∂x
∂x ∂t
+ ∂∂yf ∂y
∂t
=y s2 + (x + 2y)(1)
= (s + t)s2 + s2 t + 2s + 2t
= s3 + 2s2 t + 2s + 2t .
∂f
(2) ∂s
= ∂∂xf ∂x
∂s
+ ∂∂yf ∂y
∂s
= y (2st) + (x + 2y)(1)
= (s + t)(2st) + s2 t + 2s + 2t
= 3s2 t + 2st 2 + 2s + 2t .
Implicit Differentiation 105
∂f ∂f
Example 6.10 If f (x, y, z) = x + sin(xy) + cos(xz), x = ts, y = s + t, and z = st , calculate (1) ∂t
(2) ∂s
.
Solution:
(1) ∂∂tf = ∂ f ∂x
∂x ∂t
+ ∂∂yf ∂y
∂t
+ ∂∂zf ∂z
∂t
= (1 + y cos(xy) − z sin(xz))s + x cos(xy)(1) − x sin(xz)( −st2
)
s s 2
= s + (s + t)s + ts cos(ts(s + t)) + ( t 2 )ts − ( t )s sin(s )
= s + (s2 + 2ts) cos(ts(s + t)) .
∂f ∂ f ∂x
(2) ∂s
= ∂x ∂s
+ ∂∂yf ∂y
∂s
+ ∂∂zf ∂z
∂s
= (1 + y cos(xy) − z sin(xz))t + x cos(xy)(1) − x sin(xz)( 1t )
= t + (s + t)s + ts cos(ts(s + t)) − ( st )t + ( st )t sin(s2 )
We know that F(x, y) = 0 defines y as a function of x, y = y(x). Now, differentiate both sides of F(x, y(x)) = 0 by using the chain rule.
This implies
∂F ∂F dy dy ∂F/∂x
(1) + =0⇒ =− .
∂x ∂y dx dx ∂F/∂y
Definition 6.8
1. Suppose that the equation F(x, y) = 0 defines y implicitly as a function of x, y = f (x) such that f is differentiable. Then,
dy Fx
=− .
dx Fy
2. Suppose that the equation F(x, y, z) = 0 defines z implicitly as a function of x and y, z = f (x, y) such that f is differentiable.
Then,
∂z Fx ∂z Fy
= − and =− .
∂x Fz ∂x Fz
dy
Example 6.11 Let y2 − xy + 3x2 = 0, find dx .
Solution:
Let F(x, y) = y2 − xy + 3x2 = 0, then Fx = −y + 6x and Fy = 2y − x. Hence,
dy Fx −y + 6x y − 6x
=− =− = .
dx Fy 2y − x 2y − x
106 PARTIAL DERIVATIVES
∂z Fy x2 + xz cos(xyz)
(2) =− =− .
dy Fz 2z + xy cos(xyz)
Exercises 107
Exercises
3 f (x, y) = 3x + 4y 16 f (x, y) = x2 + xy − y2
4 f (x, y) = xy3 + x2 y2 17 f (x, y) = ln(x2 − y)
5 f (x, y) = x3 y + ex 18 f (x, y) = x cos y + yex
6 f (x, y) = xe2x+3y 19 f (x, y) = y sin xy
x−y
7 f (x, y) = x+y 20 f (x, y) = 4x2 − 8xy4 + 7y3 − 3
∂f ∂f
42 - 57 Find ∂s
and ∂t
at the given point.
58 x3 − 3xy2 + y3 = 5 67 y sin y + x = 1
√ 68 x2 + y2 − 4 = 0
59 x − xy + 3y = 4
√
60 4x + 6y = 5 69 xy − y2 + 2x = 2
1
61 x2 + y2 = 1 70 y 2 − 2x2 + 5y = 1
62 4y + 2x = 8 71 x2 y3 + x = 2
63 y2 + x2 = 16 72 ln x + ln y = 4
64 5y3 + 4x5 = 20 73 sin−1 x − y = 0
65 x2 + y3 = 2
p
74 1 + x2 y2 = 2xy
66 y2 − x3 (2 − x) = 0 75 2x3 + x2 y + y3 = 1
109
Chapter 7
DIFFERENTIAL EQUATIONS
Definition 7.1 An equation that involves x, y, y0 , y00 , y000 , y(4) , ..., y(n) for a function y(x) with nth derivative of y with respect to
x is an ordinary differential equation of order n.
Example 7.1
(1) y0 = x2 + 5 is a differential equation of order 1.
Remark 7.2 y = y(x) is called a solution of a differential equation if it satisfies that differential equation.
Solution:
We want to see whether the function y satisfies the equation. By taking the derivative, we obtain y0 = 6x + 4 and this is the given
differential equation.
Note:
1. y = y(x) + c is the general solution of the differential equation.
2. If an initial condition was added to the differential equation to assign a certain value for c, then y = y(x) is called the particular solution
of the differential equation.
Example 7.3 Verify that y = 4x3 + 2x2 + x is a solution of the differential equation y0 = 12x2 + 4x + 1. Then, with the initial condition
y(0) = 2, find the particular solution of the equation.
Solution:
110 DIFFERENTIAL EQUATIONS
First, we want to check whether the function y satisfies the differential equation. By taking the derivative, we have y0 = 12x2 + 4x + 1
and this is the given differential equation. Hence, y = 4x3 + 2x2 + x + c is the general solution of the given differential equation. Now,
since y(0) = 2, then
y(0) = 403 + 202 + 0 + c = 2 ⇒ c = 2 .
Therefore, y = 4x3 + 2x2 + x + 2 is the particular solution of the differential equation y0 = 12x2 + 4x + 1.
y−1
⇒ = ex + c
−1
1
⇒y = − . solve for y
ex + c
dy
Example 7.5 Solve the differential equation dx = y2 ex , with y(0) = 1.
Solution: Write the differential equation in the form N(y)dy = −M(x)dx.
dy dy
= yx ⇒ = x dx
dx y
1
Z Z
⇒ dy = x dx integrate both sides
y
⇒ ln |y| = x2 + c
2
+c
⇒y = ex . solve for y
2
With y(0) = 1, we have 1 = ec . This implies c = ln(1) = 0. Hence, the particular solution is y = ex .
dy
Example 7.7 Solve the differential equation dx − 12 y = 32 , with y(0) = 4.
Solution: Manipulate the differential equation to become N(y)dy = −M(x)dx.
dy 1 3 dy dy
− y = ⇒ 2 −y = 3 ⇒ = 2 dx
dx 2 2 dx 3+y
1
Z Z
⇒ dy = 2 dx integrate both sides
3+y
⇒ ln |3 + y| = 2x + c
⇒y = e2x+c − 3 . solve for y
With y(0) = 4, we have 4 = ec − 3. Hence, ec = 7 and this implies c = ln(7). Therefore, the particular solution is y = 7e2x − 3.
cos2 x dy sin x
e−y sin x − y0 cos2 x = 0 ⇒ e−y − = 0 ⇒ey = 2x
dx
sin x dx Z
cos Z
⇒ ey dy = tan x sec x dx integrate both sides
⇒ey = sec x + c
⇒y = ln | sec x + c| . solve for y by taking ln for both sides
y0 + P(x)y = Q(x) ,
dy
Example 7.10 Solve the differential equation x dx + y = x2 + 1.
Solution: Write the differential equation in the form y0 + P(x)y = Q(x).
dy 1 x2 + 1
x + y = x2 + 1 ⇒ y0 + y = .
dx x x
2
R 1
From this, we have P(x) = 1x and Q(x) = x x+1 . Hence, the integrating factor is µ(x) = e x dx = eln |x| = x.
The general solution of the first-order linear differential equation is
1 x2 + 1
Z
y(x) = x dx
x x
1
Z
= (x2 + 1) dx
x
1 x3
= ( + x) + c
x 3
x2 c
= +1+ .
3 x
= x2 ex + c .
With y(1) = e, we have e = 1 + c and this implies c = e − 1. The particular solution is y = x2 (ex + e − 1).
Solution: The
R differential equation takes the form y0 + P(x)y = Q(x) where P(x) = 1 and Q(x) = cos(ex ). Hence, the integrating factor
is µ(x) = e 1 dx = ex .
The general solution of the first-order linear differential equation is
Z
y(x) = e−x ex cos(ex ) dx Use integration by substitution
with u = ex and du = ex dx
= e−x sin(ex ) + c .
Exercises
4 y0 = 1 + y 12 y0 = y cos x
6 dy
dx = (1 + y2 ) sin x 14 y0 − y = e−x
7 y0 + y = e2x 15 y0 + 1+x
2x
2y=1
8 xy0 − y = x3 ex 16 dy
+ y − ex 1+1 = 0
dx
17 - 40 Solve the differential equation with the given initial condition.
1+3x2
17 y0 + 2y = x, y(0) = 1 29 y0 = 3y2
, y(0) = 1
√
18 dy
dx + 2y = e−x , y(0) = 3
4
30 cos xy0 + sin xy = 2 cos3 x sin x, y( π4 ) = 3 2, 0 ≤ x < π
2
dy dy 1
19 dx − 2xy = x, y(0) = 0 31 dx = 5y2 x, y(1) = 25
21 xy0 + y = sin x, y( π3 ) = 2 33 y0 =
2
√xy , y(0) = −1
1+x2
22 y0 − 13 y = e−x , y(0) = a 34 y0 = e−y (2x − 4), y(5) = 0
23 xy0 + 2y = 4x2 , y(1) = 2 dy y2
35 dx = x, y(1) = 2
dy 3x2 +4x+2
24 = , y(0) = −1
dx 2y 36 y0 = ex−sin y sec(y), y(0) = 0
25 y0 = 3x2 + 3x2 y, y(0) = 0 37 xy0 + y = x3 , y(1) = −3
26 xy0 + y = x3 , y(−1) = 3 38 y0 = 1+2x
y(0) = 0
tan y ,
27 y0 − 1x y = x, y(1) = 2
39 y0 = 2x2 + 2x2 y2 , y(0) = 0
28 y0 − y2 = 0, y(0) = 1 40 xy0 + 2y = x2 − x + 1, y(1) = 1
2
115
Chapter 8
Definition 8.1 The polar coordinate system is a two-dimensional system consisted of a pole and polar axis (half line). Each
point P in a polar plane is determined by a distance r from a fixed point O called the pole (or origin) and an angle θ from a fixed
direction.
Figure 8.1: The Cartesian coordinate system (on the left) and the polar coordinate system (on the right).
Note:
1. The point P in the polar coordinate system is represented by the ordered pair (r, θ) where r and θ are called polar coordinates.
2. The angle θ takes positive numbers if it is measured counterclockwise from the polar axis, but if the angle is measured clockwise, it
takes negative numbers.
3. In the polar coordinate system, if r > 0, the point P(r, θ) will be in the same quadrant as θ. However, if r < 0, the point will be in the
quadrant on the opposite side of the pole. That is, the points P(r, θ) and P(−r, θ) lie in the same line through the pole O, but on opposite
sides.
4. In the Cartesian coordinate system, every point has only one representation while in the polar coordinate system, each point has many
representations. The following formula gives all representations of the point P(r, θ) in the polar coordinate system
116 Polar Coordinates System
Figure 8.2
Example 8.1 Plot the points whose polar coordinates are given.
(1) (2, 5π/4)
(2) (2, −3π/4)
(3) (2, 13π/4)
(4) (−2, π/4)
The Relationship between Cartesian and Polar Coordinates 117
Solution:
(1) (2)
(3) (4)
Figure 8.3
Let (x, y) be a Cartesian coordinate and (r, θ) be a polar coordinate of the same point P. Let the pole be at the origin of the Cartesian
coordinates system, and let the polar axis lies on the positive x-axis and the line θ = π2 lies on the positive y-axis as shown in Figure 8.4.
x = r cos θ, y = r sin θ
y
tan θ = for x 6= 0
x
x2 + y2 = r2
(2) From the polar point (2, π), we have r = 2 and θ = π. Hence,
x = r cos θ = 2 cos π = −2 ,
y = r sin θ = 2 sin π = 0.
Hence, the polar point (2, π) is (−2, 0) in the Cartesian coordinates.
−2π
(3) From the polar point (2, −2π/3), we have r = 2 and θ = 3 . Hence,
−2π
x = r cos θ = 2 cos = −1 ,
3
−2π √
y = r sin θ = 2 sin = − 3.
3
√
Therefore, the Cartesian coordinate (−1, − 3) is the point corresponding to the polar point (2, −2π/3).
3π
(4) From the polar point (4, 3π/4), we have r = 4 and θ = 4 . Hence,
3π √
x = r cos θ = 4 cos = −2 2 ,
4
3π √
y = r sin θ = 4 sin = 2 2.
4
√ √
In the Cartesian coordinates, the point (4, 3π/4) is represented by (−2 2, 2 2).
Example 8.3 For the given Cartesian point, find one representation in the polar coordinates.
(1) (1, −1) (3) (−2, 2)
√
(2) (2 3, −2) (4) (1, 1)
Solution:
(1) From the given Cartesian point, we have x = 1 and y = −1. Hence,
√
x2 + y2 = r2 ⇒ r = 2,
y π
tan θ = = −1 ⇒ θ = − .
x 4
√
In the polar coordinates, the Cartesian point (1, −1) can be represented by ( 2, − π4 ).
Remember, there are infinitely polar representations of the point (x, y) (see Note 4 on page 115).
The Relationship between Cartesian and Polar Coordinates 119
√
(2) From the Cartesian point, we have x = 2 3 and y = −2. Hence,
x2 + y2 = r2 ⇒ r = 4,
y −1 5π
tan θ = = √ ⇒θ= .
x 3 6
√
Therefore, the polar point (4, 5π
6 ) is one representation of the Cartesian point (2 3, −2).
In the Cartesian coordinates, the function y = f (x) is a dependent relation which can be represented by a curve in the Cartesian plane. In
polar coordinates, the function r = f (θ) is a dependent relation between coordinates r and θ which also can be represented by a curve
called a polar curve. For example, r = 2 cos θ is a polar equation represents the dependent relation between coordinates r and θ:
π π π π
θ 0 √
6 4 √ 3 2
r 2 3 2/ 2 1 0
Table 8.1
Example 8.4 Find a polar equation that has the same graph as the equation in x and y.
(1) x = −5 (3) x2 + y2 = 2
(2) y = 3 (4) y2 = 9x
Solution:
(1) x = 7 ⇒ r cos θ = −5 ⇒ r = −5 sec θ.
Example 8.5 Find an equation in x and y that has the same graph as the polar equation.
(1) r = 4 (3) r = 6 cos θ
(2) r = 3 sin θ ⇒ r = 3 yr ⇒ r2 = 3y ⇒ x2 + y2 = 3y ⇒ x2 + y2 − 3y = 0.
120 Polar Coordinates System
x
(3) r = 6 cos θ ⇒ r = 6 r ⇒ r2 = 6x ⇒ x2 + y2 − 6x = 0.
1
(4) r = sec θ ⇒ r = cos θ ⇒ r cos θ = 1 ⇒ x = 1.
Theorem 8.2
1. Symmetry about the polar axis.
The graph of r = f (θ) is symmetric with respect to the polar axis if replacing (r, θ) with (r, −θ) or with (−r, π − θ) does
not change the polar equation.
A B C
Figure 8.5: Symmetry of the curves in the polar coordinates system. (A) symmetry about the polar axis, (B) symmetry about the vertical line θ = π2 ,
and (C) symmetry about the pole θ = 0.
Example 8.6
(1) The graph of r = 4 cos θ is symmetric about the polar axis. By replacing (r, θ) with (r, −θ), we have
4 cos (−θ) = 4 cos θ = r, thus (r, θ) = (r, −θ) .
Also, by replacing (r, θ) with (−r, π − θ), we have
−4 cos(π − θ) = 4 cos θ = r, thus (r, θ) = (−r, π − θ) .
(2) The graph of r = 2 sin θ is symmetric about the vertical line θ = π2 . By replacing (r, θ) with (r, π − θ), we obtain
2 sin (π − θ) = 2 sin θ = r, so (r, θ) = (r, π − θ) .
Also, by replacing (r, θ) with (−r, −θ), we have
−2 sin (−θ) = 2 sin θ = r, so (r, θ) = (−r, −θ) .
Polar Curves 121
(3) The graph of r2 = a2 sin 2θ is symmetric about the pole. If we replace (r, θ) with (−r, θ), we have
(−r)2 = a2 sin 2θ and this implies r2 = a2 sin 2θ, thus (r, θ) = (−r, θ) .
r2 = a2 sin 2(π + θ) = a2 sin (2π + 2θ) = a2 sin 2θ, thus (r, θ) = (r, θ + π) .
c
ax + by = c ⇒ r(a cos θ + b sin θ) = c ⇒ r =
(a cos θ + b sin θ)
4. The polar equation of a line that passes the origin point and makes an angle θ0 with the positive x-axis is θ = θ0 .
Example 8.7 Sketch the graph of θ = π4 .
Solution:
{(r, θ) |, r ∈ R}
π π π π π π
θ 4 4 4 4 4 4
r −6 −3 −1 1 3 6
Table 8.2
Figure 8.6
π π π π
θ 0 6 4 √ 3√ 2
r 0 2 4/ 2 2 3 4
Table 8.3
Cardioid curves
1. r = a(1 ± cos θ) 2. r = a(1 ± sin θ)
r = a(1 + cos θ) r = a(1 − cos θ)
Polar Curves 123
π π 2π
θ 0 3 2 3 π
r 0 a/2 a 3a/2 2a
Table 8.4
Limaçons curves
1. r = a ± b cos θ 2. r = a ± b sin θ
r = a + b cos θ
124 Polar Coordinates System
r = a − b cos θ
r = a + b sin θ
Polar Curves 125
r = a − b sin θ
r = a sin (nθ)
126 Polar Coordinates System
Spiral of Archimedes
r=aθ
To find the area of R, let P = {θ1 , θ2 , ..., θn } be a regular partition of the interval [α, β]. Consider the interval [θi−1 , θi ] where
Area in Polar Coordinates 127
∆θi = θi − θi−1 . By choosing ωi ∈ [θi−1 , θi ], we have a circular sector where its angle and radius are ∆θi and f (ωi ), respectively. The
area between θi−1 and θi can be approximated by the area of a circular sector.
Figure 8.16
By summing from i = 1 to i = n, we obtain
n n n
1 2 1 2
∑2 f (ui ) ∆θi f (ui ) ≤ ∑ ∆Ai ≤ ∑ 2 f (vi ) ∆θi f (vi )
i=1 i=1 i=1
| {z }
=A
n n
1 1 β 1
2 2 Z 2
lim ∑2 f (ui ) ∆θi f (ui ) = lim ∑2 f (ui ) ∆θi f (vi ) = f (θ) dθ .
||P||→0 i=1 ||P||→0 i=1 α 2
Therefore,
1
Z β 2
A= f (θ) dθ
2 α
Similarly, assume f and g are continuous on the interval [α, β] such that f (θ) ≥ g(θ). The area of the region bounded by the polar graphs
of f and g on the interval [α, β] is
1
Z βh 2 2 i
A= f (θ) − g(θ) dθ
2 α
Example 8.10 Find the area of the region bounded by the graph of the polar equation.
(1) r = 3
(2) r = 2 cos θ
(3) r = 4 sin θ
(4) r = 6 − 6 sin θ
Solution:
128 Polar Coordinates System
Note that we can evaluate the area in the first quadrant and multiply
the result by 4 to find the area of the whole region i.e.
1 Z π2 Z π
2
h iπ
2
A=4 32 dθ = 2 9 dθ = 18 θ = 9π.
2 0 0 0
Figure 8.17
(2) We find the area of the upper half circle and multiply the result by 2
as follows:
1 Z π2 Z π2
A=2 (2 cos θ)2 dθ = 4 cos2 θ dθ
2 0 0
Z π
2
=2 (1 + cos 2θ) dθ
0
h sin 2θ i π2
= 2 θ+
2 0
hπ i
= 2 −0
2
= π.
Figure 8.18
= 4 π−0
= 4π.
Figure 8.19
Area in Polar Coordinates 129
= 18 (2π + 2 + π) − 2
= 54π.
Figure 8.20
Example 8.11 Find the area of the region that is between the curves r = 2 and r = 3 in the first quadrant.
Solution: The region bounded by the two curves r1 = 2 and r2 = 3 is displayed in the figure.
π
1 2 2
Z
A= (r − r12 ) dθ
2 0 2
Z π
1 2
= 5 dθ
2 0
5 h i π2
= θ
2 0
5 π
h i
= −0
2 2
5π
= .
4
Figure 8.21
√
Example 8.12 Find the area of the region that is inside the graphs of the equations r = sin θ and r = 3 cos θ.
130 Polar Coordinates System
Solution:
First, we find the intersection points of the two curves
√ √ π
sin θ = 3 cos θ ⇒ tan θ = 3 ⇒ θ = .
3
The origin O is in each circle, but it cannot be found by solving the
equations. Therefore, when looking for the intersection points of
the polar graphs, we sometimes take under consideration the graphs.
The region is divided into two small regions: below and above the
line π3 .
Figure 8.22
Figure 8.23
1
Z π
2 √ Z π
3 2
A2 = ( 3 cos θ)2 dθ = (1 + cos 2θ) dθ
2 π
3
4 π3
3h sin 2θ i π2
= θ+
4 2 π
3
√ i
3h π π 3
= −0 − +
4 2 3 4
√ i
3hπ 3
= − .
4 6 4
Figure 8.24
√
5π 3
Total area A = A1 + A2 = 24 − 4 .
Area in Polar Coordinates 131
Example 8.13 Find the area of the region that is outside the graph of r = 3 and inside the graph of r = 2 + 2 cos θ.
Solution: The intersection point of the two curves in the first quadrant is
1 π
2 + 2 cos θ = 3 ⇒ cos θ = ⇒θ= .
2 3
As shown in the figure, we find the area in the first quadrant, then we double the result to find the area of the whole region.
1 Z π3
A=2 4(1 + cos θ)2 − 9 dθ
2 0
Z π
3
4(1 + 2 cos θ + cos2 θ) − 9 dθ
=
0
Z π
3
= (8 cos θ + 4 cos2 θ − 5) dθ
0
h iπ
3
= 8 sin θ + sin 2θ − 3θ
0
9√
= 3 − π.
2
Figure 8.25
Exercises 132
Exercises
1-8 Find the corresponding Cartesian coordinates for the given polar coordinates.
1 (2, π3 ) 5 ( 12 , 3π
2 )
3 (−2, π6 ) 7 (−7, 3π
4 )
4 (3, π) 8 (3, π3 )
9 - 16 Find the corresponding polar coordinates for the given Cartesian coordinates for r ≥ 0 and 0 ≤ θ ≤ π.
√
9 (1, 1) 13 ( 3, 1)
√
10 (1, 3) 14 ( 21 , 21 )
√
11 (−1, 1) 15 (−1, 3)
√
12 ( 3, 3) 16 (3, 0)
17 - 24 Find a polar equation that has the same graph as the equation in x and y and vice versa.
17 x = 4 21 x2 = 2y
18 x2 + y2 = 5 22 x2 − y2 = 9x
19 r = csc θ 3
23 r = 1−sin θ
20 r = 6 cos θ 24 r = 2 − 3 sin θ
25 - 28 Sketch the curve of the polar equations.
25 r = 3 sec θ 27 r = 2 + 2 sin θ
26 r = 4 cos θ 28 r = 3 + 2 cos θ
29 - 34 Find the area of the region bounded by the graph of the polar equation.
29 r = 3 sin θ 32 r = 4 cos θ
31 r = 2 34 r = 2(1 − sin θ)
35 - 41 Find the area of the region bounded by the graph of the polar equations.
35 inside r = 2
36 between r = 2 and r = 3
Integration Rules: Z
Z Z Z
f 0 g(x) g0 (x) dx = f g(x) + c
f (x) ± g(x) dx = f (x) dx ± g(x) dx
Z Z Z b
k f (x) dx = k f (x) dx f 0 (x) dx = f (b) − f (a)
a
Elementary Integrals: Z
xr+1 sec x tan x dx = sec x
Z
xr dx = if r 6= −1
r+1 Z
csc x cot x dx = − csc x
Z
sin x dx = cos x
1 x
Z
dx = sin−1
Z
cos x dx = − sin x
p
a2 − x2 a
1 1 x
Z Z
sec2 x dx = tan x dx = tan−1
a2 + x 2 a a
1 1 x
Z Z
csc2 x dx = − cot x p dx = sec−1 | |
x x 2 − a2 a a
Z p
sin−1 x dx = x sin−1 x + 1 − x2 + c
1
Z
tan−1 x dx = x tan−1 x − ln(1 + x2 ) + c
2
Z p
sec−1 x dx = x sec−1 x − ln | x + x2 − 1 | +c
xn+1 1 xn+1
Z Z
xn sin−1 x dx = sin−1 x − p dx + c if n 6= −1
n+1 n+1 1 − x2
xn+1 1 xn+1
Z Z
xn tan−1 x dx = tan−1 x − dx + c if n 6= −1
n+1 n+1 1 + x2
xn+1 1 xn
Z Z
xn sec−1 x dx = sec−1 x − p dx + c if n 6= −1
n+1 n+1 x2 − 1
Appendix 134
Trigonometric Integrals:
x sin 2x
Z
sin2 x dx = − +c
2 4
x sin 2x
Z
cos2 x dx = + +c
2 4
Z
tan2 x dx = tan x − x + c
Z
cot2 x dx = − cot x − x + c
1 1
Z
sec3 x dx = sec x tan x + ln | sec x + tan x | +c
2 2
1 1
Z
sec3 x dx = csc x cot x + ln | csc x − cot x | +c
2 2
1 n−1
Z Z
sinn x dx = − sinn−1 x cos x + sinn−2 x dx + c
n n
1 n−1
Z Z
cosn x dx = cosn−1 x sin x + cosn−2 x dx + c
n n
tann−1 x
Z Z
tann x dx = − tann−2 x dx + c if n 6= 1
n−1
cotn−1 x
Z Z
cotn x dx = − − cotn−2 x dx + c if n 6= 1
n−1
1 n−2
Z Z
secn x dx = secn−2 x tan x + secn−2 x dx + c if n 6= 1
n−1 n−1
1 n−2
Z Z
cscn x dx = − cscn−2 x cot x + cscn−2 x dx + c if n 6= 1
n−1 n−1
sinn−1 x cosm+1 x n−1
Z Z
sinn x cosm x dx = − + sinn−2 x cosm x dx + c if n 6= m
n+m n+m
sinn+1 x cosm−1 x m − 1
Z Z
sinn x cosm x dx = + sinn x cosm−2 x dx + c if m 6= n
n+m n+m
Z Z
xn sin x dx = −xn cos x + n xn−1 cos x dx + c
Z Z
xn cos x dx = xn sin x − n xn−1 sin x dx + c
Miscellaneous Integrals:
x b
Z
x(ax + b)−1 dx = − ln |ax + b| + c
a a2
1 b
Z
x(ax + b)−2 dx = 2 ln | ax + b | + +c
a ax + b
(ax + b)n+1 ax + b b
Z
x(ax + b)n dx = − +c
a2 n+2 n−1
a 1 x 1
Z Z
dx = 2 + (2n − 3) dx if n 6= −1
(a2 ± x2 )n 2a (n − 1) (a2 ± x2 )n−1 (a2 ± x2 )n−1
Z √ 2
x ax + b dx = (3ax − 2b)(ax + b)3/2 + c
15a2
Z √ 2
Z √
xn ax + b dx = xn (ax + b)3/2 − nb xn−1 ax + b dx
a(2n + 3)
Z
x 2 √
√ dx = 2 (ax − 2b) ax + b + c
ax + b 3a
Z
xn 2 √ Z
xn−1
xn ax + b − nb √
√ dx = dx
ax + b a(2n + 1) ax + b
√ √
1 1 ax + b − b
Z
√ dx = √ ln | √ √ | +c if b > 0
x ax + b b ax + b + b
r
1 1 ax + b
Z
√ dx = √ tan−1 + c if b < 0
x ax + b −b −b
√
1 ax + b (2n − 3)a 1
Z Z
√ dx = − − √ dx if n 6= 1
xn ax + b b(n − 1)xn−1 2(n − 1)b xn−1 ax + b
x−a p a2 a−x
Z p
2ax − x2 dx = 2ax − x2 + cos−1 ( )+c
2 2 a
2x2 − ax − 3a3 p a3 a−x
Z p
x 2ax − x2 dx = 2ax − x2 + cos−1 ( )+c
6 2 a
Z p
2ax − x2 p a−x
dx = 2ax − x2 + a cos−1 ( )+c
x a
p p
2ax − x2 2 2ax − x2 a−x
Z
dx = − − cos−1 ( )+c
x2 x a
dx a−x
Z
p = cos−1 ( )+c
2ax − x2 a
Appendix 135
x a−x
Z p
p dx = − 2ax − x2 + a cos−1 ( )+c
2ax − x2 a
79
0
5
15 X = −10 20 X = 3
96 1
5
Chapter 2:
− 16
6
3 0 173
1 1 11 −16 5 16 X = 4 21 X = 3
13 3 −3 2 10
15 14 3
−4 −13 7
−1
4 0 14 Not possible 2
17 X = − 27 22 X = 13
2 2 29
4
−5 3 1
35 31 7 − 54
15 46 27 16
65 39 23 3
19
−7 0 10
18 X = 45 10
23 X = 75
3 −1 3 7 1 1
7
−5 −8 16 57 36 21 10 − 21
10
37 22 13 7
4 Not possible 2 6
17 Not possible
19 X = −1 24 X = 8
28 34 −2 3 1 1
10
18 2
5 81 50 1 6 4
29 58
8 2 4
19
6 Not possible −2 10 14
Chapter 4:
4 −1
1 5 0 1
20 1 5 1 tan(3x − 5) + c
7 3 −4 9 3
2 7
2 6 2
21 −21 2 sin−1 4x + c
3 15 0 22 −6 1 2x 1 2x
8 9 −12 27 23 −3
3 2 xe − 4 e + c
6 18 6
24 −6 4 x sin 2x + 41 cos 2x + c
2
9 −2 25 11 p
10 −16 26 −43 5 x sin−1 x + 1 − x2 + c
27 29
− 31 ln |x + 1| + 13 ln |x − 2| + c
2 0 6
11 4 11 28 638
3 11 29 5 1 2 9
7 36 (2x − 3) + c
30 −104 √
20 −5 8 3 sin 3 x + c
31 73
12 5 25
10 35 32 12 21
9 2
10 8
3
11 1
Chapter 3:
12 1
5 6 √
1 X = 6 8 X =8 13 1 + √1 − 2
3
7 10
14 12
0 5
15 5
2 X = 3 9 X = 6 2
1 7
16 275
4 6
3
3
3 X = −3 10 17 0
10 X = 54
− 11 7
3 10 18 − 11
20
8
7
7
1
√ √
4 X = − 97
3
5
19 3 5 5−2 2
11 X = 11
− 97 29 17
33 20 2
1
3 13 1 x2 ln |x| − 1 x2 + c
5 X = 3 543 21 4 8
7 12 X = − 10
3 29 32
10 22 3
3
10 2 23 ln(3) − ln(e − 2)
6 X = 45 13 X = 0
7
10 1 24 2 ln(7) − ln(4)
9 25 17
−4 25 6
2
7 X = 21 14 X = 2 √
5 − 11 26 40 10 + 10
4 2
√
27 4( 2 − 1)
28 1
138
29 10
3
30 √1
3
9 1
44 373 π
31 9 3 14
2 10 2
45 10π
2 11 32
32 3 46 2
3 15 π
12 12
33 2 47 4387
1 4480 π
13 3
34 1 512 π
2 11 48 15
14 6
35 −2 (e4 −1)
15 1 49 2 π
2
36 1
sin x2 + c √ 1944 π
2 23−8 2 50 5
16 6
√
37 −2 cot x + c 51 2
17 1 15 π
6
38 tan x + sec x + c 4 52 16 ln2 (2) − 16 ln(2) +
18 3
3 3 6 π
39 2x2 ln |x| − 4x2 +c 9
3 9 19 2 53 π
2
40 x tan x + ln | cos x| + c 20 4
3 54 π2
4
1 −4x + c 21 1
41 − 4x − 16 e 6 55 3π
22 1 64
56 15 π
42 x(ln x)2 − 2x ln |x| + 2x + c 3
23 1 3
3 57 5π
43 − 14 ln |x − 1| + 2(x−1)
1 + 41 ln |x − 3| + c √
24 2 2 − 2 58 8π
1 72 π
44 2 ln |x + 2| − 12 ln |x + 4| + c 25 3 59 5
26 3 60 (π − 2)π
45 1 ln |x + 1| + 23 ln |x − 2| + c 2
3 √π
e3 −1
61 −2 π
27 2
e2
46 − 14 ln | x+1 x+1
2 + 1| − ln | 2 − 1| + c √ 62 16 π
3−1 15
28 2 10 π
47 4 + 23 ln(50) − ln(2) 29 2 ln(2) − 1
63 3
64 6π
2
−10 + 6(x − 2 + 2 ln |x − 2|) + c 30 1
48 − 3xx−2 65 9π
31 1 2
2
(x−1)2 √ 66 4
49 + 2(x − 1) − 8 ln |x − 1| + c 2−1
√ 21 π
2 32
2 16 π
67 3
50 x2 − x − 32 ln |x + 1| + x+1
1 + 1 ln |x − 3| + c
2 33 2
√ 68 π
6
34 ln( 2)
51 x − ln(ex + 1) + c 3π
4 69 2
35 3
52 53 1 ln |x − 2| + 55 +
ln |x + 1| + 27 1 +c
27 9(x−2) 3(x−2)2
70 8π
36 e2 + 1 π
4
71 2
53 −2 37 3π
72 4π
38 16π 3
54 4
√ 73 (d)
39 64 2 π
55 7 3
74 (a)
40 2π
56 −1 3 75 (a)
41 3
57 65 10 π 76 (a)
3 √
42 2 2π 77 (a)
58 0 3
2 78 (d)
65 (a) 43 35 π
66 (c)
67 (d) Chapter 6:
68 (a) 1 fx = 8x3 y3 − y2
69 (d) fy = 6x4 y2 − 2xy + 3
fxx = 24x2 y3
70 (c)
fyy = 12x4 y − 2x
71 (b)
2 3
2 fx = 8xy3 ex y
2 3
fy = 12x2 y2 ex y
Chapter 5: 2 3 2 3
fxx = 8y3 ex y + 16x2 y6 ex y
9 5 2
1 3 2 3 2 3
2 fyy = 24x2 yex y + +36x4 y4 ex y
2 4 6 44
3 32 7 1
3
8 11
4 4 4
139
3 fx = 3 16 fx = 2x + y
fy = 4 fy = x − 2y
fxx = 0 fxx = 2
fyy = 0 fyy = −2
4 fx = y3 + 2xy2 17 fx = 22x
x −y
fy = 3xy2 + 2x2 y
fy = − 21
x −y
fxx = 2y2
−2(x2 +y)
fyy = 6xy + 2x2 fxx =
(x2 −y)2
fyy = −1
5 fx = 3x2 y + ex (x2 −y)2
f y = x3
fxx = 6xy + ex 18 fx = cos y + yex
fyy = 0 fy = −x sin y + ex
fxx = yex
6 fx = e2x+3y + 2xe2x+3y fyy = −x cos y
fy = 3xe2x+3y
19 fx = y2 cos(xy)
fxx = 4e2x+3y + 4xe2x+3y
fy = sin(xy) + xy cos(xy)
fyy = 9xe2x+3y
fxx = −y3 sin(xy)
2y
7 fx = fyy = 2x cos(xy) − x2 y sin(xy)
(x+y)2
fy = −2x 2 20 fx = 8x − 8y4
(x+y)
−4y fy = −32xy3 + 21y2
fxx =
(x+y)3 fxx = 8
4x
fyy =
(x+y)3 fyy = −96xy2 + 42y
21 fx = y cos(xy)
8 fx = 2 sin(x2 y) + 4x2 y cos(x2 y)
fy = x cos(xy)
fy = 2x3 cos(x2 y)
fxx = −y2 sin(xy)
fxx = 12xy cos(x2 y) − 8x3 y2 sin(x2 y)
fyy = −x2 sin(xy)
fyy = −2x5 sin(x2 y)
22 fx = 3x2 + 6xy + 4
9 fx = 2x sin y − y2 sin x
fy = 3x2 + 2y
fy = x2 cos y + 2y cos x
fxx = 6x + 6y
fxx = 2 sin y − y2 cos x
fyy = 2
fyy = x2 sin y + 2 cos x
23 fx = 2xy + 4y3
10 fx = 3x2 + y2
fy = x2 + 12xy2
fy = 2xy + 1
fxx = 2y
fxx = 6x
fyy = 24xy
fyy = 2x
24 fx = 2x tan y
11 fx = 2xy2 + y2 fy = x2 sec2 y + 2y
fy = 2x2 y + 2xy fxx = 2 tan y
fxx = 2y2 fyy = 2x2 sec2 y tan y + 2
fyy = 2x(x + 1)
25 fx = 3x2 ln y + y4
12 fx = 3x2 + 1 3
fy = xy + 4xy3
fy = 4y + 1
fxx = 6x ln y
fxx = 6x 3
fyy = 4 fyy = − x2 + 12xy2
y
13 fx = 3yx2 + y4 − 3 26 fx = 3x2 y + y3
fy = x3 + 4xy3 − 3 fy = x3 − 3xy2
fxx = 6yx fxx = 6xy
fyy = 12xy2 fyy = −6xy
14 fx = − y2 ln x + y2 27 6, 2, 3, 1
x x
fy = lnxx 28 0, −1, 0, 0
fxx = 2y3 ln x − 3y3
x x 29 3, − 14 , − 21 , − 14
4
fyy = 0
30 1, 0, −6, 6
15 fx = −2x
(x2 +y2 )2
31 1, 0, 1, 0
−2y
fy =
(x2 +y2 )2
32 0, 3, −4, 0
2 2
fxx = 6x2 −2y
2 3
(x +y )
33 −6, 3, 0, −18
2 2
fyy = 6y2 −2x
2 3
(x +y ) 34 0, 0, 0, 0
140
Chapter 7:
35 0, 0, 1, 0
x
1 y = − 1+cx
36 0, 0, 0, 4
2 y = − tan1x+c
37 12, 1, 0, 0
38 3, 0, 0, 2 3 y = x2 (sec x + c)
39 2, 0, 6, 6 4 y = ex (−e−x + c)
7 y = e−x 31 e3x + c
42 12, 12
43 3, −7 8 y = x(xex − ex + c)
44 −2, 2 9 y = x(−e−x + c)
x −x
45 0, 0 10 y = e 2 4e 2 + c
52 1, 0 16 y = e−x ln(ex + 1) + c
53 3, 10
17 y = e−2x 2x e2x − 14 e2x + 54
54 0, 4
18 y = e−2x (ex − 41 )
55 3 − 3 sin(3), 6 − 6 sin(3) + 3 cos(3)
2 2
19 y = ex − 12 e−x + 12
56 3e6 , 5e6
20 y = 1 1 ln(1 + x2 ) + 1
57 10 sin2 (1), 4 sin(1) 2 cos(1) + sin(1)
1+x2 2
21 y = 1x − cos x + 4π+3
y2 −x2 6
58
2xy+y2
x −4x
√ 22 y = e 3 − 34 e 3 + a + 34
y−2 xy
59 √
6 xy−x
23 y = 12 x4 + 1
60 − 23 x
− xy
p
61 24 y = x3 + 2x2 + 2x + 1
3 3
62 − 12 25 y = ex − e−x + 1
63 − xy 4
26 y = 1x x4 − 13
4
4
64 − 4x2 27 y = x(x + 1)
3y
1
28 y = − x−1
65 − 2x2
3y
3 3
p
29 y = x +x+1
66 3x2 −2x3
y
30 y = cos x sin2 x + 11
2
67 1
− sin y+y cos y
31 y = − 5 21
2 (x −11)
68 − xy
p
√
y+4 xy
32 y = x3 + 2x2 − 4x + 10
69 √
4y xy−x
33 y = − √ 1
1+x2
70 8x
−1
y 2 +10 34 y = ln |x2 − 4x − 4|
3 1
71 − 2xy2 +1 35 y = −
2
3x y ln |x|− 1
2
72 − xy 36 y = sin−1 x
√1 4
73 37 y = 1x x4 − 13
1−x2 4
2
38 y = π2 cos−1 e−(x +x)
q
2y 1+x2 y2 −xy2
74 q
x2 y−2x 1+x2 y2
39 y = tan 32 x3
75 − 6x+2xy
2 2 4 3 2
x +3y 40 y = 12 x4 − x3 + x2 + 12
1
x
141
Chapter 8: 27
√
1 (1, 3)
2 (0, 1)
√
3 (− 3, −1)
4 (−3, 0)
5 (0, − 21 )
6 (3, 0)
7 ( √7 , − √7 ) 28
2 2
√
8 ( 32 , 3 2 3 )
√
9 ( 2, π4 )
10 (2, π3 )
√
11 ( 2, 3π
4 )
√
12 (2 3, π3 )
13 (2, π6 )
29 9π
4
14 ( √1 , π4 )
2 30 3π
2
15 (2, 5π )
6 31 4π
16 (3, 0) 32 4π
17 r = 4 sec θ 33 36( 3π
4 + 2)
√
18 r = 5 34 6π
19 y = 1 35 4π
20 x2 + y2 − 6x = 0 36 5π
21 r = 2 tan θ sec θ 37 2π
9 cos θ √
22 r = 9 3
cos2 θ−sin2 θ 38 2 −π
p
23 x2 + y2 − y − 3 = 0 39 10π
p
40 5π−8
24 x2 + y2 − 2 x2 + y2 + 3y = 0 4
41 π−2
25 2
26
Basic Mathematical Concepts
In this part of the book, we prepared some mathematical concepts that hopefully help students to understand the main ideas of the book.
By taking into account the different scientific levels of the students, it is necessary to present these concepts with some examples and
figures. The necessity of this part is not limited to this course, but it is for other courses. I personally recommend the students to give this
additional part a primary attention before starting the course, where the necessity of it is not limited to this course, but it is for other
courses.
Mathematical Expressions
1. ⇒ is the symbol for implying.
2. ⇔ is the symbol for “⇒ and ⇐". Also, the expression “iff" means if and only if .
3. b > a means b is greater than a and a < b means a is less than b.
4. b ≥ a means b is greater than or equal to a.
Sets of Numbers & Notations
(2) y0 = √1x 2√
1
x
1
= 2x
Exponents
Assume n is a positive integer and a is a real number. The nth power of a is
an = a . a . . . a .
The Natural Exponential Function
y = exp x ⇔ ln y = x
Some properties: If a, b > 0 and r ∈ Q, then
• 1. ea eb = ea+b
ea y = ex
2. b = ea−b
e
3. (ea )r = ear
• Note that ex and ln x are inverse functions, so
1
ln ex = x, ∀x ∈ R , and eln x = x, ∀x ∈ (0, ∞).
Solution:
√
(1) y0 = e x 2√
1
x
y = ax
1 1
y = ax
x x
Figure A.4: The function y = ax for a > 1. Figure A.5: The function y = ax for a < 1.
2. xa xb = xa+b 5. (xy)a = xa ya
xa
3. xb
= xa−b 6. x−a = 1
xa
Example A.19
(1) 24 2−7 = 23−7 = 2−3 = 1
23
= 1
8
(3) (5x)2 = 25x2
x 6 y3 x6 y3 x 6 y3 1
= x6−5 y3−5 z15 = x
2
3
(2) 3−2
= 32−(−2) = 34 = 81 (4) (xyz) 5 = x5 y5 z5 = x5 y5 z5 y2 z5
Example A.20
√
Find the derivative of the function.
(1) y = 2 x (2) y = 3tan x
Solution: √
(1) y0 = 2 x ln 2 2√
1
x
x = ay ⇔ y = loga x.
y y
y = loga x
x
1
x
1
y = loga x
Figure A.6: The function y = loga x for a > 1. Figure A.7: The function y = loga x for a < 1.
3. loga xr = r loga x
(2) y0 = √1 1
√
x ln 3 2 x
= 1
2x ln 3
Algebraic Expressions
Let a and b be real numbers. Then,
1. (a + b)2 = a2 + 2ab + b2 5. (a − b)3 = a3 − 3a2 b + 3ab2 − b3
3. (a + b)(a − b) = a2 − b2 7. a3 − b3 = (a − b)(a2 + ab + b2 )
(3) The domain of the numerator is R, but the domain of the denominator is [−2, ∞). Also, the denominator g(x) = 0 if x = −2. Thus, the
domain is (−2, ∞).
• Operations on Functions
Let f and g be two functions, then
1. ( f ± g)(x) = f (x) ± g(x) .
2. ( f g)(x) = f (x)g(x) .
f (x)
3. ( gf )(x) = g(x) where g(x) 6= 0 .
−1
1 The −1 in f −1 is not exponent where 1
f (x)
is written as f (x) .
Basic Mathematical Concepts 6
3. f −1 f (x) = x, ∀x ∈ D( f ).
4. f f −1 (x) = x, ∀x ∈ D( f −1 ).
xy = 0 ⇒ x = 0 or y = 0 .
x2 + 2x − 8 = (x − 2)(x + 4) = 0 ⇒ x − 2 = 0 or x + 4 = 0 ⇒ x = 2 or x = −4 .
(x + 2)(x + 3) = 0 ⇒ x + 2 = 0 or x + 3 = 0 ⇒ x = −2 or x = −3 .
x2 + 2x + 1 = 8 + 1 ⇒ (x + 1)2 = 9 .
Step 4: x + 1 = ±3 ⇒ x = ±3 − 1 ⇒ x = 2 or x = −4.
Systems of Equations
A system of equations consists of two or more equations with the same set of unknowns. The equations in the system can be linear or
non-linear, but for the purpose of this book, we only consider the linear ones.
Consider a system of two equations in two unknowns x and y
ax + by = c
dx + ey = f .
To solve the system, we try to find values of the unknowns that will satisfy each equation in the system. To do this, we can use elimination
or substitution.
Example A.33 Solve the following system of equations:
x − 3y = 4 → 1
2x + y = 6 → 2
Solution:
• By using the elimination method.
Multiply equation 2 by 3, then add the result to equation 1 . This implies 7x = 22 ⇒ x = 22 7 . Substitute the value of x into the first
2
or the second equation to obtain y = − 7 .
• By using the substitution method.
From the first equation, we have x = 4 + 3y. By substituting that into the second equation, we obtain
2
2(4 + 3y) + y = 6 ⇒ 7y + 8 = 6 ⇒ y = −
7
22
Substitute value of y into x = 4 + 3y to have x = 7 .
Pythagorean Theorem
Basic Mathematical Concepts 8
x
Solution: 4
a = 3, b = 4 ⇒ x2 = 42 + 32 = 25 ⇒ x = 5
cos θ = 53 θ
sin θ = 54 3
Figure A.10
Trigonometric Functions
• If (x, y) is a point on the unit circle, and if the ray
from the origin (0, 0) to that point (x, y) makes an angle
θ with the positive x-axis, then
cos θ = x , sin θ = y ,
sin θ cos θ 1 1
tan θ = cot θ = sec θ = csc θ =
cos θ sin θ cos θ sin θ
• Trigonometric functions of negative angles
cos(−θ) = cos(θ), sin(−θ) = − sin(θ), tan(−θ) = − tan(θ)
• Double and half angle formulas
Degrees 0 30 45 60 90 120 135 150 180 210 225 240 270 300 315 330 360
π π π π 2π 3π 5π 7π 5π 4π 3π 5π 7π 11π
Radians 0 6 4 3
√ 2 √3 4 6 π 6 4 3√ 2 3√ 4 6 2π
1 √1 3 3 1 1 −1 −1 − 3 − 3 −1 −1
2 −1
sin θ 0 1 √ 0 √ √ 0
2
√ 2 2 2 2 2 √ 2 2 √ 2 2 2 √
3 √1 1 −1 −1 − 3 − 3 √
−1 −1 1 √1 3
cos θ 1 2 2 2 0 2
√
2 2 −1 2 2 2 0 2 2 2 1
Table A.1
2
x x
y =−π/4
tan x
− π4 π
4 π/4
−2
y = cos x
y = sin x −1
Figure A.12: The graphs of sin x and cos x. Figure A.13: The graph of tan x.
Distance Formula
Solution:
q √ √
D = (−3 − 1)2 + (4 − 1)2 = 16 + 9 = 25 = 5 .
Figure A.14: The distance between two points.
Differentiation of Functions
Differentiation
Rules
d 0 0
dx f (x) + g(x) = f (x) + g (x) d 1
−g0 (x)
dx g(x) = 2
g(x)
d
f (x)g(x) = f 0 (x)g(x) + f (x)g0 (x)
dx d
c f (x) = c f 0 (x)
f 0 (x)g(x)− f (x)g0 (x) dx
d f (x)
dx g(x) = 2
g(x)
Elementary Derivatives
d r d 1 1 d √ 1
dx x = rxr−1 dx x = − x2
√
dx x = 2 x
Derivative of Composite Functions (Chain Rule)
Basic Mathematical Concepts 10
Let y = f (u) and u = g(x) such that dy/du and du/dx exist. Then, the derivative of the composite function ( f ◦ g)(x) exists and
dy dy du
= f 0 (u)g0 (x) = f 0 g(x) g0 (x) .
=
dx du dx
Derivative of Inverse Functions
If a function f has an inverse function f −1 , then d −1
dx f
1
(x) = 0 −1 .
f f (x)
Graphs of Functions
• The First and Second Derivative Test
1. Let f be continuous on [a, b] and f 0 exists on (a, b).
• If f 0 (x) > 0, ∀x ∈ (a, b), then f is increasing on [a, b].
• If f 0 (x) < 0, ∀x ∈ (a, b), then f is decreasing on [a, b].
2. Let f be continuous at a critical number c and differentiable on an open interval (a, b), except possibly at c.
• f (c) is a local maximum of f if f 0 changes from positive to negative at c.
• f (c) is a local minimum of f if f 0 changes from negative to positive at c.
Figure A.15: The local maximum and minimum value of the function f .
• Lines
The general linear equation in two variables x and y can be written in the form:
ax + by + c = 0 ,
2x + y = 4
a = 2, b = −1, c = −4
To plot the line, we rewrite the equation to become
y = −2x + 4. Then, we use the following table to make
points on the plane:
The line 2x + y = 4 passes
x 0 2
through the points (0, 4)
y 4 0 and (2, 0).
3. Slope-Intercept form:
If b 6= 0, the general linear equation can be rewritten as
a c
ax + by + c = 0 ⇒ by = −ax − c ⇒ y = − x − ⇒ y = mx + d ,
b b
y2 − y1 3−1 2
m= = = .
x2 − x1 3 − (−2) 5
3. Let L1 and L2 be two lines in a plane, and let m1 and m2 be the corresponding slopes, respectively.
−1
• If L1 and L2 are parallel, then m1 = m2 . • If L1 and L2 are vertical, then m1 = m2 .
y = x2 y = x2 + a y = x2 − a
√
y = (x + a)2 y = (x − a)2 x= y
√
x = y2 − a x = (y − a)2 y= x
x = y2 y = x3 y =| x |
Figure A.21
Figure A.22
Integral Calculus and Applications