 function x = logmap(x0,r,N)
% x = logmap(x0,r,N)
% Computes 100 iterates of the nonlinear difference eq:
%      x(n+1) = r*x(n)*(1-x(n))
% Input initial condition 'x0' (between 0 and 1)
% and value of the initial gowth rate 'r'o

   x(1) = x0;

   for i =1:N
     x(i+1) = r*x(i)*(1 - x(i));
   end

   return

