% #9
A=[4,0,-7,-7;-6,1,11,9;7,-5,10,19;-1,2,3,-1]
rref(A)
% #42
x1 = rand(4,1)
b = A*x1
x2 = A\b
x1-x2