<< MB.m

<< AMBRE.m

<< KinematicsGen.m

a = {p1^2 -> m^2, p2^2 -> m^2, p3^2 -> m^2, 
          p4^2 -> 0, p5^2 -> m^2};

b = {s12, s23, s34, s45, s15};

invariants = Kinematics[a, b]

Fullintegral[{1}, {PR[p1 + k1, 0, n1]*PR[p1 + p5 + k1, m, n2]*
            PR[p1 + p5 + p4 + k1, m, n3] PR[p1 + p5 + p4 + p3 + k1, 0, 
              n4] PR[k1, m, n5]}, {k1}];

FUPolynomials[
          DS[p1 + k1, 0, n1]*DS[p1 + p5 + k1, m, n2]*
            DS[p1 + p5 + p4 + k1, m, n3] DS[p1 + p5 + p4 + p3 + k1, 0, 
              n4] DS[k1, m, n5], {k1}, invariants][[1]] // 
      Expand

IntPart[1]

SubLoop[integral];

Fauto[0]

fupc = fupc /. s45 -> s45p + m^2 /. s34 -> s34p + m^2 // 
        Expand

repr = SubLoop[integral];

BarnesLemma[repr, 1]

repr = % /. {n1 -> 1, n2 -> 1, n3 -> 1, n4 -> 1, 
          n5 -> 1}

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

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

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

MBintegrate[ser /.  s45p -> s45 - m^2  /.  
          s34p -> s34 - m^2, {s12 -> -1/2, s15 -> -3, 
          s23 -> -5, s34 -> -2, s45 -> -10, 
          m -> 1}, Verbose -> False] // AbsoluteTiming