multivariate zeta function