function df = dfunc(x,a)
% df = exp(x);
  df = sech(x)^2 + a;
