function eta = back_norm(A,b,x)
% BACK_NORM: Compute the backward error
eta = norm(b-A*x,inf)/(norm(A,inf)*norm(x,inf)+norm(b,inf));