%
% Read the comments in fom.m and matrix5p.m first ;)
%
function test_krylov(k,m)
[A,b]=matrix5p(k);
n = size(A,1);
fprintf(1,'Grid size k for PDE = %d. ',k);
fprintf(1,'So A is %d x %d matrix...\n',n,n);
fprintf(1,'Number of vectors in K_m(A,r0) is at most %d\n ',m);
x0=rand(size(b));
x=fom(A,b,x0,m);