mu=1;
x0=0.4;
x=zeros(201,1);
xcob=zeros(402,1);
ycob=zeros(402,1);
x(1)=x0;
xold=x0;
xcob(1)=xold;
fofx=mu-xold^2;
ycob(1)=fofx;
xcob(2)=fofx;
ycob(2)=fofx;
for i=2:201
   fofx=mu-xold^2;
   xnew=fofx;
   xcob(2*i-1)=xnew;
   fofx=mu-xnew^2;
   ycob(2*i-1)=fofx;
   xcob(2*i)=fofx;
   ycob(2*i)=fofx;
   xold=xnew;
   x(i)=xold;
end
axx=[-1.2;1.2];
axy=[0;0];
plot(axx,axy)
hold on
ayx=[0;0];
ayy=[-1.2;1.2];
plot(ayx,ayy)
xcurve=[-1.2:0.02:1.2]';
ycurve=mu-xcurve.^2;
plot(xcurve,xcurve)
plot(xcurve,ycurve)
plot(xcob,ycob,'r')
hold off