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

clc; 
close all; 

x = [ round( (2/3)*50 ) ];
while x(end)~=1
  x = [ x, round( (2/3)*x(end) ) ];
end
x