Matrix Fundamentals
Definition and Notation
An m × n m \times n m × n matrix A A A is a rectangular array of real numbers with m m m rows and n n n columns:
A = ( a 11 a 12 ⋯ a 1 n a 21 a 22 ⋯ a 2 n ⋮ ⋮ ⋱ ⋮ a m 1 a m 2 ⋯ a m n ) A = \begin{pmatrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{m1} & a_{m2} & \cdots & a_{mn} \end{pmatrix} A = a 11 a 21 ⋮ a m 1 a 12 a 22 ⋮ a m 2 ⋯ ⋯ ⋱ ⋯ a 1 n a 2 n ⋮ a mn
The entry in row i i i and column j j j is a i j a_{ij} a ij . When m = n m = n m = n , A A A is a square matrix of order
n n n .
The identity matrix I n I_n I n is the n × n n \times n n × n matrix with 1 1 1 s on the diagonal and 0 0 0 s
elsewhere. For any compatible matrix A A A : A I = I A = A AI = IA = A A I = I A = A .
Matrix Operations
Addition and Scalar Multiplication
If A , B ∈ ′ { ′ M ′ } ′ m × n ( ′ { ′ R ′ } ′ ) A, B \in \mathcal{'\{'}M{'\}'}_{m \times n}(\mathbb{'\{'}R{'\}'}) A , B ∈ ′ { ′ M ′ } ′ m × n ( ′ { ′ R ′ } ′ ) and λ ∈ ′ { ′ R ′ } ′ \lambda \in \mathbb{'\{'}R{'\}'} λ ∈ ′ { ′ R ′ } ′ :
( A + B ) i j = a i j + b i j , ( λ A ) i j = λ ⋅ a i j (A + B)_{ij} = a_{ij} + b_{ij}, \qquad (\lambda A)_{ij} = \lambda \cdot a_{ij} ( A + B ) ij = a ij + b ij , ( λ A ) ij = λ ⋅ a ij
Matrix Multiplication
If A A A is m × p m \times p m × p and B B B is p × n p \times n p × n , the product C = A B C = AB C = A B is m × n m \times n m × n with:
c i j = ∑ k = 1 p a i k b k j c_{ij} = \sum_{k=1}^{p} a_{ik} b_{kj} c ij = ∑ k = 1 p a ik b k j
Key properties:
Associative: ( A B ) C = A ( B C ) (AB)C = A(BC) ( A B ) C = A ( B C )
Distributive: A ( B + C ) = A B + A C A(B + C) = AB + AC A ( B + C ) = A B + A C
NOT commutative: A B ≠ B A AB \ne BA A B = B A in general
Transpose of product: ( A B ) T = B T A T (AB)^T = B^T A^T ( A B ) T = B T A T
Transpose
The transpose A T A^T A T is obtained by reflecting A A A across its main diagonal:
( A T ) i j = a j i (A^T)_{ij} = a_{ji} ( A T ) ij = a j i .
Properties:
( A T ) T = A (A^T)^T = A ( A T ) T = A
( A + B ) T = A T + B T (A + B)^T = A^T + B^T ( A + B ) T = A T + B T
( λ A ) T = λ A T (\lambda A)^T = \lambda A^T ( λ A ) T = λ A T
( A B ) T = B T A T (AB)^T = B^T A^T ( A B ) T = B T A T
Determinants
2 × 2 2 \times 2 2 × 2 Determinant
det ( a b c d ) = a d − b c \det\begin{pmatrix} a & b \\ c & d \end{pmatrix} = ad - bc det ( a c b d ) = a d − b c
3 × 3 3 \times 3 3 × 3 Determinant
Expand along any row or column. Expanding along row 1:
det ( A ) = a 11 ∣ a 22 a 23 a 32 a 33 ∣ − a 12 ∣ a 21 a 23 a 31 a 33 ∣ + a 13 ∣ a 21 a 22 a 31 a 32 ∣ \det(A) = a_{11}\begin{vmatrix} a_{22} & a_{23} \\ a_{32} & a_{33} \end{vmatrix} - a_{12}\begin{vmatrix} a_{21} & a_{23} \\ a_{31} & a_{33} \end{vmatrix} + a_{13}\begin{vmatrix} a_{21} & a_{22} \\ a_{31} & a_{32} \end{vmatrix} det ( A ) = a 11 a 22 a 32 a 23 a 33 − a 12 a 21 a 31 a 23 a 33 + a 13 a 21 a 31 a 22 a 32
The signs follow the checkerboard pattern + − + + - + + − + along row 1, − + − - + - − + − along row 2, etc.
Properties of Determinants
det ( A B ) = det ( A ) det ( B ) \det(AB) = \det(A)\det(B) det ( A B ) = det ( A ) det ( B )
det ( A T ) = det ( A ) \det(A^T) = \det(A) det ( A T ) = det ( A )
det ( λ A ) = λ n det ( A ) \det(\lambda A) = \lambda^n \det(A) det ( λ A ) = λ n det ( A ) for an n × n n \times n n × n matrix
Swapping two rows multiplies the determinant by − 1 -1 − 1
A matrix with a zero row (or column) has determinant 0 0 0
A triangular matrix has determinant equal to the product of its diagonal entries
Singular and Non-Singular Matrices
A square matrix A A A is singular if det ( A ) = 0 \det(A) = 0 det ( A ) = 0 and non-singular (invertible) if
det ( A ) ≠ 0 \det(A) \ne 0 det ( A ) = 0 .
Inverse Matrices
Definition
The inverse of a square matrix A A A is the matrix A − 1 A^{-1} A − 1 satisfying:
A A − 1 = A − 1 A = I AA^{-1} = A^{-1}A = I A A − 1 = A − 1 A = I
A A A is invertible if and only if det ( A ) ≠ 0 \det(A) \ne 0 det ( A ) = 0 .
2 × 2 2 \times 2 2 × 2 Inverse
A − 1 = 1 det ( A ) ( d − b − c a ) A^{-1} = \frac{1}{\det(A)}\begin{pmatrix} d & -b \\ -c & a \end{pmatrix} A − 1 = d e t ( A ) 1 ( d − c − b a )
Adjugate Method (3 × 3 3 \times 3 3 × 3 )
A − 1 = 1 det ( A ) ⋅ a d j ( A ) A^{-1} = \frac{1}{\det(A)} \cdot \mathrm{adj}(A) A − 1 = d e t ( A ) 1 ⋅ adj ( A )
where the adjugate (classical adjoint) is the transpose of the cofactor matrix. The
( i , j ) (i, j) ( i , j ) -cofactor is C i j = ( − 1 ) i + j M i j C_{ij} = (-1)^{i+j} M_{ij} C ij = ( − 1 ) i + j M ij , where M i j M_{ij} M ij is the ( i , j ) (i, j) ( i , j ) -minor (determinant
of the submatrix obtained by deleting row i i i and column j j j ).
Properties of Inverses
( A − 1 ) − 1 = A (A^{-1})^{-1} = A ( A − 1 ) − 1 = A
( A B ) − 1 = B − 1 A − 1 (AB)^{-1} = B^{-1}A^{-1} ( A B ) − 1 = B − 1 A − 1
( A T ) − 1 = ( A − 1 ) T (A^T)^{-1} = (A^{-1})^T ( A T ) − 1 = ( A − 1 ) T
( λ A ) − 1 = 1 λ A − 1 (\lambda A)^{-1} = \dfrac{1}{\lambda} A^{-1} ( λ A ) − 1 = λ 1 A − 1 , λ ≠ 0 \lambda \ne 0 λ = 0
Solving Systems of Linear Equations
A system of n n n linear equations in n n n unknowns:
A ′ { ′ x ′ } ′ = ′ { ′ b ′ } ′ A\mathbf{'\{'}x{'\}'} = \mathbf{'\{'}b{'\}'} A ′ { ′ x ′ } ′ = ′ { ′ b ′ } ′
has the unique solution ′ { ′ x ′ } ′ = A − 1 ′ { ′ b ′ } ′ \mathbf{'\{'}x{'\}'} = A^{-1}\mathbf{'\{'}b{'\}'} ′ { ′ x ′ } ′ = A − 1 ′ { ′ b ′ } ′ when det ( A ) ≠ 0 \det(A) \ne 0 det ( A ) = 0 .
Cramer's Rule
For the system A ′ { ′ x ′ } ′ = ′ { ′ b ′ } ′ A\mathbf{'\{'}x{'\}'} = \mathbf{'\{'}b{'\}'} A ′ { ′ x ′ } ′ = ′ { ′ b ′ } ′ where A A A is n × n n \times n n × n and det ( A ) ≠ 0 \det(A) \ne 0 det ( A ) = 0 :
x i = det ( A i ) det ( A ) x_i = \frac{\det(A_i)}{\det(A)} x i = d e t ( A ) d e t ( A i )
where A i A_i A i is the matrix A A A with column i i i replaced by ′ { ′ b ′ } ′ \mathbf{'\{'}b{'\}'} ′ { ′ b ′ } ′ .
Gaussian Elimination
For larger systems, use elementary row operations to reduce the augmented matrix
[ A ∣ ′ { ′ b ′ } ′ ] [A \mid \mathbf{'\{'}b{'\}'}] [ A ∣ ′ { ′ b ′ } ′ ] to row echelon form, then solve by back-substitution. The three elementary
operations are:
Swap two rows (R i ↔ R j R_i \leftrightarrow R_j R i ↔ R j )
Multiply a row by a nonzero scalar (R i → k R i R_i \to kR_i R i → k R i )
Add a multiple of one row to another (R i → R i + k R j R_i \to R_i + kR_j R i → R i + k R j )
Existence and Uniqueness
det ( A ) ≠ 0 \det(A) \ne 0 det ( A ) = 0 : unique solution
det ( A ) = 0 \det(A) = 0 det ( A ) = 0 and system is consistent: infinitely many solutions
det ( A ) = 0 \det(A) = 0 det ( A ) = 0 and system is inconsistent: no solution
Each transformation matrix T T T maps a point ( x , y ) (x, y) ( x , y ) to a new point by left-multiplication:
( x ′ y ′ ) = T ( x y ) \begin{pmatrix} x' \\ y' \end{pmatrix} = T\begin{pmatrix} x \\ y \end{pmatrix} ( x ′ y ′ ) = T ( x y )
Transformation Matrix Reflection in x x x -axis ( 1 0 0 − 1 ) \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} ( 1 0 0 − 1 ) Reflection in y y y -axis ( − 1 0 0 1 ) \begin{pmatrix} -1 & 0 \\ 0 & 1 \end{pmatrix} ( − 1 0 0 1 ) Reflection in y = x y = x y = x ( 0 1 1 0 ) \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} ( 0 1 1 0 ) Rotation by θ \theta θ about origin ( cos θ − sin θ sin θ cos θ ) \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} ( cos θ sin θ − sin θ cos θ ) Enlargement by factor k k k ( k 0 0 k ) \begin{pmatrix} k & 0 \\ 0 & k \end{pmatrix} ( k 0 0 k ) Horizontal stretch by factor k k k ( k 0 0 1 ) \begin{pmatrix} k & 0 \\ 0 & 1 \end{pmatrix} ( k 0 0 1 ) Vertical stretch by factor k k k ( 1 0 0 k ) \begin{pmatrix} 1 & 0 \\ 0 & k \end{pmatrix} ( 1 0 0 k ) Shear parallel to x x x -axis, factor k k k ( 1 k 0 1 ) \begin{pmatrix} 1 & k \\ 0 & 1 \end{pmatrix} ( 1 0 k 1 )
Determinant and Area
For a 2 × 2 2 \times 2 2 × 2 transformation matrix T T T :
∣ det ( T ) ∣ |\det(T)| ∣ det ( T ) ∣ is the area scale factor
det ( T ) > 0 \det(T) \gt 0 det ( T ) > 0 : orientation preserved
det ( T ) < 0 \det(T) \lt 0 det ( T ) < 0 : orientation reversed
Applying transformation A A A followed by B B B corresponds to the matrix product B A BA B A (right to left
order):
′ { ′ v ′ } ′ ′ = B ( A ′ { ′ v ′ } ′ ) = ( B A ) ′ { ′ v ′ } ′ \mathbf{'\{'}v{'\}'}' = B(A\mathbf{'\{'}v{'\}'}) = (BA)\mathbf{'\{'}v{'\}'} ′ { ′ v ′ } ′ ′ = B ( A ′ { ′ v ′ } ′ ) = ( B A ) ′ { ′ v ′ } ′
Invariant Points and Lines
An invariant point satisfies T ′ { ′ v ′ } ′ = ′ { ′ v ′ } ′ T\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}v{'\}'} T ′ { ′ v ′ } ′ = ′ { ′ v ′ } ′ , i.e. ( T − I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ (T - I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} ( T − I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ .
An invariant line is mapped to itself (points may move along the line).
Eigenvalues and Eigenvectors (HL)
Definition
Let A A A be an n × n n \times n n × n matrix. A scalar λ \lambda λ is an eigenvalue of A A A if there exists a
nonzero vector ′ { ′ v ′ } ′ \mathbf{'\{'}v{'\}'} ′ { ′ v ′ } ′ (an eigenvector ) such that:
A ′ { ′ v ′ } ′ = λ ′ { ′ v ′ } ′ A\mathbf{'\{'}v{'\}'} = \lambda\mathbf{'\{'}v{'\}'} A ′ { ′ v ′ } ′ = λ ′ { ′ v ′ } ′
Geometrically, A A A stretches or compresses the eigenvector ′ { ′ v ′ } ′ \mathbf{'\{'}v{'\}'} ′ { ′ v ′ } ′ by factor λ \lambda λ without
changing its direction.
Characteristic Equation
Rearranging A ′ { ′ v ′ } ′ = λ ′ { ′ v ′ } ′ A\mathbf{'\{'}v{'\}'} = \lambda\mathbf{'\{'}v{'\}'} A ′ { ′ v ′ } ′ = λ ′ { ′ v ′ } ′ gives ( A − λ I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ (A - \lambda I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} ( A − λ I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ . For a
nonzero solution ′ { ′ v ′ } ′ \mathbf{'\{'}v{'\}'} ′ { ′ v ′ } ′ , we require:
det ( A − λ I ) = 0 \det(A - \lambda I) = 0 det ( A − λ I ) = 0
This is the characteristic equation . Its roots are the eigenvalues of A A A .
For a 2 × 2 2 \times 2 2 × 2 matrix A = ( a b c d ) A = \begin{pmatrix} a & b \\ c & d \end{pmatrix} A = ( a c b d ) :
det ( a − λ b c d − λ ) = ( a − λ ) ( d − λ ) − b c = 0 \det\begin{pmatrix} a - \lambda & b \\ c & d - \lambda \end{pmatrix} = (a - \lambda)(d - \lambda) - bc = 0 det ( a − λ c b d − λ ) = ( a − λ ) ( d − λ ) − b c = 0
λ 2 − ( a + d ) λ + ( a d − b c ) = 0 \lambda^2 - (a + d)\lambda + (ad - bc) = 0 λ 2 − ( a + d ) λ + ( a d − b c ) = 0
Note: a + d = t r ( A ) a + d = \mathrm{tr}(A) a + d = tr ( A ) (the trace) and a d − b c = det ( A ) ad - bc = \det(A) a d − b c = det ( A ) , so:
λ 2 − t r ( A ) λ + det ( A ) = 0 \lambda^2 - \mathrm{tr}(A)\,\lambda + \det(A) = 0 λ 2 − tr ( A ) λ + det ( A ) = 0
Finding Eigenvectors
For each eigenvalue λ i \lambda_i λ i , substitute into ( A − λ i I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ (A - \lambda_i I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} ( A − λ i I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ and
solve for ′ { ′ v ′ } ′ \mathbf{'\{'}v{'\}'} ′ { ′ v ′ } ′ . The solution space is the eigenspace corresponding to λ i \lambda_i λ i .
Properties
The trace equals the sum of eigenvalues: t r ( A ) = ∑ λ i \mathrm{tr}(A) = \sum \lambda_i tr ( A ) = ∑ λ i
The determinant equals the product of eigenvalues: det ( A ) = ∏ λ i \det(A) = \prod \lambda_i det ( A ) = ∏ λ i
If A A A is symmetric (A = A T A = A^T A = A T ), all eigenvalues are real and eigenvectors corresponding to
distinct eigenvalues are orthogonal.
If A A A has n n n linearly independent eigenvectors, it is diagonalisable : A = P D P − 1 A = PDP^{-1} A = P D P − 1 where
D D D is diagonal with eigenvalues on the diagonal.
Diagonalisation (HL Extension)
If A A A has eigenvalues λ 1 , … , λ n \lambda_1, \ldots, \lambda_n λ 1 , … , λ n and corresponding eigenvectors
′ { ′ v ′ } ′ 1 , … , ′ { ′ v ′ } ′ n \mathbf{'\{'}v{'\}'}_1, \ldots, \mathbf{'\{'}v{'\}'}_n ′ { ′ v ′ } ′ 1 , … , ′ { ′ v ′ } ′ n forming a linearly independent set, then:
A = P D P − 1 A = PDP^{-1} A = P D P − 1
where P = ( ′ { ′ v ′ } ′ 1 ⋯ ′ { ′ v ′ } ′ n ) P = \begin{pmatrix} \mathbf{'\{'}v{'\}'}_1 & \cdots & \mathbf{'\{'}v{'\}'}_n \end{pmatrix} P = ( ′ { ′ v ′ } ′ 1 ⋯ ′ { ′ v ′ } ′ n ) and
D = d i a g ( λ 1 , … , λ n ) D = \mathrm{diag}(\lambda_1, \ldots, \lambda_n) D = diag ( λ 1 , … , λ n ) .
This is useful for computing large powers: A n = P D n P − 1 A^n = PD^nP^{-1} A n = P D n P − 1 .
Common Pitfalls
Matrix multiplication order. A B ≠ B A AB \ne BA A B = B A in general. When composing transformations, the
rightmost matrix is applied first.
Transpose of a product. ( A B ) T = B T A T (AB)^T = B^T A^T ( A B ) T = B T A T , not A T B T A^T B^T A T B T . The order reverses.
Determinant of a scalar multiple. det ( k A ) = k n det ( A ) \det(kA) = k^n\det(A) det ( k A ) = k n det ( A ) for an n × n n \times n n × n matrix, not
k det ( A ) k\det(A) k det ( A ) .
Forgetting to check invertibility. Only non-singular matrices have inverses. Always verify
det ( A ) ≠ 0 \det(A) \ne 0 det ( A ) = 0 before attempting to compute A − 1 A^{-1} A − 1 .
Eigenvector scaling. Eigenvectors are determined only up to a nonzero scalar multiple. Any
nonzero multiple of an eigenvector is also an eigenvector.
Confusing invariant points with fixed points of composition. A point fixed by A A A is not
necessarily fixed by B A BA B A .
Cofactor sign errors. The ( i , j ) (i, j) ( i , j ) -cofactor sign is ( − 1 ) i + j (-1)^{i+j} ( − 1 ) i + j . Losing track of the sign
pattern is the most common error in computing adjugates.
Practice Problems
Problem 1 Let A = ( 2 3 1 4 ) A = \begin{pmatrix} 2 & 3 \\ 1 & 4 \end{pmatrix} A = ( 2 1 3 4 ) and
B = ( 1 − 1 2 0 ) B = \begin{pmatrix} 1 & -1 \\ 2 & 0 \end{pmatrix} B = ( 1 2 − 1 0 ) . Find A B AB A B and B A BA B A . Is A B = B A AB = BA A B = B A ?
Problem 2 Find det ( A ) \det(A) det ( A ) and A − 1 A^{-1} A − 1 where
A = ( 1 2 0 3 1 − 1 2 0 1 ) A = \begin{pmatrix} 1 & 2 & 0 \\ 3 & 1 & -1 \\ 2 & 0 & 1 \end{pmatrix} A = 1 3 2 2 1 0 0 − 1 1 .
Problem 3 Solve the system using an inverse matrix: x + 2 y = 5 x + 2y = 5 x + 2 y = 5 3 x + 7 y = 18 3x + 7y = 18 3 x + 7 y = 18
Problem 4 Find the matrix representing a rotation of 90 ∘ 90^\circ 9 0 ∘ anticlockwise about the origin, and verify
that applying it twice gives a rotation of 180 ∘ 180^\circ 18 0 ∘ .
Problem 5 Find the eigenvalues and eigenvectors of A = ( 4 1 2 3 ) A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix} A = ( 4 2 1 3 ) .
Problem 6 The transformation represented by matrix T = ( 3 1 0 2 ) T = \begin{pmatrix} 3 & 1 \\ 0 & 2 \end{pmatrix} T = ( 3 0 1 2 ) is
applied to the unit square. Find the area of the image.
Problem 7 Use Cramer's rule to solve: 2 x − y + z = 3 2x - y + z = 3 2 x − y + z = 3 x + 3 y − z = 1 x + 3y - z = 1 x + 3 y − z = 1 3 x + y + 2 z = 7 3x + y + 2z = 7 3 x + y + 2 z = 7
Problem 8 For the matrix A = ( 5 − 6 2 − 2 ) A = \begin{pmatrix} 5 & -6 \\ 2 & -2 \end{pmatrix} A = ( 5 2 − 6 − 2 ) , find A 10 A^{10} A 10 by
diagonalisation.
Answers to Selected Problems Problem 1: A B = ( 8 − 2 9 − 1 ) AB = \begin{pmatrix} 8 & -2 \\ 9 & -1 \end{pmatrix} A B = ( 8 9 − 2 − 1 ) ,
B A = ( 1 − 1 4 6 ) BA = \begin{pmatrix} 1 & -1 \\ 4 & 6 \end{pmatrix} B A = ( 1 4 − 1 6 ) . A B ≠ B A AB \ne BA A B = B A , confirming non-commutativity.
Problem 2: Expanding along row 1: det ( A ) = 1 ( 1 − 0 ) − 2 ( 3 − ( − 2 ) ) + 0 = 1 − 10 = − 9 \det(A) = 1(1 - 0) - 2(3 - (-2)) + 0 = 1 - 10 = -9 det ( A ) = 1 ( 1 − 0 ) − 2 ( 3 − ( − 2 )) + 0 = 1 − 10 = − 9 .
Cofactors:
C 11 = 1 , C 12 = − 5 , C 13 = − 2 , C 21 = − 2 , C 22 = 1 , C 23 = 4 , C 31 = − 2 , C 32 = 1 , C 33 = − 5 C_{11} = 1, C_{12} = -5, C_{13} = -2, C_{21} = -2, C_{22} = 1, C_{23} = 4, C_{31} = -2, C_{32} = 1, C_{33} = -5 C 11 = 1 , C 12 = − 5 , C 13 = − 2 , C 21 = − 2 , C 22 = 1 , C 23 = 4 , C 31 = − 2 , C 32 = 1 , C 33 = − 5 .
A − 1 = 1 − 9 ( 1 2 − 2 − 5 1 1 − 2 4 − 5 ) = ( − 1 / 9 − 2 / 9 2 / 9 5 / 9 − 1 / 9 − 1 / 9 2 / 9 − 4 / 9 5 / 9 ) A^{-1} = \dfrac{1}{-9}\begin{pmatrix} 1 & 2 & -2 \\ -5 & 1 & 1 \\ -2 & 4 & -5 \end{pmatrix} = \begin{pmatrix} -1/9 & -2/9 & 2/9 \\ 5/9 & -1/9 & -1/9 \\ 2/9 & -4/9 & 5/9 \end{pmatrix} A − 1 = − 9 1 1 − 5 − 2 2 1 4 − 2 1 − 5 = − 1/9 5/9 2/9 − 2/9 − 1/9 − 4/9 2/9 − 1/9 5/9 .
Problem 3: A = ( 1 2 3 7 ) A = \begin{pmatrix} 1 & 2 \\ 3 & 7 \end{pmatrix} A = ( 1 3 2 7 ) , det ( A ) = 7 − 6 = 1 \det(A) = 7 - 6 = 1 det ( A ) = 7 − 6 = 1 .
A − 1 = ( 7 − 2 − 3 1 ) A^{-1} = \begin{pmatrix} 7 & -2 \\ -3 & 1 \end{pmatrix} A − 1 = ( 7 − 3 − 2 1 ) .
( x y ) = ( 7 − 2 − 3 1 ) ( 5 18 ) = ( 35 − 36 − 15 + 18 ) = ( − 1 3 ) \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 7 & -2 \\ -3 & 1 \end{pmatrix}\begin{pmatrix} 5 \\ 18 \end{pmatrix} = \begin{pmatrix} 35 - 36 \\ -15 + 18 \end{pmatrix} = \begin{pmatrix} -1 \\ 3 \end{pmatrix} ( x y ) = ( 7 − 3 − 2 1 ) ( 5 18 ) = ( 35 − 36 − 15 + 18 ) = ( − 1 3 ) .
So x = − 1 , y = 3 x = -1, y = 3 x = − 1 , y = 3 .
Problem 4: Rotation by 90 ∘ 90^\circ 9 0 ∘ : R = ( 0 − 1 1 0 ) R = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} R = ( 0 1 − 1 0 ) .
R 2 = ( 0 − 1 1 0 ) ( 0 − 1 1 0 ) = ( − 1 0 0 − 1 ) R^2 = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} = \begin{pmatrix} -1 & 0 \\ 0 & -1 \end{pmatrix} R 2 = ( 0 1 − 1 0 ) ( 0 1 − 1 0 ) = ( − 1 0 0 − 1 ) ,
which is a 180 ∘ 180^\circ 18 0 ∘ rotation.
Problem 5: Characteristic equation: λ 2 − 7 λ + 10 = 0 \lambda^2 - 7\lambda + 10 = 0 λ 2 − 7 λ + 10 = 0 , so λ = 2 \lambda = 2 λ = 2 or
λ = 5 \lambda = 5 λ = 5 . For λ = 2 \lambda = 2 λ = 2 :
( A − 2 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 2 1 2 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ = ( 1 − 2 ) (A - 2I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \begin{pmatrix} 2 & 1 \\ 2 & 1 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'} = \begin{pmatrix} 1 \\ -2 \end{pmatrix} ( A − 2 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 2 2 1 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ = ( 1 − 2 )
(up to scalar). For λ = 5 \lambda = 5 λ = 5 :
( A − 5 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( − 1 1 2 − 2 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ = ( 1 1 ) (A - 5I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \begin{pmatrix} -1 & 1 \\ 2 & -2 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'} = \begin{pmatrix} 1 \\ 1 \end{pmatrix} ( A − 5 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( − 1 2 1 − 2 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ = ( 1 1 ) .
Problem 6: Area = ∣ det ( T ) ∣ = ∣ 6 − 0 ∣ = 6 = |\det(T)| = |6 - 0| = 6 = ∣ det ( T ) ∣ = ∣6 − 0∣ = 6 .
Worked Examples
Worked Example: Finding the Inverse of a 3 × 3 3 \times 3 3 × 3 Matrix
Find the inverse of A = ( 2 1 0 1 3 − 1 0 1 2 ) A = \begin{pmatrix} 2 & 1 & 0 \\ 1 & 3 & -1 \\ 0 & 1 & 2 \end{pmatrix} A = 2 1 0 1 3 1 0 − 1 2 .
Solution Step 1: Compute det ( A ) \det(A) det ( A ) . Expanding along row 1:
det ( A ) = 2 ∣ 3 − 1 1 2 ∣ − 1 ∣ 1 − 1 0 2 ∣ + 0 = 2 ( 7 ) − 1 ( 2 ) = 14 − 2 = 12 \det(A) = 2\begin{vmatrix} 3 & -1 \\ 1 & 2 \end{vmatrix} - 1\begin{vmatrix} 1 & -1 \\ 0 & 2 \end{vmatrix} + 0 = 2(7) - 1(2) = 14 - 2 = 12 det ( A ) = 2 3 1 − 1 2 − 1 1 0 − 1 2 + 0 = 2 ( 7 ) − 1 ( 2 ) = 14 − 2 = 12
Step 2: Compute the cofactor matrix.
C 11 = + ∣ 3 − 1 1 2 ∣ = 7 C_{11} = +\begin{vmatrix} 3 & -1 \\ 1 & 2 \end{vmatrix} = 7 C 11 = + 3 1 − 1 2 = 7 , C 12 = − ∣ 1 − 1 0 2 ∣ = − 2 C_{12} = -\begin{vmatrix} 1 & -1 \\ 0 & 2 \end{vmatrix} = -2 C 12 = − 1 0 − 1 2 = − 2 , C 13 = + ∣ 1 3 0 1 ∣ = 1 C_{13} = +\begin{vmatrix} 1 & 3 \\ 0 & 1 \end{vmatrix} = 1 C 13 = + 1 0 3 1 = 1
C 21 = − ∣ 1 0 1 2 ∣ = − 2 C_{21} = -\begin{vmatrix} 1 & 0 \\ 1 & 2 \end{vmatrix} = -2 C 21 = − 1 1 0 2 = − 2 , C 22 = + ∣ 2 0 0 2 ∣ = 4 C_{22} = +\begin{vmatrix} 2 & 0 \\ 0 & 2 \end{vmatrix} = 4 C 22 = + 2 0 0 2 = 4 , C 23 = − ∣ 2 1 0 1 ∣ = − 2 C_{23} = -\begin{vmatrix} 2 & 1 \\ 0 & 1 \end{vmatrix} = -2 C 23 = − 2 0 1 1 = − 2
C 31 = + ∣ 1 0 3 − 1 ∣ = − 1 C_{31} = +\begin{vmatrix} 1 & 0 \\ 3 & -1 \end{vmatrix} = -1 C 31 = + 1 3 0 − 1 = − 1 , C 32 = − ∣ 2 0 1 − 1 ∣ = 2 C_{32} = -\begin{vmatrix} 2 & 0 \\ 1 & -1 \end{vmatrix} = 2 C 32 = − 2 1 0 − 1 = 2 , C 33 = + ∣ 2 1 1 3 ∣ = 5 C_{33} = +\begin{vmatrix} 2 & 1 \\ 1 & 3 \end{vmatrix} = 5 C 33 = + 2 1 1 3 = 5
Cofactor matrix: ( 7 − 2 1 − 2 4 − 2 − 1 2 5 ) \begin{pmatrix} 7 & -2 & 1 \\ -2 & 4 & -2 \\ -1 & 2 & 5 \end{pmatrix} 7 − 2 − 1 − 2 4 2 1 − 2 5 .
Step 3: Transpose to get the adjugate.
a d j ( A ) = ( 7 − 2 − 1 − 2 4 2 1 − 2 5 ) \mathrm{adj}(A) = \begin{pmatrix} 7 & -2 & -1 \\ -2 & 4 & 2 \\ 1 & -2 & 5 \end{pmatrix} adj ( A ) = 7 − 2 1 − 2 4 − 2 − 1 2 5
Step 4: Divide by the determinant.
A − 1 = 1 12 ( 7 − 2 − 1 − 2 4 2 1 − 2 5 ) A^{-1} = \frac{1}{12}\begin{pmatrix} 7 & -2 & -1 \\ -2 & 4 & 2 \\ 1 & -2 & 5 \end{pmatrix} A − 1 = 12 1 7 − 2 1 − 2 4 − 2 − 1 2 5
Worked Example: Composition of Transformations
Find the matrix for the transformation that first reflects in the line y = x y = x y = x , then rotates 90 ∘ 90^\circ 9 0 ∘
anticlockwise about the origin. Verify the result by applying the combined transformation to the point
( 3 , 1 ) (3, 1) ( 3 , 1 ) .
Solution Reflection in y = x y = x y = x : R f = ( 0 1 1 0 ) R_f = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} R f = ( 0 1 1 0 ) .
Rotation 90 ∘ 90^\circ 9 0 ∘ anticlockwise: R θ = ( 0 − 1 1 0 ) R_\theta = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} R θ = ( 0 1 − 1 0 ) .
The combined transformation (reflect first, then rotate) is:
T = R θ ⋅ R f = ( 0 − 1 1 0 ) ( 0 1 1 0 ) = ( − 1 0 0 1 ) T = R_\theta \cdot R_f = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} = \begin{pmatrix} -1 & 0 \\ 0 & 1 \end{pmatrix} T = R θ ⋅ R f = ( 0 1 − 1 0 ) ( 0 1 1 0 ) = ( − 1 0 0 1 )
This is a reflection in the y y y -axis.
Verification: ( 3 , 1 ) (3, 1) ( 3 , 1 ) reflected in y = x y = x y = x gives ( 1 , 3 ) (1, 3) ( 1 , 3 ) . Rotating ( 1 , 3 ) (1, 3) ( 1 , 3 ) by 90 ∘ 90^\circ 9 0 ∘
anticlockwise: R θ ( 1 3 ) = ( − 3 1 ) R_\theta\begin{pmatrix} 1 \\ 3 \end{pmatrix} = \begin{pmatrix} -3 \\ 1 \end{pmatrix} R θ ( 1 3 ) = ( − 3 1 ) .
Direct application: T ( 3 1 ) = ( − 3 1 ) T\begin{pmatrix} 3 \\ 1 \end{pmatrix} = \begin{pmatrix} -3 \\ 1 \end{pmatrix} T ( 3 1 ) = ( − 3 1 ) . Confirmed.
Worked Example: Solving a System Using Gaussian Elimination
Solve the system using Gaussian elimination:
2 x + y − z = 3 2x + y - z = 3 2 x + y − z = 3
x − y + 2 z = 6 x - y + 2z = 6 x − y + 2 z = 6
3 x + 2 y + z = 1 3x + 2y + z = 1 3 x + 2 y + z = 1
Solution Augmented matrix:
( 2 1 − 1 ∣ 3 1 − 1 2 ∣ 6 3 2 1 ∣ 1 ) \begin{pmatrix} 2 & 1 & -1 & | & 3 \\ 1 & -1 & 2 & | & 6 \\ 3 & 2 & 1 & | & 1 \end{pmatrix} 2 1 3 1 − 1 2 − 1 2 1 ∣ ∣ ∣ 3 6 1
R 1 ↔ R 2 R_1 \leftrightarrow R_2 R 1 ↔ R 2 :
( 1 − 1 2 ∣ 6 2 1 − 1 ∣ 3 3 2 1 ∣ 1 ) \begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 2 & 1 & -1 & | & 3 \\ 3 & 2 & 1 & | & 1 \end{pmatrix} 1 2 3 − 1 1 2 2 − 1 1 ∣ ∣ ∣ 6 3 1
R 2 → R 2 − 2 R 1 R_2 \to R_2 - 2R_1 R 2 → R 2 − 2 R 1 , R 3 → R 3 − 3 R 1 R_3 \to R_3 - 3R_1 R 3 → R 3 − 3 R 1 :
( 1 − 1 2 ∣ 6 0 3 − 5 ∣ − 9 0 5 − 5 ∣ − 17 ) \begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 0 & 3 & -5 & | & -9 \\ 0 & 5 & -5 & | & -17 \end{pmatrix} 1 0 0 − 1 3 5 2 − 5 − 5 ∣ ∣ ∣ 6 − 9 − 17
R 2 → R 2 / 3 R_2 \to R_2/3 R 2 → R 2 /3 :
( 1 − 1 2 ∣ 6 0 1 − 5 / 3 ∣ − 3 0 5 − 5 ∣ − 17 ) \begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 0 & 1 & -5/3 & | & -3 \\ 0 & 5 & -5 & | & -17 \end{pmatrix} 1 0 0 − 1 1 5 2 − 5/3 − 5 ∣ ∣ ∣ 6 − 3 − 17
R 3 → R 3 − 5 R 2 R_3 \to R_3 - 5R_2 R 3 → R 3 − 5 R 2 :
( 1 − 1 2 ∣ 6 0 1 − 5 / 3 ∣ − 3 0 0 20 / 3 ∣ − 2 ) \begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 0 & 1 & -5/3 & | & -3 \\ 0 & 0 & 20/3 & | & -2 \end{pmatrix} 1 0 0 − 1 1 0 2 − 5/3 20/3 ∣ ∣ ∣ 6 − 3 − 2
Back-substitution: R 3 ⟹ z = − 2 ⋅ 3 20 = − 3 10 R_3 \implies z = -2 \cdot \dfrac{3}{20} = -\dfrac{3}{10} R 3 ⟹ z = − 2 ⋅ 20 3 = − 10 3 .
R 2 R_2 R 2 : y − 5 3 ( − 3 10 ) = − 3 ⟹ y + 1 2 = − 3 ⟹ y = − 7 2 y - \dfrac{5}{3}\!\left(-\dfrac{3}{10}\right) = -3 \implies y + \dfrac{1}{2} = -3 \implies y = -\dfrac{7}{2} y − 3 5 ( − 10 3 ) = − 3 ⟹ y + 2 1 = − 3 ⟹ y = − 2 7 .
R 1 R_1 R 1 : x − ( − 7 2 ) + 2 ( − 3 10 ) = 6 ⟹ x + 7 2 − 3 5 = 6 ⟹ x = 6 − 29 10 = 31 10 x - \left(-\dfrac{7}{2}\right) + 2\!\left(-\dfrac{3}{10}\right) = 6 \implies x + \dfrac{7}{2} - \dfrac{3}{5} = 6 \implies x = 6 - \dfrac{29}{10} = \dfrac{31}{10} x − ( − 2 7 ) + 2 ( − 10 3 ) = 6 ⟹ x + 2 7 − 5 3 = 6 ⟹ x = 6 − 10 29 = 10 31 .
Solution: x = 31 10 x = \dfrac{31}{10} x = 10 31 , y = − 7 2 y = -\dfrac{7}{2} y = − 2 7 , z = − 3 10 z = -\dfrac{3}{10} z = − 10 3 .
Worked Example: Diagonalisation and Large Powers
For A = ( 3 1 0 2 ) A = \begin{pmatrix} 3 & 1 \\ 0 & 2 \end{pmatrix} A = ( 3 0 1 2 ) , find A 5 A^5 A 5 using diagonalisation.
Solution The matrix is already upper triangular, so eigenvalues are λ 1 = 3 \lambda_1 = 3 λ 1 = 3 and λ 2 = 2 \lambda_2 = 2 λ 2 = 2 .
For λ 1 = 3 \lambda_1 = 3 λ 1 = 3 : ( A − 3 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 0 1 0 − 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 1 = ( 1 0 ) (A - 3I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \begin{pmatrix} 0 & 1 \\ 0 & -1 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'}_1 = \begin{pmatrix} 1 \\ 0 \end{pmatrix} ( A − 3 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 0 0 1 − 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 1 = ( 1 0 ) .
For λ 2 = 2 \lambda_2 = 2 λ 2 = 2 : ( A − 2 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 1 1 0 0 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 2 = ( 1 − 1 ) (A - 2I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \begin{pmatrix} 1 & 1 \\ 0 & 0 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'}_2 = \begin{pmatrix} 1 \\ -1 \end{pmatrix} ( A − 2 I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 1 0 1 0 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 2 = ( 1 − 1 ) .
P = ( 1 1 0 − 1 ) P = \begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix} P = ( 1 0 1 − 1 ) , D = ( 3 0 0 2 ) D = \begin{pmatrix} 3 & 0 \\ 0 & 2 \end{pmatrix} D = ( 3 0 0 2 ) .
P − 1 = 1 − 1 ( − 1 − 1 0 1 ) = ( 1 1 0 − 1 ) P^{-1} = \dfrac{1}{-1}\begin{pmatrix} -1 & -1 \\ 0 & 1 \end{pmatrix} = \begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix} P − 1 = − 1 1 ( − 1 0 − 1 1 ) = ( 1 0 1 − 1 ) .
A 5 = P D 5 P − 1 = ( 1 1 0 − 1 ) ( 243 0 0 32 ) ( 1 1 0 − 1 ) A^5 = PD^5P^{-1} = \begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix}\begin{pmatrix} 243 & 0 \\ 0 & 32 \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix} A 5 = P D 5 P − 1 = ( 1 0 1 − 1 ) ( 243 0 0 32 ) ( 1 0 1 − 1 )
= ( 243 32 0 − 32 ) ( 1 1 0 − 1 ) = ( 243 211 0 32 ) = \begin{pmatrix} 243 & 32 \\ 0 & -32 \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix} = \begin{pmatrix} 243 & 211 \\ 0 & 32 \end{pmatrix} = ( 243 0 32 − 32 ) ( 1 0 1 − 1 ) = ( 243 0 211 32 ) .
Additional Common Pitfalls
Inverse of a product order. ( A B ) − 1 = B − 1 A − 1 (AB)^{-1} = B^{-1}A^{-1} ( A B ) − 1 = B − 1 A − 1 , not A − 1 B − 1 A^{-1}B^{-1} A − 1 B − 1 . The order reverses,
just as with the transpose.
Gaussian elimination arithmetic errors. A single arithmetic mistake in row operations propagates
through all subsequent steps. Always verify the solution by substituting back into the original
equations.
Misidentifying eigenvectors. For each eigenvalue, solve ( A − λ I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ (A - \lambda I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} ( A − λ I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ .
Substituting the eigenvalue into the original matrix A A A instead of A − λ I A - \lambda I A − λ I is a common error.
Transformation composition order. The matrix B A BA B A represents "apply A A A first, then B B B ." Writing
A B AB A B when B B B should be applied first is the single most common transformation error.
Determinant and area relationship for reflections. A reflection matrix has det ( T ) = − 1 \det(T) = -1 det ( T ) = − 1 . The
area scale factor is ∣ det ( T ) ∣ = 1 |\det(T)| = 1 ∣ det ( T ) ∣ = 1 , but the negative sign indicates orientation reversal. The area
is always positive; the sign encodes geometric information.
Exam-Style Problems
Problem 9 Given A = ( 1 2 3 4 ) A = \begin{pmatrix} 1 & 2 \\ 3 & 4 \end{pmatrix} A = ( 1 3 2 4 ) and B = ( 0 − 1 1 2 ) B = \begin{pmatrix} 0 & -1 \\ 1 & 2 \end{pmatrix} B = ( 0 1 − 1 2 ) ,
verify that ( A B ) − 1 = B − 1 A − 1 (AB)^{-1} = B^{-1}A^{-1} ( A B ) − 1 = B − 1 A − 1 .
Problem 10 Find the invariant points and invariant lines of the transformation represented by
T = ( 3 1 0 2 ) T = \begin{pmatrix} 3 & 1 \\ 0 & 2 \end{pmatrix} T = ( 3 0 1 2 ) .
Problem 11 The matrix A = ( 1 0 2 0 1 − 1 0 0 1 ) A = \begin{pmatrix} 1 & 0 & 2 \\ 0 & 1 & -1 \\ 0 & 0 & 1 \end{pmatrix} A = 1 0 0 0 1 0 2 − 1 1 represents a
transformation in 3D. Describe the geometric effect of A A A and find A n A^n A n for any positive integer n n n .
Problem 12 Use Cramer's rule to solve: 3 x + 2 y − z = 4 3x + 2y - z = 4 3 x + 2 y − z = 4 , x − y + 3 z = 5 x - y + 3z = 5 x − y + 3 z = 5 , 2 x + y + z = 3 2x + y + z = 3 2 x + y + z = 3 .
Problem 13 The transformation T T T consists of an enlargement by factor 2 2 2 followed by a rotation of 45 ∘ 45^\circ 4 5 ∘
anticlockwise. Find the matrix T T T and the area scale factor.
Problem 14 For A = ( 2 1 1 2 ) A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix} A = ( 2 1 1 2 ) , find the eigenvalues and eigenvectors. Verify
that the eigenvectors are orthogonal and diagonalise A A A .
Problem 15 Determine all values of k k k for which the system has (a) a unique solution, (b) no solution,
(c) infinitely many solutions:
x + 2 y + k z = 1 x + 2y + kz = 1 x + 2 y + k z = 1
2 x + k y + 2 z = 2 2x + ky + 2z = 2 2 x + k y + 2 z = 2
k x + 2 y + z = 3 kx + 2y + z = 3 k x + 2 y + z = 3
Answers to Additional Problems Problem 9: A B = ( 2 3 4 5 ) AB = \begin{pmatrix} 2 & 3 \\ 4 & 5 \end{pmatrix} A B = ( 2 4 3 5 ) , det ( A B ) = 10 − 12 = − 2 \det(AB) = 10 - 12 = -2 det ( A B ) = 10 − 12 = − 2 .
( A B ) − 1 = 1 − 2 ( 5 − 3 − 4 2 ) = ( − 5 / 2 3 / 2 2 − 1 ) (AB)^{-1} = \dfrac{1}{-2}\begin{pmatrix} 5 & -3 \\ -4 & 2 \end{pmatrix} = \begin{pmatrix} -5/2 & 3/2 \\ 2 & -1 \end{pmatrix} ( A B ) − 1 = − 2 1 ( 5 − 4 − 3 2 ) = ( − 5/2 2 3/2 − 1 ) .
A − 1 = 1 − 2 ( 4 − 2 − 3 1 ) = ( − 2 1 3 / 2 − 1 / 2 ) A^{-1} = \dfrac{1}{-2}\begin{pmatrix} 4 & -2 \\ -3 & 1 \end{pmatrix} = \begin{pmatrix} -2 & 1 \\ 3/2 & -1/2 \end{pmatrix} A − 1 = − 2 1 ( 4 − 3 − 2 1 ) = ( − 2 3/2 1 − 1/2 ) .
B − 1 = 1 1 ( 2 1 − 1 0 ) = ( 2 1 − 1 0 ) B^{-1} = \dfrac{1}{1}\begin{pmatrix} 2 & 1 \\ -1 & 0 \end{pmatrix} = \begin{pmatrix} 2 & 1 \\ -1 & 0 \end{pmatrix} B − 1 = 1 1 ( 2 − 1 1 0 ) = ( 2 − 1 1 0 ) .
B − 1 A − 1 = ( 2 1 − 1 0 ) ( − 2 1 3 / 2 − 1 / 2 ) = ( − 5 / 2 3 / 2 2 − 1 ) B^{-1}A^{-1} = \begin{pmatrix} 2 & 1 \\ -1 & 0 \end{pmatrix}\begin{pmatrix} -2 & 1 \\ 3/2 & -1/2 \end{pmatrix} = \begin{pmatrix} -5/2 & 3/2 \\ 2 & -1 \end{pmatrix} B − 1 A − 1 = ( 2 − 1 1 0 ) ( − 2 3/2 1 − 1/2 ) = ( − 5/2 2 3/2 − 1 ) .
Verified: ( A B ) − 1 = B − 1 A − 1 (AB)^{-1} = B^{-1}A^{-1} ( A B ) − 1 = B − 1 A − 1 .
Problem 10: Invariant points: ( T − I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 2 1 0 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ 2 x + y = 0 (T - I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \begin{pmatrix} 2 & 1 \\ 0 & 1 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies 2x + y = 0 ( T − I ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ( 2 0 1 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ 2 x + y = 0 and y = 0 y = 0 y = 0 , so x = 0 x = 0 x = 0 , y = 0 y = 0 y = 0 . Only the origin is invariant. For invariant lines y = m x + c y = mx + c y = m x + c : T ( x m x + c ) T\begin{pmatrix} x \\ mx + c \end{pmatrix} T ( x m x + c ) must lie on y = m x + c y = mx + c y = m x + c . This gives conditions leading to the invariant line y = 0 y = 0 y = 0 (the x x x -axis is mapped to itself).
Problem 11: A A A is an upper triangular matrix with 1 1 1 s on the diagonal. This represents a shear
transformation. For any n n n , A n = ( 1 0 2 n 0 1 − n 0 0 1 ) A^n = \begin{pmatrix} 1 & 0 & 2n \\ 0 & 1 & -n \\ 0 & 0 & 1 \end{pmatrix} A n = 1 0 0 0 1 0 2 n − n 1 , which
can be verified by induction.
Problem 12: det ( A ) = 3 ( − 1 − 3 ) − 2 ( 1 − 6 ) + ( − 1 ) ( 1 − ( − 2 k ) ) = − 12 + 10 − 1 + 2 k = 2 k − 3 \det(A) = 3(-1 - 3) - 2(1 - 6) + (-1)(1 - (-2k)) = -12 + 10 - 1 + 2k = 2k - 3 det ( A ) = 3 ( − 1 − 3 ) − 2 ( 1 − 6 ) + ( − 1 ) ( 1 − ( − 2 k )) = − 12 + 10 − 1 + 2 k = 2 k − 3 .
A x = ( 4 2 − 1 5 − 1 3 3 1 1 ) A_x = \begin{pmatrix} 4 & 2 & -1 \\ 5 & -1 & 3 \\ 3 & 1 & 1 \end{pmatrix} A x = 4 5 3 2 − 1 1 − 1 3 1 , det ( A x ) = 4 ( − 1 − 3 ) − 2 ( 5 − 9 ) + ( − 1 ) ( 5 + 3 ) = − 16 + 8 − 8 = − 16 \det(A_x) = 4(-1-3) - 2(5-9) + (-1)(5+3) = -16 + 8 - 8 = -16 det ( A x ) = 4 ( − 1 − 3 ) − 2 ( 5 − 9 ) + ( − 1 ) ( 5 + 3 ) = − 16 + 8 − 8 = − 16 .
A y = ( 3 4 − 1 1 5 3 2 3 1 ) A_y = \begin{pmatrix} 3 & 4 & -1 \\ 1 & 5 & 3 \\ 2 & 3 & 1 \end{pmatrix} A y = 3 1 2 4 5 3 − 1 3 1 , det ( A y ) = 3 ( 5 − 9 ) − 4 ( 1 − 6 ) + ( − 1 ) ( 3 − 10 ) = − 12 + 20 + 7 = 15 \det(A_y) = 3(5-9) - 4(1-6) + (-1)(3-10) = -12 + 20 + 7 = 15 det ( A y ) = 3 ( 5 − 9 ) − 4 ( 1 − 6 ) + ( − 1 ) ( 3 − 10 ) = − 12 + 20 + 7 = 15 .
A z = ( 3 2 4 1 − 1 5 2 1 3 ) A_z = \begin{pmatrix} 3 & 2 & 4 \\ 1 & -1 & 5 \\ 2 & 1 & 3 \end{pmatrix} A z = 3 1 2 2 − 1 1 4 5 3 , det ( A z ) = 3 ( − 3 − 5 ) − 2 ( 3 − 10 ) + 4 ( 1 + 2 ) = − 24 + 14 + 12 = 2 \det(A_z) = 3(-3-5) - 2(3-10) + 4(1+2) = -24 + 14 + 12 = 2 det ( A z ) = 3 ( − 3 − 5 ) − 2 ( 3 − 10 ) + 4 ( 1 + 2 ) = − 24 + 14 + 12 = 2 .
x = − 16 / ( 2 k − 3 ) x = -16/(2k-3) x = − 16/ ( 2 k − 3 ) , y = 15 / ( 2 k − 3 ) y = 15/(2k-3) y = 15/ ( 2 k − 3 ) , z = 2 / ( 2 k − 3 ) z = 2/(2k-3) z = 2/ ( 2 k − 3 ) .
Problem 13: Enlargement by 2 2 2 : E = ( 2 0 0 2 ) E = \begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix} E = ( 2 0 0 2 ) .
Rotation by 45 ∘ 45^\circ 4 5 ∘ : R = ( cos 45 ∘ − sin 45 ∘ sin 45 ∘ cos 45 ∘ ) = 1 2 ( 1 − 1 1 1 ) R = \begin{pmatrix} \cos 45^\circ & -\sin 45^\circ \\ \sin 45^\circ & \cos 45^\circ \end{pmatrix} = \dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 & -1 \\ 1 & 1 \end{pmatrix} R = ( cos 4 5 ∘ sin 4 5 ∘ − sin 4 5 ∘ cos 4 5 ∘ ) = 2 1 ( 1 1 − 1 1 ) .
T = R ⋅ E = 1 2 ( 1 − 1 1 1 ) ( 2 0 0 2 ) = 2 2 ( 1 − 1 1 1 ) = 2 ( 1 − 1 1 1 ) T = R \cdot E = \dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 & -1 \\ 1 & 1 \end{pmatrix}\begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix} = \dfrac{2}{\sqrt{2}}\begin{pmatrix} 1 & -1 \\ 1 & 1 \end{pmatrix} = \sqrt{2}\begin{pmatrix} 1 & -1 \\ 1 & 1 \end{pmatrix} T = R ⋅ E = 2 1 ( 1 1 − 1 1 ) ( 2 0 0 2 ) = 2 2 ( 1 1 − 1 1 ) = 2 ( 1 1 − 1 1 ) .
Area scale factor = ∣ det ( T ) ∣ = ∣ 2 − ( − 2 ) ∣ = 4 = |\det(T)| = |2 - (-2)| = 4 = ∣ det ( T ) ∣ = ∣2 − ( − 2 ) ∣ = 4 .
Problem 14: Characteristic equation: λ 2 − 4 λ + 3 = 0 ⟹ λ = 1 \lambda^2 - 4\lambda + 3 = 0 \implies \lambda = 1 λ 2 − 4 λ + 3 = 0 ⟹ λ = 1 or λ = 3 \lambda = 3 λ = 3 .
For λ = 1 \lambda = 1 λ = 1 : ( 1 1 1 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 1 = ( 1 − 1 ) \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'}_1 = \begin{pmatrix} 1 \\ -1 \end{pmatrix} ( 1 1 1 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 1 = ( 1 − 1 ) .
For λ = 3 \lambda = 3 λ = 3 : ( − 1 1 1 − 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 2 = ( 1 1 ) \begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'}_2 = \begin{pmatrix} 1 \\ 1 \end{pmatrix} ( − 1 1 1 − 1 ) ′ { ′ v ′ } ′ = ′ { ′ 0 ′ } ′ ⟹ ′ { ′ v ′ } ′ 2 = ( 1 1 ) .
Orthogonality check: ′ { ′ v ′ } ′ 1 ⋅ ′ { ′ v ′ } ′ 2 = 1 − 1 = 0 \mathbf{'\{'}v{'\}'}_1 \cdot \mathbf{'\{'}v{'\}'}_2 = 1 - 1 = 0 ′ { ′ v ′ } ′ 1 ⋅ ′ { ′ v ′ } ′ 2 = 1 − 1 = 0 . Confirmed.
A = P D P − 1 A = PDP^{-1} A = P D P − 1 where P = ( 1 1 − 1 1 ) P = \begin{pmatrix} 1 & 1 \\ -1 & 1 \end{pmatrix} P = ( 1 − 1 1 1 ) , D = ( 1 0 0 3 ) D = \begin{pmatrix} 1 & 0 \\ 0 & 3 \end{pmatrix} D = ( 1 0 0 3 ) .
Problem 15: det ( A ) = 1 ( k − 4 ) − 2 ( 2 − 2 k ) + k ( 4 − k 2 ) = k − 4 − 4 + 4 k + 4 k − k 3 = − k 3 + 9 k − 8 \det(A) = 1(k - 4) - 2(2 - 2k) + k(4 - k^2) = k - 4 - 4 + 4k + 4k - k^3 = -k^3 + 9k - 8 det ( A ) = 1 ( k − 4 ) − 2 ( 2 − 2 k ) + k ( 4 − k 2 ) = k − 4 − 4 + 4 k + 4 k − k 3 = − k 3 + 9 k − 8 .
Factor: − k 3 + 9 k − 8 = − ( k − 1 ) ( k 2 + k − 8 ) -k^3 + 9k - 8 = -(k-1)(k^2 + k - 8) − k 3 + 9 k − 8 = − ( k − 1 ) ( k 2 + k − 8 ) . The discriminant of k 2 + k − 8 k^2 + k - 8 k 2 + k − 8 is 1 + 32 = 33 1 + 32 = 33 1 + 32 = 33 .
(a) Unique solution when det ( A ) ≠ 0 \det(A) \ne 0 det ( A ) = 0 : k ≠ 1 k \ne 1 k = 1 and k ≠ − 1 ± 33 2 k \ne \dfrac{-1 \pm \sqrt{33}}{2} k = 2 − 1 ± 33 .
(b) No solution when det ( A ) = 0 \det(A) = 0 det ( A ) = 0 and the system is inconsistent.
(c) Infinitely many solutions when det ( A ) = 0 \det(A) = 0 det ( A ) = 0 and the system is consistent.
If You Get These Wrong, Revise:
Solving simultaneous equations → Review ./matrices (sections on Solving Systems and Cramer's Rule)
Trigonometry for rotation matrices → Review geometry and trigonometry topics
Algebraic manipulation and factorisation → Review algebra fundamentals
Vectors and linear independence → Review ./vectors for related concepts
Determinants and their geometric meaning → Review ./matrices (section on Determinants and Area)
For the A-Level Further Maths treatment of this topic, see Matrices .