function y = df(x)
  y = 1 + 2*x*sin(1/x) - cos(1/x); 
