The math optional, made finite. Daily Practice

Cayley-Hamilton theorem

At a Glance

Why This Chapter Matters

Cayley-Hamilton delivers 10–20 marks across three different question styles: compute a high matrix power, find the inverse via the characteristic equation, or evaluate a matrix polynomial by polynomial division. All three reduce to the same workflow — compute the characteristic polynomial, apply Cayley-Hamilton to reduce powers, then evaluate numerically. The 3×33\times 3 matrix with repeated eigenvalue (appearing in both 2019 and 2023) adds one step: the double-root derivative condition. A student who has practised this template can solve any variant in under 15 minutes.

Minimum Theory

Cayley-Hamilton theorem. Every n×nn\times n matrix AA satisfies its own characteristic equation: if p(λ)=det(λIA)p(\lambda)=\det(\lambda I-A) then p(A)=Op(A)=O. This means every matrix power AkA^k (for knk\ge n) can be expressed as a polynomial in AA of degree n1\le n-1.

Power reduction via Lagrange interpolation. To find ANA^N for large NN, write λN=q(λ)p(λ)+r(λ)\lambda^N=q(\lambda)p(\lambda)+r(\lambda) where degr<n\deg r<n. Since p(A)=Op(A)=O, AN=r(A)A^N=r(A). Determine the nn coefficients of rr by substituting each eigenvalue (with multiplicity): a simple root λ0\lambda_0 gives λ0N=r(λ0)\lambda_0^N=r(\lambda_0); a double root λ1\lambda_1 gives both λ1N=r(λ1)\lambda_1^N=r(\lambda_1) and Nλ1N1=r(λ1)N\lambda_1^{N-1}=r'(\lambda_1) (differentiate both sides of λN=q(λ)p(λ)+r(λ)\lambda^N=q(\lambda)p(\lambda)+r(\lambda) and use p(λ1)=p(λ1)=0p(\lambda_1)=p'(\lambda_1)=0).

Inverse via Cayley-Hamilton. For 2×22\times 2: from A2(trA)A+(detA)I=OA^2-(\text{tr}\,A)A+(\det A)I=O, rearrange to isolate II on one side and factor out AA: A1=((trA)IA)/detAA^{-1}=\bigl((\text{tr}\,A)I-A\bigr)/\det A.

Pipeline: from matrix A to A^n via the Cayley-Hamilton power-reduction technique

Question Archetypes

ArchetypeYou are seeing this when…
cayley-hamilton-power”Use the Cayley-Hamilton theorem to find ANA^N” for large NN
cayley-hamilton-application”Verify CH; find A1A^{-1}; evaluate a matrix polynomial p(A)p(A)

cayley-hamilton-power (2 question(s); 2019, 2023)

Recognition Cues

Solution Template

  1. State the theorem. Every n×nn\times n matrix satisfies its characteristic equation p(A)=Op(A)=O.
  2. Compute p(λ)=det(λIA)p(\lambda)=\det(\lambda I-A). Expand and factorise.
  3. Write λN=q(λ)p(λ)+r(λ)\lambda^N=q(\lambda)p(\lambda)+r(\lambda) with rr of degree 2\le 2 (i.e.\ r=aλ2+bλ+cr=a\lambda^2+b\lambda+c). So AN=aA2+bA+cIA^N=aA^2+bA+cI.
  4. Find a,b,ca,b,c from eigenvalue conditions. Simple root: substitute; double root: substitute and differentiate.
  5. Compute A2A^2 entry-by-entry. Assemble AN=aA2+bA+cIA^N=aA^2+bA+cI.

Worked Example

2019 Paper 1, 2019-P1-Q4a (15 marks)

State the Cayley-Hamilton theorem. Use it to find A100A^{100} for

A=(100101010).A=\begin{pmatrix}1&0&0\\1&0&1\\0&1&0\end{pmatrix}.

Step 1 — State theorem. Every square matrix satisfies its own characteristic equation: p(A)=Op(A)=O.

Step 2 — Characteristic polynomial.

p(λ)=det(λ1001λ101λ)=(λ1)(λ21)=(λ1)2(λ+1).p(\lambda)=\det\begin{pmatrix}\lambda-1&0&0\\-1&\lambda&-1\\0&-1&\lambda\end{pmatrix}=(\lambda-1)(\lambda^2-1)=(\lambda-1)^2(\lambda+1).

So p(λ)=λ3λ2λ+1p(\lambda)=\lambda^3-\lambda^2-\lambda+1 and by CH, A3=A2+AIA^3=A^2+A-I.

