clear all

d=2.3;

x=linspace(0,2.5,201);


f=x.*exp((1-x)*d);
f2=f.*exp((1-f)*d);
f3=f2.*exp((1-f2)*d);
f4=f3.*exp((1-f3)*d);



figure(10)
plot(x,f,'b',x,x,'r',x,f2,'k',x,f4,'g')
xlabel('x')
ylim([0 1.1*max(f2)])
%hold on

%primitive way for constructing a web diagram
y(1)=0.6;
x1(1)=y(1);x2(1)=x1(1);
gh=linspace(0,1,5);
for ia=2:20;
    y(ia)=y(ia-1)*exp((1-y(ia-1))*d);
    g1=y(ia-1)*exp(0*gh);g2=linspace(y(ia-1),y(ia),5);
    g3=y(ia)*exp(0*gh);
    if ia==2;
    figure(12)
    plot(x,f,'b',x,x,'r',x,f2,'g',g1,g2,'--k',g2,g3,'--k')
    ylim([0 1.1*max(f2)])
    xlabel('$x_n$','Interpreter','Latex','Fontsize',18)
    ylabel('$x_{n+1}$','Interpreter','Latex','Fontsize',18)
    drawnow
    hold on
    else;
    figure(12)     
    plot(x,f,'b',x,x,'r',g1,g2,'--k',g2,g3,'--k')
    drawnow
    hold on  
    end;
end;
hold off