Consider the reduced row echelon matrix
\begin{equation*}
\left[
\begin{array}{rrr|r}
1 \amp 0 \amp 2 \amp -1 \\
0 \amp 1 \amp 1 \amp 2 \\
\end{array}
\right]
\end{equation*}
and its corresponding linear system as
\begin{equation*}
\begin{alignedat}{4}
x \amp \amp \amp {}+{} \amp 2z \amp {}={} \amp -1 \\
\amp \amp y \amp {}+{} \amp z \amp {}={} \amp 2. \\
\end{alignedat}
\end{equation*}
Let’s rewrite the equations as
\begin{equation*}
\begin{alignedat}{2}
x \amp {}={} \amp -1 -2z\\
y \amp {}={} \amp 2-z. \\
\end{alignedat}
\end{equation*}
From this description, it is clear that we obtain a solution for any value of the variable \(z\text{.}\) For instance, if \(z=2\text{,}\) then \(x =
-5\) and \(y=0\) so that \((x,y,z) = (-5,0,2)\) is a solution. Similarly, if \(z=0\text{,}\) we see that \((x,y,z) = (-1,2,0)\) is also a solution.
We will call this description of the solution space, in which the basic variables are written in terms of the free variables, a parametric description of the solution space.