Automated Computation of Multiple Mellin–Barnes Integrals Having Polygamma Functions in Their Integrand

The conic hulls/triangulation method is extended to the analytic computation of multiple Mellin–Barnes having polygamma functions of arbitrary order in their integrand. This computational approach is automatized in the <span class="tt">MBConicHulls.wl</span> <span class="sf">Mathematica</span> package, which we give a brief description of the new corresponding functionalities, basing our discussion on the example of a simple two-loop sunset integral. Abstract Published by the Jagiellonian University 2026 authors