  We present an analytical method to calculate the three-loop massive Feynman
integral in arbitrary dimensions. The method is based on the Mellin-Barnes
representation of the Feynman integral. The Meijer theorem and its corollary
are used to perform the integration over the Gamma functions, exponential
functions, and hypergeometric functions. We also discuss the application of the
method in other multi-loop Feynman integrals.
