Solve the system of equations using matrices. Use Gaussian elimination with back- substitution.

x+y+z = -5
x-y+3z = -1
4x+y+z = -2