function []=prob_9_3_17()
%
% Written by:
% -- 
% John L. Weatherwax                2006-12-20
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----
  
clc; close all; 

A = [ 2 -1; -1 2 ]; [L,U] = lu(A); 
fprintf('The eigenvalues of A are...\n'); 
e = eig(A)
fprintf('Thus the eigenvalues of A^(-1) are...\n'); 
1./e

M = 12; u_k = [ 1 0 ].'; 

res=[]; 
for k=1:M,
  u_k = U \ ( L \ u_k );
  res = [ res, u_k ]; 
end
fprintf('the sequence of eigenvectors generated is given by ...\n'); 
res