Chaos them: explore the sequence xn+1 = f(xn) where f(x) =r x(1-x),...

function cobweb(r,x0,N, to_wait)
%Inputs
% r -- the coefficient of f(x)=rx(1-x); r should be in the range of 0 to 4
% x0 -- initial condition (.5 default)
% N -- the number of iterates for which to draw the cobweb (100 default)
% to_wait -- 0 or 1 (0 default)
%       0 -- draw the cobweb straight from the beginning
%       1 -- iterate the point x0 1000 times BEFORE drawing the cobweb N
%       times