Application of the Meijer theorem in calculation of three-loop massive vacuum Feynman integrals and beyond

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.