% MATLAB program to check results of grid7 and grid8
% use bmat = 0 for Identity, bmat = 1 for BMAT file...

bmat = 1;
MATA
A = a;
if(bmat)
  MATB
  B= a;
else
  B = speye(size(A,1));
end
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-shift*B)*x(:,i)-b(:,i))/norm(b(:,i));
end
format long e
resid