function f = fun(xx)
  f(1,1)=xx(1)^2+xx(2)^2+xx(3)^2-100;
f(2,1)=xx(1)*xx(2)*xx(3)-1;
f(3,1)=xx(1)-xx(2)-sin(xx(3));
end