About

2.2 The Inverse of a Matrix

Theory: The Inverse of a Matrix

The Inverse of a Matrix (Lay §2.2)

An n×nn \times n matrix AA is invertible if there is an n×nn \times n matrix CC with

CA=IandAC=ICA = I \quad \text{and} \quad AC = I

This CC is unique, is called the inverse of AA, and is written A−1A^{-1}:

A−1A=IandAA−1=IA^{-1}A = I \quad \text{and} \quad AA^{-1} = I

Theorem 4 (Lay §2.2): Let A=[abcd]A = \begin{bmatrix} a & b \\ c & d \end{bmatrix}. If ad−bc≠0ad - bc \ne 0, then AA is invertible and

A−1=1ad−bc[d−b−ca]A^{-1} = \frac{1}{ad - bc}\begin{bmatrix} d & -b \\ -c & a \end{bmatrix}

If ad−bc=0ad - bc = 0, then AA is not invertible.

The number ad−bcad - bc is the determinant of AA: det⁡A=ad−bc\det A = ad - bc. So a 2×22 \times 2 matrix is invertible if and only if det⁡A≠0\det A \ne 0.

Recipe for the 2×22 \times 2 formula: swap the diagonal entries a↔da \leftrightarrow d, negate the off-diagonal entries b,cb, c, divide by ad−bcad - bc.

Mini-example. A=[3456]A = \begin{bmatrix} 3 & 4 \\ 5 & 6 \end{bmatrix}: det⁡A=18−20=−2\det A = 18 - 20 = -2, so

A−1=1−2[6−4−53]=[−325/2−3/2]A^{-1} = \frac{1}{-2}\begin{bmatrix} 6 & -4 \\ -5 & 3 \end{bmatrix} = \begin{bmatrix} -3 & 2 \\ 5/2 & -3/2 \end{bmatrix}

Solving Ax=bA\mathbf{x} = \mathbf{b} with A−1A^{-1}

Theorem 5 (Lay §2.2): If AA is an invertible n×nn \times n matrix, then for each b\mathbf{b} in Rn\mathbb{R}^n, the equation Ax=bA\mathbf{x} = \mathbf{b} has the unique solution

x=A−1b\mathbf{x} = A^{-1}\mathbf{b}

Why: A(A−1b)=Ib=bA(A^{-1}\mathbf{b}) = I\mathbf{b} = \mathbf{b}, so it is a solution; and if Au=bA\mathbf{u} = \mathbf{b}, multiply on the left by A−1A^{-1} to get u=A−1b\mathbf{u} = A^{-1}\mathbf{b}, so it is the only one.

In practice, row reducing [Ab]\left[\begin{array}{c|c} A & \mathbf{b} \end{array}\right] is faster than computing A−1A^{-1} first. Use x=A−1b\mathbf{x} = A^{-1}\mathbf{b} when the question gives you A−1A^{-1}, when AA is 2×22 \times 2, or when you must solve many systems with the same AA.

Properties of Inverses

Theorem 6 (Lay §2.2):

  1. If AA is invertible, then A−1A^{-1} is invertible and (A−1)−1=A(A^{-1})^{-1} = A.
  2. If AA and BB are n×nn \times n invertible matrices, then so is ABAB, and (AB)−1=B−1A−1(AB)^{-1} = B^{-1}A^{-1}.
  3. If AA is invertible, then so is ATA^T, and (AT)−1=(A−1)T(A^T)^{-1} = (A^{-1})^T.

Generalisation: a product of invertible matrices is invertible, and its inverse is the product of the inverses in reverse order: (ABC)−1=C−1B−1A−1(ABC)^{-1} = C^{-1}B^{-1}A^{-1}.

Common mistake: (AB)−1≠A−1B−1(AB)^{-1} \ne A^{-1}B^{-1} in general. Check: (AB)(B−1A−1)=A(BB−1)A−1=AIA−1=AA−1=I(AB)(B^{-1}A^{-1}) = A(BB^{-1})A^{-1} = AIA^{-1} = AA^{-1} = I. The reverse order is what makes the middle cancel (“socks and shoes”: you take off the last thing you put on first).

Elementary Matrices (brief)

An elementary matrix is obtained by performing one elementary row operation (switch / scale / replace) on the identity InI_n.

