clear;
close all;

a=1.4;b=0.3;
xn=0;yn=0;

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