function y = df(x)
  
  y = -7.*(1-x).^6;
% y =  7 + 42.*x - 105.*x.^2 + 140.*x.^3 - 105.*x.^4 + 42.*x.^5 - 7.*x.^6;
