poly2_newton_system

poly2_newton_system:
  Use Newton's method for a root of the poly2 function.
it = 2: alpha = 0.177797, log(alpha)= -1.72711
it = 3: alpha = 0.0329564, log(alpha)= -3.41257, r = 1.97588
it = 4: alpha = 0.00110876, log(alpha)= -6.80451, r = 1.99396
it = 5: alpha = 1.21189e-06, log(alpha)= -13.6233, r = 2.0021

  Number of steps = 5
  x = (   1.203166963347774,   -1.374080534239942)
  f(x) = (-1.77636e-15, 0)
diary off
poly2_newton_system2

poly2_newton_system2:
  Use Newton's method for a root of the poly2 function.
it = 2: alpha = 0.177797, log(alpha)= -1.72711
it = 3: alpha = 0.0329564, log(alpha)= -3.41257, r = 1.97588
it = 4: alpha = 0.00110876, log(alpha)= -6.80451, r = 1.99396
it = 5: alpha = 1.21189e-06, log(alpha)= -13.6233, r = 2.0021

  Number of steps = 5
  x = (   1.203166963347774,   -1.374080534239942)
  f(x) = (-1.77636e-15, 0)
diary off