% Script to run tests on RayleighQuotientF.m
% K. Ming Leung, 09/17/03

clear all;
% Example 4.13: Heath, p.178
A=[3 1; 1 3]; [evec,eval]=PowerRayQuotF(A,[0;1],100,1e-5);
% 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]=PowerRayQuotF(A,[1;1;1;1],100,1e-5);
% % From MATLAB's eigenfunction
% [Mevec,Meval]=eig(A),