% To show a bug in Matlab's gmres.m routine!
% From Alison Ramage (alison@maths.strath.ac.uk)

clear all
load matdata
K=n^2
x0 = zeros(K,1);
maxits = K;
tol = 1e-6;
[x_mat,flag,relres,iter,resvec] = gmres(A,b,[],tol,maxits);
iter
plot(resvec)