A = rand(10,10);
b = rand(10,1);
x = inv(A)*b;
A*x-b