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