Step 3 — Write remainder. λ100=q(λ)(λ1)2(λ+1)+(aλ2+bλ+c)\lambda^{100}=q(\lambda)(\lambda-1)^2(\lambda+1)+(a\lambda^2+b\lambda+c), so A100=aA2+bA+cIA^{100}=aA^2+bA+cI.

Step 4 — Eigenvalue conditions.

λ=1\lambda=1 (value): 1=a+b+c.(1)1=a+b+c. \qquad(1)

λ=1\lambda=1 (derivative, double root): 100=2a+b.(2)100=2a+b. \qquad(2)

λ=1\lambda=-1 (simple root): (1)100=1=ab+c.(3)(-1)^{100}=1=a-b+c. \qquad(3)

From (1)(3)(1)-(3): 2b=02b=0, so b=0b=0. From (2)(2): a=50a=50. From (1)(1): c=150=49c=1-50=-49.

Step 5 — Compute A2A^2 and assemble.

A2=(100110101).A^2=\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix}.

A100=50A249I=50(100110101)49(100010001)=(10050105001).A^{100}=50A^2-49I=50\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix}-49\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}=\begin{pmatrix}1&0&0\\50&1&0\\50&0&1\end{pmatrix}.

  A100=(10050105001).  \boxed{\;A^{100}=\begin{pmatrix}1&0&0\\50&1&0\\50&0&1\end{pmatrix}.\;}


2023 Paper 1, 2023-P1-Q3a (20 marks)

(i) Verify CH for A=(100101010)A=\begin{pmatrix}1&0&0\\1&0&1\\0&1&0\end{pmatrix}. (ii) Show An=An2+A2IA^n=A^{n-2}+A^2-I for n3n\ge 3 and find A40A^{40}.

Part (i) — Verify CH.

Characteristic polynomial: p(λ)=(λ1)2(λ+1)=λ3λ2λ+1p(\lambda)=(\lambda-1)^2(\lambda+1)=\lambda^3-\lambda^2-\lambda+1.

Compute A2A^2 and A3A^3 entry-by-entry:

A2=(100110101),A3=(100201110).A^2=\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix},\qquad A^3=\begin{pmatrix}1&0&0\\2&0&1\\1&1&0\end{pmatrix}.

A3A2A+I=(111+100211+0010+1101+0110+0101+0010+1)=(000000000).  A^3-A^2-A+I=\begin{pmatrix}1-1-1+1&0&0\\2-1-1+0&0-1-0+1&1-0-1+0\\1-1-0+0&1-0-1+0&0-1-0+1\end{pmatrix}=\begin{pmatrix}0&0&0\\0&0&0\\0&0&0\end{pmatrix}.\;\checkmark

Part (ii) — Recurrence and A40A^{40}.

From CH: A3=A2+AIA^3=A^2+A-I, i.e.\ A3A=A2IA^3-A=A^2-I, i.e.\ A3A1=A2IA^3-A^1=A^2-I.

Claim: AnAn2=A2IA^n-A^{n-2}=A^2-I for all n3n\ge 3.

Proof by induction. Base n=3n=3: A3A=A2IA^3-A=A^2-I follows directly from CH. Inductive step: assume AkAk2=A2IA^k-A^{k-2}=A^2-I. Multiply both sides on the left by AA: Ak+1Ak1=A3A=A2IA^{k+1}-A^{k-1}=A^3-A=A^2-I (using the base case again). This proves Ak+1Ak1=A2IA^{k+1}-A^{k-1}=A^2-I.

Telescope to A40A^{40}: sum the identity AnAn2=A2IA^n-A^{n-2}=A^2-I for n=4,6,8,,40n=4,6,8,\ldots,40 (19 terms, all even steps):

A40A2=19(A2I)A40=20A219I.A^{40}-A^2=19(A^2-I)\qquad\Rightarrow\qquad A^{40}=20A^2-19I.

A40=20(100110101)19I=(10020102001).A^{40}=20\begin{pmatrix}1&0&0\\1&1&0\\1&0&1\end{pmatrix}-19I=\begin{pmatrix}1&0&0\\20&1&0\\20&0&1\end{pmatrix}.

  A40=(10020102001).  \boxed{\;A^{40}=\begin{pmatrix}1&0&0\\20&1&0\\20&0&1\end{pmatrix}.\;}

Common Traps


cayley-hamilton-application (1 question(s); 2014)

Recognition Cues

