<< MBv1.2.m

<<AMBRE.m

<<MBnum.m

<<barnesroutines.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}];

SetOptions[MBnum, Analytical -> True, Numerical -> False, ShowMBrep -> True];

MBanalytic = MBnum[repr, 1, {s -> -11}, {n1 -> 1, n2 -> 1, n3 -> 1, n4 -> 1, n5 -> 1}, 2];

Print["before lemmas"]

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