function mappa_logistica(x0,rmin,rmax,ntot)
% programma dinamica_popolazione

if nargin < 3,
    fprintf(' inserire x0, r, ntot\n')
    return
end

x=zeros(1,ntot);
x(1)=x0;

if rmin == rmax
 r = rmin;
    for i = 1:ntot
        x(i+1)=r*x(i)*(1-x(i));
        if x(i+1)<0
            x(i+1)=0;
        end
    end
figure(1);
plot(x);
ylabel('xn')
xlabel('n');
s_title=sprintf('Mappa logistica con r=%.2f',r);
title(s_title);
ylim([0 1]); 
return 
end

iterazioni=200;
dr=(rmax-rmin)/iterazioni;

for r=rmin:dr:rmax
    for i = 1:ntot
        x(i+1)=r*x(i)*(1-x(i));
        if x(i+1)<0
            x(i+1)=0;
        end
    end
figure(1);
plot(x);
ylabel('xn')
xlabel('n');
s_title=sprintf('Mappa logistica con r=%.2f',r);
title(s_title);
ylim([0 1]); 
pause(0.2);
    
end

return