function [g,Jg,domerr] = mpecdcon(x,y)
% [g,Jg,domerr] = mpecdcon(x,y)

[h,Jh,g,Jg,domerr] = mpecfunc(6,x,y);
g = [h; g];
Jg = [Jh sparse(length(h),length(y)); Jg];

return;
