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

clc; 

t = [ (1-sqrt(5))/2, (1+sqrt(5))/2 ] 

et = 8*t 

rho = t(2); 
p_0 = ((1-rho^3)/(1-rho))^(-1);
p_0