We need to start by rearranging the equations into the form: ax+by+cz=d
Equation 1:

Equation 2:

Equation 3:

Now, rewrite the system as an augmented matrix:

1. Subtract the first row from the second row, and the result will be the second row. It is written as:

The new matrix is:

The next operation is:

The matrix is:

Next step:

Next step:

Then, we have the following system of equations:

Replace equation 3 into equation 2 and solve for y:

Replace y and z into equation 1 and solve for x:

The solution set is { 1 , 5 , 1 }