function y = f(x)

% same as y = (x+1)*(x-1)^2;
y = 1-x-x^2+x^3;