function z=objf(X)
x=X(1); y=X(2);
z=100*(y-x.^2).^2 + (1-x).^2;
