<<MB/MB.m";
<<MB/MBresolve.m";

F1=-((S^(-5-4 ep-z7)) (T^z7) Gamma[1+z1] Gamma[-1-ep-z1-z2] Gamma[-z2] Gamma[-1-ep-z1-z3] Gamma[-z3] Gamma[1+z1+z2+z3] Gamma[2+ep+z1+z2+z3] Gamma[z10-z4] Gamma[-z1+z4] Gamma[-ep+z1+z2-z4-z5] Gamma[-z5] Gamma[-ep+z1+z3-z4-z6] Gamma[-z6] Gamma[1+z4+z5+z6] Gamma[1+ep-z1-z2-z3+z4+z5+z6] Gamma[-z7] Gamma[1+z7] Gamma[-z10+z7] Gamma[-ep-z10+z4+z5-z8] Gamma[-z8] Gamma[-ep+z10-z7+z8] Gamma[-ep-z10+z4+z6-z9] Gamma[1+ep-z10+z7-z8-z9] Gamma[-z9] Gamma[-ep+z10-z7+z9] Gamma[1+z10+z8+z9] Gamma[1+ep+z10-z4-z5-z6+z8+z9])/(Gamma[-2 ep] Gamma[1-z2] Gamma[1-z3] Gamma[1-2 ep+z1+z2+z3] Gamma[1-z5] Gamma[1-z6] Gamma[1-2 ep+z4+z5+z6] Gamma[1-z8] Gamma[1-z9] Gamma[1-2 ep+z10+z8+z9]);

F1cont=MBresolve[F1,ep]

Length[F1cont]

F1exp=MBexpand[F1select,E^(4EulerGamma ep),{ep,0,0}]//Timing

F1rules=MBoptimizedRules[F1,ep->0,{},{ep}]

F1rules=MBcorrectContours[ MBoptimizedRules[F1,ep->0,{},{ep}] ,10000]

M1cont=MBcontinue[F1,ep->0,F1rules]

2200" integral(s) found"