MATH 133 Linear Algebra and Geometry • McGill University, Montreal

Corrected exercises: determinants, inverses and Cramer's rule (MATH 133)

This is the corrected exercise set for the chapter of MATH 133, Linear Algebra and Geometry at McGill University, that turns the determinant from a computation into a tool (Nicholson, section 3.2). Computing a determinant was the previous chapter. Here you use det⁡(AB)=det⁡A det⁡B\det(AB) = \det A \, \det B, det⁡(AT)=det⁡A\det(A^{T}) = \det A, det⁡(A−1)=1/det⁡A\det(A^{-1}) = 1/\det A and det⁡(kA)=kndet⁡A\det(kA) = k^{n}\det A to evaluate a determinant without its entries, decide for which values of a parameter a matrix is invertible, build the inverse from the adjugate, and extract one unknown of a system by Cramer's rule. No calculator, as on the exam: every number reduces by hand.

The thread running through the whole set: the determinant is a MULTIPLICATIVE volume factor. Products multiply, inverses divide, the transpose changes nothing, and a scalar kk in front of an n×nn \times n matrix scales each of its nn edges, hence knk^{n}. The same number sits in the denominator of every entry of A−1A^{-1} and of every unknown given by Cramer's rule, which is why det⁡A=0\det A = 0 is the single place where all of them stop at once, without that meaning the system has no solution.

The traps named explicitly in the solutions: writing kdet⁡Ak\det A for det⁡(kA)\det(kA), confusing (3A)−1(3A)^{-1} with 3A−13A^{-1}, keeping one branch when det⁡(A)2\det(A)^{2} is known, dividing by det⁡A\det A before knowing it is not zero, forgetting the transpose in the adjugate and so swapping the indices of a single entry of A−1A^{-1}, proving a statement about (−1)n(-1)^{n} without ever mentioning nn, and reading det⁡A=0\det A = 0 as no solution when Cramer's rule has simply stopped applying.

10 corrected exercises • 100 points • 150 minutes

Revision sheet for this chapter → Every MATH 133 chapter →

Course recap

  • • det⁡(AB)=det⁡A det⁡B\det(AB) = \det A \, \det B in any order, det⁡(AT)=det⁡A\det(A^{T}) = \det A, det⁡(Am)=(det⁡A)m\det(A^{m}) = (\det A)^{m}.
  • • AA is invertible if and only if det⁡A≠0\det A \neq 0, and then det⁡(A−1)=1det⁡A\det(A^{-1}) = \dfrac{1}{\det A}.
  • • det⁡(kA)=kndet⁡A\det(kA) = k^{n}\det A for an n×nn \times n matrix; ONE row times kk gives only kdet⁡Ak\det A.
  • • Cofactor Cij=(−1)i+jMijC_{ij} = (-1)^{i+j}M_{ij}; adj(A)=[Cij]T\text{adj}(A) = [C_{ij}]^{T} and A adj(A)=det⁡(A) IA\,\text{adj}(A) = \det(A)\,I.
  • • A−1=1det⁡A adj(A)A^{-1} = \dfrac{1}{\det A}\,\text{adj}(A), so (A−1)ij=Cjidet⁡A(A^{-1})_{ij} = \dfrac{C_{ji}}{\det A}, indices swapped.
  • • Cramer: if det⁡A≠0\det A \neq 0, the solution of Ax=bA\mathbf{x} = \mathbf{b} has xj=det⁡Aj(b)det⁡Ax_j = \dfrac{\det A_j(\mathbf{b})}{\det A}, column jj replaced by b\mathbf{b}.
  • • det⁡(adj A)=(det⁡A)n−1\det(\text{adj}\,A) = (\det A)^{n-1}; an integer matrix has an integer inverse if and only if det⁡A=±1\det A = \pm 1.

Part A: the basics (/50)

Exercise 1: One number per matrix: a composite determinant without a single entry

The determinant turns every square matrix into one number, and it turns products into products: det⁡(AB)=det⁡A det⁡B\det(AB) = \det A \, \det B. Three companions follow from it: det⁡(AT)=det⁡A\det(A^{T}) = \det A, det⁡(A−1)=1/det⁡A\det(A^{-1}) = 1/\det A, and det⁡(kA)=kndet⁡A\det(kA) = k^{n}\det A for an n×nn \times n matrix. With those four rules, a determinant built from five factors is five small numbers multiplied together, and the entries of the matrices are never needed.

The figure shows the product rule in the plane: PP doubles every area, QQ triples every area, so doing PP and then QQ multiplies areas by 66, whatever the shapes look like.

In parts a) to e), AA and BB are 3×33 \times 3 matrices with det⁡A=2\det A = 2 and det⁡B=−3\det B = -3.

Pdet 2Qdet 3area 1area 2area 6 = 2 × 3
  • a) Compute det⁡(AB)\det(AB) and det⁡(BA)\det(BA). Does it follow that AB=BAAB = BA?
  • b) Compute det⁡(ATB−1)\det(A^{T}B^{-1}).
  • c) Compute det⁡(2A−1BTA3)\det(2A^{-1}B^{T}A^{3}).
  • d) Compute det⁡(B−1AB)\det(B^{-1}AB) and explain why the answer does not depend on BB.
  • e) Compute det⁡((3A)−1)\det\big((3A)^{-1}\big) and det⁡(3A−1)\det(3A^{-1}), and say why they differ.
Show the solution

Answers

  • a) det⁡(AB)=det⁡(BA)=−6\det(AB) = \det(BA) = -6; no, equal determinants do not make equal matrices.
  • b) −23-\frac{2}{3}
  • c) −96-96
  • d) 22, for every invertible BB
  • e) 154\frac{1}{54} and 272\frac{27}{2}

a) det⁡(AB)=det⁡A det⁡B=2×(−3)=−6\det(AB) = \det A \, \det B = 2 \times (-3) = -6, and det⁡(BA)=det⁡B det⁡A=−6\det(BA) = \det B \, \det A = -6 as well, because the two determinants are real NUMBERS and numbers commute even when matrices do not. Nothing follows about ABAB and BABA themselves: two different matrices can share a determinant exactly as two different parallelograms can share an area. Writing AB=BAAB = BA here costs the part. The determinant forgets almost everything about a matrix and keeps one number, which is precisely why it is so easy to compute with and so dangerous to over-read.

b) Take the factors one at a time and write one line for each. det⁡(AT)=det⁡A=2\det(A^{T}) = \det A = 2: transposing swaps rows and columns, and expanding ATA^{T} along a row is expanding AA along a column. det⁡(B−1)=1det⁡B=−13\det(B^{-1}) = \dfrac{1}{\det B} = -\dfrac{1}{3}: from BB−1=IBB^{-1} = I you get det⁡B det⁡(B−1)=det⁡I=1\det B \, \det(B^{-1}) = \det I = 1. So det⁡(ATB−1)=2×(−13)=−23\det(A^{T}B^{-1}) = 2 \times \left(-\dfrac{1}{3}\right) = -\dfrac{2}{3}. The inverse DIVIDES by det⁡B\det B, it does not change its sign: writing det⁡(B−1)=3\det(B^{-1}) = 3 is the second most common slip of the chapter.

c) One line per factor before multiplying, it is what earns the method marks. The scalar 22 multiplies a 3×33 \times 3 matrix, so it contributes 23=82^{3} = 8, not 22. Then det⁡(A−1)=12\det(A^{-1}) = \dfrac{1}{2}, det⁡(BT)=−3\det(B^{T}) = -3 and det⁡(A3)=(det⁡A)3=8\det(A^{3}) = (\det A)^{3} = 8. Multiply: 8×12×(−3)×8=−968 \times \dfrac{1}{2} \times (-3) \times 8 = -96. Two traps sit in this single line. Writing 22 instead of 232^{3} gives −24-24. Reading A3A^{3} as 3A3A gives det⁡(3A)=27×2=54\det(3A) = 27 \times 2 = 54 instead of 88. Notice also what you did NOT need: the order of the factors. The matrix A−1BTA3A^{-1}B^{T}A^{3} cannot be simplified to BTA2B^{T}A^{2}, since the factors do not commute, but its determinant is a product of four numbers and those do.

