Skip to main content

Matrices

Matrix Fundamentals

Definition and Notation

An m×nm \times n matrix AA is a rectangular array of real numbers with mm rows and nn columns:

A=(a11a12a1na21a22a2nam1am2amn)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}

The entry in row ii and column jj is aija_{ij}. When m=nm = n, AA is a square matrix of order nn.

The identity matrix InI_n is the n×nn \times n matrix with 11s on the diagonal and 00s elsewhere. For any compatible matrix AA: AI=IA=AAI = IA = A.


Matrix Operations

Addition and Scalar Multiplication

If A,B{M}m×n({R})A, B \in \mathcal{'\{'}M{'\}'}_{m \times n}(\mathbb{'\{'}R{'\}'}) and λ{R}\lambda \in \mathbb{'\{'}R{'\}'}:

(A+B)ij=aij+bij,(λA)ij=λaij(A + B)_{ij} = a_{ij} + b_{ij}, \qquad (\lambda A)_{ij} = \lambda \cdot a_{ij}

Matrix Multiplication

If AA is m×pm \times p and BB is p×np \times n, the product C=ABC = AB is m×nm \times n with:

cij=k=1paikbkjc_{ij} = \sum_{k=1}^{p} a_{ik} b_{kj}

Key properties:

  • Associative: (AB)C=A(BC)(AB)C = A(BC)
  • Distributive: A(B+C)=AB+ACA(B + C) = AB + AC
  • NOT commutative: ABBAAB \ne BA in general
  • Transpose of product: (AB)T=BTAT(AB)^T = B^T A^T

Transpose

The transpose ATA^T is obtained by reflecting AA across its main diagonal: (AT)ij=aji(A^T)_{ij} = a_{ji}.

Properties:

  • (AT)T=A(A^T)^T = A
  • (A+B)T=AT+BT(A + B)^T = A^T + B^T
  • (λA)T=λAT(\lambda A)^T = \lambda A^T
  • (AB)T=BTAT(AB)^T = B^T A^T

Determinants

2×22 \times 2 Determinant

det(abcd)=adbc\det\begin{pmatrix} a & b \\ c & d \end{pmatrix} = ad - bc

3×33 \times 3 Determinant

Expand along any row or column. Expanding along row 1:

det(A)=a11a22a23a32a33a12a21a23a31a33+a13a21a22a31a32\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}

The signs follow the checkerboard pattern +++ - + along row 1, +- + - along row 2, etc.

Properties of Determinants

  • det(AB)=det(A)det(B)\det(AB) = \det(A)\det(B)
  • det(AT)=det(A)\det(A^T) = \det(A)
  • det(λA)=λndet(A)\det(\lambda A) = \lambda^n \det(A) for an n×nn \times n matrix
  • Swapping two rows multiplies the determinant by 1-1
  • A matrix with a zero row (or column) has determinant 00
  • A triangular matrix has determinant equal to the product of its diagonal entries

Singular and Non-Singular Matrices

A square matrix AA is singular if det(A)=0\det(A) = 0 and non-singular (invertible) if det(A)0\det(A) \ne 0.


Inverse Matrices

Definition

The inverse of a square matrix AA is the matrix A1A^{-1} satisfying:

AA1=A1A=IAA^{-1} = A^{-1}A = I

AA is invertible if and only if det(A)0\det(A) \ne 0.

2×22 \times 2 Inverse

A1=1det(A)(dbca)A^{-1} = \frac{1}{\det(A)}\begin{pmatrix} d & -b \\ -c & a \end{pmatrix}

Adjugate Method (3×33 \times 3)

A1=1det(A)adj(A)A^{-1} = \frac{1}{\det(A)} \cdot \mathrm{adj}(A)

where the adjugate (classical adjoint) is the transpose of the cofactor matrix. The (i,j)(i, j)-cofactor is Cij=(1)i+jMijC_{ij} = (-1)^{i+j} M_{ij}, where MijM_{ij} is the (i,j)(i, j)-minor (determinant of the submatrix obtained by deleting row ii and column jj).

Properties of Inverses

  • (A1)1=A(A^{-1})^{-1} = A
  • (AB)1=B1A1(AB)^{-1} = B^{-1}A^{-1}
  • (AT)1=(A1)T(A^T)^{-1} = (A^{-1})^T
  • (λA)1=1λA1(\lambda A)^{-1} = \dfrac{1}{\lambda} A^{-1}, λ0\lambda \ne 0

