
clear all
global a e1 Tp e2 e3
e1=0.;
e2=0.5;
e3=0;

a=0.2;

x0=0.1;
Tp=1;
Nt=100;
tt=linspace(0,1,Nt);
dt=tt(2)-tt(1);


XX=linspace(0.1,1,51);
XX=0.4;

for ix=1:length(XX)
    x0=XX(ix);

xt(1)=x0;

for it=2:length(tt)
    xa=xt(it-1);ta=tt(it-1);
     %f1=xa*(1-xa)-a*(1+e*sin(2*pi*ta/Tp));
     f2=logisticdot(xa,ta);
    
     xt(it)=xa+dt*f2;
     % if xt(it)<0
     %     xt(it)=NaN;
     % end
end

figure(1)
plot(tt,xt,'b','linewidth',2)
xlabel('t');ylabel('x(t)')
set(gca,'Fontsize',24)
hold on
P(ix)=xt(end);

end

figure(20)
plot(XX,P,'b',XX,XX,'--k','linewidth',2)
xlabel('$x_n$','interpreter','latex');
ylabel('$x_{n+1}=P(x_n)$','interpreter','latex')
set(gca,'Fontsize',24)
hold on







