clear all;

a = 1.4;
b = 0.3;
x(1) = 1;
y(1) = 1;
numpts = 1000000;
for n=1:numpts
    x(n+1) = 1-a*x(n)^2 + y(n);
    y(n+1) = b*x(n);
end


figure(1)
subplot(2,1,1)
plot(x(1,2000:numpts),y(1,2000:numpts),'.')
axis([-1.5,1.5,-0.5,0.5])
subplot(2,1,2)
plot(x,'.')