Solving Systems of Linear Equations

Matrix Form

A system of nn linear equations in nn unknowns:

A{x}={b}A\mathbf{'\{'}x{'\}'} = \mathbf{'\{'}b{'\}'}

has the unique solution {x}=A1{b}\mathbf{'\{'}x{'\}'} = A^{-1}\mathbf{'\{'}b{'\}'} when det(A)0\det(A) \ne 0.

Cramer's Rule

For the system A{x}={b}A\mathbf{'\{'}x{'\}'} = \mathbf{'\{'}b{'\}'} where AA is n×nn \times n and det(A)0\det(A) \ne 0:

xi=det(Ai)det(A)x_i = \frac{\det(A_i)}{\det(A)}

where AiA_i is the matrix AA with column ii replaced by {b}\mathbf{'\{'}b{'\}'}.

Gaussian Elimination

For larger systems, use elementary row operations to reduce the augmented matrix [A{b}][A \mid \mathbf{'\{'}b{'\}'}] to row echelon form, then solve by back-substitution. The three elementary operations are:

  1. Swap two rows (RiRjR_i \leftrightarrow R_j)
  2. Multiply a row by a nonzero scalar (RikRiR_i \to kR_i)
  3. Add a multiple of one row to another (RiRi+kRjR_i \to R_i + kR_j)

Existence and Uniqueness

  • det(A)0\det(A) \ne 0: unique solution
  • det(A)=0\det(A) = 0 and system is consistent: infinitely many solutions
  • det(A)=0\det(A) = 0 and system is inconsistent: no solution

Geometric Transformations

2×22 \times 2 Transformation Matrices

Each transformation matrix TT maps a point (x,y)(x, y) to a new point by left-multiplication:

(xy)=T(xy)\begin{pmatrix} x' \\ y' \end{pmatrix} = T\begin{pmatrix} x \\ y \end{pmatrix}

TransformationMatrix
Reflection in xx-axis(1001)\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}
Reflection in yy-axis(1001)\begin{pmatrix} -1 & 0 \\ 0 & 1 \end{pmatrix}
Reflection in y=xy = x(0110)\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}
Rotation by θ\theta about origin(cosθsinθsinθcosθ)\begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}
Enlargement by factor kk(k00k)\begin{pmatrix} k & 0 \\ 0 & k \end{pmatrix}
Horizontal stretch by factor kk(k001)\begin{pmatrix} k & 0 \\ 0 & 1 \end{pmatrix}
Vertical stretch by factor kk(100k)\begin{pmatrix} 1 & 0 \\ 0 & k \end{pmatrix}
Shear parallel to xx-axis, factor kk(1k01)\begin{pmatrix} 1 & k \\ 0 & 1 \end{pmatrix}

Determinant and Area

For a 2×22 \times 2 transformation matrix TT:

  • det(T)|\det(T)| is the area scale factor
  • det(T)>0\det(T) \gt 0: orientation preserved
  • det(T)<0\det(T) \lt 0: orientation reversed

Composition of Transformations

Applying transformation AA followed by BB corresponds to the matrix product BABA (right to left order):

{v}=B(A{v})=(BA){v}\mathbf{'\{'}v{'\}'}' = B(A\mathbf{'\{'}v{'\}'}) = (BA)\mathbf{'\{'}v{'\}'}

Invariant Points and Lines

