function x=logistica2(lambda, trans, nmax, x0)
% x=logistica2(lambda, trans, nout, x0)
% esegue nmax iterazioni dell'equazione x(k+1)=lambda*x(k)(1-x(k)) a
% partire da ognuno dei valori di lambda
% esegue un plot dei risultati al variare di lambda

n=length(lambda); % numero parametri lambda
xk=x0*ones(1, n); % inizializza il valore iniziale delle n successioni
x=zeros(n, nmax); % matrice che conterra' i risultati: ogni riga si riferisce ad un diverso lambda
for i=1: trans
    xk=lambda.*xk.*(1-xk);
end
for i=1:nmax
    xk=lambda.*xk.*(1-xk);
    x(:, i)=xk';
end
plot(lambda, x, '.');

    
 