Пресс-релиз популярных книг
.
Авторы: 111 А Б В Г Д Е Ж З И Й К Л М Н О П Р С Т У Ф Х Ц Ч Ш Щ Э Ю Я
Книги: 164 А Б В Г Д Е Ж З И Й К Л М Н О П Р С Т У Ф Х Ц Ч Ш Щ Э Ю Я
На сайте 111 авторов, 92 книг, 72 статей, 5913 глав.
CONTENTS
PROBLEMS 1.6 ............................................ 1
PROBLEMS 2.4 ............................................ 12
PROBLEMS 2.7 ............................................ 18
PROBLEMS 3.6 ............................................ 32
PROBLEMS 4.1 ............................................ 45
PROBLEMS 5.8 ............................................ 58
PROBLEMS 6.3 ............................................ 69
PROBLEMS 7.3 ............................................ 83
PROBLEMS 8.8 ............................................ 91
i
SECTION 1:6
2. (i)
·
0 0 0
2 4 0
¸
R1 $ R2
·
2 4 0
0 0 0
¸
R1 ! 1
2R1
·
1 2 0
0 0 0
¸
;
(ii)
·
0 1 3
1 2 4
¸
R1 $ R2
·
1 2 4
0 1 3
¸
R1 ! R1 ¡ 2R2
·
1 0 ¡2
0 1 3
¸
;
(iii)
2
4
1 1 1
1 1 0
1 0 0
3
5 R2 ! R2 ¡ R1
R3 ! R3 ¡ R1
2
4
1 1 0
0 0 ¡1
0 ¡1 ¡1
3
5
R1 ! R1 + R3
R3 ! ¡R3
R2 $ R3
2
4
1 0 0
0 1 1
0 0 ¡1
3
5 R2 ! R2 + R3
R3 ! ¡R3
2
4
1 0 0
0 1 0
0 0 1
3
5;
(iv)
2
4
2 0 0
0 0 0
¡4 0 0
3
5 R3 ! R3 + 2R1
R1 ! 1
2R1
2
4
1 0 0
0 0 0
0 0 0
3
5.
3. (a)
2
4
1 1 1 2
2 3 ¡1 8
1 ¡1 ¡1 ¡8
3
5 R2 ! R2 ¡ 2R1
R3 ! R3 ¡ R1
2
4
1 1 1 2
0 1 ¡3 4
0 ¡2 ¡2 ¡10
3
5
R1 ! R1 ¡ R2
R3 ! R3 + 2R2
2
4
1 0 4 ¡2
0 1 ¡3 4
0 0 ¡8 ¡2
3
5R3 ! ¡1
8 R3
2
4
1 0 4 2
0 1 ¡3 4
0 0 1 1
4
3
5
R1 ! R1 ¡ 4R3
R2 ! R2 + 3R3
2
4
1 0 0 ¡3
0 1 0 19
4
0 0 1 1
4
3
5.
The augmented matrix has been converted to reduced row{echelon form
and we read o® the unique solution x = ¡3; y = 19
4 ; z = 1
4 .
(b)
2
4
1 1 ¡1 2 10
3 ¡1 7 4 1
¡5 3 ¡15 ¡6 9
3
5 R2 ! R2 ¡ 3R1
R3 ! R3 + 5R1
2
4
1 1 ¡1 2 10
0 ¡4 10 ¡2 ¡29
0 8 ¡20 4 59
3
5
R3 ! R3 + 2R2
2
4
1 1 ¡1 2 10
0 ¡4 10 ¡2 ¡29
0 0 0 0 1
3
5.
From the last matrix we see that the original system is inconsistent.
1
(c)
2
664
3 ¡1 7 0
2 ¡1 4 1
2
1 ¡1 1 1
6 ¡4 10 3
3
775
R1 $ R3
2
664
1 ¡1 1 1
2 ¡1 4 1
2
3 ¡1 7 0
6 ¡4 10 3
3
775
R2 ! R2 ¡ 2R1
R3 ! R3 ¡ 3R1
R4 ! R4 ¡ 6R1
2
664
1 ¡1 1 1
0 1 2 ¡3
2
0 2 4 ¡3
0 2 4 ¡3
3
775
R1 ! R1 + R2
R4 ! R4 ¡ R3
R3 ! R3 ¡ 2R2
2
664
1 0 3 ¡1
2
0 1 2 ¡3
2
0 0 0 0
0 0 0 0
3
775
.
The augmented matrix has been converted to reduced row{echelon form
and we read o® the complete solution x = ¡1
2 ¡ 3z; y = ¡3
2 ¡ 2z, with z
arbitrary.
4.
2
4
2 ¡1 3 a
3 1 ¡5 b
¡5 ¡5 21 c
3
5R2 ! R2 ¡ R1
2
4
2 ¡1 3 a
1 2 ¡8 b ¡ a
¡5 ¡5 21 c
3
5
R1 $ R2
2
4
1 2 ¡8 b ¡ a
2 ¡1 3 a
¡5 ¡5 21 c
3
5 R2 ! R2 ¡ 2R1
R3 ! R3 + 5R1
2
4
1 2 ¡8 b ¡ a
0 ¡5 19 ¡2b + 3a
0 5 ¡19 5b ¡ 5a + c
3
5
R3 ! R3 + R2
R2 ! ¡1
5 R2
2
4
1 2 ¡8 b ¡ a
0 1 ¡19
5
2b¡3a
5
0 0 0 3b ¡ 2a + c
3
5
R1 ! R1 ¡ 2R2
2
4
1 0 ¡2
5
(b+a)
5
0 1 ¡19
5
2b¡3a
5
0 0 0 3b ¡ 2a + c
3
5.
From the last matrix we see that the original system is inconsistent if
3b¡2a+c 6= 0. If 3b¡2a+c = 0, the system is consistent and the solution
is
x =
(b + a)
5
+
2
5
z; y =
(2b ¡ 3a)
5
+
19
5
z;
where z is arbitrary.
5.
2
4
1 1 1
t 1 t
1 + t 2 3
3
5 R2 ! R2 ¡ tR1
R3 ! R3 ¡ (1 + t)R1
2
4
1 1 1
0 1 ¡ t 0
0 1 ¡ t 2 ¡ t
3
5
R3 ! R3 ¡ R2
2
4
1 1 1
0 1 ¡ t 0
0 0 2 ¡ t
3
5 = B:
Case 1. t 6= 2. No solution.
2
Case 2. t = 2. B =
2
4
1 0 1
0 ¡1 0
0 0 0
3
5 !
2
4
1 0 1
0 1 0
0 0 0
3
5 :
We read o® the unique solution x = 1; y = 0.
6. M2ethod 1.
664
¡3 1 1 1
1 ¡3 1 1
1 1 ¡3 1
1 1 1 ¡3
3
775
R1 ! R1 ¡ R4
R2 ! R2 ¡ R4
R3 ! R3 ¡ R4
2
664
¡4 0 0 4
0 ¡4 0 4
0 0 ¡4 4
1 1 1 ¡3
3
775
!
2
664
1 0 0 ¡1
0 1 0 ¡1
0 0 1 ¡1
1 1 1 ¡3
3
775
R4 ! R4 ¡ R3 ¡ R2 ¡ R1
2
664
1 0 0 ¡1
0 1 0 ¡1
0 0 1 ¡1
0 0 0 0
3
775
:
Hence the given homogeneous system has complete solution
x1 = x4; x2 = x4; x3 = x4;
with x4 arbitrary.
Method 2. Write the system as
x1 + x2 + x3 + x4 = 4x1
x1 + x2 + x3 + x4 = 4x2
x1 + x2 + x3 + x4 = 4x3
x1 + x2 + x3 + x4 = 4x4:
Then it is immediate that any solution must satisfy x1 = x2 = x3 = x4.
Conversely, if x1; x2; x3; x4 satisfy x1 = x2 = x3 = x4, we get a solution.
7.
·
¸ ¡ 3 1
1 ¸ ¡ 3
¸
R1 $ R2
·
1 ¸ ¡ 3
¸ ¡ 3 1
¸
R2 ! R2 ¡ (¸ ¡ 3)R1
·
1 ¸ ¡ 3
0 ¡¸2 + 6¸ ¡ 8
¸
= B:
Case 1: ¡¸2 + 6¸ ¡ 8 6= 0. That is ¡(¸ ¡ 2)(¸ ¡ 4) 6= 0 or ¸ 6= 2; 4. Here B is
row equivalent to
·
1 0
0 1
¸
:
R2 ! 1
¡¸2+6¸¡8R2
·
1 ¸ ¡ 3
0 1
¸
R1 ! R1 ¡ (¸ ¡ 3)R2
·
1 0
0 1
¸
:
Hence we get the trivial solution x = 0; y = 0.
3
Case 2: ¸ = 2. Then B =
·
1 ¡1
0 0
¸
and the solution is x = y, with y
arbitrary.
Case 3: ¸ = 4. Then B =
·
1 1
0 0
¸
and the solution is x = ¡y, with y
arbitrary.
8.
·
3 1 1 1
5 ¡1 1 ¡1
¸
R1 !
1
3
R1
·
1 1
3
1
3
1
3
5 ¡1 1 ¡1
¸
R2 ! R2 ¡ 5R1
·
1 1
3
1
3
1
3
0 ¡8
3 ¡2
3 ¡8
3
¸
R2 ! ¡3
8
R2
·
1 1
3
1
3
1
3
0 1 1
4 1
¸
R1 ! R1 ¡
1
3
R2
·
1 0 1
4 0
0 1 1
4 1
¸
:
Hence the solution of the associated homogeneous system is
x1 = ¡
1
4
x3; x2 = ¡
1
4
x3 ¡ x4;
with x3 and x4 arbitrary.
9.
A =
2
6664
1 ¡ n 1 ¢ ¢ ¢ 1
1 1 ¡ n ¢ ¢ ¢ 1
...
...
¢ ¢ ¢
...
1 1 ¢ ¢ ¢ 1 ¡ n
3
7775
R1 ! R1 ¡ Rn
R2 ! R2 ¡ Rn
...
Rn¡1 ! Rn¡1 ¡ Rn
2
6664
¡n 0 ¢ ¢ ¢ n
0 ¡n ¢ ¢ ¢ n
...
...
¢ ¢ ¢
...
1 1 ¢ ¢ ¢ 1 ¡ n
3
7775
!
2
6664
1 0 ¢ ¢ ¢ ¡1
0 1 ¢ ¢ ¢ ¡1
...
...
¢ ¢ ¢
...
1 1 ¢ ¢ ¢ 1 ¡ n
3
7775
Rn ! Rn ¡ Rn¡1 ¢ ¢ ¢ ¡ R1
2
6664
1 0 ¢ ¢ ¢ ¡1
0 1 ¢ ¢ ¢ ¡1
...
...
¢ ¢ ¢
...
0 0 ¢ ¢ ¢ 0
3
7775
:
The last matrix is in reduced row{echelon form.
Consequently the homogeneous system with coe±cient matrix A has the
solution
x1 = xn; x2 = xn; : : : ; xn¡1 = xn;
4
with xn arbitrary.
Alternatively, writing the system in the form
x1 + ¢ ¢ ¢ + xn = nx1
x1 + ¢ ¢ ¢ + xn = nx2
...
x1 + ¢ ¢ ¢ + xn = nxn
shows that any solution must satisfy nx1 = nx2 = ¢ ¢ ¢ = nxn, so x1 = x2 =
¢ ¢ ¢ = xn. Conversely if x1 = xn; : : : ; xn¡1 = xn, we see that x1; : : : ; xn is a
solution.
10. Let A =
·
a b
c d
¸
and assume that ad ¡ bc 6= 0.
Case 1: a 6= 0.
·
a b
c d
¸
R1 ! 1
aR1
·
1 b
a
c d
¸
R2 ! R2 ¡ cR1
·
1 b
a
0 ad¡bc
a
¸
R2 ! a
ad¡bcR2
·
1 b
a
0 1
¸
R1 ! R1 ¡ b
aR2
·
1 0
0 1
¸
:
Case 2: a = 0. Then bc 6= 0 and hence c 6= 0.
A =
·
0 b
c d
¸
R1 $ R2
·
c d
0 b
¸
!
·
1 d
c
0 1
¸
!
·
1 0
0 1
¸
.
So in both cases, A has reduced row{echelon form equal to
·
1 0
0 1
¸
.
11. We simplify the augmented matrix of the system using row operations:
2
4
1 2 ¡3 4
3 ¡1 5 2
4 1 a2 ¡ 14 a + 2
3
5 R2 ! R2 ¡ 3R1
R3 ! R3 ¡ 4R1
2
4
1 2 ¡3 4
0 ¡7 14 ¡10
0 ¡7 a2 ¡ 2 a ¡ 14
3
5
R3 ! R3 ¡ R2
R2 ! ¡1
7 R2
R1 ! R1 ¡ 2R2
2
4
1 2 ¡3 4
0 1 ¡2 10
7
0 0 a2 ¡ 16 a ¡ 4
3
5 R1 ! R1 ¡ 2R2
2
4
1 0 1 8
7
0 1 ¡2 10
7
0 0 a2 ¡ 16 a ¡ 4
3
5 :
Denote the last matrix by B.
5
Case 1: a2 ¡ 16 6= 0. i.e. a 6= §4. Then
R3 ! 1
a2¡16R3
R1 ! R1 ¡ R3
R2 ! R2 + 2R3
2
64
1 0 0 8a+25
7(a+4)
0 1 0 10a+54
7(a+4)
0 0 1 1
a+4
3
75
and we get the unique solution
x =
8a + 25
7(a + 4)
; y =
10a + 54
7(a + 4)
; z =
1
a + 4
:
Case 2: a = ¡4. Then B =
2
4
1 0 1 8
7
0 1 ¡2 10
7
0 0 0 ¡8
3
5, so our system is inconsistent.
Case 3: a = 4. Then B =
2
4
1 0 1 8
7
0 1 ¡2 10
7
0 0 0 0
3
5. We read o® that the system is
consistent, with complete solution x = 8
7 ¡ z; y = 10
7 + 2z, where z is
arbitrary.
12. We reduce the augmented array of the system to reduced row{echelon
form:
2
664
1 0 1 0 1
0 1 0 1 1
1 1 1 1 0
0 0 1 1 0
3
775
R3 ! R3 + R1
2
664
1 0 1 0 1
0 1 0 1 1
0 1 0 1 1
0 0 1 1 0
3
775
R3 ! R3 + R2
2
664
1 0 1 0 1
0 1 0 1 1
0 0 0 0 0
0 0 1 1 0
3
775
R1 ! R1 + R4
R3 $ R4
2
664
1 0 0 1 1
0 1 0 1 1
0 0 1 1 0
0 0 0 0 0
3
775
:
The last matrix is in reduced row{echelon form and we read o® the solution
of the corresponding homogeneous system:
x1 = ¡x4 ¡ x5 = x4 + x5
x2 = ¡x4 ¡ x5 = x4 + x5
x3 = ¡x4 = x4;
6
where x4 and x5 are arbitrary elements of Z2. Hence there are four solutions:
x1 x2 x3 x4 x5
0 0 0 0 0
1 1 0 0 1
1 1 1 1 0
0 0 1 1 1
:
13. (a) We reduce the augmented matrix to reduced row{echelon form:
2
4
2 1 3 4
4 1 4 1
3 1 2 0
3
5 R1 ! 3R1
2
4
1 3 4 2
4 1 4 1
3 1 2 0
3
5
R2 ! R2 + R1
R3 ! R3 + 2R1
2
4
1 3 4 2
0 4 3 3
0 2 0 4
3
5 R2 ! 4R2
2
4
1 3 4 2
0 1 2 2
0 2 0 4
3
5
R1 ! R1 + 2R2
R3 ! R3 + 3R2
2
4
1 0 3 1
0 1 2 2
0 0 1 0
3
5 R1 ! R1 + 2R3
R2 ! R2 + 3R3
2
4
1 0 0 1
0 1 0 2
0 0 1 0
3
5 :
Consequently the system has the unique solution x = 1; y = 2; z = 0.
(b) Again we reduce the augmented matrix to reduced row{echelon form:
2
4
2 1 3 4
4 1 4 1
1 1 0 3
3
5 R1 $ R3
2
4
1 1 0 3
4 1 4 1
2 1 3 4
3
5
R2 ! R2 + R1
R3 ! R3 + 3R1
2
4
1 1 0 3
0 2 4 4
0 4 3 3
3
5 R2 ! 3R2
2
4
1 1 0 3
0 1 2 2
0 4 3 3
3
5
R1 ! R1 + 4R2
R3 ! R3 + R2
2
4
1 0 3 1
0 1 2 2
0 0 0 0
3
5 :
We read o® the complete solution
x = 1 ¡ 3z = 1 + 2z
y = 2 ¡ 2z = 2 + 3z;
where z is an arbitrary element of Z5.
7
14. Suppose that (®1; : : : ; ®n) and (¯1; : : : ; ¯n) are solutions of the system
of linear equations
Xn
j=1
aijxj = bi; 1 · i · m:
Then
Xn
j=1
aij®j = bi and
Xn
j=1
aij¯j = bi
for 1 · i · m.
Let °i = (1 ¡t)®i +t¯i for 1 · i · m. Then (°1; : : : ; °n) is a solution of
the given system. For
Xn
j=1
aij°j =
Xn
j=1
aijf(1 ¡ t)®j + t¯jg
=
Xn
j=1
aij(1 ¡ t)®j +
Xn
j=1
aijt¯j
= (1 ¡ t)bi + tbi
= bi:
15. Suppose that (®1; : : : ; ®n) is a solution of the system of linear equations
Xn
j=1
aijxj = bi; 1 · i · m: (1)
Then the system can be rewritten as
Xn
j=1
aijxj =
Xn
j=1
aij®j ; 1 · i · m;
or equivalently
Xn
j=1
aij(xj ¡ ®j) = 0; 1 · i · m:
So we have
Xn
j=1
aijyj = 0; 1 · i · m:
where xj ¡ ®j = yj . Hence xj = ®j + yj ; 1 · j · n, where (y1; : : : ; yn) is
a solution of the associated homogeneous system. Conversely if (y1; : : : ; yn)
8
is a solution of the associated homogeneous system and xj = ®j + yj ; 1 ·
j · n, then reversing the argument shows that (x1; : : : ; xn) is a solution of
the system 1 .
16. We simplify the augmented matrix using row operations, working to-
wards row{echelon form:
2
4
1 1 ¡1 1 1
a 1 1 1 b
3 2 0 a 1 + a
3
5 R2 ! R2 ¡ aR1
R3 ! R3 ¡ 3R1
2
4
1 1 ¡1 1 1
0 1 ¡ a 1 + a 1 ¡ a b ¡ a
0 ¡1 3 a ¡ 3 a ¡ 2
3
5
R2 $ R3
R2 ! ¡R2
2
4
1 1 ¡1 1 1
0 1 ¡3 3 ¡ a 2 ¡ a
0 1 ¡ a 1 + a 1 ¡ a b ¡ a
3
5
R3 ! R3 + (a ¡ 1)R2
2
4
1 1 ¡1 1 1
0 1 ¡3 3 ¡ a 2 ¡ a
0 0 4 ¡ 2a (1 ¡ a)(a ¡ 2) ¡a2 + 2a + b ¡ 2
3
5 = B:
Case 1: a 6= 2. Then 4 ¡ 2a 6= 0 and
B !
2
4
1 1 ¡1 1 1
0 1 ¡3 3 ¡ a 2 ¡ a
0 0 1 a¡1
2 ¡a2+2a+b¡2
4¡2a
3
5 :
Hence we can solve for x; y and z in terms of the arbitrary variable w.
Case 2: a = 2. Then
B =
2
4
1 1 ¡1 1 1
0 1 ¡3 1 0
0 0 0 0 b ¡ 2
3
5 :
Hence there is no solution if b 6= 2. However if b = 2, then
B =
2
4
1 1 ¡1 1 1
0 1 ¡3 1 0
0 0 0 0 0
3
5 !
2
4
1 0 2 0 1
0 1 ¡3 1 0
0 0 0 0 0
3
5
and we get the solution x = 1 ¡ 2z; y = 3z ¡ w, where w is arbitrary.
17. (a) We ¯rst prove that 1 + 1 + 1 + 1 = 0. Observe that the elements
1 + 0; 1 + 1; 1 + a; 1 + b
9
are distinct elements of F by virtue of the cancellation law for addition. For
this law states that 1+x = 1+y ) x = y and hence x 6= y ) 1+x 6= 1+y.
Hence the above four elements are just the elements 0; 1; a; b in some
order. Consequently
(1 + 0) + (1 + 1) + (1 + a) + (1 + b) = 0 + 1 + a + b
(1 + 1 + 1 + 1) + (0 + 1 + a + b) = 0 + (0 + 1 + a + b);
so 1 + 1 + 1 + 1 = 0 after cancellation.
Now 1 + 1 + 1 + 1 = (1 + 1)(1 + 1), so we have x2 = 0, where x = 1 + 1.
Hence x = 0. Then a + a = a(1 + 1) = a ¢ 0 = 0.
Next a + b = 1. For a + b must be one of 0; 1; a; b. Clearly we can't
have a + b = a or b; also if a + b = 0, then a + b = a + a and hence b = a;
hence a + b = 1. Then
a + 1 = a + (a + b) = (a + a) + b = 0 + b = b:
Similarly b + 1 = a. Consequently the addition table for F is
+ 0 1 a b
0 0 1 a b
1 1 0 b a
a a b 0 1
b b a 1 0
.
We now ¯nd the multiplication table. First, ab must be one of 1; a; b;
however we can't have ab = a or b, so this leaves ab = 1.
Next a2 = b. For a2 must be one of 1; a; b; however a2 = a ) a = 0 or
a = 1; also
a2 = 1 ) a2 ¡ 1 = 0 ) (a ¡ 1)(a + 1) = 0 ) (a ¡ 1)2 = 0 ) a = 1;
hence a2 = b. Similarly b2 = a. Consequently the multiplication table for F
is
£ 0 1 a b
0 0 0 0 0
1 0 1 a b
a 0 a b 1
b 0 b 1 a
.
(b) We use the addition and multiplication tables for F:
A =
2
4
1 a b a
a b b 1
1 1 1 a
3
5 R2 ! R2 + aR1
R3 ! R3 + R1
2
4
1 a b a
0 0 a a
0 b a 0
3
5
10
R2 $ R3
2
4
1 a b a
0 b a 0
0 0 a a
3
5 R2 ! aR2
R3 ! bR3
2
4
1 a b a
0 1 b 0
0 0 1 1
3
5
R1 $ R1 + aR2
2
4
1 0 a a
0 1 b 0
0 0 1 1
3
5 R1 ! R1 + aR3
R2 ! R2 + bR3
2
4
1 0 0 0
0 1 0 b
0 0 1 1
3
5 :
The last matrix is in reduced row{echelon form.
11
Section 2:4
2. Suppose B =
2
4
a b
c d
e f
3
5 and that AB = I2. Then
·
¡1 0 1
0 1 0
¸ 2
4
a b
c d
e f
3
5 =
·
1 0
0 1
¸
=
·
¡a + e ¡b + f
c + e d + f
¸
:
Hence
¡a + e = 1
c + e = 0
; ¡b + f = 0
d + f = 1
;
e = a + 1
c = ¡e = ¡(a + 1)
;
f = b
d = 1 ¡ f = 1 ¡ b
;
B =
2
4
a b
¡a ¡ 1 1 ¡ b
a + 1 b
3
5 :
Next,
(BA)2B = (BA)(BA)B = B(AB)(AB) = BI2I2 = BI2 = B
4. Let pn denote the statement
An = (3n¡1)
2 A + (3¡3n)
2 I2:
Then p1 asserts that A = (3¡1)
2 A + (3¡3)
2 I2; which is true. So let n ¸ 1 and
assume pn. Then from (1),
An+1 = A ¢ An = A
n
(3n¡1)
2 A + (3¡3n)
2 I2
o
= (3n¡1)
2 A2 + (3¡3n)
2 A
= (3n¡1)
2 (4A ¡ 3I2) + (3¡3n)
2 A = (3n¡1)4+(3¡3n)
2 A + (3n¡1)(¡3)
2 I2
= (4¢3n¡3n)¡1
2 A + (3¡3n+1)
2 I2
= (3n+1¡1)
2 A + (3¡3n+1)
2 I2:
Hence pn+1 is true and the induction proceeds.
5. The equation xn+1 = axn + bxn¡1 is seen to be equivalent to
·
xn+1
xn
¸
=
·
a b
1 0
¸ ·
xn
xn¡1
¸
12
or
Xn = AXn¡1;
where Xn =
·
xn+1
xn
¸
and A =
·
a b
1 0
¸
. Then
Xn = AnX0
if n ¸ 1. Hence by Question 3,
·
xn+1
xn
¸
=
½
(3n ¡ 1)
2
A +
(3 ¡ 3n)
2
I2
¾·
x1
x0
¸
=
½
(3n ¡ 1)
2
·
4 ¡3
1 0
¸
+
· 3¡3n
2 0
0 3¡3n
2
¸¾·
x1
x0
¸
=
2
4
(3n ¡ 1)2 + 3¡3n
2 (3n ¡ 1)(¡3)
3n¡1
2
3¡3n
2
3
5
·
x1
x0
¸
Hence, equating the (2; 1) elements gives
xn =
(3n ¡ 1)
2
x1 +
(3 ¡ 3n)
2
x0 if n ¸ 1
7. Note: ¸1 + ¸2 = a + d and ¸1¸2 = ad ¡ bc.
Then
(¸1 + ¸2)kn ¡ ¸1¸2kn¡1 = (¸1 + ¸2)(¸n¡1
1 + ¸n¡2
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡2
2 + ¸n¡1
2 )
¡¸1¸2(¸n¡2
1 + ¸n¡3
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡3
2 + ¸n¡2
2 )
= (¸n
1 + ¸n¡1
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡1
2 )
+(¸n¡1
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡1
2 + ¸n
2 )
¡(¸n¡1
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡1
2 )
= ¸n
1 + ¸n¡1
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡1
2 + ¸n
2 = kn+1
If ¸1 = ¸2, we see
kn = ¸n¡1
1 + ¸n¡2
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡2
2 + ¸n¡1
2
= ¸n¡1
1 + ¸n¡2
1 ¸1 + ¢ ¢ ¢ + ¸1¸n¡2
1 + ¸n¡1
1
= n¸n¡1
1
13
If ¸1 6= ¸2, we see that
(¸1 ¡ ¸2)kn = (¸1 ¡ ¸2)(¸n¡1
1 + ¸n¡2
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡2
2 + ¸n¡1
2 )
= ¸n
1 + ¸n¡1
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡1
2
¡(¸n¡1
1 ¸2 + ¢ ¢ ¢ + ¸1¸n¡1
2 + ¸n
2 )
= ¸n
1 ¡ ¸n
2 :
Hence kn = ¸n
1 ¡¸n
2
¸1¡¸2
:
We have to prove
An = knA ¡ ¸1¸2kn¡1I2: ¤
n=1:
A1 = A; also k1A ¡ ¸1¸2k0I2 = k1A ¡ ¸1¸20I2
= A:
Let n ¸ 1 and assume equation ¤ holds. Then
An+1 = An ¢ A = (knA ¡ ¸1¸2kn¡1I2)A
= knA2 ¡ ¸1¸2kn¡1A:
Now A2 = (a + d)A ¡ (ad ¡ bc)I2 = (¸1 + ¸2)A ¡ ¸1¸2I2: Hence
An+1 = kn(¸1 + ¸2)A ¡ ¸1¸2I2 ¡ ¸1¸2kn¡1A
= fkn(¸1 + ¸2) ¡ ¸1¸2kn¡1gA ¡ ¸1¸2knI2
= kn+1A ¡ ¸1¸2knI2;
and the induction goes through.
8. Here ¸1; ¸2 are the roots of the polynomial x2 ¡ 2x ¡ 3 = (x ¡ 3)(x + 1).
So we can take ¸1 = 3; ¸2 = ¡1: Then
kn =
3n ¡ (¡1)n
3 ¡ (¡1)
=
3n + (¡1)n+1
4
:
Hence
An =
½
3n + (¡1)n+1
4
¾
A ¡ (¡3)
½
3n¡1 + (¡1)n
4
¾
I2
=
3n + (¡1)n+1
4
·
1 2
2 1
¸
+ 3
½
3n¡1 + (¡1)n
4
¾·
1 0
0 1
¸
;
14
which is equivalent to the stated result.
9. In terms of matrices, we have
·
Fn+1
Fn
¸
=
·
1 1
1 0
¸ ·
Fn
Fn¡1
¸
for n ¸ 1:
·
Fn+1
Fn
¸
=
·
1 1
1 0
¸n ·
F1
F0
¸
=
·
1 1
1 0
¸n ·
1
0
¸
:
Now ¸1; ¸2 are the roots of the polynomial x2 ¡ x ¡ 1 here.
Hence ¸1 = 1+p5
2 and ¸2 = 1¡
p5
2 and
kn =
³
1+p5
2
´n¡1
¡
³
1¡
p5
2
´n¡1
1+p5
2 ¡
³
1¡
p5
2
´
=
³
1+p5
2
´n¡1
¡
³
1¡
p5
2
´n¡1
p5
:
Hence
An = knA ¡ ¸1¸2kn¡1I2
= knA + kn¡1I2
So
·
Fn+1
Fn
¸
= (knA + kn¡1I2)
·
1
0
¸
= kn
·
1
1
¸
+ kn¡1
·
1
0
¸
=
·
kn + kn¡1
kn
¸
:
Hence
Fn = kn =
³
1+p5
2
´n¡1
¡
³
1¡
p5
2
´n¡1
p5
:
10. From Question 5, we know that
·
xn
yn
¸
=
·
1 r
1 1
¸n ·
a
b
¸
:
15
Now by Question 7, with A =
·
1 r
1 1
¸
,
An = knA ¡ ¸1¸2kn¡1I2
= knA ¡ (1 ¡ r)kn¡1I2;
where ¸1 = 1 + pr and ¸2 = 1 ¡ pr are the roots of the polynomial
x2 ¡ 2x + (1 ¡ r) and
kn =
¸n
1 ¡ ¸n
2
2pr
:
Hence
·
xn
yn
¸
= (knA ¡ (1 ¡ r)kn¡1I2)
·
a
b
¸
=
µ·
kn knr
kn kn
¸
¡
·
(1 ¡ r)kn¡1 0
0 (1 ¡ r)kn¡1
¸¶·
a
b
¸
=
·
kn ¡ (1 ¡ r)kn¡1 knr
kn kn ¡ (1 ¡ r)kn¡1
¸ ·
a
b
¸
=
·
a(kn ¡ (1 ¡ r)kn¡1) + bknr
akn + b(kn ¡ (1 ¡ r)kn¡1)
¸
:
Hence, in view of the fact that
kn
kn¡1
=
¸n
1 ¡ ¸n
2
¸n¡1
1 ¡ ¸n¡1
2
=
¸n
1 (1 ¡ f¸2
¸1 gn)
¸n¡1
1 (1 ¡ f¸2
¸1 gn¡1) ! ¸1; as n ! 1;
we have
·
xn
yn
¸
=
a(kn ¡ (1 ¡ r)kn¡1) + bknr
akn + b(kn ¡ (1 ¡ r)kn¡1)
=
a( kn
kn¡1 ¡ (1 ¡ r)) + b kn
kn¡1
r
a kn
kn¡1
+ b( kn
kn¡1 ¡ (1 ¡ r))
!
a(¸1 ¡ (1 ¡ r)) + b¸1r
a¸1 + b(¸1 ¡ (1 ¡ r))
=
a(pr + r) + b(1 + pr)r
a(1 + pr) + b(pr + r)
=
prfa(1 + pr) + b(1 + pr)prg
a(1 + pr) + b(pr + r)
= pr:
16
Section 2:7
1. [AjI2] =
·
1 4
¡3 1
¯¯¯¯
1 0
0 1
¸
R2 ! R2 + 3R1
·
1 4
0 13
¯¯¯¯
1 0
3 1
¸
R2 ! 1
13R2
·
1 4
0 1
¯¯¯¯
1 0
3=13 1=13
¸
R1 ! R1 ¡ 4R2
·
1 0
0 1
¯¯¯¯
1=13 ¡4=13
3=13 1=13
¸
.
Hence A is non{singular and A¡1 =
·
1=13 ¡4=13
3=13 1=13
¸
.
Moreover
E12(¡4)E2(1=13)E21(3)A = I2;
so
A¡1 = E12(¡4)E2(1=13)E21(3):
Hence
A = fE21(3)g¡1fE2(1=13)g¡1fE12(¡4)g¡1 = E21(¡3)E2(13)E12(4):
2. Let D = [dij ] be an m£m diagonal matrix and let A = [ajk] be an m£n
matrix. Then
(DA)ik =
Xn
j=1
dijajk = diiaik;
as dij = 0 if i 6= j. It follows that the ith row of DA is obtained by
multiplying the ith row of A by dii.
Similarly, post{multiplication of a matrix by a diagonal matrix D results
in a matrix whose columns are those of A, multiplied by the respective
diagonal elements of D.
In particular,
diag (a1; : : : ; an)diag (b1; : : : ; bn) = diag (a1b1; : : : ; anbn);
as the left{hand side can be regarded as pre{multiplication of the matrix
diag (b1; : : : ; bn) by the diagonal matrix diag (a1; : : : ; an).
Finally, suppose that each of a1; : : : ; an is non{zero. Then a¡1
1 ; : : : ; a¡1
n
all exist and we have
diag (a1; : : : ; an)diag (a¡1
1 ; : : : ; a¡1
n ) = diag (a1a¡1
1 ; : : : ; ana¡1
n )
= diag (1; : : : ; 1) = In:
Hence diag (a1; : : : ; an) is non{singular and its inverse is diag (a¡1
1 ; : : : ; a¡1
n ).
17
Next suppose that ai = 0. Then diag (a1; : : : ; an) is row{equivalent to a
matix containing a zero row and is hence singular.
3. [AjI3] =
2
4
0 0 2
1 2 6
3 7 9
¯¯¯¯¯¯
1 0 0
0 1 0
0 0 1
3
5 R1 $ R2
2
4
1 2 6 0 1 0
0 0 2 1 0 0
3 7 9 0 0 1
3
5
R3 ! R3 ¡ 3R1
2
4
1 2 6 0 1 0
0 0 2 1 0 0
0 1 ¡9 0 ¡3 1
3
5 R2 $ R3
2
4
1 2 6 0 1 0
0 1 ¡9 0 ¡3 1
0 0 2 1 0 0
3
5
R3 ! 1
2R3
2
4
1 2 6 0 1 0
0 1 ¡9 0 ¡3 1
0 0 1 1=2 0 0
3
5 R1 ! R1 ¡ 2R2
2
4
1 0 24 0 7 ¡2
0 1 ¡9 0 ¡3 1
0 0 1 1=2 0 0
3
5
R1 ! R1 ¡ 24R3
R2 ! R2 + 9R3
2
4
1 0 0 ¡12 7 ¡2
0 1 0 9=2 ¡3 1
0 0 1 1=2 0 0
3
5.
Hence A is non{singular and A¡1 =
2
4 ¡12 7 ¡2
9=2 ¡3 1
1=2 0 0
3
5.
Also
E23(9)E13(¡24)E12(¡2)E3(1=2)E23E31(¡3)E12A = I3:
Hence
A¡1 = E23(9)E13(¡24)E12(¡2)E3(1=2)E23E31(¡3)E12;
so
A = E12E31(3)E23E3(2)E12(2)E13(24)E23(¡9):
4.
A =
2
4
1 2 k
3 ¡1 1
5 3 ¡5
3
5 !
2
4
1 2 k
0 ¡7 1 ¡ 3k
0 ¡7 ¡5 ¡ 5k
3
5 !
2
4
1 2 k
0 ¡7 1 ¡ 3k
0 0 ¡6 ¡ 2k
3
5 = B:
Hence if ¡6¡2k 6= 0, i.e. if k 6= ¡3, we see that B can be reduced to I3
and hence A is non{singular.
If k = ¡3, then B =
2
4
1 2 ¡3
0 ¡7 10
0 0 0
3
5 = B and consequently A is singu-
lar, as it is row{equivalent to a matrix containing a zero row.
18
5. E21(2)
·
1 2
¡2 ¡4
¸
=
·
1 2
0 0
¸
. Hence, as in the previous question,
·
1 2
¡2 ¡4
¸
is singular.
6. Starting from the equation A2 ¡ 2A + 13I2 = 0, we deduce
A(A ¡ 2I2) = ¡13I2 = (A ¡ 2I2)A:
Hence AB = BA = I2, where B = ¡1
13 (A ¡ 2I2). Consequently A is non{
singular and A¡1 = B.
7. We assume the equation A3 = 3A2 ¡ 3A + I3.
(ii) A4 = A3A = (3A2 ¡ 3A + I3)A = 3A3 ¡ 3A2 + A
= 3(3A2 ¡ 3A + I3) ¡ 3A2 + A = 6A2 ¡ 8A + 3I3:
(iii) A3 ¡ 3A2 + 3A = I3. Hence
A(A2 ¡ 3A + 3I3) = I3 = (A2 ¡ 3A + 3I3)A:
Hence A is non{singular and
A¡1 = A2 ¡ 3A + 3I3
=
2
4 ¡1 ¡3 1
2 4 ¡1
0 1 0
3
5 :
8. (i) If B3 = 0 then
(In ¡ B)(In + B + B2) = In(In + B + B2) ¡ B(In + B + B2)
= (In + B + B2) ¡ (B + B2 + B3)
= In ¡ B3 = In ¡ 0 = In:
Similarly (In + B + B2)(In ¡ B) = In.
Hence A = In ¡ B is non{singular and A¡1 = In + B + B2.
It follows that the system AX = b has the unique solution
X = A¡1b = (In + B + B2)b = b + Bb + B2b:
19
(ii) Let B =
2
4
0 r s
0 0 t
0 0 0
3
5. Then B2 =
2
4
0 0 rt
0 0 0
0 0 0
3
5 and B3 = 0. Hence
from the preceding question
(I3 ¡ B)¡1 = I3 + B + B2
=
2
4
1 0 0
0 1 0
0 0 1
3
5 +
2
4
0 r s
0 0 t
0 0 0
3
5 +
2
4
0 0 rt
0 0 0
0 0 0
3
5
=
2
4
1 r s + rt
0 1 t
0 0 1
3
5 :
9. (i) Suppose that A2 = 0. Then if A¡1 exists, we deduce that A¡1(AA) =
A¡10, which gives A = 0 and this is a contradiction, as the zero matrix is
singular. We conclude that A does not have an inverse.
(ii). Suppose that A2 = A and that A¡1 exists. Then
A¡1(AA) = A¡1A;
which gives A = In. Equivalently, if A2 = A and A 6= In, then A does not
have an inverse.
10. The system of linear equations
x + y ¡ z = a
z = b
2x + y + 2z = c
is equivalent to the matrix equation AX = B, where
A =
2
4
1 1 ¡1
0 0 1
2 1 2
3
5 ; X =
2
4
x
y
z
3
5 ; B =
2
4
a
b
c
3
5 :
By Question 7, A¡1 exists and hence the system has the unique solution
X =
2
4 ¡1 ¡3 1
2 4 ¡1
0 1 0
3
5
2
4
a
b
c
3
5 =
2
4 ¡a ¡ 3b + c
2a + 4b ¡ c
b
3
5 :
Hence x = ¡a ¡ 3b + c; y = 2a + 4b ¡ c; z = b.
20
12.
A = E3(2)E14E42(3) = E3(2)E14
2
664
1 0 0 0
0 1 0 0
0 0 1 0
0 3 0 1
3
775
= E3(2)
2
664
0 3 0 1
0 1 0 0
0 0 1 0
1 0 0 0
3
775
=
2
664
0 3 0 1
0 1 0 0
0 0 2 0
1 0 0 0
3
775
:
Also
A¡1 = (E3(2)E14E42(3))¡1
= (E42(3))¡1E¡1
14 (E3(2))¡1
= E42(¡3)E14E3(1=2)
= E42(¡3)E14
2
664
1 0 0 0
0 1 0 0
0 0 1=2 0
0 0 0 1
3
775
= E42(¡3)
2
664
0 0 0 1
0 1 0 0
0 0 1=2 0
1 0 0 0
3
775
=
2
664
0 0 0 1
0 1 0 0
0 0 1=2 0
1 ¡3 0 0
3
775
:
13. (All matrices in this question are over Z2.)
(a)
2
664
1 1 0 1
0 0 1 1
1 1 1 1
1 0 0 1
¯¯¯¯¯¯¯¯
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
3
775
!
2
664
1 1 0 1
0 0 1 1
0 0 1 0
0 1 0 0
¯¯¯¯¯¯¯¯
1 0 0 0
0 1 0 0
1 0 1 0
1 0 0 1
3
775
!
2
664
1 1 0 1
0 1 0 0
0 0 1 0
0 0 1 1
¯¯¯¯¯¯¯¯
1 0 0 0
1 0 0 1
1 0 1 0
0 1 0 0
3
775
!
2
664
1 0 0 1
0 1 0 0
0 0 1 0
0 0 0 1
¯¯¯¯¯¯¯¯
0 0 0 1
1 0 0 1
1 0 1 0
1 1 1 0
3
775
21
!
2
664
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
¯¯¯¯¯¯¯¯
1 1 1 1
1 0 0 1
1 0 1 0
1 1 1 0
3
775
:
Hence A is non{singular and
A¡1 =
2
664
1 1 1 1
1 0 0 1
1 0 1 0
1 1 1 0
3
775
:
(b) A =
2
664
1 1 0 1
0 1 1 1
1 0 1 0
1 1 0 1
3
775
R4 ! R4 + R1
2
664
1 1 0 1
0 1 1 1
1 0 1 0
0 0 0 0
3
775
, so A is singular.
14.
(a)
2
4
1 1 1
¡1 1 0
2 0 0
¯¯¯¯¯¯
1 0 0
0 1 0
0 0 1
3
5
R3 ! 1
2R3
R1 ! R1 ¡ R3
R2 ! R2 + R3
R1 $ R3
2
4
1 0 0
0 1 0
0 1 1
¯¯¯¯¯¯
0 0 1=2
0 1 1=2
1 0 ¡1=2
3
5
R3 ! R3 ¡ R2
2
4
1 0 0
0 1 0
0 0 1
¯¯¯¯¯¯
0 0 1=2
0 1 1=2
1 ¡1 ¡1
3
5 :
Hence A¡1 exists and
A¡1 =
2
4
0 0 1=2
0 1 1=2
1 ¡1 ¡1
3
5 :
(b)
2
4
2 2 4
1 0 1
0 1 0
¯¯¯¯¯¯
1 0 0
0 1 0
0 0 1
3
5
R1 ! R1 ¡ 2R2
R1 $ R2
R2 $ R3
2
4
1 0 1
0 1 0
0 2 2
¯¯¯¯¯¯
0 1 0
0 0 1
1 ¡2 0
3
5
R3 ! R3 ¡ 2R2
2
4
1 0 1
0 1 0
0 0 2
¯¯¯¯¯¯
0 1 0
0 0 1
1 ¡2 ¡2
3
5
R3 ! 1
2R3
2
4
1 0 1
0 1 0
0 0 1
¯¯¯¯¯¯
0 1 0
0 0 1
1=2 ¡1 ¡1
3
5
22
R1 ! R1 ¡ R3
2
4
1 0 0
0 1 0
0 0 1
¯¯¯¯¯¯
¡1=2 2 1
0 0 1
1=2 ¡1 ¡1
3
5 :
Hence A¡1 exists and
A¡1 =
2
4 ¡1=2 2 1
0 0 1
1=2 ¡1 ¡1
3
5 :
(c)
2
4
4 6 ¡3
0 0 7
0 0 5
3
5 R2 ! 1
7R2
R3 ! 1
5R3
2
4
4 6 ¡3
0 0 1
0 0 1
3
5 R3 ! R3 ¡ R2
2
4
4 6 ¡3
0 0 1
0 0 0
3
5 :
Hence A is singular by virtue of the zero row.
(d)
2
4
2 0 0
0 ¡5 0
0 0 7
¯¯¯¯¯¯
1 0 0
0 1 0
0 0 1
3
5
R1 ! 1
2R1
R2 ! ¡1
5 R2
R3 ! 1
7R3
2
4
1 0 0
0 1 0
0 0 1
¯¯¯¯¯¯
1=2 0 0
0 ¡1=5 0
0 0 1=7
3
5 :
Hence A¡1 exists and A¡1 = diag (1=2; ¡1=5; 1=7).
(Of course this was also immediate from Question 2.)
(e)
2
664
1 2 4 6
0 1 2 0
0 0 1 2
0 0 0 2
¯¯¯¯¯¯¯¯
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
3
775
R1 ! R1 ¡ 2R2
2
664
1 0 0 6
0 1 2 0
0 0 1 2
0 0 0 2
¯¯¯¯¯¯¯¯
1 ¡2 0 0
0 1 0 0
0 0 1 0
0 0 0 1
3
775
R2 ! R2 ¡ 2R3
2
664
1 0 0 6
0 1 0 ¡4
0 0 1 2
0 0 0 2
¯¯¯¯¯¯¯¯
1 ¡2 0 0
0 1 ¡2 0
0 0 1 0
0 0 0 1
3
775
R1 ! R1 ¡ 3R4
R2 ! R2 + 2R4
R3 ! R3 ¡ R4
R4 ! 1
2R4
2
664
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
¯¯¯¯¯¯¯¯
1 ¡2 0 ¡3
0 1 ¡2 2
0 0 1 ¡1
0 0 0 1=2
3
775
:
Hence A¡1 exists and
A¡1 =
2
664
1 ¡2 0 ¡3
0 1 ¡2 2
0 0 1 ¡1
0 0 0 1=2
3
775
:
23
(f)
2
4
1 2 3
4 5 6
5 7 9
3
5 R2 ! R2 ¡ 4R1
R3 ! R3 ¡ 5R1
2
4
1 2 3
0 ¡3 ¡6
0 ¡3 ¡6
3
5 R3 ! R3 ¡ R2
2
4
1 2 3
0 ¡3 ¡6
0 0 0
3
5 :
Hence A is singular by virtue of the zero row.
15. Suppose that A is non{singular. Then
AA¡1 = In = A¡1A:
Taking transposes throughout gives
(AA¡1)t = It
n = (A¡1A)t
(A¡1)tAt = In = At(A¡1)t;
so At is non{singular and (At)¡1 = (A¡1)t.
16. Let A =
·
a b
c d
¸
, where ad ¡ bc = 0. Then the equation
A2 ¡ (a + d)A + (ad ¡ bc)I2 = 0
reduces to A2 ¡ (a + d)A = 0 and hence A2 = (a + d)A. From the last
equation, if A¡1 exists, we deduce that A = (a + d)I2, or
·
a b
c d
¸
=
·
a + d 0
0 a + d
¸
:
Hence a = a + d; b = 0; c = 0; d = a + d and a = b = c = d = 0, which
contradicts the assumption that A is non{singular.
17.
A =
2
4
1 a b
¡a 1 c
¡b ¡c 1
3
5 R2 ! R2 + aR1
R3 ! R3 + bR1
2
4
1 a b
0 1 + a2 c + ab
0 ab ¡ c 1 + b2
3
5
R2 ! 1
1+a2R2
2
4
1 a b
0 1 c+ab
1+a2
0 ab ¡ c 1 + b2
3
5
R3 ! R3 ¡ (ab ¡ c)R2
2
64
1 a b
0 1 c+ab
1+a2
0 0 1 + b2 + (c¡ab)(c+ab)
1+a2
3
75
= B:
24
Now
1 + b2 +
(c ¡ ab)(c + ab)
1 + a2 = 1 + b2 +
c2 ¡ (ab)2
1 + a2
=
1 + a2 + b2 + c2
1 + a2 6= 0:
Hence B can be reduced to I3 using four more row operations and conse-
quently A is non{singular.
18. The proposition is clearly true when n = 1. So let n ¸ 1 and assume
(P¡1AP)n = P¡1AnP. Then
(P¡1AP)n+1 = (P¡1AP)n(P¡1AP)
= (P¡1AnP)(P¡1AP)
= P¡1An(PP¡1)AP
= P¡1AnIAP
= P¡1(AnA)P
= P¡1An+1P
and the induction goes through.
19. Let A =
·
2=3 1=4
1=3 3=4
¸
and P =
·
1 3
¡1 4
¸
. Then P¡1 = 1
7
·
4 ¡3
1 1
¸
.
We then verify that P¡1AP =
·
5=12 0
0 1
¸
. Then from the previous ques-
tion,
P¡1AnP = (P¡1AP)n =
·
5=12 0
0 1
¸n
=
·
(5=12)n 0
0 1n
¸
=
·
(5=12)n 0
0 1
¸
:
Hence
An = P
·
(5=12)n 0
0 1
¸
P¡1 =
·
1 3
¡1 4
¸ ·
(5=12)n 0
0 1
¸
1
7
·
4 ¡3
1 1
¸
=
1
7
·
(5=12)n 3
¡(5=12)n 4
¸ ·
4 ¡3
1 1
¸
=
1
7
·
4(5=12)n + 3 (¡3)(5=12)n + 3
¡4(5=12)n + 4 3(5=12)n + 4
¸
=
1
7
·
3 3
4 4
¸
+
1
7
(5=12)n
·
4 ¡3
¡4 3
¸
:
25
Notice that An ! 1
7
·
3 3
4 4
¸
as n ! 1. This problem is a special case of
a more general result about Markov matrices.
20. Let A =
·
a b
c d
¸
be a matrix whose elements are non{negative real
numbers satisfying
a ¸ 0; b ¸ 0; c ¸ 0; d ¸ 0; a + c = 1 = b + d:
Also let P =
·
b 1
c ¡1
¸
and suppose that A 6= I2.
(i) det P = ¡b ¡ c = ¡(b + c). Now b + c ¸ 0. Also if b + c = 0, then we
would have b = c = 0 and hence d = a = 1, resulting in A = I2. Hence
det P < 0 and P is non{singular.
Next,
P¡1AP = ¡1
b + c
·
¡1 ¡1
¡c b
¸ ·
a b
c d
¸ ·
b 1
c ¡1
¸
= ¡1
b + c
·
¡a ¡ c ¡b ¡ d
¡ac + bc ¡cb + bd
¸ ·
b 1
c ¡1
¸
= ¡1
b + c
·
¡1 ¡1
¡ac + bc ¡cb + bd
¸ ·
b 1
c ¡1
¸
= ¡1
b + c
·
¡b ¡ c 0
(¡ac + bc)b + (¡cb + bd)c ¡ac + bc + cb ¡ bd
¸
:
Now
¡acb + b2c ¡ c2b + bdc = ¡cb(a + c) + bc(b + d)
= ¡cb + bc = 0:
Also
¡(a + d ¡ 1)(b + c) = ¡ab ¡ ac ¡ db ¡ dc + b + c
= ¡ac + b(1 ¡ a) + c(1 ¡ d) ¡ bd
= ¡ac + bc + cb ¡ bd:
Hence
P¡1AP = ¡1
b + c
·
¡(b + c) 0
0 ¡(a + d ¡ 1)(b + c)
¸
=
·
1 0
0 a + d ¡ 1
¸
:
26
(ii) We next prove that if we impose the extra restriction that A 6=
·
0 1
1 0
¸
,
then ja + d ¡ 1j < 1. This will then have the following consequence:
A = P
·
1 0
0 a + d ¡ 1
¸
P¡1
An = P
·
1 0
0 a + d ¡ 1
¸n
P¡1
= P
·
1 0
0 (a + d ¡ 1)n
¸
P¡1
! P
·
1 0
0 0
¸
P¡1
=
·
b 1
c ¡1
¸ ·
1 0
0 0
¸
¡1
b + c
·
¡1 ¡1
¡c b
¸
= ¡1
b + c
·
b 0
c 0
¸ ·
¡1 ¡1
¡c b
¸
= ¡1
b + c
·
¡b ¡b
¡c ¡c
¸
=
1
b + c
·
b b
c c
¸
;
where we have used the fact that (a + d ¡ 1)n ! 0 as n ! 1.
We ¯rst prove the inequality ja + d ¡ 1j · 1:
a + d ¡ 1 · 1 + d ¡ 1 = d · 1
a + d ¡ 1 ¸ 0 + 0 ¡ 1 = ¡1:
Next, if a+d¡1 = 1, we have a+d = 2; so a = 1 = d and hence c = 0 = b,
contradicting our assumption that A 6= I2. Also if a + d ¡ 1 = ¡1, then
a + d = 0; so a = 0 = d and hence c = 1 = b and hence A =
·
0 1
1 0
¸
.
22. The system is inconsistent: We work towards reducing the augmented
matrix:
2
4
1 2
1 1
3 5
¯¯¯¯¯¯
4
5
12
3
5 R2 ! R2 ¡ R1
R3 ! R3 ¡ 3R1
2
4
1 2
0 ¡1
0 ¡1
¯¯¯¯¯¯
4
1
0
3
5
R3 ! R3 ¡ R2
2
4
1 2
0 ¡1
0 0
¯¯¯¯¯¯
4
1
¡1
3
5 :
27
The last row reveals inconsistency.
The system in matrix form is AX = B, where
A =
2
4
1 2
1 1
3 5
3
5 ; X =
·
x
y
¸
; B =
2
4
4
5
12
3
5 :
The normal equations are given by the matrix equation
AtAX = AtB:
Now
AtA =
·
1 1 3
2 1 5
¸ 2
4
1 2
1 1
3 5
3
5 =
·
11 18
18 30
¸
AtB =
·
1 1 3
2 1 5
¸ 2
4
4
5
12
3
5 =
·
45
73
¸
:
Hence the normal equations are
11x + 18y = 45
18x + 30y = 73:
These may be solved, for example, by Cramer's rule:
x =
¯¯¯¯
45 18
73 30
¯¯¯¯
¯¯¯¯
11 18
18 30
¯¯¯¯
=
36
6
= 6
y =
¯¯¯¯
11 45
18 73
¯¯¯¯
¯¯¯¯
11 18
18 30
¯¯¯¯
= ¡7
6
:
23. Substituting the coordinates of the ¯ve points into the parabola equation
gives the following equations:
a = 0
a + b + c = 0
a + 2b + 4c = ¡1
a + 3b + 9c = 4
a + 4b + 16c = 8:
28
The associated normal equations are given by
2
4
5 10 30
10 30 100
30 100 354
3
5
2
4
a
b
c
3
5 =
2
4
11
42
160
3
5 ;
which have the solution a = 1=5; b = ¡2; c = 1.
24. Suppose that A is symmetric, i.e. At = A and that AB is de¯ned. Then
(BtAB)t = BtAt(Bt)t = BtAB;
so BtAB is also symmetric.
25. Let A be m£n and B be n£m, where m > n. Then the homogeneous
system BX = 0 has a non{trivial solution X0, as the number of unknowns
is greater than the number of equations. Then
(AB)X0 = A(BX0) = A0 = 0
and the m £ m matrix AB is therefore singular, as X0 6= 0.
26. (i) Let B be a singular n £ n matrix. Then BX = 0 for some non{zero
column vector X. Then (AB)X = A(BX) = A0 = 0 and hence AB is also
singular.
(ii) Suppose A is a singular n £ n matrix. Then At is also singular and
hence by (i) so is BtAt = (AB)t. Consequently AB is also singular
29
Section 3.6
1. (a) Let S be the set of vectors [x; y] satisfying x = 2y. Then S is a vector
subspace of R2. For
(i) [0; 0] 2 S as x = 2y holds with x = 0 and y = 0.
(ii) S is closed under addition. For let [x1; y1] and [x2; y2] belong to S.
Then x1 = 2y1 and x2 = 2y2. Hence
x1 + x2 = 2y1 + 2y2 = 2(y1 + y2)
and hence
[x1 + x2; y1 + y2] = [x1; y1] + [x2; y2]
belongs to S.
(iii) S is closed under scalar multiplication. For let [x; y] 2 S and t 2 R.
Then x = 2y and hence tx = 2(ty). Consequently
[tx; ty] = t[x; y] 2 S:
(b) Let S be the set of vectors [x; y] satisfying x = 2y and 2x = y. Then S is
a subspace of R2. This can be proved in the same way as (a), or alternatively
we see that x = 2y and 2x = y imply x = 4x and hence x = 0 = y. Hence
S = f[0; 0]g, the set consisting of the zero vector. This is always a subspace.
(c) Let S be the set of vectors [x; y] satisfying x = 2y + 1. Then S doesn't
contain the zero vector and consequently fails to be a vector subspace.
(d) Let S be the set of vectors [x; y] satisfying xy = 0. Then S is not
closed under addition of vectors. For example [1; 0] 2 S and [0; 1] 2 S, but
[1; 0] + [0; 1] = [1; 1] 62 S.
(e) Let S be the set of vectors [x; y] satisfying x ¸ 0 and y ¸ 0. Then S is
not closed under scalar multiplication. For example [1; 0] 2 S and ¡1 2 R,
but (¡1)[1; 0] = [¡1; 0] 62 S.
2. Let X; Y; Z be vectors in Rn. Then by Lemma 3.2.1
hX + Y; X + Z; Y + Zi µ hX; Y; Zi;
as each of X + Y; X + Z; Y + Z is a linear combination of X; Y; Z.
30
Also
X =
1
2
(X + Y ) +
1
2
(X + Z) ¡
1
2
(Y + Z);
Y =
1
2
(X + Y ) ¡
1
2
(X + Z) +
1
2
(Y + Z);
Z = ¡1
2
(X + Y ) +
1
2
(X + Z) +
1
2
(Y + Z);
so
hX; Y; Zi µ hX + Y; X + Z; Y + Zi:
Hence
hX; Y; Zi = hX + Y; X + Z; Y + Zi:
3. Let X1 =
2
664
1
0
1
2
3
775
; X2 =
2
664
0
1
1
2
3
775
and X3 =
2
664
1
1
1
3
3
775
. We have to decide if
X1; X2; X3 are linearly independent, that is if the equation xX1 + yX2 +
zX3 = 0 has only the trivial solution. This equation is equivalent to the
folowing homogeneous system
x + 0y + z = 0
0x + y + z = 0
x + y + z = 0
2x + 2y + 3z = 0:
We reduce the coe±cient matrix to reduced row{echelon form:
2
664
1 0 1
0 1 1
1 1 1
2 2 3
3
775
!
2
664
1 0 0
0 1 0
0 0 1
0 0 0
3
775
and consequently the system has only the trivial solution x = 0; y = 0; z =
0. Hence the given vectors are linearly independent.
4. The vectors
X1 =
2
4
¸
¡1
¡1
3
5 ; X2 =
2
4 ¡1
¸
¡1
3
5 ; X3 =
2
4 ¡1
¡1
¸
3
5
31
are linearly dependent for precisely those values of ¸ for which the equation
xX1+yX2+zX3 = 0 has a non{trivial solution. This equation is equivalent
to the system of homogeneous equations
¸x ¡ y ¡ z = 0
¡x + ¸y ¡ z = 0
¡x ¡ y + ¸z = 0:
Now the coe±cient determinant of this system is
¯¯¯¯¯¯
¸ ¡1 ¡1
¡1 ¸ ¡1
¡1 ¡1 ¸
¯¯¯¯¯¯
= (¸ + 1)2(¸ ¡ 2):
So the values of ¸ which make X1; X2; X3 linearly independent are those ¸
satisfying ¸ 6= ¡1 and ¸ 6= 2.
5. Let A be the following matrix of rationals:
A =
2
664
1 1 2 0 1
2 2 5 0 3
0 0 0 1 3
8 11 19 0 11
3
775
:
Then A has reduced row{echelon form
B =
2
664
1 0 0 0 ¡1
0 1 0 0 0
0 0 1 0 1
0 0 0 1 3
3
775
:
From B we read o® the following:
(a) The rows of B form a basis for R(A). (Consequently the rows of A
also form a basis for R(A).)
(b) The ¯rst four columns of A form a basis for C(A).
(c) To ¯nd a basis for N(A), we solve AX = 0 and equivalently BX = 0.
From B we see that the solution is
x1 = x5
x2 = 0
x3 = ¡x5
x4 = ¡3x5;
32
with x5 arbitrary. Then
X =
2
66664
x5
0
¡x5
¡3x5
x5
3
77775
= x5
2
66664
1
0
¡1
¡3
1
3
77775
;
so [1; 0; ¡1; ¡3; 1]t is a basis for N(A).
6. In Section 1.6, problem 12, we found that the matrix
A =
2
664
1 0 1 0 1
0 1 0 1 1
1 1 1 1 0
0 0 1 1 0
3
775
has reduced row{echelon form
B =
2
664
1 0 0 1 1
0 1 0 1 1
0 0 1 1 0
0 0 0 0 0
3
775
:
From B we read o® the following:
(a) The three non{zero rows of B form a basis for R(A).
(b) The ¯rst three columns of A form a basis for C(A).
(c) To ¯nd a basis for N(A), we solve AX = 0 and equivalently BX = 0.
From B we see that the solution is
x1 = ¡x4 ¡ x5 = x4 + x5
x2 = ¡x4 ¡ x5 = x4 + x5
x3 = ¡x4 = x4;
with x4 and x5 arbitrary elements of Z2. Hence
X =
2
66664
x4 + x5
x4 + x5
x4
x4
x5
3
77775
= x4
2
66664
1
1
1
1
0
3
77775
+
2
66664
1
1
0
0
1
3
77775
:
Hence [1; 1; 1; 1; 0]t and [1; 1; 0; 0; 1]t form a basis for N(A).
33
7. Let A be the following matrix over Z5:
A =
2
664
1 1 2 0 1 3
2 1 4 0 3 2
0 0 0 1 3 0
3 0 2 4 3 2
3
775
:
We ¯nd that A has reduced row{echelon form B:
B =
2
664
1 0 0 0 2 4
0 1 0 0 4 4
0 0 1 0 0 0
0 0 0 1 3 0
3
775
:
From B we read o® the following:
(a) The four rows of B form a basis for R(A). (Consequently the rows of
A also form a basis for R(A).
(b) The ¯rst four columns of A form a basis for C(A).
(c) To ¯nd a basis for N(A), we solve AX = 0 and equivalently BX = 0.
From B we see that the solution is
x1 = ¡2x5 ¡ 4x6 = 3x5 + x6
x2 = ¡4x5 ¡ 4x6 = x5 + x6
x3 = 0
x4 = ¡3x5 = 2x5;
where x5 and x6 are arbitrary elements of Z5. Hence
X = x5
2
6666664
3
1
0
2
1
0
3
7777775 +
x
6
2
6666664 1
1
0
0
0
1
3
7777775
;
so [3; 1; 0; 2; 1; 0]t and [1; 1; 0; 0; 0; 1]t form a basis for R(A).
8. Let F = f0; 1; a; bg be a ¯eld and let A be the following matrix over F:
A =
2
4
1 a b a
a b b 1
1 1 1 a
3
5 :
34
In Section 1.6, problem 17, we found that A had reduced row{echelon form
B =
2
4
1 0 0 0
0 1 0 b
0 0 1 1
3
5 :
From B we read o® the following:
(a) The rows of B form a basis for R(A). (Consequently the rows of A
also form a basis for R(A).
(b) The ¯rst three columns of A form a basis for C(A).
(c) To ¯nd a basis for N(A), we solve AX = 0 and equivalently BX = 0.
From B we see that the solution is
x1 = 0
x2 = ¡bx4 = bx4
x3 = ¡x4 = x4;
where x4 is an arbitrary element of F. Hence
X = x4
2
664
0
b
1
1
3
775
;
so [0; b; 1; 1]t is a basis for N(A).
9. Suppose that X1; : : : ;Xm form a basis for a subspace S. We have to
prove that
X1; X1 + X2; : : : ;X1 + ¢ ¢ ¢ + Xm
also form a basis for S.
First we prove the independence of the family: Suppose
x1X1 + x2(X1 + X2) + ¢ ¢ ¢ + xm(X1 + ¢ ¢ ¢ + Xm) = 0:
Then
(x1 + x2 + ¢ ¢ ¢ + xm)X1 + ¢ ¢ ¢ + xmXm = 0:
Then the linear independence of X1; : : : ;Xm gives
x1 + x2 + ¢ ¢ ¢ + xm = 0; : : : ; xm = 0;
35
form which we deduce that x1 = 0; : : : ; xm = 0.
Secondly we have to prove that every vector of S is expressible as a linear
combination of X1; X1 + X2; : : : ;X1 + ¢ ¢ ¢ + Xm. Suppose X 2 S. Then
X = a1X1 + ¢ ¢ ¢ + amXm:
We have to ¯nd x1; : : : ; xm such that
X = x1X1 + x2(X1 + X2) + ¢ ¢ ¢ + xm(X1 + ¢ ¢ ¢ + Xm)
= (x1 + x2 + ¢ ¢ ¢ + xm)X1 + ¢ ¢ ¢ + xmXm:
Then
a1X1 + ¢ ¢ ¢ + amXm = (x1 + x2 + ¢ ¢ ¢ + xm)X1 + ¢ ¢ ¢ + xmXm:
So if we can solve the system
x1 + x2 + ¢ ¢ ¢ + xm = a1; : : : ; xm = am;
we are ¯nished. Clearly these equations have the unique solution
x1 = a1 ¡ a2; : : : ; xm¡1 = am ¡ am¡1; xm = am:
10. Let A =
·
a b c
1 1 1
¸
. If [a; b; c] is a multiple of [1; 1; 1], (that is,
a = b = c), then rankA = 1. For if
[a; b; c] = t[1; 1; 1];
then
R(A) = h[a; b; c]; [1; 1; 1]i = ht[1; 1; 1]; [1; 1; 1]i = h[1; 1; 1]i;
so [1; 1; 1] is a basis for R(A).
However if [a; b; c] is not a multiple of [1; 1; 1], (that is at least two
of a; b; c are distinct), then the left{to{right test shows that [a; b; c] and
[1; 1; 1] are linearly independent and hence form a basis for R(A). Conse-
quently rankA = 2 in this case.
11. Let S be a subspace of Fn with dim S = m. Also suppose that
X1; : : : ;Xm are vectors in S such that S = hX1; : : : ;Xmi. We have to
prove that X1; : : : ;Xm form a basis for S; in other words, we must prove
that X1; : : : ;Xm are linearly independent.
36
However if X1; : : : ;Xm were linearly dependent, then one of these vec-
tors would be a linear combination of the remaining vectors. Consequently
S would be spanned by m ¡ 1 vectors. But there exist a family of m lin-
early independent vectors in S. Then by Theorem 3.3.2, we would have the
contradiction m · m ¡ 1.
12. Let [x; y; z]t 2 S. Then x + 2y + 3z = 0. Hence x = ¡2y ¡ 3z and
2
4
x
y
z
3
5 =
2
4 ¡2y ¡ 3z
y
z
3
5 = y
2
4 ¡2
1
0
3
5 + z
2
4 ¡3
0
1
3
5 :
Hence [¡2; 1; 0]t and [¡3; 0; 1]t form a basis for S.
Next (¡1) + 2(¡1) + 3(1) = 0, so [¡1; ¡1; 1]t 2 S.
To ¯nd a basis for S which includes [¡1; ¡1; 1]t, we note that [¡2; 1; 0]t
is not a multiple of [¡1; ¡1; 1]t. Hence we have found a linearly independent
family of two vectors in S, a subspace of dimension equal to 2. Consequently
these two vectors form a basis for S.
13. Without loss of generality, suppose that X1 = X2. Then we have the
non{trivial dependency relation:
1X1 + (¡1)X2 + 0X3 + ¢ ¢ ¢ + 0Xm = 0:
14. (a) Suppose that Xm+1 is a linear combination of X1; : : : ;Xm. Then
hX1; : : : ;Xm; Xm+1i = hX1; : : : ;Xmi
and hence
dim hX1; : : : ;Xm; Xm+1i = dim hX1; : : : ;Xmi:
(b) Suppose that Xm+1 is not a linear combination of X1; : : : ;Xm. If not
all of X1; : : : ;Xm are zero, there will be a subfamily Xc1 ; : : : ;Xcr which is
a basis for hX1; : : : ;Xmi.
Then as Xm+1 is not a linear combination of Xc1 ; : : : ;Xcr , it follows that
Xc1 ; : : : ;Xcr ; Xm+1 are linearly independent. Also
hX1; : : : ;Xm; Xm+1i = hXc1 ; : : : ;Xcr ; Xm+1i:
Consequently
dim hX1; : : : ;Xm; Xm+1i = r + 1 = dim hX1; : : : ;Xmi + 1:
37
Our result can be rephrased in a form suitable for the second part of the
problem:
dim hX1; : : : ;Xm; Xm+1i = dim hX1; : : : ;Xmi
if and only if Xm+1 is a linear combination of X1; : : : ;Xm.
If X = [x1; : : : ; xn]t, then AX = B is equivalent to
B = x1A¤1 + ¢ ¢ ¢ + xnA¤n:
So AX = B is soluble for X if and only if B is a linear combination of the
columns of A, that is B 2 C(A). However by the ¯rst part of this question,
B 2 C(A) if and only if dimC([AjB]) = dimC(A), that is, rank [AjB] =
rankA.
15. Let a1; : : : ; an be elements of F, not all zero. Let S denote the set of
vectors [x1; : : : ; xn]t, where x1; : : : ; xn satisfy
a1x1 + ¢ ¢ ¢ + anxn = 0:
Then S = N(A), where A is the row matrix [a1; : : : ; an]. Now rankA = 1
as A 6= 0. So by the \rank + nullity" theorem, noting that the number of
columns of A equals n, we have
dimN(A) = nullity (A) = n ¡ rankA = n ¡ 1:
16. (a) (Proof of Lemma 3.2.1) Suppose that each of X1; : : : ;Xr is a linear
combination of Y1; : : : ; Ys. Then
Xi =
Xs
j=1
aijYj ; (1 · i · r):
Now let X =
P r
i=1 xiXi be a linear combination of X1; : : : ;Xr. Then
X = x1(a11Y1 + ¢ ¢ ¢ + a1sYs)
+ ¢ ¢ ¢
+ xr(ar1Y1 + ¢ ¢ ¢ + arsYs)
= y1Y1 + ¢ ¢ ¢ + ysYs;
where yj = a1jx1+¢ ¢ ¢+arjxr. Hence X is a linear combination of Y1; : : : ; Ys.
Another way of stating Lemma 3.2.1 is
hX1; : : : ;Xri µ hY1; : : : ; Ysi; (1)
38
if each of X1; : : : ;Xr is a linear combination of Y1; : : : ; Ys.
(b) (Proof of Theorem 3.2.1) Suppose that each of X1; : : : ;Xr is a linear
combination of Y1; : : : ; Ys and that each of Y1; : : : ; Ys is a linear combination
of X1; : : : ;Xr. Then by (a) equation (1) above
hX1; : : : ;Xri µ hY1; : : : ; Ysi
and
hY1; : : : ; Ysi µ hX1; : : : ;Xri:
Hence
hX1; : : : ;Xri = hY1; : : : ; Ysi:
(c) (Proof of Corollary 3.2.1) Suppose that each of Z1; : : : ;Zt is a linear
combination of X1; : : : ;Xr. Then each of X1; : : : ;Xr; Z1; : : : ;Zt is a linear
combination of X1; : : : ;Xr.
Also each of X1; : : : ;Xr is a linear combination of X1; : : : ;Xr; Z1; : : : ;Zt,
so by Theorem 3.2.1
hX1; : : : ;Xr; Z1; : : : ;Zti = hX1; : : : ;Xri:
(d) (Proof of Theorem 3.3.2) Let Y1; : : : ; Ys be vectors in hX1; : : : ;Xri
and assume that s > r. We have to prove that Y1; : : : ; Ys are linearly
dependent. So we consider the equation
x1Y1 + ¢ ¢ ¢ + xsYs = 0:
Now Yi =
P r
j=1 aijXj , for 1 · i · s. Hence
x1Y1 + ¢ ¢ ¢ + xsYs = x1(a11X1 + ¢ ¢ ¢ + a1rXr)
+ ¢ ¢ ¢
+ xr(as1X1 + ¢ ¢ ¢ + asrXr):
= y1X1 + ¢ ¢ ¢ + yrXr; (1)
where yj = a1jx1 + ¢ ¢ ¢ + asjxs. However the homogeneous system
y1 = 0; ¢ ¢ ¢ ; yr = 0
has a non{trivial solution x1; : : : ; xs, as s > r and from (1), this results in a
non{trivial solution of the equation
x1Y1 + ¢ ¢ ¢ + xsYs = 0:
39
Hence Y1; : : : ; Ys are linearly dependent.
17. Let R and S be subspaces of Fn, with R µ S. We ¯rst prove
dimR · dim S:
Let X1; : : : ;Xr be a basis for R. Now by Theorem 3.5.2, because X1; : : : ;Xr
form a linearly independent family lying in S, this family can be extended
to a basis X1; : : : ;Xr; : : : ;Xs for S. Then
dim S = s ¸ r = dimR:
Next suppose that dimR = dim S. Let X1; : : : ;Xr be a basis for R. Then
because X1; : : : ;Xr form a linearly independent family in S and S is a sub-
space whose dimension is r, it follows from Theorem 3.4.3 that X1; : : : ;Xr
form a basis for S. Then
S = hX1; : : : ;Xri = R:
18. Suppose that R and S are subspaces of Fn with the property that R[S
is also a subspace of Fn. We have to prove that R µ S or S µ R. We argue
by contradiction: Suppose that R 6µ S and S 6µ R. Then there exist vectors
u and v such that
u 2 R and v 62 S; v 2 S and v 62 R:
Consider the vector u+v. As we are assuming R[S is a subspace, R[S is
closed under addition. Hence u + v 2 R [ S and so u + v 2 R or u + v 2 S.
However if u + v 2 R, then v = (u + v) ¡ u 2 R, which is a contradiction;
similarly if u + v 2 S.
Hence we have derived a contradiction on the asumption that R 6µ S and
S 6µ R. Consequently at least one of these must be false. In other words
R µ S or S µ R.
19. Let X1; : : : ;Xr be a basis for S.
(i) First let
Y1 = a11X1 + ¢ ¢ ¢ + a1rXr
...
(2)
Yr = ar1X1 + ¢ ¢ ¢ + arrXr;
40
where A = [aij ] is non{singular. Then the above system of equations can
be solved for X1; : : : ;Xr in terms of Y1; : : : ; Yr. Consequently by Theorem
3.2.1
hY1; : : : ; Yri = hX1; : : : ;Xri = S:
It follows from problem 11 that Y1; : : : ; Yr is a basis for S.
(ii) We show that all bases for S are given by equations 2. So suppose
that Y1; : : : ; Yr forms a basis for S. Then because X1; : : : ;Xr form a basis
for S, we can express Y1; : : : ; Yr in terms of X1; : : : ;Xr as in 2, for some
matrix A = [aij ]. We show A is non{singular by demonstrating that the
linear independence of Y1; : : : ; Yr implies that the rows of A are linearly
independent.
So assume
x1[a11; : : : ; a1r] + ¢ ¢ ¢ + xr[ar1; : : : ; arr] = [0; : : : ; 0]:
Then on equating components, we have
a11x1 + ¢ ¢ ¢ + ar1xr = 0
...
a1rx1 + ¢ ¢ ¢ + arrxr = 0:
Hence
x1Y1 + ¢ ¢ ¢ + xrYr = x1(a11X1 + ¢ ¢ ¢ + a1rXr) + ¢ ¢ ¢ + xr(ar1X1 + ¢ ¢ ¢ + arrXr)
= (a11x1 + ¢ ¢ ¢ + ar1xr)X1 + ¢ ¢ ¢ + (a1rx1 + ¢ ¢ ¢ + arrxr)Xr
= 0X1 + ¢ ¢ ¢ + 0Xr = 0:
Then the linear independence of Y1; : : : ; Yr implies x1 = 0; : : : ; xr = 0.
(We mention that the last argument is reversible and provides an alter-
native proof of part (i).)
41
P2
P1
P3
O
¡
¡
¡
¡
¡
¡
¡
¡
¡
¡
¡
¡
©©©©©©©©
¢
¢
¢
¢
¢
¢
¢
¢@
@
@
@
©©©©©©©©
¢
¢
¢
¢
¢
¢
¢
¢
Section 4.1
1. We ¯rst prove that the area of a triangle P1P2P3, where the points
are in anti{clockwise orientation, is given by the formula
1
2
½¯¯¯¯
x1 x2
y1 y2
¯¯¯¯
+
¯¯¯¯
x2 x3
y2 y3
¯¯¯¯
+
¯¯¯¯
x3 x1
y3 y1
¯¯¯¯
¾
:
Referring to the above diagram, we have
Area P1P2P3 = AreaOP1P2 + AreaOP2P3 ¡ AreaOP1P3
=
1
2
¯¯¯¯
x1 x2
y1 y2
¯¯¯¯
+
1
2
¯¯¯¯
x2 x3
y2 y3
¯¯¯¯
¡
1
2
¯¯¯¯
x1 x3
y1 y3
¯¯¯¯
;
which gives the desired formula.
We now turn to the area of a quadrilateral. One possible con¯guration
occurs when the quadrilateral is convex as in ¯gure (a) below. The interior
diagonal breaks the quadrilateral into two triangles P1P2P3 and P1P3P4.
Then
Area P1P2P3P4 = Area P1P2P3 + Area P1P3P4
=
1
2
½¯¯¯¯
x1 x2
y1 y2
¯¯¯¯
+
¯¯¯¯
x2 x3
y2 y3
¯¯¯¯
+
¯¯¯¯
x3 x1
y3 y1
¯¯¯¯
¾
42
¶
¶
¶
¶
¶
¶
¶
¶
` ` ` ` ` ` ` ` ` `
"
"
"
"
"
"
"
"
"
"
HH HH
P1
P2
P3
P4
(a)
` ` ` ` ` ` ` ` ` `
©©©©©©
L
L
L
L
L
L
L
L
L
L
\
\
\
\
\
P1
P2
P3
(b) P4
+
1
2
½¯¯¯¯
x1 x3
y1 y3
¯¯¯¯
+
¯¯¯¯
x3 x4
y3 y4
¯¯¯¯
+
¯¯¯¯
x4 x1
y4 y1
¯¯¯¯
¾
=
1
2
½¯¯¯¯
x1 x2
y1 y2
¯¯¯¯
+
¯¯¯¯
x2 x3
y2 y3
¯¯¯¯
+
¯¯¯¯
x3 x4
y3 y4
¯¯¯¯
+
¯¯¯¯
x4 x1
y4 y1
¯¯¯¯
¾
;
after cancellation.
Another possible con¯guration for the quadrilateral occurs when it is not
convex, as in ¯gure (b). The interior diagonal P2P4 then gives two triangles
P1P2P4 and P2P3P4 and we can proceed similarly as before.
2.
¢ =
¯¯¯¯¯¯
a + x b + y c + z
x + u y + v z + w
u + a v + b w + c
¯¯¯¯¯¯
=
¯¯¯¯¯¯
a b c
x + u y + v z + w
u + a v + b w + c
¯¯¯¯¯¯
+
¯¯¯¯¯¯
x y z
x + u y + v z + w
u + a v + b w + c
¯¯¯¯¯¯
:
Now
¯¯¯¯¯¯
a b c
x + u y + v z + w
u + a v + b w + c
¯¯¯¯¯¯
=
¯¯¯¯¯¯
a b c
x y z
u + a v + b w + c
¯¯¯¯¯¯+
¯¯¯¯¯¯ a
b
c
u v w
u + a v + b w + c
¯¯¯¯¯¯
=
¯¯¯¯¯¯
a b c
x y z
u v w
¯¯¯¯¯¯
+
¯¯¯¯¯¯
a b c
x y z
a b c
¯¯¯¯¯¯+
¯¯¯¯¯¯
a b c
u v w
u v w
¯¯¯¯¯¯
¯¯¯¯¯¯
a b c
u v w
a b c
¯¯¯¯¯¯
=
¯¯¯¯¯¯
a b c
x y z
u v w
¯¯¯¯¯¯
:
Similarly
¯¯¯¯¯¯
x y z
x + u y + v z + w
u + a v + b w + c
¯¯¯¯¯¯
=
¯¯¯¯¯¯
x y z
u v w
a b c
¯¯¯¯¯¯
= ¡
¯¯¯¯¯¯
x y z
a b c
u v w
¯¯¯¯¯¯
=
¯¯¯¯¯¯
a b c
x y z
u v w
¯¯¯¯¯¯
:
43
Hence ¢ = 2
¯¯¯¯¯¯
a b c
x y z
u v w
¯¯¯¯¯¯
.
3:
¯¯¯¯¯¯
n2 (n + 1)2 (n + 2)2
(n + 1)2 (n + 2)2 (n + 3)2
(n + 2)2 (n + 3)2 (n + 4)2
¯¯¯¯¯¯
C3 ! C3 ¡ C2
C2 ! C2 ¡ C1
=
¯¯¯¯¯¯
n2 2n + 1 2n + 3
(n + 1)2 2n + 3 2n + 5
(n + 2)2 2n + 5 2n + 7
¯¯¯¯¯¯
C3 ! C3 ¡ C2
=
¯¯¯¯¯¯
n2 2n + 1 2
(n + 1)2 2n + 3 2
(n + 2)2 2n + 5 2
¯¯¯¯¯¯
R3 ! R3 ¡ R2
R2 ! R2 ¡ R1
=
¯¯¯¯¯¯
n2 2n + 1 2
2n + 1 2 0
2n + 3 2 0
¯¯¯¯¯¯
= ¡8:
4. (a)
¯¯¯¯¯¯
246 427 327
1014 543 443
¡342 721 621
¯¯¯¯¯¯
=
¯¯¯¯¯¯
246 100 327
1014 100 443
¡342 100 621
¯¯¯¯¯¯
= 100
¯¯¯¯¯¯
246 1 327
1014 1 443
¡342 1 621
¯¯¯¯¯¯
= 100
¯¯¯¯¯¯
246 1 327
768 0 116
¡588 0 294
¯¯¯¯¯¯
= 100(¡1)
¯¯¯¯
768 116
¡588 294
¯¯¯¯
= ¡29400000:
(b)
¯¯¯¯¯¯¯¯
1 2 3 4
¡2 1 ¡4 3
3 ¡4 ¡1 2
4 3 ¡2 ¡1
¯¯¯¯¯¯¯¯
=
¯¯¯¯¯¯¯¯
1 2 3 4
0 5 2 11
0 ¡10 ¡10 ¡10
0 ¡5 ¡14 ¡17
¯¯¯¯¯¯¯¯
=
¯¯¯¯¯¯
5 2 11
¡10 ¡10 ¡10
¡5 ¡14 ¡17
¯¯¯¯¯¯
= ¡10
¯¯¯¯¯¯
5 2 11
1 1 1
¡5 ¡14 ¡17
¯¯¯¯¯¯
= ¡10
¯¯¯¯¯¯
5 ¡3 6
1 0 0
¡5 ¡9 ¡12
¯¯¯¯¯¯
= ¡10(¡1)
¯¯¯¯
¡3 6
¡9 ¡12
¯¯¯¯
= 900:
5. detA =
¯¯¯¯¯¯
1 0 ¡2
3 1 4
5 2 ¡3
¯¯¯¯¯¯
=
¯¯¯¯¯¯
1 0 0
3 1 10
5 2 7
¯¯¯¯¯¯
=
¯¯¯¯
1 10
2 7
¯¯¯¯
= ¡13.
44
Hence A is non{singular and
A¡1 =
1
¡13
adjA =
1
¡13
2
4
C11 C21 C31
C12 C22 C32
C13 C23 C33
3
5 =
1
¡13
2
4 ¡11 ¡4 2
29 7 ¡10
1 ¡2 1
3
5 :
6. (i)
¯¯¯¯¯¯
2a 2b b ¡ c
2b 2a a + c
a + b a + b b
¯¯¯¯¯¯
R1 ! R1 + R2
=
¯¯¯¯¯¯
2a + 2b 2b + 2a b + a
2b 2a a + c
a + b a + b b
¯¯¯¯¯¯
= (a+b)
¯¯¯¯¯¯
2 2 1
2b 2a a + c
a + b a + b b
¯¯¯¯¯¯
C1 ! C1 ¡ C2
=
(a+b)
¯¯¯¯¯¯
0 2 1
2(b ¡ a) 2a a + c
0 a + b b
¯¯¯¯¯¯
= 2(a + b)(a ¡ b)
¯¯¯¯
2 1
a + b b
¯¯¯¯
= ¡2(a + b)(a ¡ b)2:
(ii)
¯¯¯¯¯¯
b + c b c
c c + a a
b a a + b
¯¯¯¯¯¯
C1 ! C1 ¡ C2
=
¯¯¯¯¯¯
c b c
¡a c + a a
b ¡ a a a + b
¯¯¯¯¯¯
C3 ! C3 ¡ C1
=
¯¯¯¯¯¯
c b 0
¡a c + a 2a
b ¡ a a 2a
¯¯¯¯¯¯
= 2a
¯¯¯¯¯¯
c b 0
¡a c + a 1
b ¡ a a 1
¯¯¯¯¯¯
R3 ! R3 ¡ R2
=
2a
¯¯¯¯¯¯
c b 0
¡a c + a 1
b ¡c 0
¯¯¯¯¯¯
= ¡2a
¯¯¯¯
c b
b ¡c
¯¯¯¯
= 2a(c2 + b2):
7. Suppose that the curve y = ax2 + bx + c passes through the points
(x1; y1); (x2; y2); (x3; y3), where xi 6= xj if i 6= j. Then
ax2
1 + bx1 + c = y1
ax2
2 + bx2 + c = y2
ax2
3 + bx3 + c = y3:
The coe±cient determinant is essentially a Vandermonde determinant:
¯¯¯¯¯¯
x2
1 x1 1
x2
2 x2 1
x2
3 x3 1
¯¯¯¯¯¯
=
¯¯¯¯¯¯ x
21
x2
2 x2
3
x1 x2 x3
1 1 1
¯¯¯¯¯¯
= ¡
¯¯¯¯¯¯
1 1 1
x1 x2 x3
x2
1 x2
2 x2
3
¯¯¯¯¯¯
= ¡(x2¡x1)(x3¡x1)(x3¡x2):
45
Hence the coe±cient determinant is non{zero and by Cramer's rule, there
is a unique solution for a; b; c.
8. Let ¢ = detA =
¯¯¯¯¯¯
1 1 ¡1
2 3 k
1 k 3
¯¯¯¯¯¯
. Then
¢ =
C3 ! C3 + C1
C2 ! C2 ¡ C1
¯¯¯¯¯¯
1 0 0
2 1 k + 2
1 k ¡ 1 4
¯¯¯¯¯¯
=
¯¯¯¯
1 k + 2
k ¡ 1 4
¯¯¯¯
= 4 ¡ (k ¡ 1)(k + 2) = ¡(k2 ¡ k ¡ 6) = ¡(k + 3)(k ¡ 2):
Hence detA = 0 if and only if k = ¡3 or k = 2.
Consequently if k 6= ¡3 and k 6= 2, then detA 6= 0 and the given system
x + y ¡ z = 1
2x + 3y + kz = 3
x + ky + 3z = 2
has a unique solution. We consider the cases k = ¡3 and k = 2 separately.
k = ¡3 :
AM =
2
4
1 1 ¡1 1
2 3 ¡3 3
1 ¡3 3 2
3
5 R2 ! R2 ¡ 2R1
R3 ! R3 ¡ R1
2
4
1 1 ¡1 1
0 1 ¡1 1
0 ¡4 4 1
3
5
R3 ! R3 + 4R2
2
4
1 1 ¡1 1
0 1 ¡1 1
0 0 0 5
3
5 ;
from which we read o® inconsistency.
k = 2 :
AM =
2
4
1 1 ¡1 1
2 3 2 3
1 2 3 2
3
5 R2 ! R2 ¡ 2R1
R3 ! R3 ¡ R1
2
4
1 1 ¡1 1
0 1 4 1
0 1 4 1
3
5
R3 ! R3 ¡ R2
2
4
1 0 ¡5 0
0 1 4 1
0 0 0 0
3
5 :
We read o® the complete solution x = 5z; y = 1 ¡ 4z, where z is arbitrary.
46
Finally we have to determine the solution for which x2 +y2 +z2 is least.
x2 + y2 + z2 = (5z)2 + (1 ¡ 4z)2 + z2 = 42z2 ¡ 8z + 1
= 42(z2 ¡
4
21
z +
1
42
) = 42
(µ
z ¡
2
21
¶2
+
1
42 ¡
µ
2
21
¶2
)
= 42
(µ
z ¡
2
21
¶2
+
13
882
)
:
We see that the least value of x2+y2+z2 is 42£ 13
882 = 13
21 and this occurs when
z = 2=21, with corresponding values x = 10=21 and y = 1¡4£ 2
21 = 13=21.
9. Let ¢ =
2
4
1 ¡2 b
a 0 2
5 2 0
¯¯¯¯¯¯
be the coe±cient determinant of the given system.
Then expanding along column 2 gives
¢ = 2
¯¯¯¯
a 2
5 0
¯¯¯¯
¡ 2
¯¯¯¯
1 b
a 2
¯¯¯¯
= ¡20 ¡ 2(2 ¡ ab)
= 2ab ¡ 24 = 2(ab ¡ 12):
Hence ¢ = 0 if and only if ab = 12. Hence if ab 6= 12, the given system has
a unique solution.
If ab = 12 we must argue with care:
AM =
2
4
1 ¡2 b 3
a 0 2 2
5 2 0 1
3
5 !
2
4
1 ¡2 b 3
0 2a 2 ¡ ab 2 ¡ 3a
0 12 ¡5b ¡14
3
5
!
2
4
1 ¡2 b 3
0 1 ¡5b
12 ¡7
6
0 2a 2 ¡ ab 2 ¡ 3a
3
5 !
2
4
1 ¡2 b 3
0 1 ¡5b
12 ¡7
6
0 0 12¡ab
6
6¡2a
3
3
5
=
2
4
1 ¡2 b 3
0 1 ¡5b
12 ¡7
6
0 0 0 6¡2a
3
3
5 = B:
Hence if 6 ¡ 2a 6= 0, i.e. a 6= 3, the system has no solution.
If a = 3 (and hence b = 4), then
B =
2
4
1 ¡2 4 3
0 1 ¡5
3 ¡7
6
0 0 0 0
3
5 !
2
4
1 0 ¡2=3 2=3
0 1 ¡5
3 ¡7
6
0 0 0 0
3
5 :
47
Consequently the complete solution of the system is x = 2
3+2
3z; y = ¡7
6 +5
3z,
where z is arbitrary. Hence there are in¯nitely many solutions.
10.
¢ =
¯¯¯¯¯¯¯¯
1 1 2 1
1 2 3 4
2 4 7 2t + 6
2 2 6 ¡ t t
¯¯¯¯¯¯¯¯
R4 ! R4 ¡ 2R1
R3 ! R3 ¡ 2R1
R2 ! R2 ¡ R1
=
¯¯¯¯¯¯¯¯
1 1 2 1
0 1 1 3
0 2 3 2t + 4
0 0 2 ¡ t t ¡ 2
¯¯¯¯¯¯¯¯
=
¯¯¯¯¯¯
1 1 3
2 3 2t + 4
0 2 ¡ t t ¡ 2
¯¯¯¯¯¯
R2 ! R2 ¡ 2R1
=
¯¯¯¯¯¯
1 1 3
0 1 2t ¡ 2
0 2 ¡ t t ¡ 2
¯¯¯¯¯¯
=
¯¯¯¯
1 2t ¡ 2
2 ¡ t t ¡ 2
¯¯¯¯
= (t ¡ 2)
¯¯¯¯
1 2t ¡ 2
¡1 1
¯¯¯¯
= (t ¡ 2)(2t ¡ 1):
Hence ¢ = 0 if and only if t = 2 or t = 1
2 . Consequently the given matrix
B is non{singular if and only if t 6= 2 and t 6= 1
2 .
11. Let A be a 3 £ 3 matrix with detA 6= 0. Then
(i)
AadjA = (detA)I3 (1)
(detA) det ( adjA) = det (detA ¢ I3) = (detA)3:
Hence, as detA 6= 0, dividing out by detA in the last equation gives
det ( adjA) = (detA)2:
(ii) . Also from equation (1)
µ
1
detA
A
¶
adjA = I3;
so adjA is non{singular and
( adjA)¡1 =
1
detA
A:
Finally
A¡1 adj (A¡1) = (detA¡1)I3
and multiplying both sides of the last equation by A gives
adj (A¡1) = A(detA¡1)I3 =
1
detA
A:
48
12. Let A be a real 3 £ 3 matrix satisfying AtA = I3. Then
(i) At(A ¡ I3) = AtA ¡ At = I3 ¡ At
= ¡(At ¡ I3) = ¡(At ¡ It
3 ) = ¡(A ¡ I3)t:
Taking determinants of both sides then gives
detAt det (A ¡ I3) = det (¡(A ¡ I3)t)
detAdet (A ¡ I3) = (¡1)3 det (A ¡ I3)t
= ¡det (A ¡ I3) (1):
(ii) Also detAAt = det I3, so
detAt detA = 1 = (detA)2:
Hence detA = §1.
(iii) Suppose that detA = 1. Then equation (1) gives
det (A ¡ I3) = ¡det (A ¡ I3);
so (1 + 1) det (A ¡ I3) = 0 and hence det (A ¡ I3) = 0.
13. Suppose that column 1 is a linear combination of the remaining columns:
A¤1 = x2A¤2 + ¢ ¢ ¢ + xnA¤n:
Then
detA =
¯¯¯¯¯¯¯¯¯
x2a12 + ¢ ¢ ¢ + xna1n a12 ¢ ¢ ¢ a1n
x2a22 + ¢ ¢ ¢ + xna2n a22 ¢ ¢ ¢ a2n
...
...
...
...
x2an2 + ¢ ¢ ¢ + xnann an2 ¢ ¢ ¢ ann
¯¯¯¯¯¯¯¯¯
:
Now detA is unchanged in value if we perform the operation
C1 ! C1 ¡ x2C2 ¡ ¢ ¢ ¢ ¡ xnCn :
detA =
¯¯¯¯¯¯¯¯¯
0 a12 ¢ ¢ ¢ a1n
0 a22 ¢ ¢ ¢ a2n
...
...
...
...
0 an2 ¢ ¢ ¢ ann
¯¯¯¯¯¯¯¯¯
= 0:
49
Conversely, suppose that detA = 0. Then the homogeneous system AX = 0
has a non{trivial solution X = [x1; : : : ; xn]t. So
x1A¤1 + ¢ ¢ ¢ + xnA¤n = 0:
Suppose for example that x1 6= 0. Then
A¤1 =
µ
¡
x2
x1
¶
+ ¢ ¢ ¢ +
µ
¡
xn
x1
¶
A¤n
and the ¯rst column of A is a linear combination of the remaining columns.
14. Consider the system
¡2x + 3y ¡ z = 1
x + 2y ¡ z = 4
¡2x ¡ y + z = ¡3
Let ¢ =
¯¯¯¯¯¯
¡2 3 ¡1
1 2 ¡1
¡2 ¡1 1
¯¯¯¯¯¯
=
¯¯¯¯¯¯
0 7 ¡3
1 2 ¡1
0 3 ¡1
¯¯¯¯¯¯
= ¡
¯¯¯¯
7 ¡3
3 ¡1
¯¯¯¯
= ¡2 6= 0.
Hence the system has a unique solution which can be calculated using
Cramer's rule:
x =
¢1
¢
; y =
¢2
¢
; z =
¢3
¢
;
where
¢1 =
¯¯¯¯¯¯
1 3 ¡1
4 2 ¡1
¡3 ¡1 1
¯¯¯¯¯¯
= ¡4;
¢2 =
¯¯¯¯¯¯
¡2 1 ¡1
1 4 ¡1
¡2 ¡3 1
¯¯¯¯¯¯
= ¡6;
¢3 =
¯¯¯¯¯¯
¡2 3 1
1 2 4
¡2 ¡1 ¡3
¯¯¯¯¯¯
= ¡8:
Hence x = ¡4
¡2 = 2; y = ¡6
¡2 = 3; z = ¡8
¡2 = 4.
15. In Remark 4.0.4, take A = In. Then we deduce
(a) detEij = ¡1;
(b) detEi(t) = t;
50
(c) detEij(t) = 1.
Now suppose that B is a non{singular n£n matrix. Then we know that B
is a product of elementary row matrices:
B = E1 ¢ ¢ ¢Em:
Consequently we have to prove that
detE1 ¢ ¢ ¢EmA = detE1 ¢ ¢ ¢Em det A:
We prove this by induction on m.
First the case m = 1. We have to prove detE1A = detE1 detA if E1 is
an elementary row matrix. This follows form Remark 4.0.4:
(a) detEijA = ¡detA = detEij detA;
(b) detEi(t)A = t detA = detEi(t) detA;
(c) detEij(t)A = detA = detEij(t) detA.
Let m ¸ 1 and assume the proposition holds for products of m elementary
row matrices. Then
detE1 ¢ ¢ ¢EmEm+1A = det (E1 ¢ ¢ ¢Em)(Em+1A)
= det (E1 ¢ ¢ ¢Em) det (Em+1A)
= det (E1 ¢ ¢ ¢Em) detEm+1 detA
= det ((E1 ¢ ¢ ¢Em)Em+1) detA
and the induction goes through.
Hence detBA = detB detA if B is non{singular.
If B is singular, problem 26, Chapter 2.7 tells us that BA is also singlular.
However singular matrices have zero determinant, so
detB = 0 detBA = 0;
so the equation detBA = detB detA holds trivially in this case.
16. ¯¯¯¯¯¯¯¯
a + b + c a + b a a
a + b a + b + c a a
a a a + b + c a + b
a a a + b a + b + c
¯¯¯¯¯¯¯¯
51
R1 ! R1 ¡ R2
R2 ! R2 ¡ R3
R3 ! R3 ¡ R4
=
¯¯¯¯¯¯¯¯
c ¡c 0 0
b b + c ¡b ¡ c ¡b
0 0 c ¡c
a a a + b a + b + c
¯¯¯¯¯¯¯¯
C2 ! C2 + C1
=
¯¯¯¯¯¯¯¯
c 0 0 0
b 2b + c ¡b ¡ c ¡b
0 0 c ¡c
a 2a a + b a + b + c
¯¯¯¯¯¯¯¯
= c
¯¯¯¯¯¯
2b + c ¡b ¡ c ¡b
0 c ¡c
2a a + b a + b + c
¯¯¯¯¯¯
C3 ! C3 + C2
=
c
¯¯¯¯¯¯
2b + c ¡b ¡ c ¡2b ¡ c
0 c 0
2a a + b 2a + 2b + c
¯¯¯¯¯¯
= c2
¯¯¯¯
2b + c ¡2b ¡ c
2a 2a + 2b + c
¯¯¯¯
= c2(2b + c)
¯¯¯¯
1 ¡1
2a 2a + 2b + c
¯¯¯¯
= c2(2b + c)(4a + 2b + c):
17. Let ¢ =
¯¯¯¯¯¯¯¯
1 + u1 u1 u1 u1
u2 1 + u2 u2 u2
u3 u3 1 + u3 u3
u4 u4 u4 1 + u4
¯¯¯¯¯¯¯¯
: Then using the operation
R1 ! R1 + R2 + R3 + R4
we have
¢ =
¯¯¯¯¯¯¯¯
t t t t
u2 1 + u2 u2 u2
u3 u3 1 + u3 u3
u4 u4 u4 1 + u4
¯¯¯¯¯¯¯¯
(where t = 1 + u1 + u2 + u3 + u4)
= (1 + u1 + u2 + u3 + u4)
¯¯¯¯¯¯¯¯
1 1 1 1
u2 1 + u2 u2 u2
u3 u3 1 + u3 u3
u4 u4 u4 1 + u4
¯¯¯¯¯¯¯¯
The last determinant equals
C2 ! C2 ¡ C1
C3 ! C3 ¡ C1
C4 ! C4 ¡ C1
¯¯¯¯¯¯¯¯
1 0 0 0
u2 1 0 0
u3 0 1 0
u4 0 0 1
¯¯¯¯¯¯¯¯
= 1:
52
18. Suppose that At = ¡A, that A 2 Mn£n(F), where n is odd. Then
detAt = det(¡A)
detA = (¡1)n detA = ¡det A:
Hence (1 + 1) detA = 0 and consequently detA = 0 if 1 + 1 6= 0 in F.
19.
¯¯¯¯¯¯¯¯
1 1 1 1
r 1 1 1
r r 1 1
r r r 1
¯¯¯¯¯¯¯¯
=
C4 ! C4 ¡ C3
C3 ! C3 ¡ C2
C2 ! C2 ¡ C1
=
¯¯¯¯¯¯¯¯
1 0 0 0
r 1 ¡ r 0 0
r 0 1 ¡ r 0
r 0 0 1 ¡ r
¯¯¯¯¯¯¯¯
= (1 ¡ r)3:
20.
¯¯¯¯¯¯
1 a2 ¡ bc a4
1 b2 ¡ ca b4
1 c2 ¡ ab c4
¯¯¯¯¯¯
R2 ! R2 ¡ R1
R3 ! R3 ¡ R1
=
¯¯¯¯¯¯
1 a2 ¡ bc a4
0 b2 ¡ ca ¡ a2 + bc b4 ¡ a4
0 c2 ¡ ab ¡ a2 + bc c4 ¡ a4
¯¯¯¯¯¯
=
¯¯¯¯
b2 ¡ ca ¡ a2 + bc b4 ¡ a4
c2 ¡ ab ¡ a2 + bc c4 ¡ a4
¯¯¯¯
=
¯¯¯¯
(b ¡ a)(b + a) + c(b ¡ a) (b ¡ a)(b + a)(b2 + a2)
(c ¡ a)(c + a) + b(c ¡ a) (c ¡ a)(c + a)(c2 + a2)
¯¯¯¯
=
¯¯¯¯
(b ¡ a)(b + a + c) (b ¡ a)(b + a)(b2 + a2)
(c ¡ a)(c + a + b) (c ¡ a)(c + a)(c2 + a2)
¯¯¯¯
= (b ¡ a)(c ¡ a)
¯¯¯¯
b + a + c (b + a)(b2 + a2)
c + a + b (c + a)(c2 + a2)
¯¯¯¯
= (b ¡ a)(c ¡ a)(a + b + c)
¯¯¯¯
1 (b + a)(b2 + a2)
1 (c + a)(c2 + a2)
¯¯¯¯
:
Finally
¯¯¯¯
1 (b + a)(b2 + a2)
1 (c + a)(c2 + a2)
¯¯¯¯
= (c3 + ac2 + ca2 + a3) ¡ (b3 + ab2 + ba2 + a3)
= (c3 ¡ b3) + a(c2 ¡ b2) + a2(c ¡ b)
= (c ¡ b)(c2 + cb + b2 + a(c + b) + a2)
= (c ¡ b)(c2 + cb + b2 + ac + ab + a2):
53
Section 5.8
1.
(i) (¡3 + i)(14 ¡ 2i) = (¡3)(14 ¡ 2i) + i(14 ¡ 2i)
= f(¡3)14 ¡ (¡3)(2i)g + i(14) ¡ i(2i)
= (¡42 + 6i) + (14i + 2) = ¡40 + 20i:
(ii)
2 + 3i
1 ¡ 4i
=
(2 + 3i)(1 + 4i)
(1 ¡ 4i)(1 + 4i)
=
((2 + 3i) + (2 + 3i)(4i)
12 + 42
= ¡10 + 11i
17
= ¡10
17
+
11
17
i:
(iii)
(1 + 2i)2
1 ¡ i
=
1 + 4i + (2i)2
1 ¡ i
=
1 + 4i ¡ 4
1 ¡ i
= ¡3 + 4i
1 ¡ i
=
(¡3 + 4i)(1 + i)
2
= ¡7 + i
2
= ¡
7
2
+
1
2
i:
2. (i)
iz + (2 ¡ 10i)z = 3z + 2i , z(i + 2 ¡ 10i ¡ 3) = 2i
=, z(¡1 ¡ 9i) = 2i , z = ¡2i
1 + 9i
= ¡2i(1 ¡ 9i)
1 + 81
= ¡18 ¡ 2i
82
= ¡9 ¡ i
41
:
(ii) The coe±cient determinant is
¯¯¯¯
1 + i 2 ¡ i
1 + 2i 3 + i
¯¯¯¯
= (1 + i)(3 + i) ¡ (2 ¡ i)(1 + 2i) = ¡2 + i 6= 0:
Hence Cramer's rule applies: there is a unique solution given by
z =
¯¯¯¯
¡3i 2 ¡ i
2 + 2i 3 + i
¯¯¯¯
¡2 + i
= ¡3 ¡ 11i
¡2 + i
= ¡1 + 5i
w =
¯¯¯¯
1 + i ¡3i
1 + 2i 2 + 2i
¯¯¯¯
¡2 + i
= ¡6 + 7i
¡2 + i
=
19 ¡ 8i
5
:
54
3.
1 + (1 + i) + ¢ ¢ ¢ + (1 + i)99 =
(1 + i)100 ¡ 1
(1 + i) ¡ 1
=
(1 + i)100 ¡ 1
i
= ¡i
©
(1 + i)100 ¡ 1
ª
:
Now (1 + i)2 = 2i. Hence
(1 + i)100 = (2i)50 = 250i50 = 250(¡1)25 = ¡250:
Hence ¡i
©
(1 + i)100 ¡ 1
ª
= ¡i(¡250 ¡ 1) = (250 + 1)i.
4. (i) Let z2 = ¡8 ¡ 6i and write z=x+iy, where x and y are real. Then
z2 = x2 ¡ y2 + 2xyi = ¡8 ¡ 6i;
so x2 ¡ y2 = ¡8 and 2xy = ¡6. Hence
y = ¡3=x; x2 ¡
µ
¡3
x
¶2
= ¡8;
so x4 + 8x2 ¡ 9 = 0. This is a quadratic in x2. Hence x2 = 1 or ¡9 and
consequently x2 = 1. Hence x = 1; y = ¡3 or x = ¡1 and y = 3. Hence
z = 1 ¡ 3i or z = ¡1 + 3i.
(ii) z2 ¡ (3 + i)z + 4 + 3i = 0 has the solutions z = (3 + i § d)=2, where d is
any complex number satisfying
d2 = (3 + i)2 ¡ 4(4 + 3i) = ¡8 ¡ 6i:
Hence by part (i) we can take d = 1 ¡ 3i. Consequently
z =
3 + i § (1 ¡ 3i)
2
= 2 ¡ i or 1 + 2i:
(i) The number lies in the ¯rst quadrant of
the complex plane.
j4 + ij =
p
42 + 12 = p17:
Also Arg (4 + i) = ®, where tan ® = 1=4
and 0 < ® < ¼=2. Hence ® = tan ¡1(1=4).
¾ -
6
?
x
y
4 + i
»»»®»:
55
(ii) The number lies in the third quadrant of
the complex plane.
¯¯¯¯
¡3 ¡ i
2
¯¯¯¯
= j¡3 ¡ ij
2
=
1
2
p
(¡3)2 + (¡1)2 =
1
2
p9 + 1 =
p10
2
:
Also Arg (¡3¡i
2 ) = ¡¼ +®, where tan ® =
1
2=3
2 = 1=3 and 0 < ® < ¼=2. Hence ® =
tan ¡1(1=3).
¾ -
6
?
x
y
¡3¡i
2
)³®³³
(iii) The number lies in the second quadrant of
the complex plane.
j ¡ 1 + 2ij =
p
(¡1)2 + 22 = p5:
Also Arg (¡1+2i) = ¼¡®, where tan ® =
2 and 0 < ® < ¼=2. Hence ® = tan ¡12.
¾ -
6
?
x
y ¡1 + 2i
® A
A
A
A
A
AK
(iv) The number lies in the second quadrant of
the complex plane.
¯¯¯¯¯
¡1 + ip3
2
¯¯¯¯¯
= j ¡ 1 + ip3j
2
=
1
2
q
(¡1)2 + (p3)2 =
1
2
p1 + 3 = 1:
Also Arg (¡1
2 +
p3
2 i) = ¼ ¡ ®, where
tan ® =
p3
2 =1
2 = p3 and 0 < ® < ¼=2.
Hence ® = ¼=3.
¾ -
6
?
x
y
¡1
2 +
p3
2 i
® J
J
J
J
J
J]
6. (i) Let z = (1 + i)(1 + p3i)(p3 ¡ i). Then
jzj = j1 + ijj1 + p3ijj
p3 ¡ ij
=
p
12 + 12
q
12 + (p3)2
q
(p3)2 + (¡1)2
= p2p4p4 = 4p2:
Arg z ´ Arg (1 + i) + Arg (1 + p3) + Arg (p3 ¡ i) (mod 2¼)
56
´
¼
4
+
¼
3 ¡
¼
6 ´
5
12
:
Hence Arg z = 5
12 and the polar decomposition of z is
z = 4p2
µ
cos
5¼
12
+ i sin
5¼
12
¶
:
(ii) Let z = (1+i)5(1¡ip3)5
(p3+i)4 . Then
jzj = j(1 + i)j5j(1 ¡ ip3)j5
j(p3 + i)j4
=
¡p2
¢5
25
24 = 27=2:
Arg z ´ Arg (1 + i)5 + Arg (1 ¡
p3i)5 ¡ Arg (p3 + i)4 (mod 2¼)
´ 5Arg (1 + i) + 5Arg (1 ¡
p3i) ¡ 4Arg (p3 + i)
´ 5
¼
4
+ 5
µ
¡¼
3
¶
¡ 4
¼
6 ´ ¡13¼
12 ´
11¼
12
:
Hence Arg z = 11¼
12 and the polar decomposition of z is
z = 27=2
µ
cos
11¼
12
+ i sin
11¼
12
¶
:
7. (i) Let z = 2(cos ¼
4 + i sin ¼
4 ) and w = 3(cos ¼
6 + i sin ¼
6 ). (Both of these
numbers are already in polar form.)
(a) zw = 6(cos ( ¼
4 + ¼
6 ) + i sin ( ¼
4 + ¼
6 ))
= 6(cos 5¼
12 + i sin 5¼
12 ).
(b) z
w = 2
3 (cos ( ¼
4 ¡ ¼
6 ) + i sin ( ¼
4 ¡ ¼
6 ))
= 2
3 (cos ¼
12 + i sin ¼
12 ).
(c) w
z = 3
2 (cos ( ¼
6 ¡ ¼
4 ) + i sin ( ¼
6 ¡ ¼
4 ))
= 3
2 (cos (¡¼
12 ) + i sin (¡¼
12 )).
(d) z5
w2 = 25
32 (cos ( 5¼
4 ¡ 2¼
6 ) + i sin ( 5¼
4 ¡ 2¼
6 ))
= 32
9 (cos 11¼
12 + i sin 11¼
12 ).
(a) (1 + i)2 = 2i, so
(1 + i)12 = (2i)6 = 26i6 = 64(i2)3 = 64(¡1)3 = ¡64:
57
(b) ( 1¡i p2
)2 = ¡i, so
µ
1 ¡ i
p2
¶
¡6
=
õ
1 ¡ i
p2
¶2
!
¡3
= (¡i)¡3 = ¡1
i3 = ¡1
¡i
=
1
i
= ¡i:
8. (i) To solve the equation z2 = 1 + p3i, we write 1 + p3i in modulus{
argument form:
1 + p3i = 2(cos
¼
3
+ i sin
¼
3
):
Then the solutions are
zk = p2
µ
cos
µ ¼
3 + 2k¼
2
¶
+ i sin
µ ¼
3 + 2k¼
2
¶¶
; k = 0; 1:
Now k = 0 gives the solution
z0 = p2(cos
¼
6
+ i sin
¼
6
) = p2
Ãp3
2
+
i
2
!
=
p3
p2
+
i
p2
:
Clearly z1 = ¡z0.
(ii) To solve the equation z4 = i, we write i in modulus{argument form:
i = cos
¼
2
+ i sin
¼
2
:
Then the solutions are
zk = cos
µ ¼
2 + 2k¼
4
¶
+ i sin
µ ¼
2 + 2k¼
4
¶
; k = 0; 1; 2; 3:
Now cos
³ ¼
2 +2k¼
4
´
= cos
¡¼
8 + k¼
2
¢
, so
zk = cos
µ
¼
8
+
k¼
2
¶
+ sin
µ
¼
8
+
k¼
2
¶
=
³
cos
¼
2
+ i sin
¼
2
´k
(cos
¼
8
+ i sin
¼
8
)
= ik(cos
¼
8
+ i sin
¼
8
):
58
Geometrically, the solutions lie equi{spaced on the unit circle at arguments
¼
8
;
¼
8
+
¼
2
=
5¼
8
;
¼
8
+ ¼ =
9¼
8
;
¼
8
+ 3
¼
2
=
13¼
8
:
Also z2 = ¡z0 and z3 = ¡z1.
(iii) To solve the equation z3 = ¡8i, we rewrite the equation as
µ
z
¡2i
¶3
= 1:
Then µ
z
¡2i
¶
= 1; ¡1 + p3i
2
; or ¡1 ¡ p3i
2
:
Hence z = ¡2i; p3 + i or ¡p3 + i.
Geometrically, the solutions lie equi{spaced on the circle jzj = 2, at
arguments
¼
6
;
¼
6
+
2¼
3
=
5¼
6
;
¼
6
+ 2
2¼
3
=
3¼
2
:
(iv) To solve z4 = 2 ¡ 2i, we write 2 ¡ 2i in modulus{argument form:
2 ¡ 2i = 23=2
µ
cos ¡¼
4
+ i sin ¡¼
4
¶
:
Hence the solutions are
zk = 23=8 cos
µ
¡¼
4 + 2k¼
4
¶
+ i sin
µ
¡¼
4 + 2k¼
4
¶
; k = 0; 1; 2; 3:
We see the solutions can also be written as
zk = 23=8ik
µ
cos ¡¼
16
+ i sin ¡¼
16
¶
= 23=8ik
³
cos
¼
16 ¡ i sin
¼
16
´
:
Geometrically, the solutions lie equi{spaced on the circle jzj = 23=8, at ar-
guments
¡¼
16
; ¡¼
16
+
¼
2
=
7¼
16
; ¡¼
16
+ 2
¼
2
=
15¼
16
; ¡¼
16
+ 3
¼
2
=
23¼
16
:
Also z2 = ¡z0 and z3 = ¡z1.
59
9.
2
4
2 + i ¡1 + 2i 2
1 + i ¡1 + i 1
1 + 2i ¡2 + i 1 + i
3
5 R1 ! R1 ¡ R2
R3 ! R3 ¡ R2
2
4
1 i 1
1 + i ¡1 + i 1
i ¡1 i
3
5
R2 ! R2 ¡ (1 + i)R1
R3 ! R3 ¡ iR1
2
4
1 i 1
0 0 ¡i
0 0 0
3
5 R2 ! iR2
2
4
1 i 1
0 0 1
0 0 0
3
5
R1 ! R1 ¡ R2
2
4
1 i 0
0 0 1
0 0 0
3
5 :
The last matrix is in reduced row{echelon form.
10. (i) Let p = l + im and z = x + iy. Then
pz + pz = (l ¡ im)(x + iy) + (l + im)(x ¡ iy)
= (lx + liy ¡ imx + my) + (lx ¡ liy + imx + my)
= 2(lx + my):
Hence pz + pz = 2n , lx + my = n:
(ii) Let w be the complex number which results from re°ecting the com-
plex number z in the line lx + my = n. Then because p is perpendicular to
the given line, we have
w ¡ z = tp; t 2 R: (a)
Also the midpoint w+z
2 of the segment joining w and z lies on the given line,
so
p
µ
w + z
2
¶
+ p
µ
w + z
2
¶
= n;
p
µ
w + z
2
¶
+ p
µ
w + z
2
¶
= n: (b)
Taking conjugates of equation (a) gives
w ¡ z = tp: (c)
Then substituting in (b), using (a) and (c), gives
p
µ
2w ¡ tp
2
¶
+ p
µ
2z + tp
2
¶
= n
60
and hence
pw + pz = n:
(iii) Let p = b ¡ a and n = jbj2 ¡ jaj2. Then
jz ¡ aj = jz ¡ bj , jz ¡ aj2 = jz ¡ bj2
, (z ¡ a)(z ¡ a) = (z ¡ b)(z ¡ b)
, (z ¡ a)(z ¡ a) = (z ¡ b)(z ¡ b)
, zz ¡ az ¡ za + aa = zz ¡ bz ¡ zb + bb
, (b ¡ a)z + (b ¡ a)z = jbj2 ¡ jaj2
, pz + pz = n:
Suppose z lies on the circle
¯¯¯
z¡a
z¡b
¯¯¯
and let w be the re°ection of z in the
line pz + pz = n. Then by part (ii)
pw + pz = n:
Taking conjugates gives pw + pz = n and hence
z =
n ¡ pw
p
(a)
Substituting for z in the circle equation, using (a) gives
¸ =
¯¯¯¯¯
n¡pw
p ¡ a
n¡pw
p ¡ b
¯¯¯¯¯
=
¯¯¯¯
n ¡ pw ¡ pa
n ¡ pw ¡ pb
¯¯¯¯
: (b)
However
n ¡ pa = jbj2 ¡ jaj2 ¡ (b ¡ a)a
= bb ¡ aa ¡ ba + aa
= b(b ¡ a) = bp:
Similarly n ¡ pb = ap. Consequently (b) simpli¯es to
¸ =
¯¯¯¯b
p
¡
p
w
ap ¡ pw
¯¯¯¯
=
¯¯¯¯
b ¡ w
a ¡ w
¯¯¯¯
=
¯¯¯¯
w ¡ b
w ¡ a
¯¯¯¯
;
which gives
¯¯¯
w¡a
w¡b
¯¯¯
= 1
¸.
61
11. Let a and b be distinct complex numbers and 0 < ® < ¼.
(i) When z1 lies on the circular arc shown, it subtends a constant angle
®. This angle is given by Arg (z1 ¡ a) ¡ Arg (z1 ¡ b). However
Arg
µ
z1 ¡ a
z1 ¡ b
¶
= Arg (z1 ¡ a) ¡ Arg (z1 ¡ b) + 2k¼
= ® + 2k¼:
It follows that k = 0, as 0 < ® < ¼ and ¡¼ < Arg µ · ¼. Hence
Arg
µ
z1 ¡ a
z1 ¡ b
¶
= ®:
Similarly if z2 lies on the circular arc shown, then
Arg
µ
z2 ¡ a
z2 ¡ b
¶
= ¡° = ¡(¼ ¡ ®) = ® ¡ ¼:
Replacing ® by ¼ ¡ ®, we deduce that if z4 lies on the circular arc shown,
then
Arg
µ
z4 ¡ a
z4 ¡ b
¶
= ¼ ¡ ®;
while if z3 lies on the circular arc shown, then
Arg
µ
z3 ¡ a
z3 ¡ b
¶
= ¡®:
The straight line through a and b has the equation
z = (1 ¡ t)a + tb;
62
where t is real. Then 0 < t < 1 describes the segment ab. Also
z ¡ a
z ¡ b
=
t
t ¡ 1
:
Hence z¡a
z¡b is real and negative if z is on the segment a, but is real and
positive if z is on the remaining part of the line, with corresponding values
Arg
µ
z ¡ a
z ¡ b
¶
= ¼; 0;
respectively.
(ii) Case (a) Suppose z1; z2 and z3 are not collinear. Then these points
determine a circle. Now z1 and z2 partition this circle into two arcs. If z3
and z4 lie on the same arc, then
Arg
µ
z3 ¡ z1
z3 ¡ z2
¶
= Arg
µ
z4 ¡ z1
z4 ¡ z2
¶
;
whereas if z3 and z4 lie on opposite arcs, then
Arg
µ
z3 ¡ z1
z3 ¡ z2
¶
= ®
and
Arg
µ
z4 ¡ z1
z4 ¡ z2
¶
= ® ¡ ¼:
Hence in both cases
Arg
µ
z3 ¡ z1
z3 ¡ z2
=
z4 ¡ z1
z4 ¡ z2
¶
´ Arg
µ
z3 ¡ z1
z3 ¡ z2
¶
¡ Arg
µ
z4 ¡ z1
z4 ¡ z2
¶
(mod 2¼)
´ 0 or ¼:
In other words, the cross{ratio
z3 ¡ z1
z3 ¡ z2
=
z4 ¡ z1
z4 ¡ z2
is real.
(b) If z1; z2 and z3 are collinear, then again the cross{ratio is real.
The argument is reversible.
(iii) Assume that A; B; C; D are distinct points such that the cross{ratio
r =
z3 ¡ z1
z3 ¡ z2
=
z4 ¡ z1
z4 ¡ z2
is real. Now r cannot be 0 or 1. Then there are three cases:
63
(i) 0 < r < 1;
(ii) r < 0;
(iii) r > 1.
Case (i). Here jrj + j1 ¡ rj = 1. So
¯¯¯¯
z4 ¡ z1
z4 ¡ z2 ¢
z3 ¡ z2
z3 ¡ z1
¯¯¯¯
+
¯¯¯¯
1 ¡
µ
z4 ¡ z1
z4 ¡ z2 ¢
z3 ¡ z2
z3 ¡ z1
¶¯¯¯¯ = 1:
Multiplying both sides by the denominator jz4 ¡ z2jjz3 ¡ z1j gives after
simpli¯cation
jz4 ¡ z1jjz3 ¡ z2j + jz2 ¡ z1jjz4 ¡ z3j = jz4 ¡ z2jjz3 ¡ z1j;
or
(a) AD ¢ BC + AB ¢ CD = BD ¢ AC:
Case (ii). Here 1 + jrj = j1 ¡ rj. This leads to the equation
(b) BD ¢ AC + AD ¢ BC+ = AB ¢ CD:
Case (iii). Here 1 + j1 ¡ rj = jrj. This leads to the equation
(c) BD ¢ AC + AB ¢ CD = AD ¢ BC:
Conversely if (a), (b) or (c) hold, then we can reverse the argument to deduce
that r is a complex number satisfying one of the equations
jrj + j1 ¡ rj = 1; 1 + jrj = j1 ¡ rj; 1 + j1 ¡ rj = jrj;
from which we deduce that r is real.
64
Section 6.3
1. Let A =
·
4 ¡3
1 0
¸
. Then A has characteristic equation ¸2 ¡ 4¸ + 3 = 0
or (¸ ¡ 3)(¸ ¡ 1) = 0. Hence the eigenvalues of A are ¸1 = 3 and ¸2 = 1.
¸1 = 3. The corresponding eigenvectors satisfy (A ¡ ¸1I2)X = 0, or
·
1 ¡3
1 ¡3
¸
=
·
0
0
¸
;
or equivalently x ¡ 3y = 0. Hence
·
x
y
¸
=
·
3y
y
¸
= y
·
3
1
¸
and we take X1 =
·
3
1
¸
.
Similarly for ¸2 = 1 we ¯nd the eigenvector X2 =
·
1
1
¸
.
Hence if P = [X1jX2] =
·
3 1
1 1
¸
, then P is non{singular and
P¡1AP =
·
3 0
0 1
¸
:
Hence
A = P
·
3 0
0 1
¸
P¡1
and consequently
An = P
·
3n 0
0 1n
¸
P¡1
=
·
3 1
1 1
¸ ·
3n 0
0 1n
¸
1
2
·
1 ¡1
¡1 3
¸
=
1
2
·
3n+1 1
3n 1
¸ ·
1 ¡1
¡1 3
¸
=
1
2
·
3n+1 ¡ 1 ¡3n+1 + 3
3n ¡ 1 ¡3n + 3
¸
=
3n ¡ 1
2
A +
3 ¡ 3n
2
I2:
65
2. Let A =
·
3=5 4=5
2=5 1=5
¸
. Then we ¯nd that the eigenvalues are ¸1 = 1 and
¸2 = ¡1=5, with corresponding eigenvectors
X1 =
·
2
1
¸
and X2 =
·
¡1
1
¸
:
Then if P = [X1jX2], P is non{singular and
P¡1AP =
·
1 0
0 ¡1=5
¸
and A = P
·
1 0
0 ¡1=5
¸
P¡1:
Hence
An = P
·
1 0
0 (¡1=5)n
¸
P¡1
! P
·
1 0
0 0
¸
P¡1
=
·
2 ¡1
1 1
¸ ·
1 0
0 0
¸
1
3
·
1 1
¡1 2
¸
=
1
3
·
2 0
1 0
¸ ·
1 1
¡1 2
¸
=
1
3
·
2 2
1 1
¸
=
·
2=3 2=3
1=3 1=3
¸
:
3. The given system of di®erential equations is equivalent to _X = AX,
where
A =
·
3 ¡2
5 ¡4
¸
and X =
·
x
y
¸
:
The matrix P =
·
2 1
5 1
¸
is a non-singular matrix of eigenvectors corre-
sponding to eigenvalues ¸1 = ¡2 and ¸2 = 1. Then
P¡1AP =
·
¡2 0
0 1
¸
:
The substitution X = PY , where Y = [x1; y1]t, gives
_Y
=
·
¡2 0
0 1
¸
Y;
66
or equivalently x_1 = ¡2x1 and y_1 = y1.
Hence x1 = x1(0)e¡2t and y1 = y1(0)et. To determine x1(0) and y1(0),
we note that
·
x1(0)
y1(0)
¸
= P¡1
·
x(0)
y(0)
¸
= ¡
1
3
·
1 ¡1
¡5 2
¸ ·
13
22
¸
=
·
3
7
¸
:
Hence x1 = 3e¡2t and y1 = 7et. Consequently
x = 2x1 + y1 = 6e¡2t + 7et and y = 5x1 + y1 = 15e¡2t + 7et:
4. Introducing the vector Xn =
·
xn
yn
¸
, the system of recurrence relations
xn+1 = 3xn ¡ yn
yn+1 = ¡xn + 3yn;
becomes Xn+1 = AXn, where A =
·
3 ¡1
¡1 3
¸
. Hence Xn = AnX0, where
X0 =
·
1
2
¸
.
To ¯nd An we can use the eigenvalue method. We get
An =
1
2
·
2n + 4n 2n ¡ 4n
2n ¡ 4n 2n + 4n
¸
:
Hence
Xn =
1
2
·
2n + 4n 2n ¡ 4n
2n ¡ 4n 2n + 4n
¸ ·
1
2
¸
=
1
2
·
2n + 4n + 2(2n ¡ 4n)
2n ¡ 4n + 2(2n + 4n)
¸
=
1
2
·
3 £ 2n ¡ 4n
3 £ 2n + 4n
¸
=
·
(3 £ 2n ¡ 4n)=2
(3 £ 2n + 4n)=2
¸
:
Hence xn = 1
2 (3 £ 2n ¡ 4n) and yn = 1
2(3 £ 2n + 4n).
5. Let A =
·
a b
c d
¸
be a real or complex matrix with distinct eigenvalues
¸1; ¸2 and corresponding eigenvectors X1; X2. Also let P = [X1jX2].
(a) The system of recurrence relations
xn+1 = axn + byn
yn+1 = cxn + dyn
67
has the solution
·
xn
yn
¸
= An
·
x0
y0
¸
=
µ
P
·
¸1 0
0 ¸2
¸
P¡1
¶n ·
x0
y0
¸
= P
·
¸n
1 0
0 ¸n
2
¸
P¡1
·
x0
y0
¸
= [X1jX2]
·
¸n
1 0
0 ¸n
2
¸ ·
®
¯
¸
= [X1jX2]
·
¸n
1 ®
¸n
2 ¯
¸
= ¸n
1 ®X1 + ¸n
2 ¯X2;
where ·
®
¯
¸
= P¡1
·
x0
y0
¸
:
(b) In matrix form, the system is _X = AX, where X =
·
x
y
¸
. We substitute
X = PY , where Y = [x1; y1]t. Then
_X
= P _Y = AX = A(PY );
so
_Y
= (P¡1AP)Y =
·
¸1 0
0 ¸2
¸ ·
x1
y1
¸
:
Hence x_1 = ¸1x1 and y_1 = ¸2y1. Then
x1 = x1(0)e¸1t and y1 = y1(0)e¸2t:
But ·
x(0)
y(0)
¸
= P
·
x1(0)
y1(0)
¸
;
so ·
x1(0)
y1(0)
¸
= P¡1
·
x(0)
y(0)
¸
=
·
®
¯
¸
:
Consequently x1(0) = ® and y1(0) = ¯ and
·
x
y
¸
= P
·
x1
y1
¸
= [X1jX2]
·
®e¸1t
¯e¸2t
¸
= ®e¸1tX1 + ¯e¸2tX2:
68
6. Let A =
·
a b
c d
¸
be a real matrix with non{real eigenvalues ¸ = a + ib
and ¸ = a¡ib, with corresponding eigenvectors X = U+iV and X = U¡iV ,
where U and V are real vectors. Also let P be the real matrix de¯ned by
P = [UjV ]. Finally let a + ib = reiµ, where r > 0 and µ is real.
(a) As X is an eigenvector corresponding to the eigenvalue ¸, we have AX =
¸X and hence
A(U + iV ) = (a + ib)(U + iV )
AU + iAV = aU ¡ bV + i(bU + aV ):
Equating real and imaginary parts then gives
AU = aU ¡ bV
AV = bU + aV:
(b)
AP = A[UjV ] = [AUjAV ] = [aU¡bV jbU+aV ] = [UjV ]
·
a b
¡b a
¸
= P
·
a b
¡b a
¸
:
Hence, as P can be shown to be non{singular,
P¡1AP =
·
a b
¡b a
¸
:
(The fact that P is non{singular is easily proved by showing the columns of
P are linearly independent: Assume xU + yV = 0, where x and y are real.
Then we ¯nd
(x + iy)(U ¡ iV ) + (x ¡ iy)(U + iV ) = 0:
Consequently x+iy = 0 as U¡iV and U+iV are eigenvectors corresponding
to distinct eigenvalues a¡ib and a+ib and are hence linearly independent.
Hence x = 0 and y = 0.)
(c) The system of recurrence relations
xn+1 = axn + byn
yn+1 = cxn + dyn
69
has solution
·
xn
yn
¸
= An
·
x0
y0
¸
= P
·
a b
¡b a
¸n
P¡1
·
x0
y0
¸
= P
·
r cos µ r sin µ
¡r sin µ r cos µ
¸n ·
®
¯
¸
= Prn
·
cos µ sin µ
¡sin µ cos µ
¸n ·
®
¯
¸
= rn[UjV ]
·
cos nµ sin nµ
¡sin nµ cos nµ
¸ ·
®
¯
¸
= rn[UjV ]
·
® cos nµ + ¯ sin nµ
¡® sin nµ + ¯ cos nµ
¸
= rn f(® cos nµ + ¯ sin nµ)U + (¡® sin nµ + ¯ cos nµ)V g
= rn f(cos nµ)(®U + ¯V ) + (sin nµ)(¯U ¡ ®V )g :
(d) The system of di®erential equations
dx
dt
= ax + by
dy
dt
= cx + dy
is attacked using the substitution X = PY , where Y = [x1; y1]t. Then
_Y
= (P¡1AP)Y;
so ·
x_1
y_1
¸
=
·
a b
¡b a
¸ ·
x1
y1
¸
:
Equating components gives
x_1 = ax1 + by1
y_1 = ¡bx1 + ay1:
Now let z = x1 + iy1. Then
z_ = x_1 + iy_1 = (ax1 + by1) + i(¡bx1 + ay1)
= (a ¡ ib)(x1 + iy1) = (a ¡ ib)z:
70
Hence
z = z(0)e(a¡ib)t
x1 + iy1 = (x1(0) + iy1(0))eat(cos bt ¡ i sin bt):
Equating real and imaginary parts gives
x1 = eat fx1(0) cos bt + y1(0) sin btg
y1 = eat fy1(0) cos bt ¡ x1(0) sin btg :
Now if we de¯ne ® and ¯ by
·
®
¯
¸
= P¡1
·
x(0)
y(0)
¸
;
we see that ® = x1(0) and ¯ = y1(0). Then
·
x
y
¸
= P
·
x1
y1
¸
= [UjV ]
·
eat(® cos bt + ¯ sin bt)
eat(¯ cos bt ¡ ® sin bt)
¸
= eatf(® cos bt + ¯ sin bt)U + (¯ cos bt ¡ ® sin bt)V g
= eatfcos bt(®U + ¯V ) + sin bt(¯U ¡ ®V )g:
7. (The case of repeated eigenvalues.) Let A =
·
a b
c d
¸
and suppose that
the characteristic polynomial of A, ¸2 ¡(a+d)¸+(ad¡bc), has a repeated
root ®. Also assume that A 6= ®I2.
(i)
¸2 ¡ (a + d)¸ + (ad ¡ bc) = (¸ ¡ ®)2
= ¸2 ¡ 2®¸ + ®2:
Hence a + d = 2® and ad ¡ bc = ®2 and
(a + d)2 = 4(ad ¡ bc);
a2 + 2ad + d2 = 4ad ¡ 4bc;
a2 ¡ 2ad + d2 + 4bc = 0;
(a ¡ d)2 + 4bc = 0:
71
(ii) Let B ¡ A ¡ ®I2. Then
B2 = (A ¡ ®I2)2 = A2 ¡ 2®A + ®2I2
= A2 ¡ (a + d)A + (ad ¡ bc)I2;
But by problem 3, chapter 2.4, A2 ¡ (a + d)A + (ad ¡ bc)I2 = 0, so
B2 = 0.
(iii) Now suppose that B 6= 0. Then BE1 6= 0 or BE2 6= 0, as BEi is the
i{th column of B. Hence BX2 6= 0, where X2 = E1 or X2 = E2.
(iv) Let X1 = BX2 and P = [X1jX2]. We prove P is non{singular by
demonstrating that X1 and X2 are linearly independent.
Assume xX1 + yX2 = 0. Then
xBX2 + yX2 = 0
B(xBX2 + yX2) = B0 = 0
xB2X2 + yBX2 = 0
x0X2 + yBX2 = 0
yBX2 = 0:
Hence y = 0 as BX2 6= 0. Hence xBX2 = 0 and so x = 0.
Finally, BX1 = B(BX2) = B2X2 = 0, so (A ¡ ®I2)X1 = 0 and
AX1 = ®X1: (2)
Also
X1 = BX2 = (A ¡ ®I2)X2 = AX2 ¡ ®X2:
Hence
AX2 = X1 + ®X2: (3)
Then, using (2) and (3), we have
AP = A[X1jX2] = [AX1jAX2]
= [®X1jX1 + ®X2]
= [X1jX2]
·
® 1
0 ®
¸
:
Hence
AP = P
·
® 1
0 ®
¸
and hence
72
P¡1AP =
·
® 1
0 ®
¸
:
8. The system of di®erential equations is equivalent to the single matrix
equation _X = AX, where A =
·
4 ¡1
4 8
¸
.
The characteristic polynomial of A is ¸2 ¡ 12¸ + 36 = (¸ ¡ 6)2, so we
can use the previous question with ® = 6. Let
B = A ¡ 6I2 =
·
¡2 ¡1
4 2
¸
:
Then BX2 =
·
¡2
4
¸
6=
·
0
0
¸
, if X2 =
·
1
0
¸
. Also let X1 = BX2. Then if
P = [X1jX2], we have
P¡1AP =
·
6 1
0 6
¸
:
Now make the change of variables X = PY , where Y =
·
x1
y1
¸
. Then
_Y
= (P¡1AP)Y =
·
6 1
0 6
¸
Y;
or equivalently x_1 = 6x1 + y1 and y_1 = 6y1.
Solving for y1 gives y1 = y1(0)e6t. Consequently
x_1 = 6x1 + y1(0)e6t:
Multiplying both side of this equation by e¡6t gives
d
dt
(e¡6tx1) = e¡6tx_1 ¡ 6e¡6tx1 = y1(0)
e¡6tx1 = y1(0)t + c;
where c is a constant. Substituting t = 0 gives c = x1(0). Hence
e¡6tx1 = y1(0)t + x1(0)
and hence
x1 = e6t(y1(0)t + x1(0)):
73
However, since we are assuming x(0) = 1 = y(0), we have
·
x1(0)
y1(0)
¸
= P¡1
·
x(0)
y(0)
¸
=
1
¡4
·
0 ¡1
¡4 ¡2
¸ ·
1
1
¸
=
1
¡4
·
¡1
¡6
¸
=
·
1=4
3=2
¸
:
Hence x1 = e6t( 3
2 t + 1
4 ) and y1 = 3
2e6t.
Finally, solving for x and y,
·
x
y
¸
=
·
¡2 1
4 0
¸ ·
x1
y1
¸
=
·
¡2 1
4 0
¸ 2
4
e6t( 3
2 t + 1
4 )
3
2e6t
3
5
=
2
4
(¡2)e6t( 3
2 t + 1
4 ) + 3
2e6t
4e6t( 3
2 t + 1
4 )
3
5
=
·
e6t(1 ¡ 3t)
e6t(6t + 1)
¸
:
Hence x = e6t(1 ¡ 3t) and y = e6t(6t + 1).
9. Let
A =
2
4
1=2 1=2 0
1=4 1=4 1=2
1=4 1=4 1=2
3
5 :
(a) We ¯rst determine the characteristic polynomial chA(¸).
chA(¸) = det (¸I3 ¡ A) =
¯¯¯¯¯¯
¸ ¡ 1=2 ¡1=2 0
1=4 ¸ ¡ 1=4 ¡1=2
¡1=4 ¡1=4 ¸ ¡ 1=2
¯¯¯¯¯¯
=
µ
¸ ¡
1
2
¶¯¯¯¯
¸ ¡ 1=4 ¡1=2
¡1=4 ¸ ¡ 1=2
¯¯¯¯
+
1
2
¯¯¯¯
1=4 ¡1=2
¡1=4 ¸ ¡ 1=2
¯¯¯¯
=
µ
¸ ¡
1
2
¶½µ
¸ ¡
1
4
¶µ
¸ ¡
1
2
¶
¡
1
8
¾
+
1
2
½
¡1
4
µ
¸ ¡
1
2
¶
¡
1
8
¾
=
µ
¸ ¡
1
2
¶µ
¸2 ¡
3¸
4
¶
¡
¸
8
= ¸
½µ
¸ ¡
1
2
¶µ
¸ ¡
3
4
¶
¡
1
8
¾
74
= ¸
µ
¸2 ¡
5¸
4
+
1
4
¶
= ¸(¸ ¡ 1)
µ
¸ ¡
1
4
¶
:
(b) Hence the characteristic polynomial has no repeated roots and we can
use Theorem 6.2.2 to ¯nd a non{singular matrix P such that
P¡1AP = diag(1; 0;
1
4
):
We take P = [X1jX2jX3], where X1; X2; X3 are eigenvectors corresponding
to the respective eigenvalues 1; 0; 1
4 .
Finding X1: We have to solve (A ¡ I3)X = 0. we have
A ¡ I3 =
2
4 ¡1=2 1=2 0
1=4 ¡3=4 1=2
1=4 1=4 ¡1=2
3
5 !
2
4
1 0 ¡1
0 1 ¡1
0 0 0
3
5 :
Hence the eigenspace consists of vectors X = [x; y; z]t satisfying x = z and
y = z, with z arbitrary. Hence
X =
2
4
z
z
z
3
5 = z
2
4
1
1
1
3
5
and we can take X1 = [1; 1; 1]t.
Finding X2: We solve AX = 0. We have
A =
2
4
1=2 1=2 0
1=4 1=4 1=2
1=4 1=4 1=2
3
5 !
2
4
1 1 0
0 0 1
0 0 0
3
5 :
Hence the eigenspace consists of vectors X = [x; y; z]t satisfying x = ¡y
and z = 0, with y arbitrary. Hence
X =
2
4 ¡y
y
0
3
5 = y
2
4 ¡1
1
0
3
5
and we can take X2 = [¡1; 1; 0]t.
Finding X3: We solve (A ¡ 1
4 I3)X = 0. We have
A ¡
1
4
I3 =
2
4
1=4 1=2 0
1=4 0 1=2
1=4 1=4 1=4
3
5 !
2
4
1 0 2
0 1 ¡1
0 0 0
3
5 :
75
Hence the eigenspace consists of vectors X = [x; y; z]t satisfying x = ¡2z
and y = z, with z arbitrary. Hence
X =
2
4 ¡2z
z
0
3
5 = z
2
4 ¡2
1
0
3
5
and we can take X3 = [¡2; 1; 1]t.
Hence we can take P =
2
4
1 ¡1 ¡2
1 1 1
1 0 1
3
5.
(c) A = Pdiag(1; 0; 1
4 )P¡1 so An = Pdiag(1; 0; 1
4n )P¡1.
Hence
An =
2
4
1 ¡1 ¡2
1 1 1
1 0 1
3
5
2
4
1 0 0
0 0 0
0 0 1
4n
3
5 1
3
2
4
1 1 1
0 3 ¡3
¡1 ¡1 2
3
5
=
1
3
2
4
1 0 ¡ 2
4n
1 0 1
4n
1 0 1
4n
3
5
2
4
1 1 1
0 3 ¡3
¡1 ¡1 2
3
5
=
1
3
2
4
1 + 2
4n 1 + 2
4n 1 ¡ 4
4n
1 ¡ 1
4n 1 ¡ 1
4n 1 + 2
4n
1 ¡ 1
4n 1 ¡ 1
4n 1 + 2
4n
3
5
=
1
3
2
4
1 1 1
1 1 1
1 1 1
3
5 +
1
3 ¢ 4n
2
4
2 2 ¡4
¡1 ¡1 2
¡1 ¡1 2
3
5 :
10. Let
A =
2
4
5 2 ¡2
2 5 ¡2
¡2 ¡2 5
3
5 :
(a) We ¯rst determine the characteristic polynomial chA(¸).
chA(¸) =
¯¯¯¯¯¯
¸ ¡ 5 ¡2 2
¡2 ¸ ¡ 5 2
2 2 ¸ ¡ 5
3
5 R3 ! R3 + R2
=
¯¯¯¯¯¯
¸ ¡ 5 ¡2 2
¡2 ¸ ¡ 5 2
0 ¸ ¡ 3 ¸ ¡ 3
¯¯¯¯¯¯
= (¸ ¡ 3)
¯¯¯¯¯¯
¸ ¡ 5 ¡2 2
¡2 ¸ ¡ 5 2
0 1 1
¯¯¯¯¯¯
76
C3 ! C3 ¡ C2 = (¸ ¡ 3)
¯¯¯¯¯¯
¸ ¡ 5 ¡2 4
¡2 ¸ ¡ 5 ¡¸ + 7
0 1 0
¯¯¯¯¯¯
= ¡(¸ ¡ 3)
¯¯¯¯
¸ ¡ 5 4
¡2 ¡¸ + 7
¯¯¯¯
= ¡(¸ ¡ 3) f(¸ ¡ 5)(¡¸ + 7) + 8g
= ¡(¸ ¡ 3)(¡¸2 + 5¸ + 7¸ ¡ 35 + 8)
= ¡(¸ ¡ 3)(¡¸2 + 12¸ ¡ 27)
= ¡(¸ ¡ 3)(¡1)(¸ ¡ 3)(¸ ¡ 9)
= (¸ ¡ 3)2(¸ ¡ 9):
We have to ¯nd bases for each of the eigenspaces N(A¡9I3) and N(A¡3I3).
First we solve (A ¡ 3I3)X = 0. We have
A ¡ 3I3 =
2
4
2 2 ¡2
2 2 ¡2
¡2 ¡2 2
3
5 !
2
4
1 1 ¡1
0 0 0
0 0 0
3
5 :
Hence the eigenspace consists of vectors X = [x; y; z]t satisfying x = ¡y+z,
with y and z arbitrary. Hence
X =
2
4 ¡y + z
y
z
3
5 = y
2
4 ¡1
1
0
3
5 + z
2
4
1
0
1
3
5 ;
so X1 = [¡1; 1; 0]t and X2 = [1; 0; 1]t form a basis for the eigenspace
corresponding to the eigenvalue 3.
Next we solve (A ¡ 9I3)X = 0. We have
A ¡ 9I3 =
2
4 ¡4 2 ¡2
2 ¡4 ¡2
¡2 ¡2 ¡4
3
5 !
2
4
1 0 1
0 1 1
0 0 0
3
5 :
Hence the eigenspace consists of vectors X = [x; y; z]t satisfying x = ¡z
and y = ¡z, with z arbitrary. Hence
X =
2
4 ¡z
¡z
z
3
5 = z
2
4 ¡1
¡1
1
3
5
and we can take X3 = [¡1; ¡1; 1]t as a basis for the eigenspace correspond-
ing to the eigenvalue 9.
77
Then Theorem 6.2.3 assures us that P = [X1jX2jX3] is non{singular and
P¡1AP =
2
4
3 0 0
0 3 0
0 0 9
3
5 :
78
x
1
y
1
-9 -4.5 4.5 9 13.5
4.5
9
13.5
-4.5
-9
x
y
x
1
y
1
-8 -4 4 8
4
8
-4
-8
x
y
Figure 1: (a): x2 ¡ 8x + 8y + 8 = 0; (b): y2 ¡ 12x + 2y + 25 = 0
Section 7.3
1. (i) x2¡8x+8y+8 = (x¡4)2+8(y¡1). So the equation x2¡8x+8y+8 = 0
becomes
x2
1 + 8y1 = 0 (1)
if we make a translation of axes x ¡ 4 = x1; y ¡ 1 = y1.
However equation (1) can be written as a standard form
y1 = ¡
1
8
x2
1 ;
which represents a parabola with vertex at (4; 1). (See Figure 1(a).)
(ii) y2 ¡12x+2y +25 = (y +1)2 ¡12(x¡2). Hence y2 ¡12x+2y +25 = 0
becomes
y2
1 ¡ 12x1 = 0 (2)
if we make a translation of axes x ¡ 2 = x1; y + 1 = y1.
However equation (2) can be written as a standard form
y2
1 = 12x1;
which represents a parabola with vertex at (2; ¡1). (See Figure 1(b).)
2. 4xy ¡ 3y2 = XtAX, where A =
·
0 2
2 ¡3
¸
and X =
·
x
y
¸
. The
eigenvalues of A are the roots of ¸2 + 3¸ ¡ 4 = 0, namely ¸1 = ¡4 and
¸2 = 1.
79
The eigenvectors corresponding to an eigenvalue ¸ are the non{zero vec-
tors [x; y]t satisfying
·
0 ¡ ¸ 2
2 ¡3 ¡ ¸
¸ ·
x
y
¸
=
·
0
0
¸
:
¸1 = ¡4 gives equations
4x + 2y = 0
2x + y = 0
which has the solution y = ¡2x. Hence
·
x
y
¸
=
·
x
¡2x
¸
= x
·
1
¡2
¸
:
A corresponding unit eigenvector is [1=p5; ¡2=p5]t.
¸2 = 1 gives equations
¡x + 2y = 0
2x ¡ 4y = 0
which has the solution x = 2y. Hence
·
x
y
¸
=
·
2y
y
¸
= y
·
2
1
¸
:
A corresponding unit eigenvector is [2=p5; 1=p5]t.
Hence if
P =
"
1 p5
2 p5
¡2 p5
1 p5
#
;
then P is an orthogonal matrix. Also as det P = 1, P is a proper orthogonal
matrix and the equation
·
x
y
¸
= P
·
x1
y1
¸
represents a rotation to new x1; y1 axes whose positive directions are given
by the respective columns of P. Also
PtAP =
·
¡4 0
0 1
¸
:
80
Then XtAX = ¡4x2
1 +y2
1 and the original equation 4xy ¡3y2 = 8 becomes
¡4x2
1 + y2
1 = 8, or the standard form
¡x2
1
2
+
y2
1
8
= 1;
which represents an hyperbola.
The asymptotes assist in drawing the curve. They are given by the
equations
¡x2
1
2
+
y2
1
8
= 0; or y1 = §2x1:
Now ·
x1
y1
¸
= Pt
·
x
y
¸
=
"
1 p5 ¡2 p5
2 p5
1 p5
#·
x
y
¸
;
so
x1 =
x ¡ 2y
p5
; y1 =
2x + y
p5
:
Hence the asymptotes are
2x + y
p5
= §2
µ
x ¡ 2y
p5
¶
;
which reduces to y = 0 and y = 4x=3. (See Figure 2(a).)
3. 8x2 ¡ 4xy + 5y2 = XtAX, where A =
·
8 ¡2
¡2 5
¸
and X =
·
x
y
¸
. The
eigenvalues of A are the roots of ¸2 ¡ 13¸ + 36 = 0, namely ¸1 = 4 and
¸2 = 9. Corresponding unit eigenvectors turn out to be [1=p5; 2=p5]t and
[¡2=p5; 1=p5]t. Hence if
P =
"
1 p5 ¡2 p5
2 p5
1 p5
#
;
then P is an orthogonal matrix. Also as det P = 1, P is a proper orthogonal
matrix and the equation
·
x
y
¸
= P
·
x1
y1
¸
represents a rotation to new x1; y1 axes whose positive directions are given
by the respective columns of P. Also
PtAP =
·
4 0
0 9
¸
:
81
x
2
y
2
-16 -8 8 16
8
16
-8
-16
x
y x
2
y
2
-2.85 -1.9 -0.95 0.95 1.9 2.85
0.95
1.9
2.85
-0.95
-1.9
-2.85
x
y
Figure 2: (a): 4xy ¡ 3y2 = 8; (b): 8x2 ¡ 4xy + 5y2 = 36
Then XtAX = 4x2
1 + 9y2
1 and the original equation 8x2 ¡ 4xy + 5y2 = 36
becomes 4x2
1 + 9y2
1 = 36, or the standard form
x2
1
9
+
y2
1
4
= 1;
which represents an ellipse as in Figure 2(b).
The axes of symmetry turn out to be y = 2x and x = ¡2y.
4. We give the sketch only for parts (i), (iii) and (iv). We give the working
for (ii) only. See Figures 3(a) and 4(a) and 4(b), respectively.
(ii) We have to investigate the equation
5x2 ¡ 4xy + 8y2 + 4p5x ¡ 16p5y + 4 = 0: (3)
Here 5x2 ¡ 4xy + 8y2 = XtAX, where A =
·
5 ¡2
¡2 8
¸
and X =
·
x
y
¸
.
The eigenvalues of A are the roots of ¸2 ¡13¸+36 = 0, namely ¸1 = 9 and
¸2 = 4. Corresponding unit eigenvectors turn out to be [1=p5; ¡2=p5]t and
[2=p5; 1=p5]t. Hence if
P =
"
1 p5
2 p5
¡2 p5
1 p5
#
;
then P is an orthogonal matrix. Also as det P = 1, P is a proper orthogonal
matrix and the equation
·
x
y
¸
= P
·
x1
y1
¸
82
x
1
y
1
-6 -3 3 6
3
6
-3
-6
x
y x
2
y
2
-4.5 -3 -1.5 1.5 3 4.5
1.5
3
4.5
-1.5
-3
-4.5
x
y
Figure 3: (a): 4x2 ¡ 9y2 ¡ 24x ¡ 36y ¡ 36 = 0; (b): 5x2 ¡ 4xy + 8y2 + p5x ¡ 16p5y + 4 = 0
x
2
y
2
-9 -4.5 4.5 9
4.5
9
-4.5
-9
x
y
x
2
y
2
-9 -4.5 4.5 9
4.5
9
-4.5
-9
x
y
Figure 4: (a): 4x2 +y2 ¡4xy ¡10y ¡19 = 0; (b): 77x2 +78xy ¡27y2 +
70x ¡ 30y + 29 = 0
83
represents a rotation to new x1; y1 axes whose positive directions are given
by the respective columns of P. Also
PtAP =
·
9 0
0 4
¸
:
Moreover
5x2 ¡ 4xy + 8y2 = 9x2
1 + 4y2
1:
To get the coe±cients of x1 and y1 in the transformed form of equation (3),
we have to use the rotation equations
x =
1
p5
(x1 + 2y1); y =
1
p5
(¡2x1 + y1):
Then equation (3) transforms to
9x2
1 + 4y2
1 + 36x1 ¡ 8y1 + 4 = 0;
or, on completing the square,
9(x1 + 2)2 + 4(y1 ¡ 1)2 = 36;
or in standard form
x2
2
4
+
y2
2
9
= 1;
where x2 = x1 + 2 and y2 = y1 ¡ 1. Thus we have an ellipse, centre
(x2; y2) = (0; 0), or (x1; y1) = (¡2; 1), or (x; y) = (0; p5).
The axes of symmetry are given by x2 = 0 and y2 = 0, or x1 + 2 = 0
and y1 ¡ 1 = 0, or
1
p5
(x ¡ 2y) + 2 = 0 and
1
p5
(2x + y) ¡ 1 = 0;
which reduce to x ¡ 2y + 2p5 = 0 and 2x + y ¡ p5 = 0. See Figure 3(b).
5. (i) Consider the equation
2x2 + y2 + 3xy ¡ 5x ¡ 4y + 3 = 0: (4)
¢ =
¯¯¯¯¯¯
2 3=2 ¡5=2
3=2 1 ¡2
¡5=2 ¡2 3
¯¯¯¯¯¯
= 8
¯¯¯¯¯¯
4 3 ¡5
3 2 ¡4
¡5 ¡4 6
¯¯¯¯¯¯
= 8
¯¯¯¯¯¯
1 1 ¡1
3 2 ¡4
¡2 ¡2 2
¯¯¯¯¯¯
= 0:
84
Let x = x1 + ®; y = y1 + ¯ and substitute in equation (4) to get
2(x1+®)2+(y1+¯)2+3(x1+®)(y1+¯)¡5(x1+®)¡4(y1+¯)+3 = 0 (5):
Then equating the coe±cients of x1 and y1 to 0 gives
4® + 3¯ ¡ 5 = 0
3® + 2¯ ¡ 4 = 0;
which has the unique solution ® = 2; ¯ = ¡1. Then equation (5) simpli¯es
to
2x2
1 + y2
1 + 3x1y1 = 0 = (2x1 + y1)(x1 + y1):
So relative to the x1; y1 coordinates, equation (4) describes two lines: 2x1+
y1 = 0 and x1 +y1 = 0. In terms of the original x; y coordinates, these lines
become 2(x ¡ 2) + (y +1) = 0 and (x ¡ 2) + (y + 1) = 0, i.e. 2x + y ¡ 3 = 0
and x + y ¡ 1 = 0, which intersect in the point
(x; y) = (®; ¯) = (2; ¡1):
(ii) Consider the equation
9x2 + y2 ¡ 6xy + 6x ¡ 2y + 1 = 0: (6)
Here
¢ =
¯¯¯¯¯¯
9 ¡3 3
3 1 ¡1
3 ¡1 1
¯¯¯¯¯¯
= 0;
as column 3 = ¡column 2.
Let x = x1 + ®; y = y1 + ¯ and substitute in equation (6) to get
9(x1 + ®)2 + (y1 + ¯)2 ¡ 6(x1 + ®)(y1 + ¯) + 6(x1 + ®) ¡ 2(y1 + ¯) + 1 = 0:
Then equating the coe±cients of x1 and y1 to 0 gives
18® ¡ 6¯ + 6 = 0
¡6® + 2¯ ¡ 2 = 0;
or equivalently ¡3®+¯ ¡1 = 0. Take ® = 0 and ¯ = 1. Then equation (6)
simpli¯es to
9x2
1 + y2
1 ¡ 6x1y1 = 0 = (3x1 ¡ y1)2: (7)
85
In terms of x; y coordinates, equation (7) becomes
(3x ¡ (y ¡ 1))2 = 0; or 3x ¡ y + 1 = 0:
(iii) Consider the equation
x2 + 4xy + 4y2 ¡ x ¡ 2y ¡ 2 = 0: (8)
Arguing as in the previous examples, we ¯nd that any translation
x = x1 + ®; y = y1 + ¯
where 2® + 4¯ ¡ 1 = 0 has the property that the coe±cients of x1 and y1
will be zero in the transformed version of equation (8). Take ¯ = 0 and
® = 1=2. Then (8) reduces to
x2
1 + 4x1y1 + 4y2
1 ¡
9
4
= 0;
or (x1 +2y1)2 = 3=2. Hence x1 +2y1 = §3=2, with corresponding equations
x + 2y = 2 and x + 2y = ¡1:
86
Section 8.8
1. The given line has equations
x = 3 + t(13 ¡ 3) = 3 + 10t;
y = ¡2 + t(3 + 2) = ¡2 + 5t;
z = 7 + t(¡8 ¡ 7) = 7 ¡ 15t:
The line meets the plane y = 0 in the point (x; 0; z), where 0 = ¡2+5t, or
t = 2=5. The corresponding values for x and z are 7 and 1, respectively.
2. E = 1
2 (B + C), F = (1 ¡ t)A + tE, where
t =
AF
AE
=
AF
AF + FE
=
AF=FE
(AF=FE) + 1
=
2
3
:
Hence
F =
1
3
A +
2
3
µ
1
2
(B + C)
¶
=
1
3
A +
1
3
(B + C)
=
1
3
(A + B + C):
3. Let A = (2; 1; 4); B = (1; ¡1; 2); C = (3; 3; 6). Then we prove
-
AC=
t
-
AB for some real t. We have
-
AC=
2
4
1
2
2
3
5 ;
-
AB=
2
4 ¡1
¡2
¡2
3
5 :
Hence
-
AC= (¡1)
-
AB and consequently C is on the line AB. In fact A is
between C and B, with AC = AB.
4. The points P on the line AB which satisfy AP = 2
5PB are given by
P = A + t
-
AB, where jt=(1 ¡ t)j = 2=5. Hence t=(1 ¡ t) = §2=5.
The equation t=(1 ¡ t) = 2=5 gives t = 2=7 and hence
P =
2
4
2
3
¡1
3
5 +
2
7
2
4
1
4
5
3
5 =
2
4
16=7
29=7
3=7
3
5 :
87
Hence P = (16=7; 29=7; 3=7).
The equation t=(1 ¡ t) = ¡2=5 gives t = ¡2=3 and hence
P =
2
4
2
3
¡1
3
5 ¡
2
3
2
4
1
4
5
3
5 =
2
4
4=3
1=3
¡13=3
3
5 :
Hence P = (4=3; 1=3; ¡13=3).
5. An equation for M is P = A + t
-
BC, which reduces to
x = 1 + 6t
y = 2 ¡ 3t
z = 3 + 7t:
An equation for N is Q = E + s
-
EF, which reduces to
x = 1 + 9s
y = ¡1
z = 8 + 3s:
To ¯nd if and whereMand N intersect, we set P = Q and attempt to solve
for s and t. We ¯nd the unique solution t = 1; s = 2=3, proving that the
lines meet in the point
(x; y; z) = (1 + 6; 2 ¡ 3; 3 + 7) = (7; ¡1; 10):
6. Let A = (3; 5; 6); B = (¡2; 7; 9); C = (2; 1; 7). Then
(i)
cos\ABC = (
-
BA ¢
-
BC)=(BA ¢ BC);
where
-
BA= [¡1; ¡2; ¡3]t and
-
BC= [4; ¡6; ¡2]t. Hence
cos\ABC = ¡4 + 12 + 6
p14p56
=
14
p14p56
=
1
2
:
Hence \ABC = ¼=3 radians or 60±.
88
(ii)
cos\BAC = (
-
AB ¢
-
AC)=(AB ¢ AC);
where
-
AB= [1; 2; 3]t and
-
AC= [5; ¡4; 1]t. Hence
cos\BAC =
5 ¡ 8 + 3
p14p42
= 0:
Hence \ABC = ¼=2 radians or 90±.
(iii)
cos\ACB = (
-
CA ¢
-
CB)=(CA ¢ CB);
where
-
CA= [¡5; 4; ¡1]t and
-
CB= [¡4; 6; 2]t. Hence
cos\ACB =
20 + 24 ¡ 2
p42p56
=
42
p42p56
=
p42
p56
=
p3
2
:
Hence \ACB = ¼=6 radians or 30±.
7. By Theorem 8.5.2, the closest point P on the line AB to the origin O is
given by P = A + t
-
AB, where
t =
-
AO ¢
-
AB
AB2 = ¡A¢
-
AB
AB2 :
Now
A¢
-
AB=
2
4 ¡2
1
3
3
5 ¢
2
4
3
1
1
3
5 = ¡2:
Hence t = 2=11 and
P =
2
4 ¡2
1
3
3
5 +
2
11
2
4
3
1
1
3
5 =
2
4 ¡16=11
13=11
35=11
3
5
and P = (¡16=11; 13=11; 35=11).
Consequently the shortest distance OP is given by
sµ
¡16
11
¶2
+
µ
13
11
¶2
+
µ
35
11
¶2
=
p1650
11
=
p15 £ 11 £ 10
11
=
p150
p11
:
89
Alternatively, we can calculate the distance OP 2, where P is an arbitrary
point on the line AB and then minimize OP2:
P = A + t
-
AB=
2
4 ¡2
1
3
3
5 + t
2
4
3
1
1
3
5 =
2
4 ¡2 + 3t
1 + t
3 + t
3
5 :
Hence
OP2 = (¡2 + 3t)2 + (1 + t)2 + (3 + t)2
= 11t2 ¡ 4t + 14
= 11
µ
t2 ¡
4
11
t +
14
11
¶
= 11
ý
t ¡
2
11
¾2
+
14
11 ¡
4
121
!
= 11
ý
t ¡
2
11
¾2
+
150
121
!
:
Consequently
OP2 ¸ 11 £
150
121
for all t; moreover
OP2 = 11 £
150
121
when t = 2=11.
8. We ¯rst ¯nd parametric equations for N by solving the equations
x + y ¡ 2z = 1
x + 3y ¡ z = 4:
The augmented matrix is
·
1 1 ¡2 1
1 3 ¡1 4
¸
;
which reduces to ·
1 0 ¡5=2 ¡1=2
0 1 1=2 3=2
¸
:
Hence x = ¡1
2 + 5
2z; y = 3
2 ¡ z
2 , with z arbitrary. Taking z = 0 gives a point
A = (¡1
2 ; 3
2 ; 0), while z = 1 gives a point B = (2; 1; 1).
90
Hence if C = (1; 0; 1), then the closest point on N to C is given by
P = A + t
-
AB, where t = (
-
AC ¢
-
AB)=AB2.
Now
-
AC=
2
4
3=2
¡3=2
1
3
5 and
-
AB=
2
4
5=2
¡1=2
1
3
5 ;
so
t =
3
2 £ 5
2 + ¡3
2 £ ¡1
2 + 1 £ 1
¡5
2
¢2
+
¡
¡1
2
¢2
+ 12
=
11
15
:
Hence
P =
2
4 ¡1=2
3=2
0
3
5 +
11
15
2
4
5=2
¡1=2
1
3
5 =
2
4
4=3
17=15
11=15
3
5 ;
so P = (4=3; 17=15; 11=15).
Also the shortest distance PC is given by
PC =
sµ
1 ¡
4
3
¶2
+
µ
0 ¡
17
15
¶2
+
µ
1 ¡
11
15
¶2
=
p330
15
:
9. The intersection of the planes x + y ¡ 2z = 4 and 3x ¡ 2y + z = 1 is the
line given by the equations
x =
9
5
+
3
5
z; y =
11
5
+
7
5
z;
where z is arbitrary. Hence the line L has a direction vector [3=5; 7=5; 1]t
or the simpler [3; 7; 5]t. Then any plane of the form 3x + 7y + 5z = d will
be perpendicualr to L. The required plane has to pass through the point
(6; 0; 2), so this determines d:
3 £ 6 + 7 £ 0 + 5 £ 2 = d = 28:
10. The length of the projection of the segment AB onto the line CD is
given by the formula
j
-
CD ¢
-
AB j
CD
:
Here
-
CD= [¡8; 4; ¡1]t and
-
AB= [4; ¡4; 3]t, so
j
-
CD ¢
-
AB j
CD
= j(¡8)p£ 4 + 4 £ (¡4) + (¡1) £ 3j
(¡8)2 + 42 + (¡1)2
= jp¡ 51j 81
=
51
9
=
17
3
:
91
11. A direction vector for L is given by
-
BC= [¡5; ¡2; 3]t. Hence the plane
through A perpendicular to L is given by
¡5x ¡ 2y + 3z = (¡5) £ 3 + (¡2) £ (¡1) + 3 £ 2 = ¡7:
The position vector P of an arbitrary point P on L is given by P = B+t
-
BC,
or 2
4
x
y
z
3
5 =
2
4
2
1
4
3
5 + t
2
4 ¡5
¡2
3
3
5 ;
or equivalently x = 2 ¡ 5t; y = 1 ¡ 2t; z = 4 + 3t.
To ¯nd the intersection of line L and the given plane, we substitute the
expressions for x; y; z found in terms of t into the plane equation and solve
the resulting linear equation for t:
¡5(2 ¡ 5t) ¡ 2(1 ¡ 2t) + 3(4 + 3t) = ¡7;
which gives t = ¡7=38. Hence P =
¡111
38 ; 52
38 ; 131
38
¢
and
AP =
sµ
3 ¡
111
38
¶2
+
µ
¡1 ¡
52
38
¶2
+
µ
2 ¡
131
38
¶2
=
p11134
38
=
p293 £ 38
38
=
p293
p38
:
12. Let P be a point inside the triangle ABC. Then the line through P and
parallel to AC will meet the segments AB and BC in D and E, respectively.
Then
P = (1 ¡ r)D + rE; 0 < r < 1;
D = (1 ¡ s)B + sA; 0 < s < 1;
E = (1 ¡ t)B + tC; 0 < t < 1:
Hence
P = (1 ¡ r) f(1 ¡ s)B + sAg + r f(1 ¡ t)B + tCg
= (1 ¡ r)sA + f(1 ¡ r)(1 ¡ s) + r(1 ¡ t)gB + rtC
= ®A + ¯B + °C;
92
where
® = (1 ¡ r)s; ¯ = (1 ¡ r)(1 ¡ s) + r(1 ¡ t); ° = rt:
Then 0 < ® < 1; 0 < ° < 1; 0 < ¯ < (1 ¡ r) + r = 1. Also
® + ¯ + ° = (1 ¡ r)s + (1 ¡ r)(1 ¡ s) + r(1 ¡ t) + rt = 1:
13. The line AB is given by P = A + t[3; 4; 5]t, or
x = 6 + 3t; y = ¡1 + 4t; z = 11 + 5t:
Then B is found by substituting these expressions in the plane equation
3x + 4y + 5z = 10:
We ¯nd t = ¡59=50 and consequently
B =
µ
6 ¡
177
50
; ¡1 ¡
236
50
; 11 ¡
295
50
¶
=
µ
123
50
; ¡286
50
;
255
50
¶
:
Then
AB = jj
-
AB jj = jjt
2
4
3
4
5
3
5 jj
= jtj
p
32 + 42 + 52 =
59
50 £
p50 =
59
p50
:
14. Let A = (¡3; 0; 2); B = (6; 1; 4); C = (¡5; 1; 0). Then the area of
triangle ABC is 1
2 jj
-
AB £
-
AC jj. Now
-
AB £
-
AC=
2
4
9
1
2
3
5 £
2
4 ¡2
1
¡2
3
5 =
2
4 ¡4
14
11
3
5 :
Hence jj
-
AB £
-
AC jj = p333:
15. Let A1 = (2; 1; 4); A2 = (1; ¡1; 2); A3 = (4; ¡1; 1). Then the point
P = (x; y; z) lies on the plane A1A2A3 if and only if
-
A1P ¢(
-
A1A2 £
-
A1A3) = 0;
93
or ¯¯¯¯¯¯
x ¡ 2 y ¡ 1 z ¡ 4
¡1 ¡2 ¡2
2 ¡2 ¡3
¯¯¯¯¯¯
= 2x ¡ 7y + 6z ¡ 21 = 0:
16. Non{parallel lines L and M in three dimensional space are given by
equations
P = A + sX; Q = B + tY:
(i) Suppose
-
PQ is orthogonal to both X and Y . Now
-
PQ= Q ¡ P = (B + tY ) ¡ (A + sX) =
-
AB +tY ¡ sX:
Hence
(
-
AB +tY + sX) ¢ X = 0
(
-
AB +tY + sX) ¢ Y = 0:
More explicitly
t(Y ¢ X) ¡ s(X ¢ X) = ¡
-
AB ¢X
t(Y ¢ Y ) ¡ s(X ¢ Y ) = ¡
-
AB ¢Y:
However the coe±cient determinant of this system of linear equations
in t and s is equal to
¯¯¯¯
Y ¢ X ¡X ¢ X
Y ¢ Y ¡X ¢ Y
¯¯¯¯
= ¡(X ¢ Y )2 + (X ¢ X)(Y ¢ Y )
= jjX £ Y jj2 6= 0;
as X 6= 0; Y 6= 0 and X and Y are not proportional (L and M are
not parallel).
(ii) P and Q can be viewed as the projections of C and D onto the line PQ,
where C and D are arbitrary points on the lines L andM, respectively.
Hence by equation (8.14) of Theorem 8.5.3, we have
PQ · CD:
Finally we derive a useful formula for PQ. Again by Theorem 8.5.3
PQ = j
-
AB ¢
-
PQ j
PQ
= j
-
AB ¢^nj;
94
¢
¢
¢
¢
¢
¢®
-
6
y
z
x
O
@
@
@
@
@
@
@
M
L
¡
¡
¡
¡
¡
¡
¡
¡
¡
¡
@
@
@
@
¡
¡
¡
¡ """"""""""""""""
"" "" ""
PPPPPPPPP
D
C
Q
P
@¡
""
where ^n = 1
P Q
-
PQ is a unit vector which is orthogonal to X and Y .
Hence
^n = t(X £ Y );
where t = §1=jjX £ Y jj. Hence
PQ = j
-
AB ¢(X £ Y )j
jjX £ Y jj
:
17. We use the formula of the previous question.
Line L has the equation P = A + sX, where
X =
-
AC=
2
4
2
¡3
3
3
5 :
Line M has the equation Q = B + tY , where
Y =
-
BD=
2
4
1
1
1
3
5 :
Hence X £ Y = [¡6; 1; 5]t and jjX £ Y jj = p62.
95
Hence the shortest distance between lines AC and BD is equal to
j
-
AB ¢(X £ Y )j
jjX £ Y jj
=
¯¯¯¯¯¯
2
4
0
¡2
1
3
5 ¢
2
4 ¡6
1
5
3
5
¯¯¯¯¯¯
p62
=
3
p62
:
18. Let E be the foot of the perpendicular from A4 to the plane A1A2A3.
Then
volA1A2A3A4 =
1
3
( area¢A1A2A3) ¢ A4E:
Now
area¢A1A2A3 =
1
2jj
-
A1A2 £
-
A1A3 jj:
Also A4E is the length of the projection of A1A4 onto the line A4E. See
¯gure below.)
Hence A4E = j
-
A1A4 ¢Xj, where X is a unit direction vector for the line
A4E. We can take
X =
-
A1A2 £
-
A1A3
jj
-
A1A2 £
-
A1A3 jj
:
Hence
volA1A2A3A4 =
1
6jj
-
A1A2 £
-
A1A3 jjj
-
A1A4 ¢(
-
A1A2 £
-
A1A3)j
jj
-
A1A2 £
-
A1A3 jj
=
1
6j
-
A1A4 ¢(
-
A1A2 £
-
A1A3)j
96
=
1
6j(
-
A1A2 £
-
A1A3)¢
-
A1A4 j:
19. We have
-
CB= [1; 4; ¡1]t;
-
CD= [¡3; 3; 0]t;
-
AD= [3; 0; 3]t. Hence
-
CB £
-
CD= 3i + 3j + 15k;
so the vector i + j + 5k is perpendicular to the plane BCD.
Now the plane BCD has equation x+y +5z = 9, as B = (2; 2; 1) is on
the plane.
Also the line through A normal to plane BCD has equation
2
4
x
y
z
3
5 =
2
4
1
1
5
3
5 + t
2
4
1
1
5
3
5 = (1 + t)
2
4
1
1
5
3
5 :
Hence x = 1 + t; y = 1 + t; z = 5(1 + t).
[We remark that this line meets plane BCD in a point E which is given
by a value of t found by solving
(1 + t) + (1 + t) + 5(5 + 5t) = 9:
So t = ¡2=3 and E = (1=3; 1=3; 5=3).]
The distance from A to plane BCD is
j1 £ 1 + 1 £ 1 + 5 £ 5 ¡ 9j
12 + 12 + 52 =
18
p27
= 2p3:
To ¯nd the distance between lines AD and BC, we ¯rst note that
(a) The equation of AD is
P =
2
4
1
1
5
3
5 + t
2
4
3
0
3
3
5 =
2
4
1 + 3t
1
5 + 3t
3
5 ;
(b) The equation of BC is
Q =
2
4
2
2
1
3
5 + s
2
4
1
4
¡1
3
5 =
2
4
2 + s
2 + 4s
1 ¡ s
3
5 :
97
@
@
@
@
@
¢
¢
¢
¢
¢
¢
¢
¢
¢
¢
©©©©©©©©©© @
@
@
@
@
¢
¢
¢
¢
¢
¢
¢
¢
¢
¢
C
C
C
C
C
C
C
C
C
D
B
E
Then
-
PQ= [1 + s ¡ 3t; 1 + 4s; ¡4 ¡ s ¡ 3t]t and we ¯nd s and t by solving
the equations
-
PQ ¢
-
AD= 0 and
-
PQ ¢
-
BC= 0, or
(1 + s ¡ 3t)3 + (1 + 4s)0 + (¡4 ¡ s ¡ 3t)3 = 0
(1 + s ¡ 3t) + 4(1 + 4s) ¡ (¡4 ¡ s ¡ 3t) = 0:
Hence t = ¡1=2 = s.
Correspondingly, P = (¡1=2; 1; 7=2) and Q = (3=2; 0; 3=2).
Thus we have found the closest points P and Q on the respective lines
AD and BC. Finally the shortest distance between the lines is
PQ = jj
-
PQ jj = 3:
98