% m19_newtonc.m - same as m19.m but with extended precision output.

  x = str2sym('1');
  digits(80);
  x = vpa(x)
  for k = 0:7
    pause
    s = -(x^3-2) / (3*x^2);
    x = x + s     
  end