Use Cramer's Rule To Compute The Solutions Of The System

What Is Cramer's Rule?

Here's a scenario you might recognize. That's why you've got two or three equations, each with two or three unknowns, and you need to find the values that satisfy all of them at the same time. Even so, substitution works, elimination works, but sometimes there's a cleaner path — one that turns the whole problem into a series of determinants. That path is Cramer's Rule.

Cramer's Rule is a method in linear algebra that gives you an explicit formula for the solution of a system of linear equations, provided the system has the same number of equations as unknowns and the coefficient matrix is invertible (meaning its determinant is not zero). But instead of grinding through row operations or back-substitution, you compute a handful of determinants and divide. It's elegant, it's systematic, and — when the numbers cooperate — it's surprisingly fast.

The rule is named after Gabriel Cramer, a Swiss mathematician who published it in 1750 in his book Introduction à l'Analyse des lignes Courbes Algébriques*. He wasn't the first to work with determinants, but he was the one who laid out this particular connection between determinants and solutions of linear systems in a way that stuck.

When Does Cramer's Rule Apply?

Not every system qualifies. You need three things in place. That said, first, the number of equations must equal the number of variables — so it works for a 2×2 system, a 3×3 system, and so on. Second, the determinant of the coefficient matrix must be nonzero. If it's zero, the system either has no solution or infinitely many solutions, and Cramer's Rule simply can't handle that case. Third, the equations need to be linear — no exponents, no products of variables, just straight lines (or planes, or hyperplanes, depending on the dimension).

Why Cramer's Rule Matters

You might wonder why anyone would reach for Cramer's Rule when Gaussian elimination exists and scales better for large systems. Which means that's a fair question, and honestly, for systems with dozens or hundreds of variables, elimination (or matrix factorization methods) is the practical choice. But Cramer's Rule has a different kind of value.

For small systems — 2×2 or 3×3 — it gives you a direct, formulaic answer. You write down one determinant for the denominator and one determinant per variable for the numerator. There's a transparency to it that other methods lack. You can see exactly how each coefficient and constant contributes to the final answer.

It also shows up in theoretical work. Even so, when you need to prove something about the structure of solutions, or when you're working with symbolic expressions rather than numbers, Cramer's Rule gives you a clean algebraic form. Engineers, physicists, and economists sometimes use it in derivations where the closed-form expression matters more than computational speed.

And let's be real — if you're taking a linear algebra course, there's a good chance your professor will ask you to solve a system using Cramer's Rule at some point. Knowing it isn't just about the math; it's about showing up prepared.

How to Use Cramer's Rule to Solve a System

Here's where the rubber meets the road. Let's walk through the process step by step, starting with a simple 2×2 system and then scaling up.

Step 1: Write the System in Standard Form

Before you touch a determinant, make sure every equation is in the form where variables are on the left and constants are on the right. For a 2×2 system, that looks like this:

a₁x + b₁y = c₁ a₂x + b₂y = c₂

The coefficients a₁, b₁, a₂, b₂ form the coefficient matrix. The constants c₁ and c₂ form the right-hand side vector.

Step 2: Compute the Determinant of the Coefficient Matrix

This is the denominator in every formula you'll write. For a 2×2 system, the determinant D is:

D = a₁b₂ − a₂b₁

For a 3×3 system, you'd use the standard 3×3 determinant expansion — cofactor expansion along any row or column, or the rule of Sarrus if you prefer. The key point is that D must not be zero. If D = 0, stop here. Cramer's Rule doesn't apply.

Step 3: Replace Columns to Build Numerator Determinants

To solve for x, take the coefficient matrix and replace its first column (the column of x-coefficients) with the constants column. Worth adding: compute that determinant — call it Dₓ. To solve for y, replace the second column with the constants column and call that determinant Dᵧ.

For a 3×3 system, you'd do the same thing: Dₓ replaces the x-column, D_y replaces the y-column, and D_z replaces the z-column.

Step 4: Divide to Get Each Variable

x = Dₓ / D y = Dᵧ / D

That's it. Each variable is a ratio of two determinants. For a 3×3 system, you'd have three ratios — one for each variable.

A Concrete 2×2 Example

Consider this system:

3x + 2y = 16 x − y = 1

The coefficient matrix is [[3, 2], [1, −1]]. Its determinant D = (3)(−1) − (1)(2) = −3 − 2 = −5.

For Dₓ, replace the first column with [16, 1]: Dₓ = (16)(−1) − (1)(2) = −16 − 2 = −18

For Dᵧ, replace the second column with [16, 1]: Dᵧ = (3)(1) − (1)(16) = 3 − 16 = −13

If you found this helpful, you might also enjoy find the circumference of the circle use 3.14 for π or which is the major product of the following reaction.

So x = −18 / −5 = 18/5, and y = −13 / −5 = 13/5.

You can verify by plugging back into the original equations. Which means 3(18/5) + 2(13/5) = 54/5 + 26/5 = 80/5 = 16. And 18/5 − 13/5 = 5/5 = 1. Checks out.

Extending to 3×3 Systems

The process is identical in principle for a 3×3 system, but the determinant calculations get a bit more involved. You have a 3×3 coefficient matrix, and you compute D, then Dₓ, D_y, and D_z — each time replacing one column with the constants vector.

Let's say you have

Let’s Say You Have a 3 × 3 System

To see Cramer's Rule in action with three variables, consider the following system:

[ \begin{cases} 2x + y - z = 7 \ x - 3y + 2z = -4 \ 4x + y + z = 10 \end{cases} ]

The coefficient matrix (A) and the constant vector (\mathbf{b}) are

