Let D be a polydisk in C-n and f : (D) over bar --> C-n a mapping that is analytic in D and has no zeros on the boundary of D. Then f has only a finite number of zeros in D and these zeros are all isolated. We consider the problem of computing these zeros. A multidimensional generalization of the classical logarithmic residue formula from the theory of functions of one complex variable will be our means of obtaining information about the location of these zeros. This integral formula involves the integral of a differential form, which we will transform into a sum of n Riemann integrals of dimension 2n - 1. We will show how the zeros and their multiplicities can be computed from these integrals by solving a generalized eigenvalue problem that has Hankel structure, and n Vandermonde systems. Numerical examples are included.