Solving Systems by Substitution and Elimination
A solution of a system is an ordered pair (or triple) that satisfies every equation at once. Two main hand methods:
- [leftmargin=5mm,itemsep=1pt]
- Substitution: solve one equation for a variable, then plug into the other. Best when a variable already has coefficient , or for nonlinear systems.
- Elimination: add multiples of the equations to cancel a variable. Best for linear systems with awkward coefficients.
A linear system may have one solution (lines cross), no solution (parallel lines; you reach a false statement like ), or infinitely many (same line; you reach a true statement like ). A nonlinear system (line and circle, two conics, …) can have or more solutions.
Solve
Multiply eq. 1 by and eq. 2 by to cancel :
Add: . Back-substitute into eq. 1: .
Check: . ✓ Solution: .
Solve
Substitute : .
Factor: or , giving or .
Solutions: and . Check : . ✓
Tip: After finding a variable, always back-substitute into the original equation and verify both. For nonlinear systems, each must be paired with the from the equation you substituted into.
Three-Variable Linear Systems: Gaussian Elimination
Write the system as an augmented matrix and use elementary row operations to reach row-echelon form (leading 's stepping to the right, zeros below), then back-substitute.
- [leftmargin=5mm,itemsep=1pt]
- Swap two rows: .
- Scale a row: ().
- Add a multiple of one row to another: .
Reduced row-echelon form additionally clears entries above each leading (Gauss--Jordan).
Solve
Augmented matrix, then eliminate below the first pivot with and :
Swap , then :
Back-substitute: ; ; .
Check in eq. 2: . ✓ Solution: .
Tip: A row of the form with means no solution. A row means a free variable and infinitely many solutions.
Matrices and Their Operations
For matrices of the same size, add/subtract entrywise; scalar multiply multiplies every entry.
Matrix product is defined only when is and is (inner dimensions match); the result is . Entry of is the dot product of row of with column of . In general .
Let , .
.
.
Tip: Check dimensions before multiplying: . If the inner numbers differ, the product does not exist.
The Inverse of a Matrix
satisfies . For a matrix
provided . If the matrix is singular (no inverse). To solve , use .
Solve using an inverse.
Coefficient matrix , .
Then .
Check: and . ✓ Solution .
Tip: Verify an inverse cheaply by confirming , or at least that .
Determinants and Cramer's Rule
: .
(expansion along row 1):
Cramer's Rule: for with , each variable is , where replaces column of with .
Solve for by Cramer.
.
.
So (matching the Gaussian result above).
Tip: Remember the alternating sign pattern across the top row. If , Cramer's Rule fails---the system has either no solution or infinitely many.
Partial Fraction Decomposition
For a proper fraction (numerator degree denominator degree), factor the denominator and write one term per factor:
- [leftmargin=5mm,itemsep=1pt]
- Distinct linear : term .
- Repeated linear : terms .
- Irreducible quadratic : term .
If improper, do polynomial division first. Solve for the constants by clearing denominators.
Decompose .
Factor: . Write .
Clear: . Let : . Let : .
Check: . ✓
Tip: Substituting the roots of the factors instantly isolates constants. For repeated or quadratic factors, also compare coefficients of like powers.
Systems of Inequalities and Linear Programming
Graph each inequality and shade the overlap---the feasible region. In linear programming you maximize or minimize a linear objective over that region. The optimum (if it exists) always occurs at a corner (vertex), so evaluate at every vertex and pick the best.
Maximize subject to .
Vertices: . The lines cross at .
Maximum at .
Tip: Always list the vertices explicitly---including where two boundary lines intersect---then test each. A bounded feasible region guarantees both a max and a min exist.
Going Deeper: Advanced Systems & Matrices
Every square matrix satisfies its own characteristic equation. For a matrix with trace and determinant , the Cayley--Hamilton theorem gives
Two powerful consequences:
- [leftmargin=5mm,itemsep=1pt]
- Reduce high powers: (and hence ) is always a linear combination . So for scalars found by a recurrence, never a full multiplication.
- Instant inverse: rearrange to , so if , @@BLOCK1@@
The scalar recurrence itself mirrors the matrix: if then , so and .
Let , so and .
Inverse. , matching the adjugate formula.
Power via recurrence. From :
So --- no repeated matrix multiplication needed.
Some determinants have closed forms worth memorizing:
- [leftmargin=5mm,itemsep=1pt]
- Vandermonde: , i.e. the product of all differences ().
- Identity plus all-ones (): if is the matrix of all 's, then .
- Block triangular: (the off-diagonal block does not matter).
Singular () exactly when rows/columns are linearly dependent: one row is a scalar multiple of another, a row is a sum of others, or a row (or column) is all zeros. A singular matrix has no inverse, and then has either no solution or infinitely many.
Vandermonde with :
, size : formula gives . Directly, , and expanding confirms . ✓
Block triangular: (product of the two diagonal blocks, which are themselves triangular).
When both equations are conics, eliminate the squared terms rather than a single variable. Add or subtract multiples so that (or ) cancels, leaving a simpler relation; substitute back. Two conics can meet in up to four points, so expect several ordered pairs and pair each recovered value with the correct sign from the equation you used.
Solve
Subtract eq. 2 from eq. 1 to kill : .
Substitute into eq. 2: .
Solutions (four): . Check : and . ✓
Two cases beyond the distinct-linear setup demand extra terms:
- [leftmargin=5mm,itemsep=1pt]
- Repeated linear : include one term per power, .
- Irreducible quadratic (no real roots): use a linear numerator .
Root-substitution alone will not find every constant here, so clear denominators and match coefficients of like powers of , giving a small linear system for the unknowns.
Decompose .
Clear denominators: .
Root shortcut, : .
Match coefficients: ; constant; check . ✓
The corner principle assumes a bounded feasible region. Two degenerate cases:
- [leftmargin=5mm,itemsep=1pt]
- Unbounded region: if the objective can grow without limit along an open direction, there is no maximum (though a minimum may still exist at a vertex, or vice versa). Always check whether the region is closed in the direction you are optimizing.
- Edge of optima: when the objective line is parallel to a boundary edge, every point on that edge is optimal --- infinitely many solutions, all giving the same . You detect this when two adjacent vertices tie for the best value.
Consider (region opens outward to the upper right).
Maximize : moving farther from the origin increases without bound, so no maximum exists --- the problem is unbounded.
Minimize over the same region: the constraint forces , and every point of the edge from to gives . So the minimum is , attained along an entire edge (an edge of optima), not a single corner.
Formulas, Proofs & Tips
What it means. The determinant measures how a matrix scales area; when it is the system has no unique solution.
Example. .
Why it works. Solving by elimination gives , and the numerator is the determinant of with the -column replaced by the constants. A zero determinant means the two lines are parallel or identical.
Tip. Check before dividing — that is precisely the test for a unique solution.