function x=newtone(x0,e)
  i=1
  x=x0-func(x0)/funcd(x0)
  while abs(x-x0)>e
    i=i+1
    x0=x;
    x=x0-func(x0)/funcd(x0)
  end