An invariant point satisfies T{v}={v}T\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}v{'\}'}, i.e. (TI){v}={0}(T - I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'}. An invariant line is mapped to itself (points may move along the line).


Eigenvalues and Eigenvectors (HL)

Definition

Let AA be an n×nn \times n matrix. A scalar λ\lambda is an eigenvalue of AA if there exists a nonzero vector {v}\mathbf{'\{'}v{'\}'} (an eigenvector) such that:

A{v}=λ{v}A\mathbf{'\{'}v{'\}'} = \lambda\mathbf{'\{'}v{'\}'}

Geometrically, AA stretches or compresses the eigenvector {v}\mathbf{'\{'}v{'\}'} by factor λ\lambda without changing its direction.

Characteristic Equation

Rearranging A{v}=λ{v}A\mathbf{'\{'}v{'\}'} = \lambda\mathbf{'\{'}v{'\}'} gives (AλI){v}={0}(A - \lambda I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'}. For a nonzero solution {v}\mathbf{'\{'}v{'\}'}, we require:

det(AλI)=0\det(A - \lambda I) = 0

This is the characteristic equation. Its roots are the eigenvalues of AA.

For a 2×22 \times 2 matrix A=(abcd)A = \begin{pmatrix} a & b \\ c & d \end{pmatrix}:

det(aλbcdλ)=(aλ)(dλ)bc=0\det\begin{pmatrix} a - \lambda & b \\ c & d - \lambda \end{pmatrix} = (a - \lambda)(d - \lambda) - bc = 0

λ2(a+d)λ+(adbc)=0\lambda^2 - (a + d)\lambda + (ad - bc) = 0

Note: a+d=tr(A)a + d = \mathrm{tr}(A) (the trace) and adbc=det(A)ad - bc = \det(A), so:

λ2tr(A)λ+det(A)=0\lambda^2 - \mathrm{tr}(A)\,\lambda + \det(A) = 0

Finding Eigenvectors

For each eigenvalue λi\lambda_i, substitute into (AλiI){v}={0}(A - \lambda_i I)\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} and solve for {v}\mathbf{'\{'}v{'\}'}. The solution space is the eigenspace corresponding to λi\lambda_i.

Properties

  • The trace equals the sum of eigenvalues: tr(A)=λi\mathrm{tr}(A) = \sum \lambda_i
  • The determinant equals the product of eigenvalues: det(A)=λi\det(A) = \prod \lambda_i
  • If AA is symmetric (A=ATA = A^T), all eigenvalues are real and eigenvectors corresponding to distinct eigenvalues are orthogonal.
  • If AA has nn linearly independent eigenvectors, it is diagonalisable: A=PDP1A = PDP^{-1} where DD is diagonal with eigenvalues on the diagonal.

Diagonalisation (HL Extension)

If AA has eigenvalues λ1,,λn\lambda_1, \ldots, \lambda_n and corresponding eigenvectors {v}1,,{v}n\mathbf{'\{'}v{'\}'}_1, \ldots, \mathbf{'\{'}v{'\}'}_n forming a linearly independent set, then:

A=PDP1A = PDP^{-1}

where P=({v}1{v}n)P = \begin{pmatrix} \mathbf{'\{'}v{'\}'}_1 & \cdots & \mathbf{'\{'}v{'\}'}_n \end{pmatrix} and D=diag(λ1,,λn)D = \mathrm{diag}(\lambda_1, \ldots, \lambda_n).

This is useful for computing large powers: An=PDnP1A^n = PD^nP^{-1}.


Common Pitfalls

  1. Matrix multiplication order. ABBAAB \ne BA in general. When composing transformations, the rightmost matrix is applied first.

  2. Transpose of a product. (AB)T=BTAT(AB)^T = B^T A^T, not ATBTA^T B^T. The order reverses.

  3. Determinant of a scalar multiple. det(kA)=kndet(A)\det(kA) = k^n\det(A) for an n×nn \times n matrix, not kdet(A)k\det(A).

  4. Forgetting to check invertibility. Only non-singular matrices have inverses. Always verify det(A)0\det(A) \ne 0 before attempting to compute A1A^{-1}.

  5. Eigenvector scaling. Eigenvectors are determined only up to a nonzero scalar multiple. Any nonzero multiple of an eigenvector is also an eigenvector.

  6. Confusing invariant points with fixed points of composition. A point fixed by AA is not necessarily fixed by BABA.

  7. Cofactor sign errors. The (i,j)(i, j)-cofactor sign is (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=(2314)A = \begin{pmatrix} 2 & 3 \\ 1 & 4 \end{pmatrix} and B=(1120)B = \begin{pmatrix} 1 & -1 \\ 2 & 0 \end{pmatrix}. Find ABAB and BABA. Is AB=BAAB = BA?

Problem 2

Find det(A)\det(A) and A1A^{-1} where A=(120311201)A = \begin{pmatrix} 1 & 2 & 0 \\ 3 & 1 & -1 \\ 2 & 0 & 1 \end{pmatrix}.

Problem 3

Solve the system using an inverse matrix: x+2y=5x + 2y = 5 3x+7y=183x + 7y = 18

Problem 4

Find the matrix representing a rotation of 9090^\circ anticlockwise about the origin, and verify that applying it twice gives a rotation of 180180^\circ.

Problem 5

Find the eigenvalues and eigenvectors of A=(4123)A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}.

Problem 6

The transformation represented by matrix T=(3102)T = \begin{pmatrix} 3 & 1 \\ 0 & 2 \end{pmatrix} is applied to the unit square. Find the area of the image.

Problem 7

Use Cramer's rule to solve: 2xy+z=32x - y + z = 3 x+3yz=1x + 3y - z = 1 3x+y+2z=73x + y + 2z = 7

Problem 8

For the matrix A=(5622)A = \begin{pmatrix} 5 & -6 \\ 2 & -2 \end{pmatrix}, find A10A^{10} by diagonalisation.

Answers to Selected Problems

Problem 1: AB=(8291)AB = \begin{pmatrix} 8 & -2 \\ 9 & -1 \end{pmatrix}, BA=(1146)BA = \begin{pmatrix} 1 & -1 \\ 4 & 6 \end{pmatrix}. ABBAAB \ne BA, confirming non-commutativity.

Problem 2: Expanding along row 1: det(A)=1(10)2(3(2))+0=110=9\det(A) = 1(1 - 0) - 2(3 - (-2)) + 0 = 1 - 10 = -9. Cofactors: C11=1,C12=5,C13=2,C21=2,C22=1,C23=4,C31=2,C32=1,C33=5C_{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. A1=19(122511245)=(1/92/92/95/91/91/92/94/95/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}.

Problem 3: A=(1237)A = \begin{pmatrix} 1 & 2 \\ 3 & 7 \end{pmatrix}, det(A)=76=1\det(A) = 7 - 6 = 1. A1=(7231)A^{-1} = \begin{pmatrix} 7 & -2 \\ -3 & 1 \end{pmatrix}. (xy)=(7231)(518)=(353615+18)=(13)\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}. So x=1,y=3x = -1, y = 3.

Problem 4: Rotation by 9090^\circ: R=(0110)R = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}. R2=(0110)(0110)=(1001)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}, which is a 180180^\circ rotation.