Mini-example. E=[100010−401]E = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ -4 & 0 & 1 \end{bmatrix} comes from R3←R3−4R1R_3 \leftarrow R_3 - 4R_1 on I3I_3, so EAEA subtracts 4 times row 1 of AA from row 3 of AA.

Theorem 7 (Lay §2.2): An n×nn \times n matrix AA is invertible if and only if AA is row equivalent to InI_n. In that case, any sequence of elementary row operations that reduces AA to InI_n also transforms InI_n into A−1A^{-1}.

Why: if Ep⋯E1A=IE_p \cdots E_1 A = I, then A−1=Ep⋯E1=Ep⋯E1IA^{-1} = E_p \cdots E_1 = E_p \cdots E_1 I, so the same row operations applied to II produce A−1A^{-1}.

The Algorithm for Finding A−1A^{-1}

  1. Form the n×2nn \times 2n matrix [AI]\left[\begin{array}{c|c} A & I \end{array}\right].
  2. Row reduce it.
  3. If AA reduces to II, then [AI]∼[IA−1]\left[\begin{array}{c|c} A & I \end{array}\right] \sim \left[\begin{array}{c|c} I & A^{-1} \end{array}\right] and A−1A^{-1} is read off the right half.
  4. Otherwise (a row of zeros appears on the left, i.e. fewer than nn pivots), AA has no inverse.

Example 7 (Lay §2.2). Find the inverse of A=[0121034−38]A = \begin{bmatrix} 0 & 1 & 2 \\ 1 & 0 & 3 \\ 4 & -3 & 8 \end{bmatrix}, if it exists.

[0121001030104−38001]→R1↔R2[1030100121004−38001]\left[\begin{array}{ccc|ccc} 0 & 1 & 2 & 1 & 0 & 0 \\ 1 & 0 & 3 & 0 & 1 & 0 \\ 4 & -3 & 8 & 0 & 0 & 1 \end{array}\right] \xrightarrow{R_1 \leftrightarrow R_2} \left[\begin{array}{ccc|ccc} 1 & 0 & 3 & 0 & 1 & 0 \\ 0 & 1 & 2 & 1 & 0 & 0 \\ 4 & -3 & 8 & 0 & 0 & 1 \end{array}\right]

R3←R3−4R1R_3 \leftarrow R_3 - 4R_1, then R3←R3+3R2R_3 \leftarrow R_3 + 3R_2:

∼[1030100121000−3−40−41]∼[1030100121000023−41]\sim \left[\begin{array}{ccc|ccc} 1 & 0 & 3 & 0 & 1 & 0 \\ 0 & 1 & 2 & 1 & 0 & 0 \\ 0 & -3 & -4 & 0 & -4 & 1 \end{array}\right] \sim \left[\begin{array}{ccc|ccc} 1 & 0 & 3 & 0 & 1 & 0 \\ 0 & 1 & 2 & 1 & 0 & 0 \\ 0 & 0 & 2 & 3 & -4 & 1 \end{array}\right]

R3←12R3R_3 \leftarrow \tfrac12 R_3, then R1←R1−3R3R_1 \leftarrow R_1 - 3R_3 and R2←R2−2R3R_2 \leftarrow R_2 - 2R_3:

∼[1030100121000013/2−21/2]∼[100−9/27−3/2010−24−10013/2−21/2]\sim \left[\begin{array}{ccc|ccc} 1 & 0 & 3 & 0 & 1 & 0 \\ 0 & 1 & 2 & 1 & 0 & 0 \\ 0 & 0 & 1 & 3/2 & -2 & 1/2 \end{array}\right] \sim \left[\begin{array}{ccc|ccc} 1 & 0 & 0 & -9/2 & 7 & -3/2 \\ 0 & 1 & 0 & -2 & 4 & -1 \\ 0 & 0 & 1 & 3/2 & -2 & 1/2 \end{array}\right]

The left half is II, so AA is invertible and

A−1=[−9/27−3/2−24−13/2−21/2]A^{-1} = \begin{bmatrix} -9/2 & 7 & -3/2 \\ -2 & 4 & -1 \\ 3/2 & -2 & 1/2 \end{bmatrix}

