range = -2:0.01:2;
plot(range,atan(range),'r-');
line([min(range), max(range)],[0,0]);
makelarge
hold on
x = 1;
for j = 1:5,
   xold = x;
   x = x - (1+x^2)*atan(x),
   line([xold x], [atan(xold) 0]);
   makelarge
   pause
end
hold off

pause

range = -2:0.01:2;
plot(range,atan(range),'r-');
line([min(range), max(range)],[0,0]);
makelarge
hold on
x = 1.39;
for j = 1:10,
   xold = x;
   x = x - (1+x^2)*atan(x),
   line([xold x], [atan(xold) 0]);
   makelarge
   pause
end
hold off

pause

range = -2:0.01:2;
plot(range,atan(range),'r-');
line([min(range), max(range)],[0,0]);
makelarge
hold on
x = 1.40;
for j = 1:10,
   xold = x;
   x = x - (1+x^2)*atan(x),
   line([xold x], [atan(xold) 0]);
   makelarge
   pause
end
hold off