% 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 = 300;

eval = QR_eigF(A,itMax);

