x + 2y + 3z = a
x + 3y + 8z = b
x + 2y + 2z = c
factorizing
[ 1 2 3 1 3 8 1 2 2 ] \begin{bmatrix}
1 & 2&3 \\
1 & 3&8\\
1 &2&2
\end{bmatrix} ⎣ ⎡ 1 1 1 2 3 2 3 8 2 ⎦ ⎤ [ x y z ] \begin{bmatrix}
x \\
y \\
z
\end{bmatrix} ⎣ ⎡ x y z ⎦ ⎤ =[ a b c ] \begin{bmatrix}
a \\
b\\
c
\end{bmatrix} ⎣ ⎡ a b c ⎦ ⎤
writing the argument matrix
[ 1 2 3 a 1 3 8 b 1 2 2 c ] \begin{bmatrix}
1 & 2&3&a \\
1 & 3&8&b\\
1&2&2&c
\end{bmatrix} ⎣ ⎡ 1 1 1 2 3 2 3 8 2 a b c ⎦ ⎤
reducing to row echelon
R1 to R1 , R2 -R1 , R3 -R1 we obtain
[ 1 2 3 a 0 1 5 b − a 0 0 − 1 c − a ] \begin{bmatrix}
1 & 2&3&a \\
0 & 1&5&b-a\\
0&0&-1&c-a
\end{bmatrix} ⎣ ⎡ 1 0 0 2 1 0 3 5 − 1 a b − a c − a ⎦ ⎤
replacing the variables to form system of equations
x+2y+3Z=a
y+5z=b-a
-z=c-a
taking back substitution
y+5(a-c)=b-a
y=b-a-5a+5c
x+2(b-6a+5c)+3(a-c)=a
x+2b-12a+10c+3a-3c=a
x+2b-9a+7c=a