demoAitkd2.m

%  demoAitkd2.m
%%%%%%%%%%%%%%%
global A

fprintf('\n    Highest eigenvalue and eigenvector of the Hilbert matrix\n\n')

n=inp('n   ',10,'%i');
A=hilb(n);
x=rand(n,1);
epsx=inp('epsx',1e-15,'%6.0e');
w=inp('w   ',.5,'%3.1f');
maxit=200;
ipr =inp('ipr ',1,'%i');

[x,cnt]=aitkd2('eigpwr',x,epsx,w,maxit,ipr);
smax = norm(A*x);
[V,S] = eig(hilb(n));

fprintf('\n   n iter = %i\n\n',cnt);
disp('            eig                       eigpwr')
disp(' ')
format long e
disp([S(n,n), smax])
disp(' ')
disp([V(:,n), x ])
format short