function dinamica_popolazione_r(x0,r)
% programma dinamica_popolazione

if nargin < 2, 
    fprintf(' inserire x0, r\n')
    return
end
 
time=abs(10/r);
DT=time/50;
N=floor(time/DT);
y=zeros(N,1);
y(1)=x0;
for i = 1:N
    y(i+1,1)=r*y(i)*(1-y(i));
end

t=(0:DT:time);
figure(1);
plot(t, y);
ylabel('Popolazione')
xlabel('tempo');
s_title=sprintf('Dinamica popolazione RK4 lin.');
title(s_title);
xlim([t(1) t(end)]);
ylim([0 max(y)*1.2]);

return