function [f,J] = banana(x)
f = [100*(x(2)-x(1)^2)^2; (1-x(1))^2];
if nargout > 1, J=[-20*x(1) 10; -1 0]; end
   