Always check: multiply one row of AA by one column of your A−1A^{-1}. E.g. row 1 of AA times column 1 of A−1A^{-1}: 0(−9/2)+1(−2)+2(3/2)=10(-9/2) + 1(-2) + 2(3/2) = 1 ✓, matching the (1,1)(1,1) entry of II.

Finding Just One Column of A−1A^{-1}

Since AA−1=IAA^{-1} = I, column jj of A−1A^{-1} is the solution of

Ax=ejA\mathbf{x} = \mathbf{e}_j

where ej\mathbf{e}_j is column jj of II. So to get only column jj of A−1A^{-1}, row reduce [Aej]\left[\begin{array}{c|c} A & \mathbf{e}_j \end{array}\right] instead of the full [AI]\left[\begin{array}{c|c} A & I \end{array}\right]. (The full algorithm is just this done for all jj at once.)


Solved Examples (Textbook Questions)

Exercise 1

(Adapted from Lay §2.2, Exercises 1–4)

Find the inverses of the matrices.

  1. [8654]\begin{bmatrix} 8 & 6 \\ 5 & 4 \end{bmatrix}
  2. [3274]\begin{bmatrix} 3 & 2 \\ 7 & 4 \end{bmatrix}
  3. [85−7−5]\begin{bmatrix} 8 & 5 \\ -7 & -5 \end{bmatrix}
  4. [3−47−8]\begin{bmatrix} 3 & -4 \\ 7 & -8 \end{bmatrix}
Show solution ↓
Solution

Use Theorem 4: compute ad−bcad - bc, swap the diagonal, negate the off-diagonal, divide.

1. det⁡=8(4)−6(5)=32−30=2\det = 8(4) - 6(5) = 32 - 30 = 2.

[8654]−1=12[4−6−58]=[2−3−5/24]\begin{bmatrix} 8 & 6 \\ 5 & 4 \end{bmatrix}^{-1} = \frac{1}{2}\begin{bmatrix} 4 & -6 \\ -5 & 8 \end{bmatrix} = \begin{bmatrix} 2 & -3 \\ -5/2 & 4 \end{bmatrix}

2. det⁡=3(4)−2(7)=12−14=−2\det = 3(4) - 2(7) = 12 - 14 = -2.

[3274]−1=1−2[4−2−73]=[−217/2−3/2]\begin{bmatrix} 3 & 2 \\ 7 & 4 \end{bmatrix}^{-1} = \frac{1}{-2}\begin{bmatrix} 4 & -2 \\ -7 & 3 \end{bmatrix} = \begin{bmatrix} -2 & 1 \\ 7/2 & -3/2 \end{bmatrix}

3. det⁡=8(−5)−5(−7)=−40+35=−5\det = 8(-5) - 5(-7) = -40 + 35 = -5.

[85−7−5]−1=−15[−5−578]=[11−1.4−1.6]\begin{bmatrix} 8 & 5 \\ -7 & -5 \end{bmatrix}^{-1} = -\frac{1}{5}\begin{bmatrix} -5 & -5 \\ 7 & 8 \end{bmatrix} = \begin{bmatrix} 1 & 1 \\ -1.4 & -1.6 \end{bmatrix}

4. det⁡=3(−8)−(−4)(7)=−24+28=4\det = 3(-8) - (-4)(7) = -24 + 28 = 4.

[3−47−8]−1=14[−84−73]=[−21−7/43/4]\begin{bmatrix} 3 & -4 \\ 7 & -8 \end{bmatrix}^{-1} = \frac{1}{4}\begin{bmatrix} -8 & 4 \\ -7 & 3 \end{bmatrix} = \begin{bmatrix} -2 & 1 \\ -7/4 & 3/4 \end{bmatrix}

Check (part 1): [8654][2−3−5/24]=[16−15−24+2410−10−15+16]=I\begin{bmatrix} 8 & 6 \\ 5 & 4 \end{bmatrix}\begin{bmatrix} 2 & -3 \\ -5/2 & 4 \end{bmatrix} = \begin{bmatrix} 16 - 15 & -24 + 24 \\ 10 - 10 & -15 + 16 \end{bmatrix} = I ✓