d) det⁡(B−1AB)=det⁡(B−1) det⁡A det⁡B=1det⁡B×2×det⁡B=2\det(B^{-1}AB) = \det(B^{-1}) \, \det A \, \det B = \dfrac{1}{\det B} \times 2 \times \det B = 2. The matrices B−1B^{-1} and BB cannot be cancelled inside the product, because B−1ABB^{-1}AB is in general a different matrix from AA; their DETERMINANTS are numbers and do cancel. So every matrix of the form B−1ABB^{-1}AB has the determinant of AA, whichever invertible BB is chosen. Geometrically, B−1ABB^{-1}AB is the transformation AA described in another coordinate system, and changing the coordinates does not change the factor by which volumes are scaled.

e) det⁡(3A)=33×2=54\det(3A) = 3^{3} \times 2 = 54, so det⁡((3A)−1)=154\det\big((3A)^{-1}\big) = \dfrac{1}{54}. On the other side, det⁡(3A−1)=33×det⁡(A−1)=27×12=272\det(3A^{-1}) = 3^{3} \times \det(A^{-1}) = 27 \times \dfrac{1}{2} = \dfrac{27}{2}. The two expressions look alike and differ by a factor of 729=272729 = 27^{2}. In (3A)−1(3A)^{-1} the scalar is INSIDE the inverse and gets inverted with the matrix, since (3A)−1=13A−1(3A)^{-1} = \dfrac{1}{3}A^{-1}; in 3A−13A^{-1} it is outside and is not. Reading the parentheses is part of the question, and a marker reads them too.

Exercise 2: Scaling a matrix scales every edge: the power n

Multiplying ONE row of a matrix by kk multiplies the determinant by kk. Multiplying the whole MATRIX by kk multiplies each of its nn rows by kk, so the determinant is multiplied by knk^{n}. The rule det⁡(kA)=kndet⁡A\det(kA) = k^{n}\det A is where most marks are lost in this chapter, because kdet⁡Ak\det A looks so natural.

The figure shows why the power is the dimension: doubling every edge of a square gives four copies of it, doubling every edge of a cube gives eight.

n = 2: area × 4n = 3: volume × 8every edge × 2every edge × 2
  • a) AA is 3×33 \times 3 with det⁡A=5\det A = 5. Compute det⁡(2A)\det(2A), det⁡(−A)\det(-A) and det⁡(−2AT)\det(-2A^{T}).
  • b) Redo a) with AA of size 4×44 \times 4, still with det⁡A=5\det A = 5. Which answers change, and why?
  • c) BB is 3×33 \times 3 and det⁡(2B−1)=4\det(2B^{-1}) = 4. Find det⁡B\det B, then det⁡((2B)−1)\det\big((2B)^{-1}\big).
  • d) DD is 3×33 \times 3 with det⁡D=−4\det D = -4, and EE is obtained from DD by multiplying its first row by 22. Compare det⁡E\det E with det⁡(2D)\det(2D).
  • e) Prove det⁡(kA)=kndet⁡A\det(kA) = k^{n}\det A in one line, using the product rule.
Show the solution

Answers

  • a) 4040, −5-5 and −40-40
  • b) 8080, 55 and 8080: every answer changes, because nn is an input of each rule
  • c) det⁡B=2\det B = 2 and det⁡((2B)−1)=116\det\big((2B)^{-1}\big) = \frac{1}{16}
  • d) det⁡E=−8\det E = -8 and det⁡(2D)=−32\det(2D) = -32
  • e) kA=(kI)AkA = (kI)A and det⁡(kI)=kn\det(kI) = k^{n}

a) For a 3×33 \times 3 matrix, det⁡(2A)=23det⁡A=8×5=40\det(2A) = 2^{3}\det A = 8 \times 5 = 40. For −A-A the scalar is −1-1, so det⁡(−A)=(−1)3×5=−5\det(-A) = (-1)^{3} \times 5 = -5. For −2AT-2A^{T}, transposing changes nothing and the scalar −2-2 contributes (−2)3=−8(-2)^{3} = -8, so det⁡(−2AT)=−8×5=−40\det(-2A^{T}) = -8 \times 5 = -40. The answer det⁡(2A)=10\det(2A) = 10, that is 2det⁡A2\det A, is the single most common wrong line on this topic: it treats the scalar as if it multiplied one row only.

b) With n=4n = 4: det⁡(2A)=24×5=80\det(2A) = 2^{4} \times 5 = 80, det⁡(−A)=(−1)4×5=5\det(-A) = (-1)^{4} \times 5 = 5 and det⁡(−2AT)=(−2)4×5=80\det(-2A^{T}) = (-2)^{4} \times 5 = 80. The determinant of AA is still 55 and yet every answer changed, because the size of the matrix is a hidden input of every scalar rule. The sign of det⁡(−A)\det(-A) is the sharpest symptom: negating a matrix reverses the sign of its determinant when nn is odd and leaves it alone when nn is even. So never apply a scalar rule before writing nn on the line.

c) det⁡(2B−1)=23det⁡(B−1)=8det⁡B\det(2B^{-1}) = 2^{3}\det(B^{-1}) = \dfrac{8}{\det B}. Setting this equal to 44 gives det⁡B=2\det B = 2. Then det⁡(2B)=8×2=16\det(2B) = 8 \times 2 = 16 and its inverse has determinant 116\dfrac{1}{16}. A student who forgets the cube writes 2det⁡B=4\dfrac{2}{\det B} = 4 and finds det⁡B=12\det B = \dfrac{1}{2}, wrong by a factor of 44, and the second answer inherits the error. Substituting back is a five-second check: 8/2=48/2 = 4.

d) Multiplying ONE row by 22 is an elementary row operation and it multiplies the determinant by 22: det⁡E=2×(−4)=−8\det E = 2 \times (-4) = -8. Multiplying the whole matrix multiplies all three rows: det⁡(2D)=23×(−4)=−32\det(2D) = 2^{3} \times (-4) = -32. Both operations are described as multiplying by 22, and they differ by a factor of 22=42^{2} = 4, one factor of 22 for each row the row operation did not touch. In the language of the figure: stretching a cube in one direction doubles its volume, stretching it in all three directions multiplies it by 88.

e) kA=(kI)AkA = (kI)A, where kIkI is the diagonal matrix with kk in every diagonal position. The product rule gives det⁡(kA)=det⁡(kI) det⁡A\det(kA) = \det(kI)\,\det A, and the determinant of a diagonal matrix is the product of its diagonal entries, so det⁡(kI)=k⋅k⋯k=kn\det(kI) = k \cdot k \cdots k = k^{n}. Hence det⁡(kA)=kndet⁡A\det(kA) = k^{n}\det A. The proof is worth keeping because it also says what nn is: the number of diagonal entries of kIkI, that is the SIZE of the matrix, never the number of times kk is written in the expression.

Exercise 3: A parameter that switches the determinant off

A square matrix is invertible if and only if its determinant is not zero. When the matrix depends on a parameter, its determinant is a polynomial in that parameter, and the matrix fails to be invertible exactly at the roots. The work is to compute that polynomial with as little effort as possible, keep it FACTORED, and read it.

For every real kk, let Ak=(11k−1k1112)A_k = \begin{pmatrix} 1 & 1 & k \\ -1 & k & 1 \\ 1 & 1 & 2 \end{pmatrix}.

  • a) Compute det⁡(Ak)\det(A_k) as a factored polynomial in kk, naming each operation you use.
  • b) For which values of kk is AkA_k invertible?
  • c) Find every kk for which det⁡(AkTAk)=16\det(A_k^{T}A_k) = 16.
  • d) For k=0k = 0, compute det⁡(A0−1)\det(A_0^{-1}) and det⁡(2A0−1A0T)\det(2A_0^{-1}A_0^{T}).
  • e) Show that for k=2k = 2 the product A2BA_2B is singular for every 3×33 \times 3 matrix BB.
Show the solution

Answers

  • a) det⁡(Ak)=(2−k)(k+1)\det(A_k) = (2 - k)(k + 1)
  • b) Invertible for every real kk except k=2k = 2 and k=−1k = -1
  • c) k=3k = 3 or k=−2k = -2; the branch det⁡(Ak)=4\det(A_k) = 4 has no real solution
  • d) 12\frac{1}{2} and 88
  • e) det⁡(A2B)=0×det⁡B=0\det(A_2B) = 0 \times \det B = 0

