function [c,ceq] = mycon(x)
          c = (x(1)+3)^2 + (x(2)+3)^2 - 36;
          ceq = [];
end