Problem 5: Characteristic equation: λ27λ+10=0\lambda^2 - 7\lambda + 10 = 0, so λ=2\lambda = 2 or λ=5\lambda = 5. For λ=2\lambda = 2: (A2I){v}={0}    (2121){v}={0}    {v}=(12)(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} (up to scalar). For λ=5\lambda = 5: (A5I){v}={0}    (1122){v}={0}    {v}=(11)(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}.

Problem 6: Area =det(T)=60=6= |\det(T)| = |6 - 0| = 6.


Worked Examples

Worked Example: Finding the Inverse of a 3×33 \times 3 Matrix

Find the inverse of A=(210131012)A = \begin{pmatrix} 2 & 1 & 0 \\ 1 & 3 & -1 \\ 0 & 1 & 2 \end{pmatrix}.

Solution

Step 1: Compute det(A)\det(A). Expanding along row 1:

det(A)=2311211102+0=2(7)1(2)=142=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

Step 2: Compute the cofactor matrix.

C11=+3112=7C_{11} = +\begin{vmatrix} 3 & -1 \\ 1 & 2 \end{vmatrix} = 7, C12=1102=2C_{12} = -\begin{vmatrix} 1 & -1 \\ 0 & 2 \end{vmatrix} = -2, C13=+1301=1C_{13} = +\begin{vmatrix} 1 & 3 \\ 0 & 1 \end{vmatrix} = 1

C21=1012=2C_{21} = -\begin{vmatrix} 1 & 0 \\ 1 & 2 \end{vmatrix} = -2, C22=+2002=4C_{22} = +\begin{vmatrix} 2 & 0 \\ 0 & 2 \end{vmatrix} = 4, C23=2101=2C_{23} = -\begin{vmatrix} 2 & 1 \\ 0 & 1 \end{vmatrix} = -2

C31=+1031=1C_{31} = +\begin{vmatrix} 1 & 0 \\ 3 & -1 \end{vmatrix} = -1, C32=2011=2C_{32} = -\begin{vmatrix} 2 & 0 \\ 1 & -1 \end{vmatrix} = 2, C33=+2113=5C_{33} = +\begin{vmatrix} 2 & 1 \\ 1 & 3 \end{vmatrix} = 5