a) Rows 11 and 33 share their first two entries, so R3→R3−R1R_3 \to R_3 - R_1 creates two zeros, and adding a multiple of one row to another leaves the determinant unchanged. The new third row is (0,0,2−k)(0, 0, 2 - k). Expand along that row, where only the last entry survives, with sign (−1)3+3=+1(-1)^{3+3} = +1: det⁡(Ak)=(2−k)det⁡(11−1k)=(2−k)(k+1)\det(A_k) = (2 - k)\det\begin{pmatrix} 1 & 1 \\ -1 & k \end{pmatrix} = (2 - k)(k + 1). Keeping the factored form is the whole point: expanding to −k2+k+2-k^{2} + k + 2 and factoring again wastes a line and invites a sign error. Check with k=0k = 0 on the ORIGINAL matrix, expanded along its first row: 1(0⋅2−1⋅1)−1((−1)⋅2−1⋅1)+0=−1+3=21(0 \cdot 2 - 1 \cdot 1) - 1\big((-1) \cdot 2 - 1 \cdot 1\big) + 0 = -1 + 3 = 2, and (2−0)(0+1)=2(2 - 0)(0 + 1) = 2.

b) AkA_k is invertible if and only if det⁡(Ak)≠0\det(A_k) \neq 0, that is k≠2k \neq 2 and k≠−1k \neq -1. At those two values the columns are dependent and the matrix is singular; at k=2k = 2 you can even see it, since rows 11 and 33 both become (1,1,2)(1, 1, 2). State the answer as a condition on kk, every real kk except 22 and −1-1. Listing the two bad values without saying what they mean loses the conclusion mark, because the question asked where AkA_k IS invertible.

c) det⁡(AkTAk)=det⁡(AkT) det⁡(Ak)=det⁡(Ak)2\det(A_k^{T}A_k) = \det(A_k^{T})\,\det(A_k) = \det(A_k)^{2}, so the condition is det⁡(Ak)2=16\det(A_k)^{2} = 16, that is det⁡(Ak)=4\det(A_k) = 4 OR det⁡(Ak)=−4\det(A_k) = -4. First branch: (2−k)(k+1)=4(2 - k)(k + 1) = 4 becomes k2−k+2=0k^{2} - k + 2 = 0, whose discriminant 1−8=−71 - 8 = -7 is negative, so no real kk. Second branch: (2−k)(k+1)=−4(2 - k)(k + 1) = -4 becomes k2−k−6=0k^{2} - k - 6 = 0, that is (k−3)(k+2)=0(k - 3)(k + 2) = 0, so k=3k = 3 or k=−2k = -2. The classic loss is to keep only det⁡(Ak)=4\det(A_k) = 4, find no solution and conclude that no kk works, the exact opposite of the truth. The figure of the solution shows both branches: the parabola never reaches height 44 and crosses height −4-4 twice.

d) For k=0k = 0, det⁡(A0)=2≠0\det(A_0) = 2 \neq 0, so A0−1A_0^{-1} exists, and det⁡(A0−1)=12\det(A_0^{-1}) = \dfrac{1}{2}. Then det⁡(2A0−1A0T)=23×12×2=8\det(2A_0^{-1}A_0^{T}) = 2^{3} \times \dfrac{1}{2} \times 2 = 8: the inverse and the transpose cancel in the determinant, and only the 232^{3} of the scalar is left. The sentence that says k=0k = 0 is not a forbidden value comes BEFORE any use of A0−1A_0^{-1}; without it, the inverse is used before it is known to exist.

e) For k=2k = 2, det⁡(A2)=0\det(A_2) = 0. For any 3×33 \times 3 matrix BB, det⁡(A2B)=det⁡(A2) det⁡B=0×det⁡B=0\det(A_2B) = \det(A_2)\,\det B = 0 \times \det B = 0, so A2BA_2B is singular however well behaved BB is, and the same argument covers BA2BA_2 and A25A_2^{5}. In words: once one factor has flattened space onto a plane, nothing applied before or after it can restore the lost dimension, and the determinant, being a product, records this as a factor of zero.

-3-2-11234-6-4-224height 4: never reachedk = -2 and k = 3singularsingularkdet A_k

Exercise 4: The adjugate: every entry of the inverse is a cofactor over det A

Each entry of A−1A^{-1} is a 2×22 \times 2 determinant divided by det⁡A\det A. Precisely, with Cij=(−1)i+jMijC_{ij} = (-1)^{i+j}M_{ij} the cofactor of the entry in row ii and column jj, the adjugate is the TRANSPOSE of the matrix of cofactors, adj(A)=[Cij]T\text{adj}(A) = [C_{ij}]^{T}, and A adj(A)=det⁡(A) IA\,\text{adj}(A) = \det(A)\,I. When det⁡A≠0\det A \neq 0, this gives A−1=1det⁡A adj(A)A^{-1} = \dfrac{1}{\det A}\,\text{adj}(A).

Throughout, A=(11010−112−1)A = \begin{pmatrix} 1 & 1 & 0 \\ 1 & 0 & -1 \\ 1 & 2 & -1 \end{pmatrix}.

  • a) Compute det⁡A\det A by cofactor expansion along a row you choose, and say why you chose it.
  • b) Compute the nine cofactors and write the matrix of cofactors.
  • c) Write adj(A)\text{adj}(A) and check the first row of A adj(A)=det⁡(A) IA\,\text{adj}(A) = \det(A)\,I.
  • d) Deduce A−1A^{-1}.
  • e) A student forgets the transpose and offers the matrix of cofactors divided by det⁡A\det A as the inverse. Show with a single entry of a product that it cannot be A−1A^{-1}.
Show the solution

Answers

  • a) det⁡A=2\det A = 2, along row 11 because of its zero
  • b) (2021−1−1−11−1)\begin{pmatrix} 2 & 0 & 2 \\ 1 & -1 & -1 \\ -1 & 1 & -1 \end{pmatrix}
  • c) adj(A)=(21−10−112−1−1)\text{adj}(A) = \begin{pmatrix} 2 & 1 & -1 \\ 0 & -1 & 1 \\ 2 & -1 & -1 \end{pmatrix}; first row of the product (2,0,0)(2, 0, 0)
  • d) A−1=(112−120−12121−12−12)A^{-1} = \begin{pmatrix} 1 & \frac{1}{2} & -\frac{1}{2} \\ 0 & -\frac{1}{2} & \frac{1}{2} \\ 1 & -\frac{1}{2} & -\frac{1}{2} \end{pmatrix}
  • e) Entry (1,1)(1,1) of the product is 32\frac{3}{2}, not 11

a) Row 11, (1,1,0)(1, 1, 0), has a zero, so its third term vanishes: det⁡A=1⋅C11+1⋅C12+0\det A = 1 \cdot C_{11} + 1 \cdot C_{12} + 0. With C11=+det⁡(0−12−1)=0+2=2C_{11} = +\det\begin{pmatrix} 0 & -1 \\ 2 & -1 \end{pmatrix} = 0 + 2 = 2 and C12=−det⁡(1−11−1)=−(−1+1)=0C_{12} = -\det\begin{pmatrix} 1 & -1 \\ 1 & -1 \end{pmatrix} = -(-1 + 1) = 0, we get det⁡A=2\det A = 2. Choosing the row with a zero saves one 2×22 \times 2 determinant; here it also produces two cofactors needed in b), so nothing is computed twice. The determinant is not zero, so AA is invertible and the adjugate formula applies.

b) The signs follow the checkerboard (+−+−+−+−+)\begin{pmatrix} + & - & + \\ - & + & - \\ + & - & + \end{pmatrix}, and each minor deletes the row and column of its entry. Row 11: C11=2C_{11} = 2, C12=0C_{12} = 0, C13=+det⁡(1012)=2C_{13} = +\det\begin{pmatrix} 1 & 0 \\ 1 & 2 \end{pmatrix} = 2. Row 22: C21=−det⁡(102−1)=1C_{21} = -\det\begin{pmatrix} 1 & 0 \\ 2 & -1 \end{pmatrix} = 1, C22=+det⁡(101−1)=−1C_{22} = +\det\begin{pmatrix} 1 & 0 \\ 1 & -1 \end{pmatrix} = -1, C23=−det⁡(1112)=−1C_{23} = -\det\begin{pmatrix} 1 & 1 \\ 1 & 2 \end{pmatrix} = -1. Row 33: C31=+det⁡(100−1)=−1C_{31} = +\det\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} = -1, C32=−det⁡(101−1)=1C_{32} = -\det\begin{pmatrix} 1 & 0 \\ 1 & -1 \end{pmatrix} = 1, C33=+det⁡(1110)=−1C_{33} = +\det\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} = -1. The matrix of cofactors is (2021−1−1−11−1)\begin{pmatrix} 2 & 0 & 2 \\ 1 & -1 & -1 \\ -1 & 1 & -1 \end{pmatrix}. The minus signs of the checkerboard are where the marks go: C21C_{21} is MINUS a determinant that equals −1-1.

