x(1)=0.3; y(1)=0.2; %% Initial conditions
N=500; %number of iterates
a=1.4;b=0.3; %params

for i=1:N 
    x(i+1)=y(i)+1-a*x(i)^ 2; 
    y(i+1)=b*x(i);
end

figure
set(gca,'FontSize',20); 
for i=1:N
plot(x(i),y(i),'k.'); hold on ; drawnow
end


for i=1:5000 
    x(i+1)=y(i)+1-a*x(i)^ 2; 
    y(i+1)=b*x(i);
end

figure
set(gca,'FontSize',20); 
plot(x,y,'k.');