Example#1, 2-loop planar vertex with three different masses in d=6-2*eps

The same diagram as in the AMBREv1.3_examples.nb file.
All computations are done in a fully automatic way.

<<AMBREv2.1.m

<<PlanarityTestv1.2.m

invariants={p1^2->0,p2^2->0,p1*p2->s/2};

d=6-2 eps;(* by default d=4-2 eps *)

ex=PR[k1,0,n1] PR[k1-k2,m1,n2] PR[k2,m2,n3] PR[k1+p1,m3,n4] PR[k1+p1+p2,0,n5] PR[k2+p1+p2,m2,n6];

PlanarityTest[{ex},{k1,k2},DrawGraph->True];

repr=MBrepr[{1},{ex},{k2,k1}]

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

<<MB.m

<<MBresolve.m

step1=MBresolve[#,eps]&/@finres//Flatten;

step2=MBexpand[step1,Exp[2*eps*EulerGamma],{eps,0,0}];
MBanalytic=MBmerge[step2]

MBintegrate[MBanalytic,{s->-1,m1->1, m2->2,m3->3}]

Quit[]

<<../FIESTA3.2/FIESTA3.m

UsingQLink=False;UsingC=False;

d0=6;

SDEvaluate[UF[{k1,k2},{-k1^2,-(k1-k2)^2+m1^2,-k2^2+m2^2,-(k1+p1)^2+m3^2,-(k1+p1+p2)^2,-(k2+p1+p2)^2+m2^2},{p1^2->0,p2^2->0,p1*p2->s/2,
s->-1,m1->1, m2->2,m3->3}],{1,1,1,1,1,1},0]

Quit[]