Cofactor matrix: (721242125)\begin{pmatrix} 7 & -2 & 1 \\ -2 & 4 & -2 \\ -1 & 2 & 5 \end{pmatrix}.

Step 3: Transpose to get the adjugate.

adj(A)=(721242125)\mathrm{adj}(A) = \begin{pmatrix} 7 & -2 & -1 \\ -2 & 4 & 2 \\ 1 & -2 & 5 \end{pmatrix}

Step 4: Divide by the determinant.

A1=112(721242125)A^{-1} = \frac{1}{12}\begin{pmatrix} 7 & -2 & -1 \\ -2 & 4 & 2 \\ 1 & -2 & 5 \end{pmatrix}

Worked Example: Composition of Transformations

Find the matrix for the transformation that first reflects in the line y=xy = x, then rotates 9090^\circ anticlockwise about the origin. Verify the result by applying the combined transformation to the point (3,1)(3, 1).

Solution

Reflection in y=xy = x: Rf=(0110)R_f = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}.

Rotation 9090^\circ anticlockwise: Rθ=(0110)R_\theta = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}.

The combined transformation (reflect first, then rotate) is:

T=RθRf=(0110)(0110)=(1001)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}

This is a reflection in the yy-axis.

Verification: (3,1)(3, 1) reflected in y=xy = x gives (1,3)(1, 3). Rotating (1,3)(1, 3) by 9090^\circ anticlockwise: Rθ(13)=(31)R_\theta\begin{pmatrix} 1 \\ 3 \end{pmatrix} = \begin{pmatrix} -3 \\ 1 \end{pmatrix}.

Direct application: T(31)=(31)T\begin{pmatrix} 3 \\ 1 \end{pmatrix} = \begin{pmatrix} -3 \\ 1 \end{pmatrix}. Confirmed.

Worked Example: Solving a System Using Gaussian Elimination

Solve the system using Gaussian elimination:

2x+yz=32x + y - z = 3 xy+2z=6x - y + 2z = 6 3x+2y+z=13x + 2y + z = 1

Solution

Augmented matrix:

(211311263211)\begin{pmatrix} 2 & 1 & -1 & | & 3 \\ 1 & -1 & 2 & | & 6 \\ 3 & 2 & 1 & | & 1 \end{pmatrix}

R1R2R_1 \leftrightarrow R_2:

(112621133211)\begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 2 & 1 & -1 & | & 3 \\ 3 & 2 & 1 & | & 1 \end{pmatrix}

R2R22R1R_2 \to R_2 - 2R_1, R3R33R1R_3 \to R_3 - 3R_1:

(1126035905517)\begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 0 & 3 & -5 & | & -9 \\ 0 & 5 & -5 & | & -17 \end{pmatrix}

R2R2/3R_2 \to R_2/3:

(1126015/3305517)\begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 0 & 1 & -5/3 & | & -3 \\ 0 & 5 & -5 & | & -17 \end{pmatrix}

R3R35R2R_3 \to R_3 - 5R_2:

(1126015/330020/32)\begin{pmatrix} 1 & -1 & 2 & | & 6 \\ 0 & 1 & -5/3 & | & -3 \\ 0 & 0 & 20/3 & | & -2 \end{pmatrix}

Back-substitution: R3    z=2320=310R_3 \implies z = -2 \cdot \dfrac{3}{20} = -\dfrac{3}{10}.

R2R_2: y53 ⁣(310)=3    y+12=3    y=72y - \dfrac{5}{3}\!\left(-\dfrac{3}{10}\right) = -3 \implies y + \dfrac{1}{2} = -3 \implies y = -\dfrac{7}{2}.

R1R_1: x(72)+2 ⁣(310)=6    x+7235=6    x=62910=3110x - \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}.

Solution: x=3110x = \dfrac{31}{10}, y=72y = -\dfrac{7}{2}, z=310z = -\dfrac{3}{10}.

Worked Example: Diagonalisation and Large Powers

For A=(3102)A = \begin{pmatrix} 3 & 1 \\ 0 & 2 \end{pmatrix}, find A5A^5 using diagonalisation.

Solution

The matrix is already upper triangular, so eigenvalues are λ1=3\lambda_1 = 3 and λ2=2\lambda_2 = 2.

