I describe a package written in MATHEMATICA that automatizes typical operations performed during evaluation of Feynman graphs with Mellin-Barnes (MB) techniques. The main procedure allows to analytically continue a MB integral in a given parameter without any intervention from the user and thus to resolve the singularity structure in this parameter. The package can also perform numerical integrations at specified kinematic points, as long as the integrands have satisfactory convergence properties. I demonstrate that, at least in the case of massive graphs in the physical region, the convergence may turn out to be poor, making naive numerical integration of MB integrals unusable. I present possible solutions to this problem, but argue that full automatization in such cases may not be achievable.

23 pages, 11 figures, numerical evaluation functionality added