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

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