function f = fun(x)
f = (x(1)-1)^2+10*(x(1)-x(2)^2)^2;
