multivariate probability function