c) Transpose: row ii of adj(A)\text{adj}(A) is column ii of the cofactor matrix, adj(A)=(21−10−112−1−1)\text{adj}(A) = \begin{pmatrix} 2 & 1 & -1 \\ 0 & -1 & 1 \\ 2 & -1 & -1 \end{pmatrix}. Row 11 of AA, (1,1,0)(1, 1, 0), against the three columns of adj(A)\text{adj}(A): 2+0+0=22 + 0 + 0 = 2, 1−1+0=01 - 1 + 0 = 0, −1+1+0=0-1 + 1 + 0 = 0. That is (2,0,0)(2, 0, 0), the first row of 2I2I. The zeros are not luck: a row of AA against the cofactors of ANOTHER row is the expansion of a matrix with two equal rows, hence zero. This check costs three dot products and catches almost every sign error.

d) A−1=12 adj(A)=(112−120−12121−12−12)A^{-1} = \dfrac{1}{2}\,\text{adj}(A) = \begin{pmatrix} 1 & \frac{1}{2} & -\frac{1}{2} \\ 0 & -\frac{1}{2} & \frac{1}{2} \\ 1 & -\frac{1}{2} & -\frac{1}{2} \end{pmatrix}. Every entry has denominator 2=det⁡A2 = \det A or a divisor of it, which is the adjugate formula made visible: the only division in the whole computation is by the determinant. On a copy, leave the answer as 12\frac{1}{2} times an integer matrix if the question allows it; it is shorter and easier to check.

e) With the matrix of cofactors CC instead of its transpose, entry (1,1)(1, 1) of A⋅12CA \cdot \dfrac{1}{2}C is row 11 of AA against column 11 of CC: 12(1⋅2+1⋅1+0⋅(−1))=32\dfrac{1}{2}(1 \cdot 2 + 1 \cdot 1 + 0 \cdot (-1)) = \dfrac{3}{2}, not 11. One wrong entry is enough to reject a proposed inverse. The transpose is not a convention: the product A adj(A)A\,\text{adj}(A) pairs row ii of AA with the cofactors of row ii, which is the expansion of det⁡A\det A along row ii, and that pairing only happens when those cofactors sit in COLUMN ii of the adjugate.

Exercise 5: One entry of the inverse, one unknown of the system

Two questions come back on every MATH 133 paper, and neither needs the full inverse. One entry of A−1A^{-1}: (A−1)ij=Cjidet⁡A(A^{-1})_{ij} = \dfrac{C_{ji}}{\det A}, with the indices SWAPPED because of the transpose in the adjugate. One unknown of Ax=bA\mathbf{x} = \mathbf{b}: Cramer's rule, xj=det⁡Aj(b)det⁡Ax_j = \dfrac{\det A_j(\mathbf{b})}{\det A}, where Aj(b)A_j(\mathbf{b}) is AA with column jj replaced by b\mathbf{b}, valid when det⁡A≠0\det A \neq 0.

In parts a) to c), M=(21−1132101)M = \begin{pmatrix} 2 & 1 & -1 \\ 1 & 3 & 2 \\ 1 & 0 & 1 \end{pmatrix} and the system is Mx=bM\mathbf{x} = \mathbf{b} with x=(x,y,z)\mathbf{x} = (x, y, z) and b=(−3,1,4)\mathbf{b} = (-3, 1, 4).

-1123456-112345a1a2b
  • a) Compute det⁡M\det M.
  • b) Compute only the entry in row 22, column 33 of M−1M^{-1}.
  • c) Find yy alone by Cramer's rule, then zz alone, and finish with a check by substitution.
  • d) In the plane, a1=(2,1)\mathbf{a}_1 = (2, 1), a2=(1,2)\mathbf{a}_2 = (1, 2) and b=(5,4)\mathbf{b} = (5, 4), as drawn. Find xx and yy in xa1+ya2=bx\mathbf{a}_1 + y\mathbf{a}_2 = \mathbf{b} by Cramer's rule, and explain with the figure why xx is a ratio of two determinants.
Show the solution

Answers

  • a) det⁡M=10\det M = 10
  • b) (M−1)23=C3210=−12(M^{-1})_{23} = \frac{C_{32}}{10} = -\frac{1}{2}
  • c) y=−2y = -2, z=3z = 3, then x=1x = 1
  • d) x=2x = 2, y=1y = 1: b=2a1+a2\mathbf{b} = 2\mathbf{a}_1 + \mathbf{a}_2

a) Expand along row 33, (1,0,1)(1, 0, 1), which has a zero. The signs in row 33 are +,−,++, -, +: det⁡M=+1⋅det⁡(1−132)+0+1⋅det⁡(2113)=(2+3)+(6−1)=10\det M = +1 \cdot \det\begin{pmatrix} 1 & -1 \\ 3 & 2 \end{pmatrix} + 0 + 1 \cdot \det\begin{pmatrix} 2 & 1 \\ 1 & 3 \end{pmatrix} = (2 + 3) + (6 - 1) = 10. Since det⁡M≠0\det M \neq 0, MM is invertible and Cramer's rule applies: say it in one sentence, it is the hypothesis of everything that follows.

b) The entry in row 22, column 33 of M−1M^{-1} is the cofactor of position (3,2)(3, 2) of MM, divided by det⁡M\det M. Delete row 33 and column 22: the minor is det⁡(2−112)=4+1=5\det\begin{pmatrix} 2 & -1 \\ 1 & 2 \end{pmatrix} = 4 + 1 = 5, and the sign is (−1)3+2=−1(-1)^{3+2} = -1, so C32=−5C_{32} = -5 and (M−1)23=−510=−12(M^{-1})_{23} = -\dfrac{5}{10} = -\dfrac{1}{2}. The trap is to take C23C_{23} instead: deleting row 22 and column 33 gives det⁡(2110)=−1\det\begin{pmatrix} 2 & 1 \\ 1 & 0 \end{pmatrix} = -1, so C23=1C_{23} = 1 and the wrong answer 110\dfrac{1}{10}. The two differ here, and in general they differ whenever MM is not symmetric. The swap comes from the transpose in adj(M)\text{adj}(M), nothing else.

c) Replace column 22 by b\mathbf{b}: det⁡(2−3−1112141)\det\begin{pmatrix} 2 & -3 & -1 \\ 1 & 1 & 2 \\ 1 & 4 & 1 \end{pmatrix}, expanded along row 11, is 2(1−8)+3(1−2)−1(4−1)=−14−3−3=−202(1 - 8) + 3(1 - 2) - 1(4 - 1) = -14 - 3 - 3 = -20, so y=−2010=−2y = \dfrac{-20}{10} = -2. Replace column 33 by b\mathbf{b}: det⁡(21−3131104)\det\begin{pmatrix} 2 & 1 & -3 \\ 1 & 3 & 1 \\ 1 & 0 & 4 \end{pmatrix}, expanded along row 33, is 1⋅(1+9)+4⋅(6−1)=301 \cdot (1 + 9) + 4 \cdot (6 - 1) = 30, so z=3z = 3. The third equation, x+z=4x + z = 4, then gives x=1x = 1, and all three equations check: 2−2−3=−32 - 2 - 3 = -3, 1−6+6=11 - 6 + 6 = 1, 1+3=41 + 3 = 4. The column that is replaced is the column of the unknown you want, and b\mathbf{b} goes in as a COLUMN, never as a row.

d) det⁡(2112)=3\det\begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix} = 3, with a1\mathbf{a}_1 and a2\mathbf{a}_2 as columns. Cramer: x=13det⁡(5142)=10−43=2x = \dfrac{1}{3}\det\begin{pmatrix} 5 & 1 \\ 4 & 2 \end{pmatrix} = \dfrac{10 - 4}{3} = 2 and y=13det⁡(2514)=8−53=1y = \dfrac{1}{3}\det\begin{pmatrix} 2 & 5 \\ 1 & 4 \end{pmatrix} = \dfrac{8 - 5}{3} = 1. Check: 2(2,1)+(1,2)=(5,4)2(2, 1) + (1, 2) = (5, 4). The figure of the solution says why. The numerator is the signed area of the parallelogram built on b\mathbf{b} and a2\mathbf{a}_2. Sliding b\mathbf{b} parallel to a2\mathbf{a}_2, to b−ya2=xa1\mathbf{b} - y\mathbf{a}_2 = x\mathbf{a}_1, changes neither the base nor the height, so that area is the area built on xa1x\mathbf{a}_1 and a2\mathbf{a}_2, which is xx times the shaded one. Here 6=2×36 = 2 \times 3. Cramer's rule is that sentence, written with determinants.