Common mistake: watch the sign of bcbc when bb or cc is negative. In part 3, ad−bc=−40−(5)(−7)=−40+35ad - bc = -40 - (5)(-7) = -40 + 35, not −40−35-40 - 35.

Exercise 2

(Adapted from Lay §2.2, Exercise 7)

Let A=[12512]A = \begin{bmatrix} 1 & 2 \\ 5 & 12 \end{bmatrix}, b1=[−13]\mathbf{b}_1 = \begin{bmatrix} -1 \\ 3 \end{bmatrix}, b2=[1−5]\mathbf{b}_2 = \begin{bmatrix} 1 \\ -5 \end{bmatrix}, b3=[26]\mathbf{b}_3 = \begin{bmatrix} 2 \\ 6 \end{bmatrix}, and b4=[35]\mathbf{b}_4 = \begin{bmatrix} 3 \\ 5 \end{bmatrix}.

a. Find A−1A^{-1}, and use it to solve the four equations Ax=b1A\mathbf{x} = \mathbf{b}_1, Ax=b2A\mathbf{x} = \mathbf{b}_2, Ax=b3A\mathbf{x} = \mathbf{b}_3, Ax=b4A\mathbf{x} = \mathbf{b}_4.

b. The four equations in part (a) can be solved by the same set of row operations, since the coefficient matrix is the same in each case. Solve the four equations in part (a) by row reducing the augmented matrix [Ab1b2b3b4]\begin{bmatrix} A & \mathbf{b}_1 & \mathbf{b}_2 & \mathbf{b}_3 & \mathbf{b}_4 \end{bmatrix}.

Show solution ↓
Solution

(a) det⁡A=1(12)−2(5)=2\det A = 1(12) - 2(5) = 2, so

A−1=12[12−2−51]=[6−1−2.50.5]A^{-1} = \frac{1}{2}\begin{bmatrix} 12 & -2 \\ -5 & 1 \end{bmatrix} = \begin{bmatrix} 6 & -1 \\ -2.5 & 0.5 \end{bmatrix}

By Theorem 5, x=A−1b\mathbf{x} = A^{-1}\mathbf{b}:

A−1b1=12[12−2−51][−13]=12[−188]=[−94]A^{-1}\mathbf{b}_1 = \frac{1}{2}\begin{bmatrix} 12 & -2 \\ -5 & 1 \end{bmatrix}\begin{bmatrix} -1 \\ 3 \end{bmatrix} = \frac{1}{2}\begin{bmatrix} -18 \\ 8 \end{bmatrix} = \begin{bmatrix} -9 \\ 4 \end{bmatrix}A−1b2=12[12+10−5−5]=[11−5],A−1b3=12[24−12−10+6]=[6−2],A−1b4=12[36−10−15+5]=[13−5]A^{-1}\mathbf{b}_2 = \frac12\begin{bmatrix} 12 + 10 \\ -5 - 5 \end{bmatrix} = \begin{bmatrix} 11 \\ -5 \end{bmatrix}, \quad A^{-1}\mathbf{b}_3 = \frac12\begin{bmatrix} 24 - 12 \\ -10 + 6 \end{bmatrix} = \begin{bmatrix} 6 \\ -2 \end{bmatrix}, \quad A^{-1}\mathbf{b}_4 = \frac12\begin{bmatrix} 36 - 10 \\ -15 + 5 \end{bmatrix} = \begin{bmatrix} 13 \\ -5 \end{bmatrix}

(b) Put all four right-hand sides next to AA and row reduce once.

R2←R2−5R1R_2 \leftarrow R_2 - 5R_1:

[12−11235123−565]∼[12−1123028−10−4−10]\left[\begin{array}{cc|cccc} 1 & 2 & -1 & 1 & 2 & 3 \\ 5 & 12 & 3 & -5 & 6 & 5 \end{array}\right] \sim \left[\begin{array}{cc|cccc} 1 & 2 & -1 & 1 & 2 & 3 \\ 0 & 2 & 8 & -10 & -4 & -10 \end{array}\right]

R2←12R2R_2 \leftarrow \tfrac12 R_2, then R1←R1−2R2R_1 \leftarrow R_1 - 2R_2:

