% Script to run function QR_eigFt.m that uses
% QR iteration to find the eigenvalues (& eigenvectors)
% of matrix A
% Algorithm 4.7 on p.182 of Heath
% Example 4.15 on p. 183 of Heath
% K. Ming Leung, 09/18/03

A = [2.9766 0.3945 0.4198 1.1159;
0.3945 2.7328 -0.3097 0.1129;
0.4198 -0.3097 2.5675 0.6079;
1.1159 0.1129 0.6079 1.7231];
disp(A);
itMax = 3;

eval = QR_eigF(A,itMax);