<<MBv1.2.m

<<AMBRE.m

Options[MBrepr]={Text->True, OptimizeResult->False, BarnesLemma1->False, BarnesLemma2->False};
Options[BarnesLemma]={Text->True,Shifts->False};

<<MBresolve.m

<<KinematicsGen.m

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

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

invariants = Kinematics[kin, choice]/.p1->k1/.p2->k2/.p3->k3/.p4->k4/.p5->k5

points = {s12->-1,s23->-2,s34->-3,s45->-4,s15->-5};

(*
%gluza@MoNsTeR:~/AMBRE/parallel_PentaBox_000/Equations/ints/int-normal_277.m
FUN[{k1*l2, k2*l2, k5*l1}, {PR[l1, 0, n1]*PR[-k1 + l1, 0, n2]*
   PR[-k1 - k2 + l1, 0, n3]*PR[-k1 - k2 - k3 + l1, 0, n4]*PR[l2, 0, n5]*
   PR[-k5 + l2, 0, n6]*PR[-k4 - k5 + l2, 0, n7]*PR[l1 + l2, 0, n8]}, {l1, l2}]
*)


repr1=MBrepr[{k1*l2, k2*l2, k5*l1},{PR[l1, 0, n1]*PR[-k1 + l1, 0, n2]*PR[-k1 - k2 + l1, 0, n3]*PR[-k1 - k2 - k3 + l1, 0, n4]*PR[l2, 0, n5]*PR[-k5 + l2, 0, n6]*PR[-k4 - k5 + l2, 0, n7]*PR[l1 + l2, 0, n8]},{l2,l1}];

Put[repr1,"PBox_MBrepr.m"];

SetOptions[MBintegrate, Verbose->False, PrecisionGoal->10, AccuracyGoal->24, MaxPoints -> 1000000];

epsilon=0;
powers={n1->1,n2->1,n3->1,n4->1,n5->1,n6->1,n7->1,n8->1};

step1=MBresolve[#/.powers,eps]&/@repr1;
step2=MBexpand[step1,Exp[2*eps*EulerGamma],{eps,0,epsilon}];
exp=MBmerge[step2];

DeleteFile["PBox_Exp_eps0.m"]
Put[exp,"PBox_Exp_eps0.m"]

MBintegrate[exp,points];

Expand[%]//InputForm