function [c, ceq] = constraint(x)
    c = -x(1); % f >= 0
    ceq =  [];
end