close all;
clear;

r=4;
xn0=rand(1);
yn0=xn0+0.001;
x=0:0.01:1;
n=50;
nf=10;


for(i=1:n)
        if(i==1)
          xn=xn0;
          yn=yn0;
          orbit(i)=xn0;
          orbity(i)=yn0;
        else
            %xn=xn1;
            xn=r*xn*(1-xn);
            yn=r*yn*(1-yn);
            orbit(i)=xn;
            orbity(i)=yn;
        end
end

%subplot(2,2,1:2)
plot(orbit,'b-'); hold on; plot(orbity,'r--');
xlabel('n');
ylabel('x_n azul     y_n rojo');
axis([0, max(i), 0, 1]);
hold on;

figure (2);
plot(orbity-orbit);
title('Diferencia entre las trayectorias');
xlabel('n');
ylabel('y_n - x_n');