-1123456-112345a1a2b2a136

Part B: problems and reasoning (/50)

Exercise 6: Short proofs: take the determinant of both sides

A determinant identity turns a matrix equation into an equation between NUMBERS, and number equations are easy. Every proof below has the same shape: take the determinant of both sides, use the product, transpose and scalar rules, then solve for det⁡A\det A, keeping every root.

The figure shows the two answers of part b): a quarter turn, which keeps the counterclockwise order of the columns, and a reflection, which reverses it. Both keep areas.

Qe1Qe2Re1Re2rotation, det Q = 1turn keptreflection, det R = -1turn reversed
  • a) AA is n×nn \times n with A2=AA^{2} = A. Prove that det⁡A=0\det A = 0 or det⁡A=1\det A = 1, and give a 2×22 \times 2 example, other than the zero matrix, with A2=AA^{2} = A and det⁡A=0\det A = 0.
  • b) QQ is n×nn \times n with QTQ=IQ^{T}Q = I. Prove that det⁡Q=±1\det Q = \pm 1, and check both signs on the two matrices of the figure.
  • c) SS is 3×33 \times 3 with ST=−SS^{T} = -S. Prove that det⁡S=0\det S = 0, then show with a 2×22 \times 2 example that the conclusion fails for even sizes.
  • d) AA is 3×33 \times 3 with A3=2IA^{3} = 2I. Find det⁡A\det A exactly.
  • e) AA and BB are n×nn \times n with AB=−BAAB = -BA and nn odd. Prove that AA or BB is singular.
Show the solution

Answers

  • a) det⁡(A)2=det⁡A\det(A)^{2} = \det A, so det⁡A∈{0,1}\det A \in \{0, 1\}; example (1100)\begin{pmatrix} 1 & 1 \\ 0 & 0 \end{pmatrix}
  • b) det⁡(Q)2=1\det(Q)^{2} = 1; the rotation has det⁡=1\det = 1, the reflection det⁡=−1\det = -1
  • c) det⁡S=(−1)3det⁡S\det S = (-1)^{3}\det S, so det⁡S=0\det S = 0; (01−10)\begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix} has determinant 11
  • d) det⁡A=2\det A = 2
  • e) det⁡Adet⁡B=−det⁡Adet⁡B\det A \det B = -\det A \det B, so det⁡Adet⁡B=0\det A \det B = 0

a) Take determinants: det⁡(A2)=det⁡(A)2\det(A^{2}) = \det(A)^{2} by the product rule, so det⁡(A)2=det⁡A\det(A)^{2} = \det A, that is det⁡A (det⁡A−1)=0\det A\,(\det A - 1) = 0, and det⁡A=0\det A = 0 or det⁡A=1\det A = 1. Nothing more follows: the number equation has two roots and both are reached. For the example, A=(1100)A = \begin{pmatrix} 1 & 1 \\ 0 & 0 \end{pmatrix} gives A2=(1⋅1+1⋅01⋅1+1⋅000)=AA^{2} = \begin{pmatrix} 1 \cdot 1 + 1 \cdot 0 & 1 \cdot 1 + 1 \cdot 0 \\ 0 & 0 \end{pmatrix} = A and det⁡A=0\det A = 0, while II reaches det⁡=1\det = 1. The trap is to divide by det⁡A\det A and announce det⁡A=1\det A = 1 alone: dividing by a number is allowed only once you know it is not zero, and zero is the TYPICAL case here, since a matrix with A2=AA^{2} = A projects space onto a subspace.

b) det⁡(QTQ)=det⁡(QT) det⁡Q=det⁡(Q)2\det(Q^{T}Q) = \det(Q^{T})\,\det Q = \det(Q)^{2} and det⁡I=1\det I = 1, so det⁡(Q)2=1\det(Q)^{2} = 1 and det⁡Q=1\det Q = 1 or −1-1. On the figure, the quarter turn is (0−110)\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}, with determinant 0⋅0−(−1)⋅1=10 \cdot 0 - (-1) \cdot 1 = 1, and the reflection in the horizontal axis is (100−1)\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}, with determinant −1-1. Both satisfy QTQ=IQ^{T}Q = I, since their columns are perpendicular unit vectors. The absolute value 11 says that areas are kept; the SIGN says whether the turn from the first column to the second is still counterclockwise. A matrix with QTQ=IQ^{T}Q = I never stretches: it can only turn, or flip.

c) Take determinants of ST=−SS^{T} = -S. The left side is det⁡(ST)=det⁡S\det(S^{T}) = \det S. The right side is the scalar −1-1 times a 3×33 \times 3 matrix, so det⁡(−S)=(−1)3det⁡S=−det⁡S\det(-S) = (-1)^{3}\det S = -\det S. Hence det⁡S=−det⁡S\det S = -\det S, 2det⁡S=02\det S = 0 and det⁡S=0\det S = 0. The size was used exactly once, in (−1)3=−1(-1)^{3} = -1, and that is the step that fails for even nn: with n=2n = 2, (−1)2=1(-1)^{2} = 1 and the equation reduces to det⁡S=det⁡S\det S = \det S, which says nothing. Indeed S=(01−10)S = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix} satisfies ST=−SS^{T} = -S and has determinant 0−(−1)=10 - (-1) = 1. A proof that never mentions nn is wrong, because the statement is false for half of the sizes.

d) det⁡(A3)=det⁡(A)3\det(A^{3}) = \det(A)^{3}, and det⁡(2I)=23=8\det(2I) = 2^{3} = 8 for a 3×33 \times 3 matrix. So det⁡(A)3=8\det(A)^{3} = 8 and, since a real number has exactly one real cube root, det⁡A=2\det A = 2. The trap is det⁡(2I)=2\det(2I) = 2, which leads to det⁡A=23\det A = \sqrt[3]{2}, a wrong answer that LOOKS sophisticated. Such matrices exist: A=23 IA = \sqrt[3]{2}\,I satisfies A3=2IA^{3} = 2I, its entries are cube roots, and its determinant is (23)3=2(\sqrt[3]{2})^{3} = 2, as found.

e) Take determinants of AB=−BAAB = -BA: det⁡A det⁡B=(−1)ndet⁡B det⁡A\det A\,\det B = (-1)^{n}\det B\,\det A. With nn odd, (−1)n=−1(-1)^{n} = -1, so det⁡A det⁡B=−det⁡A det⁡B\det A\,\det B = -\det A\,\det B, hence det⁡A det⁡B=0\det A\,\det B = 0 and one of the two factors is zero: AA or BB is singular. It is the shape of c) again, a sign coming from (−1)n(-1)^{n} that forces a number to equal its own opposite. For even nn the conclusion fails: A=(100−1)A = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix} and B=(0110)B = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} give AB=(01−10)=−BAAB = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix} = -BA, with both determinants equal to −1-1.

Exercise 7: The adjugate identity, used as a tool

Everything about the adjugate follows from one identity, valid for every square matrix, invertible or not: A adj(A)=adj(A) A=det⁡(A) IA\,\text{adj}(A) = \text{adj}(A)\,A = \det(A)\,I. Taking determinants of it, or multiplying it by something well chosen, answers questions that would otherwise cost nine cofactors.

  • a) Prove that det⁡(adj A)=det⁡(A)n−1\det(\text{adj}\,A) = \det(A)^{n-1} when AA is an invertible n×nn \times n matrix.
  • b) Prove that if AA is singular and n≥2n \ge 2, then adj(A)\text{adj}(A) is singular too, so the formula of a) still holds.
  • c) AA is 3×33 \times 3 with det⁡A=3\det A = 3. Compute det⁡(adj A)\det(\text{adj}\,A), det⁡(2 adj A)\det(2\,\text{adj}\,A) and det⁡(adj(2A))\det\big(\text{adj}(2A)\big).
  • d) A 3×33 \times 3 matrix BB has adj(B)=(11−21101−12)\text{adj}(B) = \begin{pmatrix} 1 & 1 & -2 \\ 1 & 1 & 0 \\ 1 & -1 & 2 \end{pmatrix}. Show that det⁡B=±2\det B = \pm 2, and explain why the adjugate alone cannot decide the sign.
  • e) You are also told that the top-left entry of BB is 11. Find det⁡B\det B, B−1B^{-1} and BB.
