A=[1 0 -1;0 2 0; 1 0 1]
b=[4 2 -1]'
x=A\b
r=A*x-b