function [F,domerr] = mpecf(x,y)
% [F,domerr] = mpecf(x,y)

[F,domerr] = mpecfunc(1,x,y);
return;