∼[12−1123014−5−2−5]∼[10−911613014−5−2−5]\sim \left[\begin{array}{cc|cccc} 1 & 2 & -1 & 1 & 2 & 3 \\ 0 & 1 & 4 & -5 & -2 & -5 \end{array}\right] \sim \left[\begin{array}{cc|cccc} 1 & 0 & -9 & 11 & 6 & 13 \\ 0 & 1 & 4 & -5 & -2 & -5 \end{array}\right]

Each column on the right is the solution of the matching system:

[−94],[11−5],[6−2],[13−5]\begin{bmatrix} -9 \\ 4 \end{bmatrix}, \quad \begin{bmatrix} 11 \\ -5 \end{bmatrix}, \quad \begin{bmatrix} 6 \\ -2 \end{bmatrix}, \quad \begin{bmatrix} 13 \\ -5 \end{bmatrix}

the same as in part (a).

This is exactly the idea behind the [AI]\left[\begin{array}{c|c} A & I \end{array}\right] algorithm: you can solve many systems with the same coefficient matrix by row reducing once with all right-hand sides attached.

Exercise 3

(Adapted from Lay §2.2, Exercises 29 and 30)

Find the inverses of the matrices, if they exist. Use the algorithm [AI]∼[IA−1]\left[\begin{array}{c|c} A & I \end{array}\right] \sim \left[\begin{array}{c|c} I & A^{-1} \end{array}\right].

  1. [1247]\begin{bmatrix} 1 & 2 \\ 4 & 7 \end{bmatrix}
  2. [51047]\begin{bmatrix} 5 & 10 \\ 4 & 7 \end{bmatrix}
Show solution ↓
Solution

Part 1 (Exercise 29). R2←R2−4R1R_2 \leftarrow R_2 - 4R_1, then R2←−R2R_2 \leftarrow -R_2, then R1←R1−2R2R_1 \leftarrow R_1 - 2R_2:

[12104701]∼[12100−1−41]∼[1210014−1]∼[10−72014−1]\left[\begin{array}{cc|cc} 1 & 2 & 1 & 0 \\ 4 & 7 & 0 & 1 \end{array}\right] \sim \left[\begin{array}{cc|cc} 1 & 2 & 1 & 0 \\ 0 & -1 & -4 & 1 \end{array}\right] \sim \left[\begin{array}{cc|cc} 1 & 2 & 1 & 0 \\ 0 & 1 & 4 & -1 \end{array}\right] \sim \left[\begin{array}{cc|cc} 1 & 0 & -7 & 2 \\ 0 & 1 & 4 & -1 \end{array}\right]A−1=[−724−1]A^{-1} = \begin{bmatrix} -7 & 2 \\ 4 & -1 \end{bmatrix}

(Check with Theorem 4: det⁡A=7−8=−1\det A = 7 - 8 = -1, so A−1=−[7−2−41]A^{-1} = -\begin{bmatrix} 7 & -2 \\ -4 & 1 \end{bmatrix} ✓.)

Part 2 (Exercise 30). R1←15R1R_1 \leftarrow \tfrac15 R_1, then R2←R2−4R1R_2 \leftarrow R_2 - 4R_1:

[510104701]∼[121/504701]∼[121/500−1−4/51]\left[\begin{array}{cc|cc} 5 & 10 & 1 & 0 \\ 4 & 7 & 0 & 1 \end{array}\right] \sim \left[\begin{array}{cc|cc} 1 & 2 & 1/5 & 0 \\ 4 & 7 & 0 & 1 \end{array}\right] \sim \left[\begin{array}{cc|cc} 1 & 2 & 1/5 & 0 \\ 0 & -1 & -4/5 & 1 \end{array}\right]

R2←−R2R_2 \leftarrow -R_2, then R1←R1−2R2R_1 \leftarrow R_1 - 2R_2:

∼[121/50014/5−1]∼[10−7/52014/5−1]\sim \left[\begin{array}{cc|cc} 1 & 2 & 1/5 & 0 \\ 0 & 1 & 4/5 & -1 \end{array}\right] \sim \left[\begin{array}{cc|cc} 1 & 0 & -7/5 & 2 \\ 0 & 1 & 4/5 & -1 \end{array}\right]A−1=[−7/524/5−1]A^{-1} = \begin{bmatrix} -7/5 & 2 \\ 4/5 & -1 \end{bmatrix}

