A = hilb(20);
b = zeros(20,1);
b(1) = 1;
x1 = A \ b
b(1) = b(1)+0.001;
x2 = A \b