function [f,J,domerr]=cpfunjac(z,jacflag) 

if jacflag == 1
  [f,J,domerr] = mcpdf(z);
else
  [f,domerr] = mcpf(z);
end

return;
