clear
fx='x1*exp(2.5*(1-x1-0.5*x2))';
fy='x2*exp(5*(1-0.5*x1-x2))';
x1=0.5;
x2=0.5;
axis([0 2 0 7]);
hold on
set(gca,'Fontsize',12)
xlabel('{\it x}_1','Fontsize',12)
ylabel('{\it x}_2','Fontsize',12)
title('Chaotic trajectory of (11.92)---Figure 11.21','Fontsize',12)
h5=plot(x1,x2,'erasemode','none');
drawnow
for i=1:30000
f1(i)=eval(fx);
f2(i)=eval(fy);
x1=f1(i);
x2=f2(i);
set(h5,'xdata',x1,'ydata',x2);
drawnow;
end
