x = 2;
delta = 1;

while abs(delta) > 10^-13
  delta = -f(x)/df(x);
  x = x+delta;
end

x
