function [c, ceq] = restr(x,X,Y,Z,V) % Nonlinear inequality constraints c = []; % Nonlinear equality constraints ceq = x(2)*x(3)+x(4);