How To Solve System Of Equations With Three Variables

How to solve systemof equations with three variables is a fundamental skill in algebra that appears in physics, engineering, economics, and many other fields. When you have three unknowns—usually denoted (x), (y), and (z)—you need three independent equations to find a unique solution (if one exists). The goal is to find the ordered triple ((x, y, z)) that satisfies all equations simultaneously. Below you will find a detailed, step‑by‑step guide covering the most reliable techniques, a worked example, common mistakes to avoid, and a FAQ section to clarify lingering doubts.

Introduction

Solving a system of three linear equations with three variables means finding values for (x), (y), and (z) that make every equation true at the same time. The system can be written in the general form

[ \begin{cases} a_1x + b_1y + c_1z = d_1\ a_2x + b_2y + c_2z = d_2\ a_3x + b_3y + c_3z = d_3 \end{cases} ]

where the coefficients (a_i, b_i, c_i) and constants (d_i) are real numbers. Depending on the relationships among the equations, the system may have a single solution, infinitely many solutions, or no solution at all. The methods described below work for any consistent system and also reveal when a system is inconsistent or dependent.

Methods Overview

There are four main algebraic techniques that are widely taught: substitution, elimination (also called addition), Gaussian elimination using matrices, and Cramer's rule. Each has its own strengths; substitution is intuitive for simple systems, elimination is efficient for hand calculations, matrix methods generalize well to larger systems, and Cramer's rule provides a direct formula when the coefficient matrix is invertible.

Substitution Method

The substitution method solves one equation for a single variable, then plugs that expression into the other equations, reducing the number of variables step by step.

1. Choose an equation that looks easiest to solve for one variable (often one with a coefficient of 1 or -1).
2. Isolate that variable, expressing it as a function of the other two.
3. Substitute the expression into the remaining two equations, giving you a system of two equations with two unknowns.
4. Repeat the substitution process to solve for the second variable, then back‑substitute to find the third.
5. Check the solution by inserting the triple into all original equations.

Elimination Method

Elimination (or addition) removes a variable by adding or subtracting equations after multiplying them by suitable constants.

1. Align the equations so like terms are in columns.
2. Pick a variable to eliminate first (commonly (x) or (y)).
3. Multiply one or both equations by constants so that the coefficients of the chosen variable are opposites.
4. Add the equations; the selected variable cancels out, leaving a new equation with two variables.
5. Repeat the elimination with another pair of original equations to obtain a second two‑variable equation.
6. Solve the resulting 2×2 system using substitution or elimination again.
7. Back‑substitute to find the third variable.
8. Verify the solution.

Matrix Method (Gaussian Elimination)

Writing the system as an augmented matrix ([A|b]) allows systematic row operations that preserve the solution set.

1. Form the augmented matrix

[ \left[\begin{array}{ccc|c} a_1 & b_1 & c_1 & d_1\ a_2 & b_2 & c_2 & d_2\ a_3 & b_3 & c_3 & d_3 \end{array}\right] ]

2. Use row operations (swap rows, multiply a row by a non‑zero scalar, add a multiple of one row to another) to achieve row‑echelon form: leading entries (pivots) move to the right as you go down, and all entries below each pivot are zero.
3. Continue to reduced row‑echelon form (RREF) if you want the solution directly; each pivot column will contain a 1 and zeros elsewhere.
4. Read off the solution from the final matrix. If a row reads ([0;0;0|k]) with (k\neq0), the system is inconsistent (no solution). If you have fewer pivots than variables, the system has infinitely many solutions (express the free variables in terms of parameters).

Cramer's Rule

Cramer's rule gives an explicit formula using determinants, but it only applies when the coefficient matrix (A) is square and (\det(A)\neq0).

1. Compute the determinant (D = \det(A)).
2. Form matrices (A_x, A_y, A_z) by replacing the respective column of (A) with the constants column (\mathbf{d} = (d_1,d_2,d_3)^T).
3. Calculate (D_x = \det(A_x), D_y = \det(A_y), D_z = \det(A_z)).
4. The solutions are

[ x = \frac{D_x}{D},\quad y = \frac{D_y}{D},\quad z = \frac{D_z}{D}. ]

If (D=0), the system either has no solution or infinitely many; Cramer's rule cannot be used.

Step‑by‑Step Example

Let's solve the following system using the elimination method, then verify with Gaussian elimination.

[ \begin{cases} 2x + 3y - z = 5 \quad\text{(1)}\ 4x - y + 2z = 6 \quad\text{(2)}\

• x + 2y + 3z = 4 \quad\text{(3)} \end{cases} ]

Elimination Phase

Step 1 – Eliminate (x) from equations (2) and (3) using (1).

Multiply (1) by 2 to match the (x)-coefficient in (2):

(2\times(1):; 4x + 6y - 2z = 10)

