A=[4 5 6 7; 8 10 5 1; 2 6 1 3; 4 1 3 6]; b=[60 47 29 39]';   
gem(A,b);