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