
figure(1)
clf
set(gca,'FontSize',14);

x=0:0.01:1;
%g = inline('2*x');
g = inline('3.84.*x.*(1-x)');

Npts=150;

x0=0.1;
X=[x0];
Y=[0];

xo=x0;

for i=1:Npts
 xn=g(xo);   
 X=[X xo xn];
 Y=[Y xn xn];
 xo=xn;
end

plot(x,0*x,'k')
hold on
plot(g(x)*0,g(x),'k')
plot(x,x)
plot(x,g(x),'r','LineWidth',2)
plot(X,Y,'m')
xlabel('x_n');
ylabel('x_{n+1}');
axis tight
