<< MB.m

<< AMBRE.m

Fullintegral[{1}, {PR[k1, 0, n1] PR[k2 + k4, m, n2] PR[k2, m, n3] 
            PR[k2 - k1, m, n4] PR[k3, 0, n5] PR[k2 - k1 + k3, 0, n6] 
            PR[k4, 0, n7] PR[k2 - k1 + k4, 0, n8]}, {k4, k3, k1, 
          k2}];
invariants = {};

IntPart[1]

SubLoop[integral]

IntPart[2]

SubLoop[integral]

IntPart[3]

SubLoop[integral]

Fauto[0]

fupc = 
      fupc /.  PR[k2, m] X[1] X[2] -> 
            PR[k2, 0] X[1] X[2] - m^2 X[1] X[2] // Expand

SubLoop[integral]

IntPart[4]

SubLoop[integral];

repr = BarnesLemma[%, 1]

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

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

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

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

MBintegrate[ser, {m -> 1}]