


clear all

% delta=input('delta = ?')
% 
% x=linspace(1.e-03,1.4,100)';
% 
% 
% f=x-delta*x.*log(x);df=1-delta-delta*log(x);
% ff=f-delta*f.*log(f);dff=(1-delta-delta*log(f)).*df;
% fff=ff-delta*ff.*log(ff);dfff=(1-delta-delta*log(ff)).*dff;
% ffff=fff-delta*fff.*log(fff);dffff=(1-delta-delta*log(fff)).*dfff;
% 
% 
% figure(20)
% plot(x,f,'r',x,ff,'b',x,ffff,'g',x,fff,'m',x,x,'--k','linewidth',2)
% xlabel('$x_n$','Interpreter','latex')
% ylabel('$x_{n+1}$','Interpreter','latex')
% set(gca,'Fontsize',24)
% 
% ghgh
% figure(21);subplot(211)
% plot(x,f,'r',x,ff,'b',x,x,'--k','linewidth',2)
% xlabel('$x_n$','Interpreter','latex')
% ylabel('$x_{n+1}$','Interpreter','latex')
% set(gca,'Fontsize',24)
% figure(21);subplot(212)
% plot(x,dff,'b','linewidth',2)
% xlabel('$x$','Interpreter','latex')
% ylabel('$D f(f(x)$','Interpreter','latex')
% ylim([-1 1])
% set(gca,'Fontsize',24)
% 
% 
% x0=linspace(0.2,1.2,7)';



D=linspace(1.8,2.6,501);
%D=linspace(2.45,2.6,501);
%D=linspace(2.5,2.51,501);
x0=sqrt(2)/2;



for id=1:length(D);
    
    delta=D(id);
    
    Nt=10000;
    xt=x0;
    ft=[];
    for it=1:Nt;
        
        xt=xt-delta*xt.*log(xt);
        
        if it>Nt-400;
            ft=[ft;xt];
        end;
    end;
    
       F(:,id)=ft;
       
end;


figure(520)
plot(D,F,'.k')
xlabel('$\delta$','Interpreter','latex')
 ylabel('$x_\infty$','Interpreter','latex')
set(gca,'Fontsize',24)
% 





    
    

