MAT
RHS
b=v;
ANS
x=v;
num_rhs = size(x,2);
resid = zeros(1,num_rhs);
for i=1:size(x,2),
	resid(i) = norm(a*x(:,i)-b(:,i))/norm(b(:,i));
end
format long e
resid