(Check: det⁡A=35−40=−5\det A = 35 - 40 = -5, and −15[7−10−45]=[−7/524/5−1]-\tfrac15\begin{bmatrix} 7 & -10 \\ -4 & 5 \end{bmatrix} = \begin{bmatrix} -7/5 & 2 \\ 4/5 & -1 \end{bmatrix} ✓.)

Exercise 4

(Adapted from Lay §2.2, Exercises 31 and 32)

Find the inverses of the matrices, if they exist. Use the algorithm introduced in this section.

  1. A=[10−2−3142−34]A = \begin{bmatrix} 1 & 0 & -2 \\ -3 & 1 & 4 \\ 2 & -3 & 4 \end{bmatrix}
  2. A=[1−214−73−26−4]A = \begin{bmatrix} 1 & -2 & 1 \\ 4 & -7 & 3 \\ -2 & 6 & -4 \end{bmatrix}
Show solution ↓
Solution

Part 1 (Exercise 31). R2←R2+3R1R_2 \leftarrow R_2 + 3R_1 and R3←R3−2R1R_3 \leftarrow R_3 - 2R_1:

[10−2100−3140102−34001]∼[10−210001−23100−38−201]\left[\begin{array}{ccc|ccc} 1 & 0 & -2 & 1 & 0 & 0 \\ -3 & 1 & 4 & 0 & 1 & 0 \\ 2 & -3 & 4 & 0 & 0 & 1 \end{array}\right] \sim \left[\begin{array}{ccc|ccc} 1 & 0 & -2 & 1 & 0 & 0 \\ 0 & 1 & -2 & 3 & 1 & 0 \\ 0 & -3 & 8 & -2 & 0 & 1 \end{array}\right]

R3←R3+3R2R_3 \leftarrow R_3 + 3R_2:

∼[10−210001−2310002731]\sim \left[\begin{array}{ccc|ccc} 1 & 0 & -2 & 1 & 0 & 0 \\ 0 & 1 & -2 & 3 & 1 & 0 \\ 0 & 0 & 2 & 7 & 3 & 1 \end{array}\right]

Three pivots, so AA is invertible. R1←R1+R3R_1 \leftarrow R_1 + R_3 and R2←R2+R3R_2 \leftarrow R_2 + R_3, then R3←12R3R_3 \leftarrow \tfrac12 R_3:

∼[1008310101041002731]∼[10083101010410017/23/21/2]\sim \left[\begin{array}{ccc|ccc} 1 & 0 & 0 & 8 & 3 & 1 \\ 0 & 1 & 0 & 10 & 4 & 1 \\ 0 & 0 & 2 & 7 & 3 & 1 \end{array}\right] \sim \left[\begin{array}{ccc|ccc} 1 & 0 & 0 & 8 & 3 & 1 \\ 0 & 1 & 0 & 10 & 4 & 1 \\ 0 & 0 & 1 & 7/2 & 3/2 & 1/2 \end{array}\right]A−1=[83110417/23/21/2]A^{-1} = \begin{bmatrix} 8 & 3 & 1 \\ 10 & 4 & 1 \\ 7/2 & 3/2 & 1/2 \end{bmatrix}

Check row 1 of AA against the columns of A−1A^{-1}: 1(8)+0−2(7/2)=11(8) + 0 - 2(7/2) = 1,   1(3)+0−2(3/2)=0\;1(3) + 0 - 2(3/2) = 0,   1(1)+0−2(1/2)=0\;1(1) + 0 - 2(1/2) = 0 ✓.

Part 2 (Exercise 32). R2←R2−4R1R_2 \leftarrow R_2 - 4R_1 and R3←R3+2R1R_3 \leftarrow R_3 + 2R_1:

[1−211004−73010−26−4001]∼[1−2110001−1−41002−2201]\left[\begin{array}{ccc|ccc} 1 & -2 & 1 & 1 & 0 & 0 \\ 4 & -7 & 3 & 0 & 1 & 0 \\ -2 & 6 & -4 & 0 & 0 & 1 \end{array}\right] \sim \left[\begin{array}{ccc|ccc} 1 & -2 & 1 & 1 & 0 & 0 \\ 0 & 1 & -1 & -4 & 1 & 0 \\ 0 & 2 & -2 & 2 & 0 & 1 \end{array}\right]

