clear;
close all;

a=1.45;b=0.3;
xn1=0;yn1=0;

for i=1:200
    xn=1-a*xn1^2+yn1;
    yn=b*xn1;
    xn1=xn;yn1=yn;
end
for i=1:10000
    xn=1-a*xn1^2+yn1;
    yn=b*xn1;
    xn1=xn;yn1=yn;
    plot(xn,yn,'*b');
    hold on;
end