function y = f(x)

  y = (1-x).^7;
% y = 1 - 7.*x + 21.*x.^2 - 35.*x.^3 + 35.*x.^4 - 21.*x.^5 + 7.*x.^6 - x.^7; 
