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

clear all; format short g;
% Example 4.13: Heath, p.178
A=[3 1; 1 3]; ; 
[evec,eval]=PowerRayQuotF(A,[0; 1],500,1e-4); 
% 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),
