  We extend the (continuous) multivariate Almkvist-Zeilberger algorithm in
order to apply it for instance to special Feynman integrals emerging in
renormalizable Quantum field Theories. We will consider multidimensional
integrals over hyperexponential integrands and try to find closed form
representations in terms of nested sums and products or iterated integrals. In
addition, if we fail to compute a closed form solution in full generality, we
may succeed in computing the first coefficients of the Laurent series
expansions of such integrals in terms of indefinite nested sums and products or
iterated integrals. In this article we present the corresponding methods and
algorithms. Our Mathematica package MultiIntegrate, can be considered as an
enhanced implementation of the (continuous) multivariate Almkvist Zeilberger
algorithm to compute recurrences or differential equations for hyperexponential
integrands and integrals. Together with the summation package Sigma and the
package HarmonicSums our package provides methods to compute closed form
representations (or coefficients of the Laurent series expansions) of
multidimensional integrals over hyperexponential integrands in terms of nested
sums or iterated integrals.
