<< MBv1.2.m

<<AMBRE.m

<<MBresolve.m

invariants = {p1^2 -> s};

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

repr = MBrepr[{k1*p1, k1*p1, k2*p1}, {PR[k1, 0, n1]*PR[k2, 0, n2]*PR[k2 - k1, 0, n3]*PR[k1 + p1, 0, n4]*PR[k2 + p1, 0, n5]}, {k2, k1}];

step1 = MBresolve[#/.{n1->1,n2->1,n3->1,n4->1,n5->1},eps]&/@repr
step2 = MBexpand[step1,Exp[2*eps*EulerGamma],{eps,0,1}];
MBanalytic=MBmerge[step2];

Print["res=",MBanalytic]

MBintegrate[MBanalytic,{s->-11}]//InputForm