Z2q=( 1 + (as*CF*(-1/2 + xi/4))/ep + 
 as^2*(((CF*Tf)/16 + CA*CF*(-5/32 + (5*xi)/64 - xi^2/128))/ep + 
   (-(CF*Tf)/16 + CA*CF*(5/32 - (5*xi)/64 + xi^2/64) + 
     CF^2*(1/8 - xi/8 + xi^2/32))/ep^2) + 
 as^3*(((CF*Tf^2)/48 + Tf*(CF^2*(-7/96 + xi/64) + 
       CA*CF*(-7/64 + (11*xi)/384)) + CA^2*CF*(1/8 - (17*xi)/256 + 
       (9*xi^2)/512 - xi^3/512) + CA*CF^2*(5/64 - (5*xi)/64 + (3*xi^2)/128 - 
       xi^3/512))/ep^2 + 
   (-(CF*Tf^2)/48 + Tf*(CF^2*(1/32 - xi/64) + CA*CF*(9/128 - xi/128)) + 
     CA^2*CF*(-55/768 + (17*xi)/512 - (5*xi^2)/512 + xi^3/768) + 
     CF^3*(-1/48 + xi/32 - xi^2/64 + xi^3/384) + 
     CA*CF^2*(-5/64 + (5*xi)/64 - (7*xi^2)/256 + xi^3/256))/ep^3 + 
   ((CF*Tf^2)/48 + Tf*(CF^2*(5/48 - z3/8) + CA*CF*(41/1536 - (23*xi)/512 + 
         z3/8)) + CF^3*(1/24 - z3/16) + CA*CF^2*(1/192 + z3/64) + 
     CA^2*CF*(-191/3072 + (89*xi)/2048 - (23*xi^2)/2048 + (5*xi^3)/3072 - 
       z3/64 + (xi*z3)/128 - (xi^2*z3)/512))/ep)
);


 Z2sq=( 1 + (as*CF*xi)/(4*ep) + 
 as^2*((-CF^2/16 - (CF*Tf)/16 + CA*CF*((5*xi)/64 - xi^2/128))/ep + 
   ((CF*Tf)/16 + (CF^2*xi^2)/32 + CA*CF*(-1/32 - (5*xi)/64 + xi^2/64))/
    ep^2) + 
 as^3*((-(CF*Tf^2)/48 + Tf*(CF^2*(1/48 - xi/64) + 
       CA*CF*(3/64 + (11*xi)/384)) - (CF^3*xi)/64 + 
     CA^2*CF*(-1/64 - (17*xi)/256 + (9*xi^2)/512 - xi^3/512) + 
     CA*CF^2*(1/32 + (5*xi^2)/256 - xi^3/512))/ep^2 + 
   ((CF*Tf^2)/48 + (CF^3*xi^3)/384 + Tf*(CA*CF*(-7/128 - xi/128) + 
       (CF^2*xi)/64) + CA^2*CF*(17/768 + (17*xi)/512 - (5*xi^2)/512 + 
       xi^3/768) + CA*CF^2*(-xi/128 - (5*xi^2)/256 + xi^3/256))/ep^3 + 
   (CF^3/16 - (CF*Tf^2)/48 + Tf*(CA*CF*(35/512 - (23*xi)/512 - z3/8) + 
       CF^2*(-1/12 + z3/8)) + CA*CF^2*(1/192 - (5*z3)/64) + 
     CA^2*CF*(-95/3072 + (275*xi)/6144 - (23*xi^2)/2048 + (5*xi^3)/3072 + 
       (5*z3)/128 + (xi*z3)/256 - (xi^2*z3)/512))/ep)
);


 Z2gl=( 1 + (as*(-Tf/2 + CA*(-1/4 + xi/4)))/ep + 
 as^2*((((3*CA)/16 - CF/4)*Tf + CA^2*(-3/32 + (5*xi)/64 - xi^2/128))/ep + 
   (CA*Tf*(5/16 - xi/8) + CA^2*(3/32 - (9*xi)/64 + (3*xi^2)/64))/ep^2) + 
 as^3*(((CA*Tf^2)/32 + CA^2*Tf*(-23/128 + (7*xi)/64 - (3*xi^2)/128) + 
     CA^3*(-11/256 + (39*xi)/512 - (21*xi^2)/512 + xi^3/128))/ep^3 + 
   (((-7*CA)/96 - CF/24)*Tf^2 + CA^3*(11/128 - (7*xi)/64 + (5*xi^2)/128 - 
       xi^3/256) + Tf*(CA*CF*(3/16 - xi/16) + 
       CA^2*(-13/192 + (7*xi)/192 + xi^2/256)))/ep^2 + 
   ((CA/24 + CF/8)*Tf^2 + Tf*(CF^2/12 + CA*CF*(1/8 - (3*z3)/8) + 
       CA^2*(-69/512 - (23*xi)/512 + (19*z3)/64)) + 
     CA^3*(-71/3072 + (89*xi)/2048 - (23*xi^2)/2048 + (5*xi^3)/3072 - 
       (3*z3)/128 + (xi*z3)/128 - (xi^2*z3)/512))/ep)
);


 Z3g=( 1 + (as*(-Tf/2 + CA*(1/4 + xi/8)))/ep + 
 as^2*(((-CA/64 - CF/4)*Tf + CA^2*(7/128 + (15*xi)/256 - xi^2/128))/ep + 
   (CA*Tf*(5/32 - xi/16) + CA^2*(-5/64 - xi/128 + xi^2/64))/ep^2) + 
 as^3*(((CA*Tf^2)/64 + CA^2*Tf*(-19/256 + (9*xi)/256 - xi^2/128) + 
     CA^3*(17/512 - (3*xi)/1024 - (3*xi^2)/512 + xi^3/512))/ep^3 + 
   (((-5*CA)/128 - CF/24)*Tf^2 + CA^3*(-167/3072 - (41*xi)/2048 + 
       (17*xi^2)/1024 - (7*xi^3)/3072) + 
     Tf*(CA*CF*(1/8 - xi/32) + CA^2*(137/1536 - (23*xi)/1536 + xi^2/256)))/
    ep^2 + (((19*CA)/384 + CF/8)*Tf^2 + 
     Tf*(CF^2/12 + CA*CF*(-5/48 - (3*z3)/16) + 
       CA^2*(-197/1536 - (35*xi)/1536 + (19*z3)/128)) + 
     CA^3*(95/1536 + (25*xi)/768 - (9*xi^2)/1024 + (7*xi^3)/6144 - 
       (3*z3)/128 + (3*xi*z3)/512 - (xi^2*z3)/1024))/ep)
);


 Z3c=( 1 + (as*CA*(1/8 + xi/16))/ep + 
 as^2*(((-7*CA*Tf)/128 + CA^2*(11/256 - xi/512))/ep + 
   ((3*CA*Tf)/64 + CA^2*(-3/64 - (3*xi)/256 + (3*xi^2)/512))/ep^2) + 
 as^3*(((-7*CA*Tf^2)/384 + Tf*((CA*CF)/32 + CA^2*(35/768 + (5*xi)/6144)) + 
     CA^3*(-13/384 - (19*xi)/4096 + (21*xi^2)/8192 - xi^3/3072))/ep^2 + 
   ((CA*Tf^2)/64 + CA^2*Tf*(-5/128 + xi/1024) + 
     CA^3*(23/1024 + (3*xi)/1024 - (3*xi^2)/1024 + (5*xi^3)/8192))/ep^3 + 
   ((-11*CA*Tf^2)/768 + Tf*(CA^2*(1/32 - xi/128 - (19*z3)/256) + 
       CA*CF*(-1/12 + (3*z3)/32)) + CA^3*(65/6144 + (7*xi)/4096 - 
       (5*xi^2)/4096 + xi^3/4096 + (3*z3)/256 - (3*xi*z3)/1024 + 
       (xi^2*z3)/2048))/ep)
);


 Z3es=( 1 + (as*(-Tf/2 + CA*(1/4 + xi/4)))/ep + 
 as^2*(((-CA/16 - CF/4)*Tf + CA^2*(1/16 + (5*xi)/64 - xi^2/128))/ep + 
   (CA*Tf*(3/16 - xi/8) + CA^2*(-3/32 - xi/64 + (3*xi^2)/64))/ep^2) + 
 as^3*(((CA*Tf^2)/96 + CA^2*Tf*(-11/128 + (5*xi)/64 - (3*xi^2)/128) + 
     CA^3*(31/768 - (5*xi)/512 - (9*xi^2)/512 + xi^3/128))/ep^3 + 
   ((-CA/32 - CF/24)*Tf^2 + CA^3*(-1/16 - xi/32 + (9*xi^2)/256 - xi^3/256) + 
     Tf*(CA*CF*(7/48 - xi/16) + CA^2*(5/48 - (5*xi)/192 + xi^2/256)))/ep^2 + 
   ((CA/24 + CF/8)*Tf^2 + Tf*(CF^2/12 + CA*CF*(-1/6 - z3/8) + 
       CA^2*(-45/512 - (23*xi)/512 + (3*z3)/32)) + 
     CA^3*(209/3072 + (275*xi)/6144 - (23*xi^2)/2048 + (5*xi^3)/3072 - 
       (3*z3)/128 + (xi*z3)/256 - (xi^2*z3)/512))/ep)
);


 Zglsqq=( 1 + (as*(CA*(-1/2 + xi/8) + CF*(-1/4 + xi/4)))/ep + 
 as^2*((-CF^2/32 + (5*CA*Tf)/32 + CA*CF*(-5/64 + (5*xi)/64 - xi^2/128) + 
     CA^2*(-9/64 + (5*xi)/128 - xi^2/256))/ep + 
   ((-3*CA*Tf)/32 + CA^2*(19/64 - (13*xi)/128 + xi^2/64) + 
     CF^2*(1/32 - xi/16 + xi^2/32) + CA*CF*(3/16 - (15*xi)/64 + (3*xi^2)/64))/
    ep^2) + 
 as^3*(((5*CA*Tf^2)/96 + Tf*(-CF^2/96 + CA^2*(-91/384 + (13*xi)/384) + 
       CA*CF*(-17/128 + (13*xi)/192)) + CF^3*(1/128 - xi/128) + 
     CA^2*CF*(33/256 - (41*xi)/256 + (43*xi^2)/1024 - xi^3/256) + 
     CA*CF^2*(13/256 - (11*xi)/256 + (11*xi^2)/512 - xi^3/512) + 
     CA^3*(25/128 - (9*xi)/128 + xi^2/64 - (3*xi^3)/2048))/ep^2 + 
   (-(CA*Tf^2)/32 + Tf*(CA*CF*(1/32 - xi/32) + CA^2*(37/256 - xi/64)) + 
     CA^3*(-97/512 + (75*xi)/1024 - (9*xi^2)/512 + xi^3/512) + 
     CF^3*(-1/384 + xi/128 - xi^2/128 + xi^3/384) + 
     CA^2*CF*(-25/192 + (23*xi)/128 - (29*xi^2)/512 + (11*xi^3)/1536) + 
     CA*CF^2*(-1/32 + (9*xi)/128 - (3*xi^2)/64 + xi^3/128))/ep^3 + 
   ((CA*Tf^2)/96 + Tf*(CF^2/96 + CA*CF*(377/1536 - (23*xi)/512 - (3*z3)/16) + 
       CA^2*(-47/3072 - (23*xi)/1024 + (19*z3)/128)) + 
     CA*CF^2*(1/192 - z3/32) + CF^3*(5/96 - z3/32) + 
     CA^2*CF*(-143/3072 + (271*xi)/6144 - (23*xi^2)/2048 + (5*xi^3)/3072 + 
       (3*z3)/256 + (3*xi*z3)/512 - (xi^2*z3)/512) + 
     CA^3*(-407/6144 + (89*xi)/4096 - (23*xi^2)/4096 + (5*xi^3)/6144 - 
       (3*z3)/256 + (xi*z3)/256 - (xi^2*z3)/1024))/ep)
);


 Zgsqsq=( 1 + (as*(CA*(-1/4 + xi/16) + (CF*xi)/4))/ep + 
 as^2*((-CF^2/16 + ((7*CA)/128 - CF/16)*Tf + 
     CA*CF*((5*xi)/64 - xi^2/128) + CA^2*(-17/256 + (15*xi)/512 - xi^2/256))/
    ep + (((-3*CA)/64 + CF/16)*Tf + (CF^2*xi^2)/32 + 
     CA^2*(15/128 - (9*xi)/256 + (3*xi^2)/512) + 
     CA*CF*(-1/32 - (9*xi)/64 + xi^2/32))/ep^2) + 
 as^3*((((7*CA)/384 - CF/48)*Tf^2 + Tf*(CF^2*(1/48 - xi/64) + 
       CA^2*(-241/3072 + (71*xi)/6144) + CA*CF*(1/32 + (59*xi)/1536)) - 
     (CF^3*xi)/64 + CA^2*CF*(-1/64 - (105*xi)/1024 + (65*xi^2)/2048 - 
       (7*xi^3)/2048) + CA*CF^2*(3/64 - xi/256 + (5*xi^2)/256 - xi^3/512) + 
     CA^3*(469/6144 - (33*xi)/1024 + (69*xi^2)/8192 - (11*xi^3)/12288))/
    ep^2 + ((-CA/64 + CF/48)*Tf^2 + (CF^3*xi^3)/384 + 
     Tf*(CA*CF*(-9/128 - xi/64) + CA^2*(31/512 - (5*xi)/1024) + 
       (CF^2*xi)/64) + CA^3*(-67/1024 + (21*xi)/1024 - (21*xi^2)/4096 + 
       (5*xi^3)/8192) + CA^2*CF*(23/768 + (41*xi)/512 - (7*xi^2)/256 + 
       (23*xi^3)/6144) + CA*CF^2*(-xi/128 - (7*xi^2)/256 + (3*xi^3)/512))/
    ep^3 + (CF^3/16 + ((11*CA)/768 - CF/48)*Tf^2 + 
     Tf*(CA*CF*(233/1536 - (23*xi)/512 - (7*z3)/32) + 
       CA^2*(-37/3072 - (35*xi)/3072 + (19*z3)/256) + CF^2*(-1/12 + z3/8)) + 
     CA*CF^2*(1/192 - (5*z3)/64) + CA^2*CF*(-95/3072 + (275*xi)/6144 - 
       (23*xi^2)/2048 + (5*xi^3)/3072 + (5*z3)/128 + (xi*z3)/256 - 
       (xi^2*z3)/512) + CA^3*(-73/3072 + (25*xi)/1536 - (9*xi^2)/2048 + 
       (7*xi^3)/12288 - (3*z3)/256 + (3*xi*z3)/1024 - (xi^2*z3)/2048))/ep)
);


 Zgqq=( 1 + (as*(CA*(-1/4 + xi/16) + CF*(-1/2 + xi/4)))/ep + 
 as^2*((((7*CA)/128 + CF/16)*Tf + CA*CF*(-5/32 + (5*xi)/64 - xi^2/128) + 
     CA^2*(-17/256 + (15*xi)/512 - xi^2/256))/ep + 
   (((-3*CA)/64 - CF/16)*Tf + CA^2*(15/128 - (9*xi)/256 + (3*xi^2)/512) + 
     CA*CF*(9/32 - (11*xi)/64 + xi^2/32) + CF^2*(1/8 - xi/8 + xi^2/32))/
    ep^2) + as^3*((((7*CA)/384 + CF/48)*Tf^2 + 
     Tf*(CA^2*(-241/3072 + (71*xi)/6144) + CF^2*(-7/96 + xi/64) + 
       CA*CF*(-47/256 + (71*xi)/1536)) + CA^2*CF*(101/512 - (65*xi)/512 + 
       (69*xi^2)/2048 - (7*xi^3)/2048) + 
     CA*CF^2*(5/64 - (5*xi)/64 + (3*xi^2)/128 - xi^3/512) + 
     CA^3*(469/6144 - (33*xi)/1024 + (69*xi^2)/8192 - (11*xi^3)/12288))/
    ep^2 + ((-CA/64 - CF/48)*Tf^2 + Tf*(CA*CF*(7/64 - (3*xi)/128) + 
       CF^2*(1/32 - xi/64) + CA^2*(31/512 - (5*xi)/1024)) + 
     CA^3*(-67/1024 + (21*xi)/1024 - (21*xi^2)/4096 + (5*xi^3)/8192) + 
     CF^3*(-1/48 + xi/32 - xi^2/64 + xi^3/384) + 
     CA^2*CF*(-65/384 + (7*xi)/64 - (31*xi^2)/1024 + (23*xi^3)/6144) + 
     CA*CF^2*(-7/64 + (15*xi)/128 - (11*xi^2)/256 + (3*xi^3)/512))/ep^3 + 
   (((11*CA)/768 + CF/48)*Tf^2 + Tf*(CF^2*(5/48 - z3/8) + 
       CA*CF*(169/1536 - (23*xi)/512 + z3/32) + 
       CA^2*(-37/3072 - (35*xi)/3072 + (19*z3)/256)) + CF^3*(1/24 - z3/16) + 
     CA*CF^2*(1/192 + z3/64) + CA^2*CF*(-191/3072 + (89*xi)/2048 - 
       (23*xi^2)/2048 + (5*xi^3)/3072 - z3/64 + (xi*z3)/128 - 
       (xi^2*z3)/512) + CA^3*(-73/3072 + (25*xi)/1536 - (9*xi^2)/2048 + 
       (7*xi^3)/12288 - (3*z3)/256 + (3*xi*z3)/1024 - (xi^2*z3)/2048))/ep)
);


 Zepqq=( 1 + (as*(CA*(-1/4 + xi/8) + CF*(-1/2 + xi/4)))/ep + 
 as^2*(((CA/32 + CF/16)*Tf + CA*CF*(-5/32 + (5*xi)/64 - xi^2/128) + 
     CA^2*(-1/16 + (5*xi)/128 - xi^2/256))/ep + 
   ((-CA/32 - CF/16)*Tf + CA^2*(7/64 - (9*xi)/128 + xi^2/64) + 
     CF^2*(1/8 - xi/8 + xi^2/32) + CA*CF*(9/32 - (13*xi)/64 + (3*xi^2)/64))/
    ep^2) + 
 as^3*(((CA/96 + CF/48)*Tf^2 + Tf*(CF^2*(-7/96 + xi/64) + 
       CA^2*(-11/192 + (7*xi)/384) + CA*CF*(-31/192 + (17*xi)/384)) + 
     CA^2*CF*(25/128 - (9*xi)/64 + (21*xi^2)/512 - xi^3/256) + 
     CA*CF^2*(5/64 - (5*xi)/64 + (3*xi^2)/128 - xi^3/512) + 
     CA^3*(9/128 - (13*xi)/256 + (15*xi^2)/1024 - (3*xi^3)/2048))/ep^2 + 
   ((-CA/96 - CF/48)*Tf^2 + Tf*(CA*CF*(13/128 - (3*xi)/128) + 
       CF^2*(1/32 - xi/64) + CA^2*(11/256 - xi/128)) + 
     CA^3*(-89/1536 + (41*xi)/1024 - (7*xi^2)/512 + xi^3/512) + 
     CF^3*(-1/48 + xi/32 - xi^2/64 + xi^3/384) + 
     CA^2*CF*(-127/768 + (69*xi)/512 - (25*xi^2)/512 + (11*xi^3)/1536) + 
     CA*CF^2*(-7/64 + xi/8 - (13*xi^2)/256 + xi^3/128))/ep^3 + 
   ((CA/96 + CF/48)*Tf^2 + Tf*(CF^2*(5/48 - z3/8) + 
       CA^2*(25/3072 - (23*xi)/1024 + (3*z3)/64) + 
       CA*CF*(121/1536 - (23*xi)/512 + z3/16)) + CF^3*(1/24 - z3/16) + 
     CA*CF^2*(1/192 + z3/64) + CA^2*CF*(-191/3072 + (89*xi)/2048 - 
       (23*xi^2)/2048 + (5*xi^3)/3072 - z3/64 + (xi*z3)/128 - 
       (xi^2*z3)/512) + CA^3*(-127/6144 + (275*xi)/12288 - (23*xi^2)/4096 + 
       (5*xi^3)/6144 - (3*z3)/256 + (xi*z3)/512 - (xi^2*z3)/1024))/ep)
);


 Zgglgl=( 1 + (as*(-Tf/2 + CA*(-1/2 + (5*xi)/16)))/ep +
 as^2*((((31*CA)/128 - CF/4)*Tf + CA^2*(-41/256 + (55*xi)/512 -
       (3*xi^2)/256))/ep + (CA*Tf*(25/64 - (5*xi)/32) +
     CA^2*(35/128 - (65*xi)/256 + (35*xi^2)/512))/ep^2) +
 as^3*(((5*CA*Tf^2)/128 + CA^2*Tf*(-125/512 + (165*xi)/1024 -
       (35*xi^2)/1024) + CA^3*(-165/1024 + (45*xi)/256 - (315*xi^2)/4096 +
       (105*xi^3)/8192))/ep^3 + (((-21*CA)/256 - CF/24)*Tf^2 +
     CA^3*(1243/6144 - (391*xi)/2048 + (513*xi^2)/8192 - (77*xi^3)/12288) +
     Tf*(CA*CF*(7/32 - (5*xi)/64) + CA^2*(-533/3072 + (361*xi)/6144 +
         (3*xi^2)/512)))/ep^2 + (((43*CA)/768 + CF/8)*Tf^2 +
     Tf*(CF^2/12 + CA*CF*(5/24 - (15*z3)/32) +
       CA^2*(-451/3072 - (173*xi)/3072 + (95*z3)/256)) +
     CA^3*(-3/64 + (367*xi)/6144 - xi^2/64 + (9*xi^3)/4096 - (9*z3)/256 +
       (11*xi*z3)/1024 - (5*xi^2*z3)/2048))/ep)
);


 Zepglgl=( 1 + (as*(-Tf/2 + CA*(-1/2 + (3*xi)/8)))/ep + 
 as^2*((((7*CA)/32 - CF/4)*Tf + CA^2*(-5/32 + (15*xi)/128 - (3*xi^2)/256))/
    ep + (CA*Tf*(13/32 - (3*xi)/16) + CA^2*(17/64 - (39*xi)/128 + 
       (3*xi^2)/32))/ep^2) + 
 as^3*(((7*CA*Tf^2)/192 + CA^2*Tf*(-67/256 + (51*xi)/256 - (3*xi^2)/64) + 
     CA^3*(-233/1536 + (213*xi)/1024 - (27*xi^2)/256 + (5*xi^3)/256))/ep^3 + 
   (((-5*CA)/64 - CF/24)*Tf^2 + CA^3*(25/128 - (111*xi)/512 + (39*xi^2)/512 - 
       (15*xi^3)/2048) + Tf*(CA*CF*(11/48 - (3*xi)/32) + 
       CA^2*(-19/128 + (17*xi)/256 + (3*xi^2)/512)))/ep^2 + 
   (((5*CA)/96 + CF/8)*Tf^2 + Tf*(CF^2/12 + CA*CF*(17/96 - (7*z3)/16) + 
       CA^2*(-389/3072 - (69*xi)/1024 + (11*z3)/32)) + 
     CA^3*(-269/6144 + (809*xi)/12288 - (69*xi^2)/4096 + (5*xi^3)/2048 - 
       (9*z3)/256 + (5*xi*z3)/512 - (3*xi^2*z3)/1024))/ep)
);


 Zggg=( 1 + (as*(-Tf/2 + (3*CA*xi)/16))/ep +
 as^2*((((5*CA)/128 - CF/4)*Tf + CA^2*(-3/256 + (45*xi)/512 -
       (3*xi^2)/256))/ep + (CA*Tf*(15/64 - (3*xi)/32) +
     CA^2*(-3/128 - (15*xi)/256 + (15*xi^2)/512))/ep^2) +
 as^3*(((3*CA*Tf^2)/128 + CA^2*Tf*(-63/512 + (69*xi)/1024 -
       (15*xi^2)/1024) + CA^3*(17/1024 + (21*xi)/1024 - (75*xi^2)/4096 +
       (35*xi^3)/8192))/ep^3 + (((-37*CA)/768 - CF/24)*Tf^2 +
     CA^3*(-17/2048 - (33*xi)/512 + (273*xi^2)/8192 - (17*xi^3)/4096) +
     Tf*(CA*CF*(5/32 - (3*xi)/64) + CA^2*(63/1024 - (25*xi)/2048 +
         (3*xi^2)/512)))/ep^2 + (((49*CA)/768 + CF/8)*Tf^2 +
     Tf*(CF^2/12 + CA*CF*(-1/48 - (9*z3)/32) +
       CA^2*(-431/3072 - (35*xi)/1024 + (57*z3)/256)) +
     CA^3*(39/1024 + (25*xi)/512 - (27*xi^2)/2048 + (7*xi^3)/4096 -
       (9*z3)/256 + (9*xi*z3)/1024 - (3*xi^2*z3)/2048))/ep)
);


 Zgcc=( 1 + (as*CA*(-1/8 + xi/8))/ep + 
 as^2*((CA^2*(-3/128 + (7*xi)/256 - xi^2/256))/ep + 
   (CA^2*(5/128 - (7*xi)/128 + xi^2/64))/ep^2) + 
 as^3*((CA^2*Tf*(-19/1536 + (19*xi)/1536) + 
     CA^3*(3/128 - (35*xi)/1024 + (25*xi^2)/2048 - (3*xi^3)/2048))/ep^2 + 
   (CA^2*Tf*(59/3072 - (59*xi)/3072) + CA^3*(-27/2048 + (221*xi)/12288 - 
       (23*xi^2)/4096 + (5*xi^3)/6144))/ep + 
   (CA^2*Tf*(1/256 - xi/256) + CA^3*(-17/1024 + (27*xi)/1024 - (3*xi^2)/256 + 
       xi^3/512))/ep^3)
);


 Zgepep=( 1 + (as*(-Tf/2 + (5*CA*xi)/16))/ep + 
 as^2*(((-CA/128 - CF/4)*Tf + CA^2*(-1/256 + (55*xi)/512 - (3*xi^2)/256))/
    ep + (CA*Tf*(17/64 - (5*xi)/32) + CA^2*(-5/128 - (25*xi)/256 + 
       (35*xi^2)/512))/ep^2) + 
 as^3*(((7*CA*Tf^2)/384 + CA^2*Tf*(-73/512 + (125*xi)/1024 - 
       (35*xi^2)/1024) + CA^3*(85/3072 + (15*xi)/512 - (175*xi^2)/4096 + 
       (105*xi^3)/8192))/ep^3 + (((-31*CA)/768 - CF/24)*Tf^2 + 
     CA^3*(-113/6144 - (181*xi)/2048 + (465*xi^2)/8192 - (77*xi^3)/12288) + 
     Tf*(CA*CF*(17/96 - (5*xi)/64) + CA^2*(271/3072 - (119*xi)/6144 + 
         (3*xi^2)/512)))/ep^2 + (((43*CA)/768 + CF/8)*Tf^2 + 
     Tf*(CF^2/12 + CA*CF*(-1/12 - (7*z3)/32) + 
       CA^2*(-307/3072 - (173*xi)/3072 + (43*z3)/256)) + 
     CA^3*(17/384 + (125*xi)/2048 - xi^2/64 + (9*xi^3)/4096 - (9*z3)/256 + 
       (7*xi*z3)/1024 - (5*xi^2*z3)/2048))/ep)
);


 Zmgl=( 1 + (as*((-3*CA)/4 + Tf/2))/ep + 
 as^2*(((-3*CA^2)/8 + (CA/4 + CF/2)*Tf)/ep + 
   ((9*CA^2)/16 - (3*CA*Tf)/4 + Tf^2/4)/ep^2) + 
 as^3*(((-21*CA^3)/64 + ((5*CA^2)/16 + (13*CA*CF)/16 - CF^2/4)*Tf + 
     (-CA/16 - (3*CF)/8)*Tf^2)/ep + 
   ((9*CA^3)/16 + ((-3*CA^2)/4 - (3*CA*CF)/4)*Tf + (CA/4 + CF/2)*Tf^2)/ep^2 + 
   ((-27*CA^3)/64 + (27*CA^2*Tf)/32 - (9*CA*Tf^2)/16 + Tf^3/8)/ep^3)
);


 Zmq=( 1 - (as*CF)/(2*ep) + as^2*(((3*CA*CF)/16 + CF^2/8 - (CF*Tf)/8)/ep^2 + 
   ((-3*CA*CF)/16 + CF^2/8 + (CF*Tf)/8)/ep) + 
 as^3*(((-3*CA^2*CF)/32 - (3*CA*CF^2)/32 - CF^3/48 + 
     ((CA*CF)/8 + CF^2/16)*Tf - (CF*Tf^2)/24)/ep^3 + 
   ((5*CA^2*CF)/32 + (CA*CF^2)/32 - CF^3/16 + (-(CA*CF)/6 - (5*CF^2)/48)*Tf + 
     (CF*Tf^2)/24)/ep^2 + (-(CA^2*CF)/16 + (CA*CF^2)/16 - CF^3/12 + 
     (CF*Tf^2)/24 + Tf*(CF^2*(1/6 - z3/4) + CA*CF*(-1/48 + z3/4)))/ep)
);


 