Subtract (2) from this result:

((4x+6y-2z) - (4x - y + 2z) = 10 - 6)

(\Rightarrow 7y - 4z = 4)

Step 2 – Eliminate (y) from equation (3) using the new equation (7y - 4z = 4).

Multiply (7y - 4z = 4) by 2/7 to get (2y - \frac{8}{7}z = \frac{8}{7}).

Multiply (3) by 2/7 to get (\frac{2}{7}x + \frac{4}{7}y + \frac{6}{7}z = \frac{8}{7}).

Subtract the modified (3) from (3):

(\left(\frac{2}{7}x + \frac{4}{7}y + \frac{6}{7}z\right) - \left(\frac{2}{7}x + \frac{4}{7}y + \frac{6}{7}z\right) = \frac{8}{7} - \frac{8}{7})

(\Rightarrow 0 = 0)

This confirms that we have successfully eliminated y.

Step 3 – Solve for z from the equation 7y - 4z = 4.

Solve for z: (4z = 7y - 4 \Rightarrow z = \frac{7y - 4}{4})

Step 4 – Substitute z into equation (1) to solve for y.

(2x + 3y - \left(\frac{7y - 4}{4}\right) = 5)

Multiply by 4 to eliminate the fraction: (8x + 12y - (7y - 4) = 20)

Simplify: (8x + 12y - 7y + 4 = 20)

(8x + 5y = 16)

Step 5 – Substitute z into equation (2) to solve for x.

(4x - y + 2\left(\frac{7y - 4}{4}\right) = 6)

Multiply by 4 to eliminate the fraction: (16x - 4y + 2(7y - 4) = 24)

Simplify: (16x - 4y + 14y - 8 = 24)

(16x + 10y = 32)

Divide by 2: (8x + 5y = 16)

Notice that we arrive at the same equation as in Step 4. This indicates that we have infinitely many solutions.

Step 6 – Express the solution in terms of a parameter.

Let y = t. Then, from (8x + 5y = 16), we have (8x + 5t = 16), so (8x = 16 - 5t) and (x = \frac{16 - 5t}{8}).

Finally, we can express z in terms of t: (z = \frac{7y - 4}{4} = \frac{7t - 4}{4}).

Therefore, the solution is:

[ x = \frac{16 - 5t}{8}, \quad y = t, \quad z = \frac{7t - 4}{4} ]

Verification using Gaussian Elimination

Let's use the augmented matrix to verify our solution.

[ \left[\begin{array}{ccc|c} 2 & 3 & -1 & 5\ 4 & -1 & 2 & 6\ -1 & 2 & 3 & 4 \end{array}\right] ]

Step 1: Swap Row 1 and Row 3

[ \left[\begin{array}{ccc|c} -1 & 2 & 3 & 4\ 4 & -1 & 2 & 6\ 2 & 3 & -1 & 5 \end{array}\right] ]

Step 2: Multiply Row 1 by -1

[ \left[\begin{array}{ccc|c} 1 & -2 & -3 & -4\ 4 & -1 & 2 & 6\ 2 & 3 & -1 & 5 \end{array}\right] ]

Step 3: Row 2 = Row 2 - 4 * Row 1

[ \left[\begin{array}{ccc|c} 1 & -2 & -3 & -4\ 0 & 9 & 14 & 22\ 2 & 3 & -1 & 5 \end{array}\right] ]

Step 4: Row 3 = Row 3 - 2 * Row 1

[ \left[\begin{array}{ccc|c} 1 & -2 & -3 & -4\ 0 & 9 & 14 & 22\ 0 & 7 & 7 & 13 \end{array}\right] ]

Step 5: Divide Row 2 by 9

[ \left[\begin{array}{ccc|c} 1 & -2 & -3 & -4\ 0 & 1 & \frac{14}{9} & \frac{

22}{9}\ 0 & 7 & 7 & 13 \end{array}\right] ]

Step 6: Row 3 = Row 3 - 7 * Row 2

[ \left[\begin{array}{ccc|c} 1 & -2 & -3 & -4\ 0 & 1 & \frac{14}{9} & \frac{22}{9}\ 0 & 0 & 0 & 0 \end{array}\right] ]

The presence of a zero row in the augmented matrix confirms that the system is dependent and has infinitely many solutions. This aligns with our previous finding that the system of equations has infinitely many solutions.

Conclusion

Through both substitution and Gaussian elimination methods, we have determined that the given system of linear equations possesses infinitely many solutions. The general solution can be expressed in terms of a parameter ( t ), where ( y = t ), ( x = \frac{16 - 5t}{8} ), and ( z = \frac{7t - 4}{4} ). This parameterization allows for a comprehensive understanding of the solution space, illustrating the flexibility and breadth of solutions available within the constraints of the original equations.