R3←R3−2R2R_3 \leftarrow R_3 - 2R_2:

∼[1−2110001−1−41000010−21]\sim \left[\begin{array}{ccc|ccc} 1 & -2 & 1 & 1 & 0 & 0 \\ 0 & 1 & -1 & -4 & 1 & 0 \\ 0 & 0 & 0 & 10 & -2 & 1 \end{array}\right]

The left block has a row of zeros: only 2 pivots, so AA is not row equivalent to I3I_3. By Theorem 7, AA is not invertible. Stop here; there is no need to continue.

Exercise 5

(Adapted from Lay §2.2, Exercise 35)

Let A=[−2−7−9256134]A = \begin{bmatrix} -2 & -7 & -9 \\ 2 & 5 & 6 \\ 1 & 3 & 4 \end{bmatrix}. Find the third column of A−1A^{-1} without computing the other columns.

Show solution ↓
Solution

Column 3 of A−1A^{-1} is the solution of Ax=e3A\mathbf{x} = \mathbf{e}_3, so row reduce [Ae3]\left[\begin{array}{c|c} A & \mathbf{e}_3 \end{array}\right].

R1↔R3R_1 \leftrightarrow R_3 (to get a convenient pivot 1):

[−2−7−9025601341]∼[13412560−2−7−90]\left[\begin{array}{ccc|c} -2 & -7 & -9 & 0 \\ 2 & 5 & 6 & 0 \\ 1 & 3 & 4 & 1 \end{array}\right] \sim \left[\begin{array}{ccc|c} 1 & 3 & 4 & 1 \\ 2 & 5 & 6 & 0 \\ -2 & -7 & -9 & 0 \end{array}\right]

R2←R2−2R1R_2 \leftarrow R_2 - 2R_1 and R3←R3+2R1R_3 \leftarrow R_3 + 2R_1, then R3←R3−R2R_3 \leftarrow R_3 - R_2:

∼[13410−1−2−20−1−12]∼[13410−1−2−20014]\sim \left[\begin{array}{ccc|c} 1 & 3 & 4 & 1 \\ 0 & -1 & -2 & -2 \\ 0 & -1 & -1 & 2 \end{array}\right] \sim \left[\begin{array}{ccc|c} 1 & 3 & 4 & 1 \\ 0 & -1 & -2 & -2 \\ 0 & 0 & 1 & 4 \end{array}\right]

R1←R1−4R3R_1 \leftarrow R_1 - 4R_3 and R2←R2+2R3R_2 \leftarrow R_2 + 2R_3:

∼[130−150−1060014]\sim \left[\begin{array}{ccc|c} 1 & 3 & 0 & -15 \\ 0 & -1 & 0 & 6 \\ 0 & 0 & 1 & 4 \end{array}\right]

R2←−R2R_2 \leftarrow -R_2, then R1←R1−3R2R_1 \leftarrow R_1 - 3R_2:

∼[130−15010−60014]∼[1003010−60014]\sim \left[\begin{array}{ccc|c} 1 & 3 & 0 & -15 \\ 0 & 1 & 0 & -6 \\ 0 & 0 & 1 & 4 \end{array}\right] \sim \left[\begin{array}{ccc|c} 1 & 0 & 0 & 3 \\ 0 & 1 & 0 & -6 \\ 0 & 0 & 1 & 4 \end{array}\right]

Answer: the third column of A−1A^{-1} is [3−64]\begin{bmatrix} 3 \\ -6 \\ 4 \end{bmatrix}.

Check Ax=e3A\mathbf{x} = \mathbf{e}_3:

  • Row 1: −2(3)−7(−6)−9(4)=−6+42−36=0-2(3) - 7(-6) - 9(4) = -6 + 42 - 36 = 0 ✓
  • Row 2: 2(3)+5(−6)+6(4)=6−30+24=02(3) + 5(-6) + 6(4) = 6 - 30 + 24 = 0 ✓
  • Row 3: 1(3)+3(−6)+4(4)=3−18+16=11(3) + 3(-6) + 4(4) = 3 - 18 + 16 = 1 ✓