% Script to run tests on PowerIterationF.m
% PowerIterationFt.m

% Example 4.13: Heath, p.178
A=[3 1; 1 3]; [evec,eval] = PowerIterationF(A,[0;1],100,1e-5), 
 
% RESULTS:
% >>        it                  m                   x
%             1            5            1          0.4
%             2          2.6            1      0.53846
%             3       3.1538            1       0.4878
%             4       2.9512            1      0.50413
%             5       3.0165            1      0.49863
%             6       2.9945            1      0.50046
%             7       3.0018            1      0.49985
%             8       2.9994            1      0.50005
%             9       3.0002            1      0.49998
%            10       2.9999            1      0.50001
% evec =
%             1
%       0.50001
% eval =
%        2.9999

% From MATLAB's eigenfunction
[Mevec,Meval]=eig(A),

% % Table 7.2: Fausett, p.196
% A=[ 4  2/3  -4/3  4/3
%    2/3   4     0    0
%   -4/3   0     6    2
%    4/3   0     2    6];
% [evec,eval]=PowerIterationF(A,[1;1;1;1],100,1e-5); 
% % From MATLAB's eigenfunction
% [Mevec,Meval]=eig(A),
