
A = mark(30);
n = size(A,1);

x = randn(n,1);

[u, val,err] = arn0(A, x, 50, 1.0) ;
semilogy(err);
val 
norm(A*u - val*u)
