fcinv( replaceZeros(abs(t)) ));
ft = @(t,dim) cos( sqrt( normsqd(yinv(t)) )) *(sqrt(pi))^dim;
fKeister = @(x) ft(x,dim); exactInteg = Keistertrue(dim);
inputArgs ={'absTol',absTol, 'relTol',relTol};
inputArgs =[inputArgs {'arbMean',true}];
[~,muhat]=cubBayesNet_g(fKeister,dim,inputArgs{:});
check = double(abs(exactInteg-muhat) < max(absTol,relTol*abs(exactInteg)))

%%
% *Example 4: Multivariate normal probability*
%
% For $\bf{X}\sim N(\bf{\mu},\Sigma)$, estimate the following
% probability:
%
% $$ P\left(\