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