[ A=\begin{pmatrix} 2 & 1 & -1\ 1 & -3 & 2\ 4 & 1 & 1 \end{pmatrix}, \qquad \mathbf{b}= \begin{pmatrix}7\-4\10\end{pmatrix}. ]

1. Compute the Main Determinant (D)

Using cofactor expansion along the first row (or any convenient row/column), we get

[ \begin{aligned} D &= 2\begin{vmatrix}-3 & 2\ 1 & 1\end{vmatrix} - 1\begin{vmatrix}1 & 2\ 4 & 1\end{vmatrix} + (-1)\begin{vmatrix}1 & -3\ 4 & 1\end{vmatrix} \[4pt] &= 2\big((-3)(1)-2(1)\big) - \big(1\cdot1-2\cdot4\big) - \big(1\cdot1-(-3)\cdot4\big) \[4pt] &= 2(-3-2) - (1-8) - (1+12) \[4pt] &= 2(-5) - (-7) - 13 \[4pt] &= -10 + 7 - 13 = -16. \end{aligned} ]

Since (D\neq0), the system has a unique solution and Cramer’s Rule can be applied.

2. Form the Numerator Determinants

For (x): Replace the first column of (A) with (\mathbf{b}):

[ A_x= \begin{pmatrix} 7 & 1 & -1\ -4 & -3 & 2\ 10 & 1 & 1 \end{pmatrix}. ]

[ \begin{aligned} D_x &= 7\begin{vmatrix}-3 & 2\ 1 & 1\end{vmatrix} - 1\begin{vmatrix}-4 & 2\ 10 & 1\end{vmatrix} + (-1)\begin{vmatrix}-4 & -3\ 10 & 1\end{vmatrix} \[4pt] &= 7(-3-2) - \big((-4)(1)-2\cdot10\big) - \big((-4)(1)-(-3)\cdot10\big) \[4pt] &= 7(-5) - (-4-20) - (-4+30) \[4pt] &= -35 - (-24) - 26 \[4pt] &= -35 + 24 - 26 = -37. \end{aligned} ]

For (y): Replace the second column:

[ A_y= \begin{pmatrix} 2 & 7 & -1\ 1 & -4 & 2\ 4 & 10 & 1 \end{pmatrix}, \qquad \begin{aligned} D_y &= 2\begin{vmatrix}-4 & 2\ 10 & 1\end{vmatrix} - 7\begin{vmatrix}1 & 2\ 4 & 1\end{vmatrix} + (-1)\begin{vmatrix}1 & -4\ 4 & 10\end{vmatrix} \[4pt] &= 2((-4)(1)-2\cdot10) - 7(1-8) - (1\cdot10-(-4)\cdot4) \[4pt] &= 2(-4-20) - 7(-7) - (10+16) \[4pt] &= 2(-24) + 49 - 26 \[4pt] &= -48 + 49 - 26 = -25. \end{aligned} ]

For (z): Replace the third column:

[ A_z= \begin{pmatrix} 2 &

[ A_z= \begin{pmatrix} 2 & 1 & 7\ 1 & -3 & -4\ 4 & 1 & 10 \end{pmatrix}, \qquad \begin{aligned} D_z &= 2\begin{vmatrix}-3 & -4\ 1 & 10\end{vmatrix} - 1\begin{vmatrix}1 & -4\ 4 & 10\end{vmatrix} + 7\begin{vmatrix}1 & -3\ 4 & 1\end{vmatrix} \[4pt] &= 2\big((-3)(10)-(-4)(1)\big) - \big(1\cdot10-(-4)\cdot4\big) + 7\big(1\cdot1-(-3)\cdot4\big) \[4pt] &= 2(-30+4) - (10+16) + 7(1+12) \[4pt] &= 2(-26) - 26 + 7\cdot13 \[4pt] &= -52 - 26 + 91 = 13 . \end{aligned} ]

Now apply Cramer’s Rule:

[ x=\frac{D_x}{D}=\frac{-37}{-16}=\frac{37}{16},\qquad y=\frac{D_y}{D}=\frac{-25}{-16}=\frac{25}{16},\qquad z=\frac{D_z}{D}=\frac{13}{-16}=-\frac{13}{16}. ]

Thus the unique solution of the system is

[ \boxed{\left(\frac{37}{16},;\frac{25}{16},;-\frac{13}{16}\right)}. ]

A quick substitution confirms the result:

[ \begin{aligned} 2x+y-z &= 2!\left(\frac{37}{16}\right)+\frac{25}{16}-!\left(-\frac{13}{16}\right) =\frac{74+25+13}{16}= \frac{112}{16}=7,\[2pt] x-3y+2z &= \frac{37}{16}-3!\left(\frac{25}{16}\right)+2!\left(-\frac{13}{16}\right) =\frac{37-75-26}{16}= \frac{-64}{16}=-4,\[2pt] 4x+y+z &= 4!

Continuing from where the verification left off:

[ 4x+y+z = 4!\left(\frac{37}{16}\right)+\frac{25}{16}-\frac{13}{16} = \frac{148+25-13}{16} = \frac{160}{16} = 10, ]

which matches the right‑hand side of the third equation. All three substitutions satisfy the original system, confirming that the vector

[ \left(\frac{37}{16},;\frac{25}{16},;-\frac{13}{16}\right) ]

is indeed the unique solution.

Boiling it down, because the determinant of the coefficient matrix is non‑zero ((D=-16\neq0)), Cramer’s Rule guarantees a unique solution, and the explicit computation yields the values above. This demonstrates both the correctness of the method and the consistency of the system.