function []=prob_9_3_14()
%
% 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 = [ 0.9 0.3; 0.1 0.7 ]; 
[L,U] = lu(A); 

M = 20; u_k = [ 0.1 0.9 ].'; 

res=[]; 
for k=1:M,
  u_k = U \ ( L \ u_k );
  res = [ res, (0.6)^k * u_k ]; 
end
res