% Script to run tests on InverseIterationF.m
% InverseIterationFt.m

%K. Ming Leung, 01/11/03

clear all;
% Example 4.13: Heath, p.178. Exact eigenvalues are 2 and 4
A=[3 1; 1 3]; [evec,eval]=InverseIterationF(A,[0;1],100,1e-5),
%A=[1 4; 1 1]; [evec,eval] = InverseIterationF(A,[1;1],100,1e-4), 
% RESULTS:
% >>        it                  m                   x
%      1     1     1     0
%             2      0.33333           -1            1
%             3       1.6667            1         -0.4
%             4      0.86667           -1      0.53846
%             5       1.0513            1      -0.4878
%             6      0.98374           -1      0.50413
%             7       1.0055            1     -0.49863
%             8      0.99817           -1      0.50046
%             9       1.0006            1     -0.49985
%            10       0.9998           -1      0.50005
%            11       1.0001            1     -0.49998
%            12      0.99998           -1      0.50001
% evec =
%            -1
%       0.50001
% eval =
%             1

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


% % Table 7.1: Fausett, p.196. Exact eigenvalues are 2 3 4
% A=[ -44  9  -6
%    -280 57 -40
%     -75 15 -13];
% [evec,eval]=InverseIterationF(A,[1;1;1],100,1e-5),
% % From MATLAB's eigenfunction
% [Mevec,Meval]=eig(A),
% 
% % Table 7.4: Fausett, p.199. Exact eigenvalues are 2 1 -3
% A=[ 66 -21  9
%    228 -73 33
%     84 -28 16];
% [evec,eval]=InverseIterationF(A,[1;1;1],100,1e-5),
% % From MATLAB's eigenfunction
% [Mevec,Meval]=eig(A),


% 
% % Table 7.2: Fausett, p.196. Exact eigenvalues are 2, 4, 6, 8
% 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]=InverseIterationF(A,[1;1;1;1],100,1e-5),
% % From MATLAB's eigenfunction
% [Mevec,Meval]=eig(A),