Show the solution

Answers

  • a) det⁡A det⁡(adj A)=det⁡(A)n\det A \, \det(\text{adj}\,A) = \det(A)^{n}, then divide by det⁡A≠0\det A \neq 0
  • b) Otherwise A=0A = 0, whose adjugate is 00: a contradiction
  • c) 99, 7272 and 576576
  • d) det⁡(B)2=4\det(B)^{2} = 4; BB and −B-B have the same adjugate
  • e) det⁡B=2\det B = 2, B−1=12 adj(B)B^{-1} = \frac{1}{2}\,\text{adj}(B), B=(101−12−1−110)B = \begin{pmatrix} 1 & 0 & 1 \\ -1 & 2 & -1 \\ -1 & 1 & 0 \end{pmatrix}

a) Take determinants of A adj(A)=det⁡(A) IA\,\text{adj}(A) = \det(A)\,I. The left side is det⁡A det⁡(adj A)\det A\,\det(\text{adj}\,A). The right side is the NUMBER det⁡A\det A times the n×nn \times n identity, so its determinant is det⁡(A)n\det(A)^{n}, by the scalar rule. Hence det⁡A det⁡(adj A)=det⁡(A)n\det A\,\det(\text{adj}\,A) = \det(A)^{n}, and dividing by det⁡A\det A, which is allowed because AA is invertible, gives det⁡(adj A)=det⁡(A)n−1\det(\text{adj}\,A) = \det(A)^{n-1}. The usual slip is to read det⁡(det⁡(A) I)\det(\det(A)\,I) as det⁡A\det A, forgetting that the scalar sits in front of all nn rows, and to conclude det⁡(adj A)=1\det(\text{adj}\,A) = 1.

b) Suppose AA is singular and adj(A)\text{adj}(A) is invertible. The identity gives A adj(A)=0⋅I=0A\,\text{adj}(A) = 0 \cdot I = 0; multiply on the right by adj(A)−1\text{adj}(A)^{-1} to get A=0A = 0. But every cofactor of the zero matrix is a determinant of a zero block, so adj(0)=0\text{adj}(0) = 0, which is not invertible: contradiction. Hence adj(A)\text{adj}(A) is singular, det⁡(adj A)=0=det⁡(A)n−1\det(\text{adj}\,A) = 0 = \det(A)^{n-1}, and the formula holds for every square matrix of size at least 22. The hypothesis n≥2n \ge 2 is used in 0n−1=00^{n-1} = 0 and in the cofactors being determinants of non-empty blocks.

c) det⁡(adj A)=32=9\det(\text{adj}\,A) = 3^{2} = 9. Then 2 adj A2\,\text{adj}\,A is a scalar times a 3×33 \times 3 matrix: det⁡(2 adj A)=23×9=72\det(2\,\text{adj}\,A) = 2^{3} \times 9 = 72. For adj(2A)\text{adj}(2A), apply a) to the matrix 2A2A, whose determinant is 23×3=242^{3} \times 3 = 24: det⁡(adj(2A))=242=576\det\big(\text{adj}(2A)\big) = 24^{2} = 576. The two expressions differ because adj(2A)=22 adj(A)\text{adj}(2A) = 2^{2}\,\text{adj}(A): each cofactor is a 2×22 \times 2 determinant and picks up 222^{2}, not 22. Check: det⁡(4 adj A)=43×9=576\det(4\,\text{adj}\,A) = 4^{3} \times 9 = 576.

d) Expand det⁡(adj B)\det(\text{adj}\,B) along row 11: 1(2−0)−1(2−0)+(−2)(−1−1)=2−2+4=41(2 - 0) - 1(2 - 0) + (-2)(-1 - 1) = 2 - 2 + 4 = 4. By a), det⁡(adj B)=det⁡(B)2\det(\text{adj}\,B) = \det(B)^{2}, so det⁡(B)2=4\det(B)^{2} = 4 and det⁡B=2\det B = 2 or −2-2. The adjugate cannot tell: each cofactor of a 3×33 \times 3 is a 2×22 \times 2 determinant, and negating the matrix multiplies each of them by (−1)2=1(-1)^{2} = 1. So adj(−B)=adj(B)\text{adj}(-B) = \text{adj}(B) while det⁡(−B)=−det⁡B\det(-B) = -\det B: the two candidates BB and −B-B have the same adjugate and opposite determinants, and no computation on adj(B)\text{adj}(B) alone can separate them.

e) For an invertible MM, adj(M)=det⁡(M) M−1\text{adj}(M) = \det(M)\,M^{-1}. Apply it to M=adj(B)M = \text{adj}(B): adj(adj B)=det⁡(B)2⋅Bdet⁡B=det⁡(B) B\text{adj}(\text{adj}\,B) = \det(B)^{2} \cdot \dfrac{B}{\det B} = \det(B)\,B, using (det⁡(B) B−1)−1=Bdet⁡B(\det(B)\,B^{-1})^{-1} = \dfrac{B}{\det B}. The nine cofactors of adj(B)\text{adj}(B), transposed, give adj(adj B)=(202−24−2−220)\text{adj}(\text{adj}\,B) = \begin{pmatrix} 2 & 0 & 2 \\ -2 & 4 & -2 \\ -2 & 2 & 0 \end{pmatrix}. Its top-left entry is 2=det⁡(B)×12 = \det(B) \times 1, so det⁡B=2\det B = 2, and B=12(202−24−2−220)=(101−12−1−110)B = \dfrac{1}{2}\begin{pmatrix} 2 & 0 & 2 \\ -2 & 4 & -2 \\ -2 & 2 & 0 \end{pmatrix} = \begin{pmatrix} 1 & 0 & 1 \\ -1 & 2 & -1 \\ -1 & 1 & 0 \end{pmatrix}. Finally B−1=12 adj(B)=(1212−11212012−121)B^{-1} = \dfrac{1}{2}\,\text{adj}(B) = \begin{pmatrix} \frac{1}{2} & \frac{1}{2} & -1 \\ \frac{1}{2} & \frac{1}{2} & 0 \\ \frac{1}{2} & -\frac{1}{2} & 1 \end{pmatrix}. Check one entry of BB−1BB^{-1}: row 11 of BB against column 11 of B−1B^{-1} gives 12+0+12=1\frac{1}{2} + 0 + \frac{1}{2} = 1.

Exercise 8: Five statements to correct

Each statement below was written by a student in a MATH 133 tutorial. Say whether it is TRUE or FALSE, correct every false statement with the smallest counterexample that settles it, and state the correct rule.

  • a) det⁡(A+B)=det⁡A+det⁡B\det(A + B) = \det A + \det B.
  • b) det⁡(−A)=−det⁡A\det(-A) = -\det A for every square matrix AA.
  • c) Since det⁡(AB)=det⁡(BA)\det(AB) = \det(BA), the matrices ABAB and BABA are equal.
  • d) When det⁡A=0\det A = 0, Cramer's rule shows that the system Ax=bA\mathbf{x} = \mathbf{b} has no solution.
  • e) det⁡(A−1BA)=det⁡B\det(A^{-1}BA) = \det B holds only when AA and BB commute.
Show the solution

Answers

  • a) False: det⁡(I2+I2)=4\det(I_2 + I_2) = 4, not 22; no rule for a sum
  • b) False for even nn: det⁡(−I2)=1\det(-I_2) = 1; the rule is (−1)ndet⁡A(-1)^{n}\det A
  • c) False: equal determinants, different matrices
  • d) False: Cramer does not apply; there may be none or infinitely many
  • e) False: it holds for every invertible AA

a) FALSE. Take A=B=I2A = B = I_2: det⁡(A+B)=det⁡(2I2)=22=4\det(A + B) = \det(2I_2) = 2^{2} = 4, while det⁡A+det⁡B=1+1=2\det A + \det B = 1 + 1 = 2. The determinant is MULTIPLICATIVE, not additive: there is no formula for det⁡(A+B)\det(A + B) in terms of det⁡A\det A and det⁡B\det B. Correct statement: to find det⁡(A+B)\det(A + B), compute the matrix A+BA + B first, then its determinant. Even when both determinants are zero, the sum can be invertible: (1000)+(0001)=I2\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix} + \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix} = I_2.

