,dim) cos( sqrt( normsqd(yinv(t)) )) *(sqrt(pi))^dim;
fKeister = @(x) ft(x,dim); exactInteg = Keistertrue(dim);
inputArgs ={<span class="string">'absTol'</span>,absTol, <span class="string">'relTol'</span>,relTol};
inputArgs =[inputArgs {<span class="string">'arbMean'</span>,true}];
[~,muhat]=cubBayesNet_g(fKeister,dim,inputArgs{:});
check = double(abs(exactInteg-muhat) &lt; max(absTol,relTol*abs(exactInteg)))
</pre><pre class="codeoutput">
check =

     1

</pre><p><b>Example 4: Multivariate normal probabi