x(1) = 0;
x(2) = 0;
a = 1.4;
b = 0.3;

for n = 2 : 500
    x(n+1) = a - (x(n))^2 + b*x(n-1);
end;


y = [];
tau = 1 ;
for i = 1 : length(x)-tau
    y(i) = x(i+tau);
end;
plot(x(1:end-tau),y, '.');
title(num2str(tau));
pause;

y = [];
tau = 2 ;
for i = 1 : length(x)-tau
    y(i) = x(i+tau);
end;
plot(x(1:end-tau),y, '.');
title(num2str(tau));
pause;


y = [];
tau = 3 ;
for i = 1 : length(x)-tau
    y(i) = x(i+tau);
end;
plot(x(1:end-tau),y, '.');
title(num2str(tau));
pause;

y = [];
tau = 4 ;
for i = 1 : length(x)-tau
    y(i) = x(i+tau);
end;
plot(x(1:end-tau),y, '.');
title(num2str(tau));
pause;


y = [];
tau = 6 ;
for i = 1 : length(x)-tau
    y(i) = x(i+tau);
end;
plot(x(1:end-tau),y, '.');
title(num2str(tau));
pause;

y = [];
tau = 8 ;
for i = 1 : length(x)-tau
    y(i) = x(i+tau);
end;
plot(x(1:end-tau),y, '.')
title(num2str(tau));
pause;

y = [];
tau = 40 ;
for i = 1 : length(x)-tau
    y(i) = x(i+tau);
end;
plot(x(1:end-tau),y, '.')
title(num2str(tau));


% 
% y = x(1:tau:end-tau);
% 
% 
% for tau = 1 : 1 :3
%     y = x(1:tau:end);