b) FALSE for even sizes. −A-A is the scalar −1-1 times an n×nn \times n matrix, so det⁡(−A)=(−1)ndet⁡A\det(-A) = (-1)^{n}\det A. For n=2n = 2, det⁡(−I2)=(−1)(−1)−0=1=det⁡I2\det(-I_2) = (-1)(-1) - 0 = 1 = \det I_2, not −1-1. Correct statement: det⁡(−A)=−det⁡A\det(-A) = -\det A when nn is odd, and det⁡(−A)=det⁡A\det(-A) = \det A when nn is even. This is the knk^{n} rule with k=−1k = -1, and it is the reason proofs about skew-symmetric matrices need nn odd.

c) FALSE. The first half is true, det⁡(AB)=det⁡A det⁡B=det⁡(BA)\det(AB) = \det A\,\det B = \det(BA), and the conclusion does not follow: equal determinants do not make equal matrices. With A=(0100)A = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} and B=(0010)B = \begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}, AB=(1000)AB = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix} and BA=(0001)BA = \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix}: both determinants are 00 and the matrices differ. Correct statement: ABAB and BABA always have the same determinant, and are in general different matrices.

d) FALSE. Cramer's rule has the hypothesis det⁡A≠0\det A \neq 0; when det⁡A=0\det A = 0 it says nothing at all, because its formula would divide by zero. The system x+y=2x + y = 2, 2x+2y=42x + 2y = 4 has determinant 1⋅2−1⋅2=01 \cdot 2 - 1 \cdot 2 = 0 and infinitely many solutions, while x+y=2x + y = 2, 2x+2y=52x + 2y = 5 has the same determinant and none. Correct statement: when det⁡A=0\det A = 0, the system has either no solution or infinitely many, and row reduction of the augmented matrix decides which.

e) FALSE: it holds for every invertible AA. det⁡(A−1BA)=1det⁡A det⁡B det⁡A=det⁡B\det(A^{-1}BA) = \dfrac{1}{\det A}\,\det B\,\det A = \det B, because the determinants are numbers and numbers commute even when the matrices do not. With A=(1101)A = \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix} and B=(1002)B = \begin{pmatrix} 1 & 0 \\ 0 & 2 \end{pmatrix}, which do not commute, A−1BA=(1−102)A^{-1}BA = \begin{pmatrix} 1 & -1 \\ 0 & 2 \end{pmatrix}, a matrix different from BB whose determinant is still 22. Correct statement: A−1BAA^{-1}BA and BB always have the same determinant, whether or not AA and BB commute.

Exercise 9: Undoing a voxel remap exactly, without rounding

A medical scanner stores a 3D image as a grid of voxels, each one at a position with INTEGER coordinates (i,j,k)(i, j, k). A registration tool realigns two scans by moving every voxel from position v\mathbf{v} to position KvK\mathbf{v}, where KK is a 3×33 \times 3 matrix with integer entries, and it must be able to undo the move EXACTLY: sending KvK\mathbf{v} back to v\mathbf{v} with no rounding, since a rounded position is a different voxel. That requires K−1K^{-1} to have integer entries too, and the adjugate formula says exactly when it does.

The remap is K=(121011122)K = \begin{pmatrix} 1 & 2 & 1 \\ 0 & 1 & 1 \\ 1 & 2 & 2 \end{pmatrix}, and two voxels flagged by a radiologist sit at positions (44,23,58)(44, 23, 58) and (25,19,43)(25, 19, 43) after the remap. The figure shows the same question in the plane.

det 1: every gridpoint is reacheddet -2: (1, 0) isnever reached
  • a) Compute det⁡K\det K, naming the row operation you use, and explain why K−1K^{-1} has integer entries before computing it.
  • b) Compute K−1K^{-1} by the adjugate.
  • c) Find the original positions of the two flagged voxels, and check one of them.
  • d) A colleague proposes the remap L=(1101−10001)L = \begin{pmatrix} 1 & 1 & 0 \\ 1 & -1 & 0 \\ 0 & 0 & 1 \end{pmatrix}. Compute det⁡L\det L and the entry in row 11, column 22 of L−1L^{-1}, and say what this remap does to the voxel grid.
  • e) Prove that if a matrix with integer entries has an inverse with integer entries, its determinant is 11 or −1-1.
Show the solution

Answers

  • a) det⁡K=1\det K = 1, so K−1=adj(K)K^{-1} = \text{adj}(K), an integer matrix
  • b) K−1=(0−2111−1−101)K^{-1} = \begin{pmatrix} 0 & -2 & 1 \\ 1 & 1 & -1 \\ -1 & 0 & 1 \end{pmatrix}
  • c) (12,9,14)(12, 9, 14) and (5,1,18)(5, 1, 18)
  • d) det⁡L=−2\det L = -2 and (L−1)12=12(L^{-1})_{12} = \frac{1}{2}: some positions come back as half-voxels
  • e) det⁡K det⁡(K−1)=1\det K \, \det(K^{-1}) = 1 with two integers, so det⁡K=±1\det K = \pm 1

a) R3→R3−R1R_3 \to R_3 - R_1 does not change the determinant and turns row 33 into (0,0,1)(0, 0, 1). Expanding along that row, with sign (−1)3+3=+1(-1)^{3+3} = +1: det⁡K=1⋅det⁡(1201)=1\det K = 1 \cdot \det\begin{pmatrix} 1 & 2 \\ 0 & 1 \end{pmatrix} = 1. Now K−1=1det⁡K adj(K)=adj(K)K^{-1} = \dfrac{1}{\det K}\,\text{adj}(K) = \text{adj}(K), and every entry of the adjugate is a signed 2×22 \times 2 determinant of INTEGERS, hence an integer. So K−1K^{-1} is an integer matrix, and we know it before computing a single entry. That is the entire reason a remap with determinant ±1\pm 1 is chosen: it can be undone in integer arithmetic, exactly.

b) Cofactors, with the checkerboard signs. Row 11: C11=+(1⋅2−1⋅2)=0C_{11} = +(1 \cdot 2 - 1 \cdot 2) = 0, C12=−(0⋅2−1⋅1)=1C_{12} = -(0 \cdot 2 - 1 \cdot 1) = 1, C13=+(0⋅2−1⋅1)=−1C_{13} = +(0 \cdot 2 - 1 \cdot 1) = -1. Row 22: C21=−(2⋅2−1⋅2)=−2C_{21} = -(2 \cdot 2 - 1 \cdot 2) = -2, C22=+(1⋅2−1⋅1)=1C_{22} = +(1 \cdot 2 - 1 \cdot 1) = 1, C23=−(1⋅2−2⋅1)=0C_{23} = -(1 \cdot 2 - 2 \cdot 1) = 0. Row 33: C31=+(2⋅1−1⋅1)=1C_{31} = +(2 \cdot 1 - 1 \cdot 1) = 1, C32=−(1⋅1−1⋅0)=−1C_{32} = -(1 \cdot 1 - 1 \cdot 0) = -1, C33=+(1⋅1−2⋅0)=1C_{33} = +(1 \cdot 1 - 2 \cdot 0) = 1. Transposing the matrix of cofactors: K−1=adj(K)=(0−2111−1−101)K^{-1} = \text{adj}(K) = \begin{pmatrix} 0 & -2 & 1 \\ 1 & 1 & -1 \\ -1 & 0 & 1 \end{pmatrix}. Check row 11 of KK against column 11 of K−1K^{-1}: 0+2−1=10 + 2 - 1 = 1.

c) K−1(44,23,58)=(0−46+58, 44+23−58, −44+0+58)=(12,9,14)K^{-1}(44, 23, 58) = (0 - 46 + 58,\ 44 + 23 - 58,\ -44 + 0 + 58) = (12, 9, 14) and K−1(25,19,43)=(0−38+43, 25+19−43, −25+0+43)=(5,1,18)K^{-1}(25, 19, 43) = (0 - 38 + 43,\ 25 + 19 - 43,\ -25 + 0 + 43) = (5, 1, 18). Both are integer positions, as a) promised, so they are genuine voxels of the original scan and no rounding was needed. Re-applying the remap is the check: K(12,9,14)=(12+18+14, 9+14, 12+18+28)=(44,23,58)K(12, 9, 14) = (12 + 18 + 14,\ 9 + 14,\ 12 + 18 + 28) = (44, 23, 58), the flagged position. Had one of the two results come out as a fraction, the error would be in K−1K^{-1}, not in the scan: an integer matrix with determinant 11 cannot produce one.

