

function [x,rip]=logisticMap(x1,r,lun)

%Calcolo di x(n+1) = r * x(n) * (1-x(n))

rip=0;
x=x1;

% continua finchè la lunghezza "lun" non è stata raggiunta o "rip=1"
while (length(x) < lun) & (rip == 0)
    
        % l'elemento "xnew" è determinato dalla formula e "x(end)" è
        % l'ultimo elemento precedentemente determinato 
        xnew=r*x(end)*(1-x(end));
        x=[x xnew];
        
        % controllo se "x(end)" e' quasi-uguale a un valore precedente
        if min(x(1:end-1)-x(end)) < 1e-4
            rip=1;
        end
end
