Linear algebra · gauss-jordan

Finding the inverse of a 3×3 matrix

Your matrix has det(A) = −15, so an inverse exists. This page walks every row operation on the augmented matrix [A | I], shows the reasoning behind each one, then verifies the answer by hand and by a second independent method. Work through it once and the method becomes mechanical.

−15
det(A)
6
row ops
6×3
augmented
Yes
invertible
01 / The problem

What you were asked to do

given
103
−41−3
−22−3

Find A⁻¹ by row reducing the augmented matrix [A | I]. The question gives two choices: give the inverse, or declare that no inverse exists.

first move

Augment A with the 3×3 identity. The identity on the right is not decoration — it is a ledger. Every row operation you perform gets applied to both halves simultaneously, so the right half quietly records the sequence of operations that undoes A.

[A | I] =
103 100
−41−3 010
−22−3 001
quick existence check

Before you touch a single row, spend ten seconds on the determinant. For a square matrix, invertible ⇔ det(A) ≠ 0. Computing det(A) with the 3×3 formula is faster than the row reduction, and it tells you which multiple-choice option to pick.

det(A) = a(ei − fh) − b(di − fg) + c(dh − eg) a = 1 b = 0 c = 3 d = −4 e = 1 f = −3 g = −2 h = 2 i = −3 det(A) = 1·((1)(−3) − (−3)(2)) − 0·((−4)(−3) − (−3)(−2)) + 3·((−4)(2) − (1)(−2)) = 1·(−3 + 6) + 3·(−8 + 2) = 3 − 18 = −15 det(A) = −15 ≠ 0 → A⁻¹ exists → choose A
02 / Row reduction

Six moves from [A | I] to [I | A⁻¹]

Start
LEFT: A → I   ·   RIGHT: I → A⁻¹
how to read the moves
GO DOWN, THEN GO UP
Eliminate below each pivot first (R₂, R₃), then above it (R₂, R₁). That order guarantees you never destroy a pivot you already made.
ELEMENTARY OPERATIONS ARE THE ONLY LEGAL MOVES
Swap two rows · scale a row by a nonzero constant · add a multiple of one row to another. Every one of these is reversible, which is exactly why the inverse can be recovered.
THE ROWS THAT CHANGE
In each frame the green entries are the ones the operation touched, and the pink entry is the pivot that was just created or protected. Trace it and the whole algorithm is obvious.
03 / The theory

Why row reducing [A | I] actually works

definition

What an inverse is

A⁻¹ is the matrix B that undoes A. Formally, A⁻¹ satisfies both A·A⁻¹ = I and A⁻¹·A = I — the identity matrix, which is the multiplicative identity: I·M = M·I = M for every matrix M.

Read it as an equation system: multiplying by A⁻¹ is the same as dividing by A. That is why every linear system Ax = b can be solved as x = A⁻¹b.

the core argument

Row operations = left-multiplying by a matrix

Every elementary row operation is itself a matrix multiplication. Swapping rows, scaling a row, adding a multiple of one row to another — each one can be written as left-multiplying the whole augmented matrix by some matrix E.

Apply operations E₁, E₂, …, Eₖ in order: Eₖ···E₂E₁ · [A | I] = [ Eₖ···E₂E₁·A | Eₖ···E₂E₁ ] The right block equals Eₖ···E₂E₁ — the product of all the row-operation matrices. So the identity recorded on the right is literally the machine that undoes A.

Now suppose the reduction reaches [I | E]. Then E·A = I, so E is a left inverse of A; and because A is square, that makes E = A⁻¹. Whatever lands on the right when the left side becomes I is the inverse. That is the whole method in one sentence.

when it fails

When does an inverse exist?

det(A) ≠ 0
The determinant is zero exactly when A collapses to a lower rank — rows or columns that are linearly dependent. Then the system Ax = b has either no solution or infinitely many.
rank(A) = n
A square matrix is invertible exactly when its rows (equivalently its columns) are linearly independent and span the whole space.
row reduction reaches I
If at any point a row of zeros appears in the left half, det(A) = 0 and no inverse exists — that is exactly when the "does not have an inverse" option is the correct choice.
Why det(A) = −15 already told you the answer

The determinant of the triangular matrix you end up with (before scaling the last row) is the product of its diagonal: 1 · 1 · (−15) = −15. That matches det(A) — and it is not a coincidence. Row addition does not change the determinant, so the determinant of the reduced form equals the determinant of the original.

det(A) = product of the diagonal of the reduced triangular matrix = 1 · 1 · (−15) = −15 ✓ consistent

One caveat worth knowing: scaling a row multiplies the determinant. Row additions preserve it; row swaps negate it; scaling by k multiplies it by k. So if you scale the last row by −1/15, the reduced matrix has determinant 1 while A has determinant −15.

