<< MB.m

<< AMBRE.m

<< KinematicsGen.m

int = 
        Fullintegral[{1}, {PR[k1, 0, n1] PR[k1 + p1, 0, n2] PR[
                k1 + p1 + p2, 0, n3] PR[k1 + p1 + p2 + p3, 0, n4] PR[
                k1 + p1 + p2 + p3 + p4, 0, n5] 
              PR[k1 + p1 + p2 + p3 + p4 + p5, 0, n6]}, {k1}];

kin = {p1^2 -> 0, p2^2 -> 0, p3^2 -> 0, 
          p4^2 -> 0, p5^2 -> 0, p6^2 -> 0};

choice = {s12, s23, s34, s45, s56, s16, s123, s234, 
          s345};

invariants = Kinematics[kin, choice]

IntPart[1]

rep = SubLoop[integral];

fin = BarnesLemma[rep, 1]

fin = fin /. n1 -> 1 /. n2 -> 1 /. 
              n3 -> 1 /. n4 -> 1 /. n5 -> 1 /. 
        n6 -> 1

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

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

ser0 = 
        MBexpand[{integrals}, 
          Exp[eps EulerGamma], {eps, 0, 0}];

AbsoluteTiming[
      MBintegrate[
        ser0, {s12 -> -1, s23 -> -2, s34 -> -3, 
          s45 -> -4, s56 -> -5, s16 -> -6, 
          s123 -> -7, s234 -> -8, 
          s345 -> -9}]]