%Sample script to produce several orbits of x^2-2

xinit=[0,0.1,0.01,0.001];
Nx(1,:)=xinit;

x=xinit;
for i=2:101
   Nx(i,:)=x.^2-2;
   x=Nx(i,:);
end

plot(Nx(:,1));  %A plot of the first orbit
figure
plot(Nx(:,2));
figure
plot(Nx(:,3));
figure
plot(Nx(:,4));

%Use the command:  close all   to close all the figures
% at once...