%   
% Written by:
% -- 
% John L. Weatherwax                2009-04-21
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

close all; drawnow; clc; clear; clear functions; 

alpha = 0.1; omega = 1 - alpha; 

j = 1:30; 

% for l=1: 
pi_1 = alpha * ( 2 - alpha * j + ( alpha / ( 1 - alpha ) ) * ( 1 - alpha * j ) * 1 ) .* omega.^( j-1 ); 
fh=figure; ph1=plot( j, pi_1, 'ko' ); 

hold on; 
pi_2 = alpha * ( 2 - alpha * j + ( alpha / ( 1 - alpha ) ) * ( 1 - alpha * j ) * 2 ) .* omega.^( j-1 ); 
ph2=plot( j, pi_2, 'bx' ); xlabel('j (index)'); ylabel( '\pi value' ); axis tight; 
legend( [ph1,ph2], {'\pi with j=1','\pi with j=2'} ); 

saveas( gcf, '../../WriteUp/Graphics/Chapter3/prob_3_8_pi', 'epsc' );