Integrates (0,inf) over multivariate normal
zeroToInfNormCDF(mu, sigma)
mu
vector of mean parameters
sigma
covariance matrix
double