Feigenbaum Diagram

feigenbaum.m

% feigenbaum diagram
r = 1:.01:4; % r values
nr = length(r); % number of r values
nI = 600; % number of iterations to perform for each r
converge = 300; 
p = nan(converge,nr);
x = zeros(1,nI);
for rx=1:nr
    x(1)=.5;
    for i=2:nI
        x(i)=r(rx)*x(i-1)*(1-x(i-1)); % logistic map
        
    end
    % pick the last 100 values. if converge to single value, they all get 
    % plotted on same point. if 2 periods, we will see 2 values
    p(:,rx) = x((end-converge+1):end);
end
plot(r,p,'k.');