function f = fun(x)
f = [x(1)^2+x(2)^2-1;...
     sin(pi*x(1)/2)+x(2)^3];