For λ1=3\lambda_1 = 3: (A3I){v}={0}    (0101){v}={0}    {v}1=(10)(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}.

For λ2=2\lambda_2 = 2: (A2I){v}={0}    (1100){v}={0}    {v}2=(11)(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}.

P=(1101)P = \begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix}, D=(3002)D = \begin{pmatrix} 3 & 0 \\ 0 & 2 \end{pmatrix}.

P1=11(1101)=(1101)P^{-1} = \dfrac{1}{-1}\begin{pmatrix} -1 & -1 \\ 0 & 1 \end{pmatrix} = \begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix}.

A5=PD5P1=(1101)(2430032)(1101)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}

=(24332032)(1101)=(243211032)= \begin{pmatrix} 243 & 32 \\ 0 & -32 \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix} = \begin{pmatrix} 243 & 211 \\ 0 & 32 \end{pmatrix}.


Additional Common Pitfalls

  • Inverse of a product order. (AB)1=B1A1(AB)^{-1} = B^{-1}A^{-1}, not A1B1A^{-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{'\}'}. Substituting the eigenvalue into the original matrix AA instead of AλIA - \lambda I is a common error.

  • Transformation composition order. The matrix BABA represents "apply AA first, then BB." Writing ABAB when BB 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. The area scale factor is 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=(1234)A = \begin{pmatrix} 1 & 2 \\ 3 & 4 \end{pmatrix} and B=(0112)B = \begin{pmatrix} 0 & -1 \\ 1 & 2 \end{pmatrix}, verify that (AB)1=B1A1(AB)^{-1} = B^{-1}A^{-1}.

Problem 10

Find the invariant points and invariant lines of the transformation represented by T=(3102)T = \begin{pmatrix} 3 & 1 \\ 0 & 2 \end{pmatrix}.

Problem 11

The matrix A=(102011001)A = \begin{pmatrix} 1 & 0 & 2 \\ 0 & 1 & -1 \\ 0 & 0 & 1 \end{pmatrix} represents a transformation in 3D. Describe the geometric effect of AA and find AnA^n for any positive integer nn.

Problem 12

Use Cramer's rule to solve: 3x+2yz=43x + 2y - z = 4, xy+3z=5x - y + 3z = 5, 2x+y+z=32x + y + z = 3.

Problem 13

The transformation TT consists of an enlargement by factor 22 followed by a rotation of 4545^\circ anticlockwise. Find the matrix TT and the area scale factor.

Problem 14

For A=(2112)A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}, find the eigenvalues and eigenvectors. Verify that the eigenvectors are orthogonal and diagonalise AA.

Problem 15

Determine all values of kk for which the system has (a) a unique solution, (b) no solution, (c) infinitely many solutions:

x+2y+kz=1x + 2y + kz = 1 2x+ky+2z=22x + ky + 2z = 2 kx+2y+z=3kx + 2y + z = 3

Answers to Additional Problems

Problem 9: AB=(2345)AB = \begin{pmatrix} 2 & 3 \\ 4 & 5 \end{pmatrix}, det(AB)=1012=2\det(AB) = 10 - 12 = -2. (AB)1=12(5342)=(5/23/221)(AB)^{-1} = \dfrac{1}{-2}\begin{pmatrix} 5 & -3 \\ -4 & 2 \end{pmatrix} = \begin{pmatrix} -5/2 & 3/2 \\ 2 & -1 \end{pmatrix}. A1=12(4231)=(213/21/2)A^{-1} = \dfrac{1}{-2}\begin{pmatrix} 4 & -2 \\ -3 & 1 \end{pmatrix} = \begin{pmatrix} -2 & 1 \\ 3/2 & -1/2 \end{pmatrix}. B1=11(2110)=(2110)B^{-1} = \dfrac{1}{1}\begin{pmatrix} 2 & 1 \\ -1 & 0 \end{pmatrix} = \begin{pmatrix} 2 & 1 \\ -1 & 0 \end{pmatrix}. B1A1=(2110)(213/21/2)=(5/23/221)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}. Verified: (AB)1=B1A1(AB)^{-1} = B^{-1}A^{-1}.

Problem 10: Invariant points: (TI){v}={0}    (2101){v}={0}    2x+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 and y=0y = 0, so x=0x = 0, y=0y = 0. Only the origin is invariant. For invariant lines y=mx+cy = mx + c: T(xmx+c)T\begin{pmatrix} x \\ mx + c \end{pmatrix} must lie on y=mx+cy = mx + c. This gives conditions leading to the invariant line y=0y = 0 (the xx-axis is mapped to itself).

