n = 100;
r = linspace(2,3.9,100);
hold on
for k = 1: length(r)
    x(1) = 0.1;
    for j = 1:n-1
        x(j+1) = r(k).*x(j).*(1-x(j));
        if (j>0.9*n)
            plot(r(k),x(j),'b.','MarkerSize', 1);
        end
    end
end
title('Logisticka funkcija');
xlabel('r');
ylabel('x');