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

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