factor = 0.01;
x = -2 + 4 * rand ();
for i = 1 : 10
  fx = f1(x);
  fpx = fp1(x);
  fprintf ( 1, '  %d  %g  %g  %g\n', i, x, fx, fpx );
  x = x - factor * fpx;
end

