%Программа для построения паутинки для нахождения неподвижных точек
% функции f=a*x.*(1-x) при различных значениях коэффициента a
% Начальное положение точки x=z0
%a=3.449499; % Значение коэффициента a
a=3.77; % Значение коэффициента a
x0=0:0.01:1;
y0=a*x0.*(1-x0);
hold on;
plot(x0,x0,'r-'); % График y=x
plot(x0,y0,'k');  % График y = ax(1-x)
%axis([-.01 1 -.01 1])

% надо ограничить значения прямоугольником 
% и надо остановить процесс при достижении периодической точки
% период - выделить жирно или другим цветом 

%draw pautinu for prityagiv. cikla.
x0=0.9; % Начальная точка x=z0
n=100;
xn=x0;
xn1=a*xn.*(1-xn); % x = ax(1-x) Первая итерация
for i=1:n
     plot([xn,xn1],[xn1,xn1]); % Гориз. Отрезок от точки графика до y=x
        pause(0.2*(1-i/(n+1)));
     yxn1=xn1; % точка на y=x 
     xn=xn1;       % переобозначим для цикла, теперь начальная точка xn на пересечении y=x 
     xn1=a*xn.*(1-xn); % точка на графике
     plot([yxn1,yxn1],[yxn1,xn1]); % вертикальная линия от y=x до графика
        pause(0.2*(1-i/(n+1)));
end
hold off;
