function y = fx(x)

y = (x-1)^2 + 2*(x+1)*(x-1);