Geometric picture: what the inverse actually does

Think of A as a linear map on R³ — a stretch, a shear, a rotation, maybe a flip. A⁻¹ is the map that undoes all of it: it stretches by the reciprocal amount along each direction and undoes every shear.

A⁻¹ · A = I (undo A, then you are back where you started) A · A⁻¹ = I (apply A⁻¹ first, then A — also identity)

This is why invertibility is really about information loss. If A squashes the cube flat along one axis (det = 0), no map can un-squash it — the data along that axis is gone, and no inverse exists. A det of −15 means A is a volume-expanding, orientation-flipping map: it multiplies volumes by 15 and reverses handedness. A⁻¹ multiplies volumes by 1/15 and flips handedness back.

Where this shows up in chemical engineering

In the classroom you compute A⁻¹ because the exercise asks for it. On the job, you almost never do — and knowing why is part of understanding the math:

SOLVE, DON'T INVERT
Real systems are huge. Computing A⁻¹ explicitly costs O(n³) and fills the matrix with noise; solving Ax = b via LU or QR factorization is faster and far more numerically stable. A⁻¹ is a theoretical object, not a computational tool.
CONDITIONING IS THE REAL WARNING
A large determinant does not mean a well-behaved matrix. A⁻¹ with huge entries means the system is ill-conditioned — small measurement errors get amplified. Watch ‖A⁻¹‖, not just det(A).
CONCRETE PLACES YOU'LL MEET IT
Material and energy balances (solve the linear system for the unknown stream rates), least-squares calibration of a reactor model, Newton–Raphson and the Jacobian in process simulation, and process-sensitivity analysis — "if I change this temperature by 1 °C, how much does my yield move?" is literally a row of A⁻¹.
04 / Verification

Check the answer: A · A⁻¹ = I

never skip this

Row reduction can fail silently if you slip on an arithmetic sign. Multiplying A by your candidate inverse and confirming you get the identity takes 60 seconds and catches every sign error. Nine dot products, each must equal 1, 0, or 0.

103
−41−3
−22−3
·
−1/5−2/51/5
2/5−1/53/5
2/52/15−1/15
=
100
010
001

Every entry lands on the identity — no rounding, no decimals. That is the signature of a correct inverse: the products cancel exactly.

05 / Second method

Cross-check with the adjugate formula

adjugate method

A⁻¹ = adj(A) / det(A), where adj(A) is the transpose of the cofactor matrix. Slower by hand than row reduction — but it is the fastest way to catch an arithmetic slip, because it uses completely different arithmetic. If both routes give the same matrix, the answer is right.

COFACTOR MATRIX C
C =
3−6−6
63−2
−3−91
adj(A) = Cᵀ =
36−3
−63−9
−6−21
A⁻¹ = adj(A) / (−15) =
−1/5−2/51/5
2/5−1/53/5
2/52/15−1/15

Identical to the row-reduction result. Two independent methods, same matrix — that is how you know.

06 / Your turn

Type it in — get it back

Enter the nine entries of A⁻¹. Fractions are fine (2/15), decimals are fine (0.1333), and a leading minus works either way. It checks each cell the moment you leave the field.

07 / Exam tips

Patterns worth memorizing

🎯
Compute the determinant first. It is 15 seconds of work and it tells you whether the answer is a matrix or "no inverse". It also lets you check the final fractions: the adjugate entries divided by det must land on the right numbers.
✂️
Work on the augmented matrix only. Never try to "solve for" the inverse by manipulating A alone. Every operation must hit all six entries of its row.
🪜
Clear below first, then above. Once R₁ has a pivot, clean R₂ and R₃ under it; once R₂ has a pivot, clean R₃ under it; then sweep upward. It is the only order that never breaks a pivot you already built.
🔢
Fractions are normal. An inverse of an integer matrix almost always has fractions. If you get all integers, suspect a mistake. Here every entry carries a denominator of 5 or 15 because det = −15.
🧊
A row of zeros on the left = no inverse. If you hit a zero row in the A-block, stop and pick the "does not have an inverse" option.
✅
Verify before you submit. Multiply A by your answer and check you get the identity. It takes a minute and catches every sign slip — the single most common way to lose a point on this question type.
🚀
Second method for free. If row reduction gives you an ugly answer, compute det(A) and adj(A) instead. Often it is faster, and it doubles as a check on your reduction.
final answer for the homework box

Select option A and enter:

−1/5−2/51/5
2/5−1/53/5
2/52/15−1/15

If the homework box wants plain ASCII, type: -1/5 -2/5 1/5 · 2/5 -1/5 3/5 · 2/5 2/15 -1/15