Solution Template

  1. Characteristic polynomial. p(λ)=λ2(trA)λ+detAp(\lambda)=\lambda^2-(\text{tr}\,A)\lambda+\det A.
  2. Verify CH. Compute A2A^2 directly and check A2(trA)A+(detA)I=OA^2-(\text{tr}\,A)A+(\det A)I=O.
  3. Inverse. Rearrange CH: A(A(trA)I)=(detA)IA(A-(\text{tr}\,A)I)=-(\det A)I, so A1=((trA)IA)/detAA^{-1}=((\text{tr}\,A)I-A)/\det A.
  4. Higher powers. Use A2=(trA)A(detA)IA^2=(\text{tr}\,A)A-(\det A)I as the recurrence; compute A3,A4,A5,A^3,A^4,A^5,\ldots iteratively as αkA+βkI\alpha_k A+\beta_k I.
  5. Sum the polynomial. Collect all AA-coefficients and all II-coefficients. Write the result as a single matrix.

Worked Example

2014 Paper 1, 2014-P1-Q2b-ii (10 marks)

Verify CH for A=(1423)A=\begin{pmatrix}1&4\\2&3\end{pmatrix}; find A1A^{-1}; evaluate A54A47A3+11A2A10IA^5-4A^4-7A^3+11A^2-A-10I.

Step 1 — Characteristic polynomial. trA=4\text{tr}\,A=4, detA=38=5\det A=3-8=-5. So p(λ)=λ24λ5p(\lambda)=\lambda^2-4\lambda-5.

Step 2 — Verify CH. A2=(916817)A^2=\begin{pmatrix}9&16\\8&17\end{pmatrix}. Check: A24A5I=(9451616088017125)=OA^2-4A-5I=\begin{pmatrix}9-4-5&16-16-0\\8-8-0&17-12-5\end{pmatrix}=O. \checkmark

Step 3 — Inverse. From A24A5I=OA^2-4A-5I=O: A(A4I)=5IA(A-4I)=5I, so

A1=A4I5=15(3421).A^{-1}=\frac{A-4I}{5}=\frac{1}{5}\begin{pmatrix}-3&4\\2&-1\end{pmatrix}.

Step 4 — Reduce higher powers. Recurrence: A2=4A+5IA^2=4A+5I.

A3=4A2+5A=4(4A+5I)+5A=21A+20I.A^3=4A^2+5A=4(4A+5I)+5A=21A+20I.

A4=AA3=21A2+20A=21(4A+5I)+20A=104A+105I.A^4=A\cdot A^3=21A^2+20A=21(4A+5I)+20A=104A+105I.

A5=AA4=104A2+105A=104(4A+5I)+105A=521A+520I.A^5=A\cdot A^4=104A^2+105A=104(4A+5I)+105A=521A+520I.

Step 5 — Evaluate polynomial. Collect:

TermAA-coeffII-coeff
A5A^5521521520520
4A4-4A^4416-416420-420
7A3-7A^3147-147140-140
11A211A^244445555
A-A1-100
10I-10I0010-10
Sum1155

Result: A+5I=(6428)A+5I=\begin{pmatrix}6&4\\2&8\end{pmatrix}.

  A54A47A3+11A2A10I=(6428).  \boxed{\;A^5-4A^4-7A^3+11A^2-A-10I=\begin{pmatrix}6&4\\2&8\end{pmatrix}.\;}

(Alternative: polynomial division of λ54λ47λ3+11λ2λ10\lambda^5-4\lambda^4-7\lambda^3+11\lambda^2-\lambda-10 by λ24λ5\lambda^2-4\lambda-5 gives remainder λ+5\lambda+5, so p(A)=A+5Ip(A)=A+5I.)

Common Traps


Marks-Aware Writing

10-mark question (2014): Verification (Steps 1–2), inverse (Step 3), and polynomial evaluation (Steps 4–5) each carry roughly 3+3+43+3+4 marks. An answer that verifies CH but miscomputes the inverse loses 3 marks; missing the polynomial evaluation entirely loses 4.

15-mark question (2019): Stating CH (2 marks), characteristic polynomial (3 marks), setting up and solving the a,b,ca,b,c system (5 marks including the double-root derivative condition), computing A2A^2 and the final matrix (5 marks).

20-mark question (2023): Verification (10 marks) and recurrence + A40A^{40} (10 marks). For verification: A2A^2 and A3A^3 computed correctly (5 marks), substitution and confirmation of zero (5 marks). For part (ii): proving the recurrence by induction (5 marks), telescoping correctly to get A40=20A219IA^{40}=20A^2-19I (3 marks), final matrix (2 marks).

Practice Set

(No additional practice items beyond the three worked examples.)

We've mapped all 13 years of this exam. Get new chapters, tools, and solutions as we release them — free.