<< MB.m

<< AMBRE.m

Fullintegral[
{1}, 
{PR[k1, m, 
              n1] PR[k1 + p1, 0, n2] PR[k1 + p1 + p2, m, n3] 
            PR[k1 - k2, 0, n4] PR[k2, m, n5] PR[k2 + p1 + p2, m, n6] 
            PR[k2 - k3, 0, n7] PR[k3, m, n8] PR[k3 + p1 + p2, m, n9] PR[
              k3 - p3, 0, n10]}, {k3, k1, k2}];

invariants = {p1^2 -> m^2, p2^2 -> m^2, p3^2 -> m^2, 
              p4^2 -> m^2, p1*p2 -> 1/2*s - m^2, 
              p3*p4 -> 1/2 s - m^2, p1*p3 -> 1/2*t - m^2, 
              p2*p4 -> 1/2*t - m^2, p2*p3 -> 1/2*u - m^2, 
              p1*p4 -> 1/2*u - m^2} /. u -> 4 m^2 - s - t // 
          Expand;

IntPart[1]

SubLoop[integral]

IntPart[2]

SubLoop[integral]

IntPart[3]

repr = SubLoop[integral];

BarnesLemma[repr, 1]

fin = % /.  {n1 -> 1, n2 -> 1, n3 -> 1, n4 -> 1, 
          n5 -> 1, n6 -> 1, n7 -> 1, n7 -> 1, 
          n8 -> 1, n9 -> 1, n10 -> 1}

rules = MBoptimizedRules[fin, eps -> 0, {}, {eps}]

integrals = 
        MBcontinue[fin, eps -> 0, rules, 
          Verbose -> False];

ser = 
        MBexpand[{integrals}, 
          Exp[3*eps*EulerGamma], {eps, 0, -1}];

AbsoluteTiming[
      MBintegrate[ser, {s -> -1/2, t -> -3, m -> 1}, 
        Verbose -> False]]

\({11.712732000000104`\ Second, {\(\(0.`\)\(\)\) - 
          0.04385500414000002`\/eps\^3 + 0.009501169729000059`\/eps\^2 + 
          0.6472366386310029`\/eps, {6.919647141789538`*^-6\/eps\^3 + 
            0.00004230482052783107`\/eps\^2 + 0.010966061249257182`\/eps, 
          0}}}\)