function [theta,g,domerr] = mpecdobj(x,y)
% [theta,g,domerr] = mpecdobj(x,y)

[theta,g,domerr] = mpecfunc(4,x,y);
return;