Problem 11: AA is an upper triangular matrix with 11s on the diagonal. This represents a shear transformation. For any nn, An=(102n01n001)A^n = \begin{pmatrix} 1 & 0 & 2n \\ 0 & 1 & -n \\ 0 & 0 & 1 \end{pmatrix}, which can be verified by induction.

Problem 12: det(A)=3(13)2(16)+(1)(1(2k))=12+101+2k=2k3\det(A) = 3(-1 - 3) - 2(1 - 6) + (-1)(1 - (-2k)) = -12 + 10 - 1 + 2k = 2k - 3. Ax=(421513311)A_x = \begin{pmatrix} 4 & 2 & -1 \\ 5 & -1 & 3 \\ 3 & 1 & 1 \end{pmatrix}, det(Ax)=4(13)2(59)+(1)(5+3)=16+88=16\det(A_x) = 4(-1-3) - 2(5-9) + (-1)(5+3) = -16 + 8 - 8 = -16. Ay=(341153231)A_y = \begin{pmatrix} 3 & 4 & -1 \\ 1 & 5 & 3 \\ 2 & 3 & 1 \end{pmatrix}, det(Ay)=3(59)4(16)+(1)(310)=12+20+7=15\det(A_y) = 3(5-9) - 4(1-6) + (-1)(3-10) = -12 + 20 + 7 = 15. Az=(324115213)A_z = \begin{pmatrix} 3 & 2 & 4 \\ 1 & -1 & 5 \\ 2 & 1 & 3 \end{pmatrix}, det(Az)=3(35)2(310)+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/(2k3)x = -16/(2k-3), y=15/(2k3)y = 15/(2k-3), z=2/(2k3)z = 2/(2k-3).

Problem 13: Enlargement by 22: E=(2002)E = \begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}. Rotation by 4545^\circ: R=(cos45sin45sin45cos45)=12(1111)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}. T=RE=12(1111)(2002)=22(1111)=2(1111)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}. Area scale factor =det(T)=2(2)=4= |\det(T)| = |2 - (-2)| = 4.

Problem 14: Characteristic equation: λ24λ+3=0    λ=1\lambda^2 - 4\lambda + 3 = 0 \implies \lambda = 1 or λ=3\lambda = 3. For λ=1\lambda = 1: (1111){v}={0}    {v}1=(11)\begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'}_1 = \begin{pmatrix} 1 \\ -1 \end{pmatrix}. For λ=3\lambda = 3: (1111){v}={0}    {v}2=(11)\begin{pmatrix} -1 & 1 \\ 1 & -1 \end{pmatrix}\mathbf{'\{'}v{'\}'} = \mathbf{'\{'}0{'\}'} \implies \mathbf{'\{'}v{'\}'}_2 = \begin{pmatrix} 1 \\ 1 \end{pmatrix}. Orthogonality check: {v}1{v}2=11=0\mathbf{'\{'}v{'\}'}_1 \cdot \mathbf{'\{'}v{'\}'}_2 = 1 - 1 = 0. Confirmed. A=PDP1A = PDP^{-1} where P=(1111)P = \begin{pmatrix} 1 & 1 \\ -1 & 1 \end{pmatrix}, D=(1003)D = \begin{pmatrix} 1 & 0 \\ 0 & 3 \end{pmatrix}.

Problem 15: det(A)=1(k4)2(22k)+k(4k2)=k44+4k+4kk3=k3+9k8\det(A) = 1(k - 4) - 2(2 - 2k) + k(4 - k^2) = k - 4 - 4 + 4k + 4k - k^3 = -k^3 + 9k - 8. Factor: k3+9k8=(k1)(k2+k8)-k^3 + 9k - 8 = -(k-1)(k^2 + k - 8). The discriminant of k2+k8k^2 + k - 8 is 1+32=331 + 32 = 33. (a) Unique solution when det(A)0\det(A) \ne 0: k1k \ne 1 and k1±332k \ne \dfrac{-1 \pm \sqrt{33}}{2}. (b) No solution when det(A)=0\det(A) = 0 and the system is inconsistent. (c) Infinitely many solutions when 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.