function y = df(x)
  y = sinh(x) - sin(x); 