function [time,x]=logis(alpha,xinit,n)

x(1)=xinit;
time(1)=0;

for k=2:n+1;
time(k)=k-1;
x(k)=alpha*x(k-1)*(1-x(k-1));
end

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