Linear systems and consistency
Before you solve a system of linear equations, ask a sharper question first: does a solution even exist, and if it does, is there exactly one or infinitely many? Guessing by trial-and-error, or blindly running an elimination algorithm and hoping for the best, wastes time and hides what is actually going on. The Rouché–Capelli theorem answers this cleanly, using nothing more than the rank of two matrices — the coefficient matrix and the augmented matrix . Comparing their ranks tells you, before you do a single back-substitution, whether the system is inconsistent, uniquely solvable, or has an infinite family of solutions. This article builds that criterion from matrix algebra and rank, works it through two contrasting examples — one consistent, one not — and sets up the numerical solvers (Gauss elimination, Gauss–Jordan, Gauss–Seidel) that come next: you must know whether a system is solvable before you spend effort finding what the solution is.
A linear system is consistent exactly when , and that common rank against the number of unknowns then decides whether the solution is unique () or an -parameter infinite family ().
The core idea
Setting up the problem. A system of linear equations in unknowns is written in matrix form as , where is an coefficient matrix, is the unknown vector, and collects the right-hand sides. The augmented matrix is formed by tacking on as one extra column; row-reducing this single matrix reveals everything about both and the system at once.
The consistency criterion. The Rouché–Capelli theorem compares with :
where: is the number of linearly independent rows (equivalently columns) of , obtained by counting nonzero rows after row-reducing to echelon form, and is the same count for the augmented matrix. Appending a column can only keep the rank the same or increase it by exactly one; it increases by one precisely when is not expressible using the row structure already present in , which shows up in row-reduced form as a row reading for some nonzero constant — an equation no values of can satisfy.
Reading off the solution count. Once consistency holds, let . If (the number of unknowns), every unknown is pinned down by a leading entry in the reduced system, so the solution is unique. If , there are unknowns left free to take any value, and the solution set is an -parameter family — a line, plane, or higher-dimensional flat, depending on .
Square systems as a shortcut. When , the determinant gives a fast route to the same conclusion. If , then automatically, the system is consistent for every , and the unique solution is . If , the determinant alone cannot tell you whether the system is inconsistent or has infinitely many solutions — you must fall back on comparing the two ranks directly.
Non-square systems. The picture below makes concrete why shape matters: rank tells the whole story, but the number of equations versus unknowns sets a ceiling on what is even possible. When (fewer equations than unknowns), automatically, so a consistent system can never have a unique solution — it is always an infinite family. When (an overdetermined system), consistency becomes the tight constraint: the extra equations generically conflict, and the system is inconsistent unless the extra rows are combinations of the others.
The goal is to establish, for a system of three equations in three unknowns given as plain algebra, whether the Rouché–Capelli ranks agree, and if they do, to recover the unique numerical solution by row reduction followed by back-substitution — this is the exact procedure every subsequent numerical method (Gauss elimination in particular) automates.
Step 1. Write the system in matrix form and build the augmented matrix by appending as a fourth column beside the three columns of . This single object carries all the information needed for both the rank test and the eventual solution.
Step 2. Row-reduce the augmented matrix using elementary row operations — replacing a row by itself minus a multiple of another row, which never changes the solution set. Track the coefficient columns and the augmented column together throughout; do not row-reduce alone and separately, since consistency depends on how the two are entangled.
Step 3. Count the nonzero rows in the reduced coefficient part to get , and count the nonzero rows of the full reduced augmented matrix to get . If a row has reduced to all zeros in the coefficient columns but a nonzero entry in the last column, that row reads with — an immediate certificate of inconsistency, and the ranks will differ by exactly one.
Step 4. If the ranks agree and equal (three unknowns here), the system is consistent with a unique solution. Recover it by back-substitution: solve the last (simplest) equation for the last variable, substitute into the equation above it to get the next variable, and continue upward until all unknowns are found. Verified against Worked example 1 below: the reduced system gives , then , then , matching .
Worked example
Determine consistency and solve, if possible:
Step 1. Write and , and form the augmented matrix .
Step 2. Eliminate below the first pivot: gives , and gives , leaving
Step 3. Eliminate below the second pivot: gives , so
Both and now have three nonzero rows, so : the system is consistent with a unique solution.
Step 4. Back-substitute. Row 3 gives . Row 2 gives . Row 1 gives . The unique solution is .
Competitive stretch
Now test a superficially similar system for consistency:
Row reduction gives : , and : . The row already signals trouble: it reads , impossible for any . Here (the coefficient rows reduce to two independent nonzero rows) but , so : the system is inconsistent — no solution exists, even though it looks just as solvable as the first at a glance.
Common traps
- Checking only for a non-square or rank-deficient square system and stopping there — you must compare with , not just inspect alone.
- Forgetting that a homogeneous system is always consistent (since always works); the real question is whether a non-trivial solution exists, decided by in the square case or by rank versus in general.
- Miscounting the free parameters in an infinite family: it is where , not — the number of equations does not directly enter the count.
- Arithmetic slips while row-reducing that create a spurious row (or erase a genuine one) — always double-check a claimed inconsistency by confirming the row truly reads , not a miscopied coefficient.
Ask the tutor: Not sure whether your system has one, none, or infinitely many solutions before you even start solving — ask the AI tutor to walk through the rank comparison with you.
Why does comparing two ranks tell you everything about a linear system, rather than just plugging numbers in and hoping?
Rank measures how much independent information the equations actually carry. If , every equation you added to the system by appending is compatible with the structure already in — no equation demands something the others forbid. If the ranks differ, appending introduced a genuinely new, unsatisfiable constraint, which always surfaces as a row like with after row reduction. Trial-and-error can miss this entirely; the rank comparison catches it every time, before you invest effort in solving.
If det A = 0 for a square system, does that automatically mean no solution?
No — only tells you the solution is not unique; it does not distinguish between "no solution" and "infinitely many." You must still compare with : if they are equal (both less than ), you get an infinite family; if they differ, the system is inconsistent. Determinant alone is a fast unique-versus-not-unique test, not a full consistency test.
How is the number of free parameters in an infinite solution family actually determined?
It is , where is the number of unknowns and is the common rank once you know the system is consistent. Each of the independent equations pins down one unknown in terms of the rest; the remaining unknowns can be assigned freely (as parameters), and every other unknown is then expressed in terms of them. It is easy to wrongly use (equations minus unknowns) instead — but redundant or dependent equations do not each remove a degree of freedom, only independent ones do, which is exactly what rank counts.