d) LL is block diagonal: det⁡L=det⁡(111−1)×1=(−1−1)×1=−2\det L = \det\begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix} \times 1 = (-1 - 1) \times 1 = -2. The entry in row 11, column 22 of L−1L^{-1} is C21det⁡L\dfrac{C_{21}}{\det L}, indices swapped: deleting row 22 and column 11 leaves det⁡(1001)=1\det\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} = 1, so C21=−1C_{21} = -1 and (L−1)12=−1−2=12(L^{-1})_{12} = \dfrac{-1}{-2} = \dfrac{1}{2}. Now undo the remap on the voxel at (1,0,0)(1, 0, 0): L−1(1,0,0)L^{-1}(1, 0, 0) is the first column of L−1L^{-1}, (12,12,0)\left(\frac{1}{2}, \frac{1}{2}, 0\right), a HALF-voxel. So no voxel of the original grid is sent to (1,0,0)(1, 0, 0): LL leaves holes in the new grid, one position out of two in every plane, since (a+b,a−b)(a + b, a - b) always has an even sum. The figure shows it in the plane: the shaded cell of the det −2-2 map has area 22 and a grid point at its centre that nothing reaches.

e) Let KK have integer entries and suppose K−1K^{-1} does too. Every determinant of an integer matrix is an integer, since the cofactor expansion only adds and multiplies integers. From KK−1=IKK^{-1} = I, det⁡K×det⁡(K−1)=1\det K \times \det(K^{-1}) = 1: two INTEGERS whose product is 11, so both are 11 or both are −1-1. Conversely, part a) showed that det⁡K=±1\det K = \pm 1 makes K−1=± adj(K)K^{-1} = \pm\,\text{adj}(K) an integer matrix. So an integer matrix has an integer inverse exactly when its determinant is ±1\pm 1.

Exercise 10: One current is all you need: Cramer's rule on a circuit

In an electrical network, the loop currents I1I_1, I2I_2, I3I_3, in amperes, satisfy one linear equation per loop. An engineer rarely needs all of them: often one branch carries a sensor or a fuse, and only its current matters. That is the situation where Cramer's rule beats elimination, since one unknown costs two determinants and nothing else.

In the network studied here, a 1212 volt battery sits in loop 11 and a second source of EE volts in loop 33. The loop equations are {2I1−I2=12−I1+3I2−I3=0−I2+4I3=−E\begin{cases} 2I_1 - I_2 = 12 \\ -I_1 + 3I_2 - I_3 = 0 \\ -I_2 + 4I_3 = -E \end{cases}, and the sensor is in loop 22.

  • a) Compute the determinant of the coefficient matrix NN, and explain why the system has exactly one solution whatever the value of EE.
  • b) For E=6E = 6, find I2I_2 alone by Cramer's rule.
  • c) Keep EE as a letter and express I2I_2 as a function of EE. For which EE does the sensor read zero?
  • d) For E=6E = 6, find I3I_3 by Cramer's rule and interpret its sign.
  • e) Using one entry of N−1N^{-1}, say by how much I2I_2 changes when EE increases by one volt, and check against c).
Show the solution

Answers

  • a) det⁡N=18≠0\det N = 18 \neq 0, and NN does not depend on EE
  • b) I2=3618=2I_2 = \frac{36}{18} = 2 A
  • c) I2=24−E9I_2 = \frac{24 - E}{9}; zero for E=24E = 24 V
  • d) I3=−1I_3 = -1 A: 11 A flowing against the chosen direction
  • e) (N−1)23=19(N^{-1})_{23} = \frac{1}{9}, so I2I_2 drops by 19\frac{1}{9} A per volt

a) N=(2−10−13−10−14)N = \begin{pmatrix} 2 & -1 & 0 \\ -1 & 3 & -1 \\ 0 & -1 & 4 \end{pmatrix}. Expand along row 11, which has a zero: det⁡N=2det⁡(3−1−14)−(−1)det⁡(−1−104)+0=2(12−1)+(−4−0)=22−4=18\det N = 2\det\begin{pmatrix} 3 & -1 \\ -1 & 4 \end{pmatrix} - (-1)\det\begin{pmatrix} -1 & -1 \\ 0 & 4 \end{pmatrix} + 0 = 2(12 - 1) + (-4 - 0) = 22 - 4 = 18. The source EE appears only on the right-hand side, never in NN, so det⁡N=18≠0\det N = 18 \neq 0 for every EE: the matrix is invertible, the system has exactly one solution, and Cramer's rule applies whatever the sources are. That is a physical statement too: a network of resistors always settles into one set of currents.

b) Replace column 22 by the right-hand side (12,0,−6)(12, 0, -6) and expand along the new column 22, whose middle entry is 00: 12 C12+(−6) C3212\,C_{12} + (-6)\,C_{32}. Here C12=−det⁡(−1−104)=−(−4)=4C_{12} = -\det\begin{pmatrix} -1 & -1 \\ 0 & 4 \end{pmatrix} = -(-4) = 4 and C32=−det⁡(20−1−1)=−(−2)=2C_{32} = -\det\begin{pmatrix} 2 & 0 \\ -1 & -1 \end{pmatrix} = -(-2) = 2, so the numerator is 48−12=3648 - 12 = 36 and I2=3618=2I_2 = \dfrac{36}{18} = 2 A. The cofactors of column 22 were taken from NN itself, not from the modified matrix, since deleting column 22 removes exactly the column that was replaced.

c) The same expansion with −E-E in place of −6-6 gives a numerator of 12×4+(−E)×2=48−2E12 \times 4 + (-E) \times 2 = 48 - 2E, so I2=48−2E18=24−E9I_2 = \dfrac{48 - 2E}{18} = \dfrac{24 - E}{9}. The sensor reads zero when E=24E = 24 volts. For E=6E = 6 this gives 189=2\dfrac{18}{9} = 2, as in b). Keeping EE as a letter costs nothing with Cramer's rule, because the letter enters only through one column; with elimination it would travel through every row operation.

d) Replace column 33 by (12,0,−6)(12, 0, -6) and expand along it: 12 C13+(−6) C3312\,C_{13} + (-6)\,C_{33}, with C13=+det⁡(−130−1)=1C_{13} = +\det\begin{pmatrix} -1 & 3 \\ 0 & -1 \end{pmatrix} = 1 and C33=+det⁡(2−1−13)=5C_{33} = +\det\begin{pmatrix} 2 & -1 \\ -1 & 3 \end{pmatrix} = 5. The numerator is 12−30=−1812 - 30 = -18, so I3=−1I_3 = -1 A. The negative sign is not an error: a current of 11 ampere flows in loop 33 in the direction OPPOSITE to the one chosen when the equations were written. Checking in the three equations with I1=7I_1 = 7: 14−2=1214 - 2 = 12, −7+6+1=0-7 + 6 + 1 = 0, −2−4=−6-2 - 4 = -6.

e) Since I=N−1(12,0,−E)\mathbf{I} = N^{-1}(12, 0, -E), the current I2I_2 is row 22 of N−1N^{-1} against the right-hand side, and EE enters through the entry in row 22, column 33: (N−1)23=C32det⁡N=218=19(N^{-1})_{23} = \dfrac{C_{32}}{\det N} = \dfrac{2}{18} = \dfrac{1}{9}, indices swapped. With −E-E in the third position, each extra volt changes I2I_2 by −19-\dfrac{1}{9} ampere. This is the slope of c), I2=249−E9I_2 = \dfrac{24}{9} - \dfrac{E}{9}, and the same row gives (N−1)21=C1218=29(N^{-1})_{21} = \dfrac{C_{12}}{18} = \dfrac{2}{9}, the 249\dfrac{24}{9} coming from 12×2912 \times \dfrac{2}{9}. One entry of the inverse is a sensitivity, and it costs one cofactor.

See also

Struggling with MATH 133?

I tutor first-year linear algebra at McGill and Concordia, in English or in French, in Montreal or online. Get in touch for a first session.

Site by Studio Squalli