% Logistic orbit for alpha varying from 2.6 to 4.

s = 1000; 
k = 700;
x0 = rand(1,s); 

A = linspace(2.6,4,s);
x = zeros(k,s); 
x(1,:) = x0;

for i = 1:k;    
    x(i+1,:) = A.*x(i,:).*(1-x(i,:));   
end

figure; 
clf;

plot(A,x(floor(2*end/3):end,:),'b.','markersize',2);
title('Logistic Orbit');
xlabel('Parameter \alpha '); 
ylabel('X');
grid


% "Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009.
% Copyright(c)