multivariate power series distribution