% Chaos control
clear
itermax=200;
mu=4;a=3.5;x=zeros(1,itermax); kontroluj=0;
p=(1-kontroluj)*mu+kontroluj*a;
x(1)=0.88;
for n=1:1:itermax
    if (n>50) & (x(n)<0.75) & (x(n)>0.7)     % kontroluj pro n>50 
    kontroluj=1;
    a=4+10*(x(n)-0.75);
    end    
 x(n+1)=((1-kontroluj)*mu+kontroluj*a)*x(n)*(1-x(n));
end
hold on
    plot(1:itermax,x(1:itermax))
    plot(1:itermax,x(1:itermax),'o')
fsize=15;
set(gca,'XTick',0:50:itermax,'FontSize',fsize)
set(gca,'YTick',[0,1],'FontSize',fsize)
xlabel('n','FontSize',fsize)
ylabel('\itx_n','FontSize',fsize)
hold off

% End of Program_6a.