(
   Gam["gamma0"] =
       + ep^(-3)*eM1^3*cf^2*ca*as0pi^3 * ( 57/512 + 93/512*xi - 9/512*
         xi^2 )

       + ep^(-3)*eM1^3*cf^3*as0pi^3 * ( 45/128 + 117/256*xi + 9/256*
         xi^2 )

       + ep^(-3)*eM1^3*nh*cf*ca*tf*as0pi^3 * (  - 1/192 + 1/192*xi )

       + ep^(-3)*eM1^3*nh*cf^2*tf*as0pi^3 * (  - 3/32 - 3/64*xi )

       + ep^(-3)*eM1^3*nl*cf^2*tf*as0pi^3 * (  - 3/64 - 3/64*xi )

       + ep^(-2)*M1^(-1)*eM1^2*cf*ca*as0pi^2*dmf * ( 19/128 + 9/128*xi
          - 3/128*xi^2 )

       + ep^(-2)*M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * ( 33/32 + 33/32*xi + 
         3/64*xi^2 )

       + ep^(-2)*M1^(-1)*eM1^2*nh*cf*tf*as0pi^2*dmf * (  - 1/8 )

       + ep^(-2)*M1^(-1)*eM1^2*nl*cf*tf*as0pi^2*dmf * (  - 1/16 )

       + ep^(-2)*eM1^2*cf^2*as0pi^2 * ( 9/32*xi )

       + ep^(-2)*eM1^3*cf*ca^2*as0pi^3 * ( 143/576 + 11/144*Zeta[2] )

       + ep^(-2)*eM1^3*cf^2*ca*as0pi^3 * ( 3221/2304 + 79/128*xi - 21/
         1024*xi^2 - 271/1152*Zeta[2] - 3/64*Zeta[2]*xi )

       + ep^(-2)*eM1^3*cf^3*as0pi^3 * ( 213/128 + 21/512*xi + 3/128*
         xi^2 + 3/16*Zeta[2] )

       + ep^(-2)*eM1^3*nh*cf*ca*tf*as0pi^3 * (  - 11/72 - 1/64*xi - 1/
         24*Zeta[2] )

       + ep^(-2)*eM1^3*nh*cf^2*tf*as0pi^3 * (  - 35/96 - 17/128*xi + 1/
         6*Zeta[2] )

       + ep^(-2)*eM1^3*nh^2*cf*tf^2*as0pi^3 * ( 1/36 )

       + ep^(-2)*eM1^3*nl*cf*ca*tf*as0pi^3 * (  - 1/6 - 1/36*Zeta[2] )

       + ep^(-2)*eM1^3*nl*cf^2*tf*as0pi^3 * (  - 4/9 - 17/128*xi + 1/9*
         Zeta[2] )

       + ep^(-2)*eM1^3*nl*nh*cf*tf^2*as0pi^3 * ( 1/18 )

       + ep^(-2)*eM1^3*nl^2*cf*tf^2*as0pi^3 * ( 1/36 )

       + ep^(-1)*M1^(-2)*eM1*cf*as0pi*dmf^2 * ( 3/4 + 3/8*xi )

       + ep^(-1)*M1^(-1)*eM1*cf*as0pi*dmf * ( 3/8*xi )

       + ep^(-1)*M1^(-1)*eM1^2*cf*ca*as0pi^2*dmf * ( 83/64 + 17/128*xi
          + 1/256*xi^2 + 3/32*Zeta[2] - 1/16*Zeta[2]*xi )

       + ep^(-1)*M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * ( 45/16 - 17/16*xi - 
         1/32*xi^2 + 1/4*Zeta[2] )

       + ep^(-1)*M1^(-1)*eM1^2*nh*cf*tf*as0pi^2*dmf * (  - 1/6 )

       + ep^(-1)*M1^(-1)*eM1^2*nl*cf*tf*as0pi^2*dmf * (  - 3/8 )

       + ep^(-1)*eM1^2*cf*ca*as0pi^2 * ( 13/64 + 1/16*Zeta[2] )

       + ep^(-1)*eM1^2*cf^2*as0pi^2 * ( 37/64 + 3/16*xi - 1/4*Zeta[2] )

       + ep^(-1)*eM1^2*nh*cf*tf*as0pi^2 * (  - 1/16 )

       + ep^(-1)*eM1^2*nl*cf*tf*as0pi^2 * (  - 1/16 )

       + ep^(-1)*eM1^3*cf*ca^2*as0pi^3 * ( 16135/6912 + 1/1536*xi - 925/
         864*Zeta[2] - 1/120*Zeta[2]^2 - 1/240*Zeta[2]^2*xi + 11/8*
         Zeta[2]*Log[2] - 281/2304*Zeta[3] + 3/512*Zeta[3]*xi )

       + ep^(-1)*eM1^3*cf^2*ca*as0pi^3 * ( 38755/6912 + 3117/1024*xi - 
         281/1024*xi^2 + 13505/27648*Zeta[2] + 47/1024*Zeta[2]*xi - 99/
         1024*Zeta[2]*xi^2 - 1/60*Zeta[2]^2 - 25/32*Zeta[2]*Log[2] + 9/
         16*Zeta[2]*Log[2]*xi - 491/1152*Zeta[3] - 15/64*Zeta[3]*xi )

       + ep^(-1)*eM1^3*cf^3*as0pi^3 * ( 3521/768 + 3247/1024*xi + 31/64
         *xi^2 + 5333/768*Zeta[2] + 999/512*Zeta[2]*xi + 99/512*Zeta[2]
         *xi^2 - 1/12*Zeta[2]^2 - 63/16*Zeta[2]*Log[2] - 9/8*Zeta[2]*
         Log[2]*xi + 81/64*Zeta[3] + 9/32*Zeta[3]*xi )

       + ep^(-1)*eM1^3*nh*cf*ca*tf*as0pi^3 * ( 307/432 + 35/576*xi - 65/
         128*Zeta[2] + 1/128*Zeta[2]*xi - 1/2*Zeta[2]*Log[2] + 1/24*
         Zeta[3] )

       + ep^(-1)*eM1^3*nh*cf^2*tf*as0pi^3 * ( 83/64 - 227/256*xi - 157/
         64*Zeta[2] + 39/128*Zeta[2]*xi + Zeta[2]*Log[2] + 1/12*Zeta[3]
          )

       + ep^(-1)*eM1^3*nh^2*cf*tf^2*as0pi^3 * (  - 16/27 + 1/3*Zeta[2]
          )

       + ep^(-1)*eM1^3*nl*cf*ca*tf*as0pi^3 * (  - 439/288 + 83/216*
         Zeta[2] - 1/2*Zeta[2]*Log[2] + 7/144*Zeta[3] )

       + ep^(-1)*eM1^3*nl*cf^2*tf*as0pi^3 * (  - 3293/1728 - 155/256*xi
          - 3899/3456*Zeta[2] - 33/128*Zeta[2]*xi + Zeta[2]*Log[2] + 1/
         18*Zeta[3] )

       + ep^(-1)*eM1^3*nl*nh*cf*tf^2*as0pi^3 * (  - 19/54 + 1/3*Zeta[2]
          )

       + ep^(-1)*eM1^3*nl^2*cf*tf^2*as0pi^3 * ( 13/54 )

       + M1^(-2)*eM1*cf*as0pi*dmf^2 * ( 1/4 - xi )

       + M1^(-1)*eM1*cf*as0pi*dmf * ( 3/4 - 1/4*xi )

       + M1^(-1)*eM1^2*cf*ca*as0pi^2*dmf * ( 125/32 + 143/128*xi - 79/
         256*xi^2 - 353/128*Zeta[2] + 13/128*Zeta[2]*xi - 15/128*Zeta[2
         ]*xi^2 + 21/8*Zeta[2]*Log[2] - 15/32*Zeta[3] - 1/8*Zeta[3]*xi
          )

       + M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * ( 87/16 + 11/2*xi + 9/16*xi^2
          + 289/32*Zeta[2] + 57/32*Zeta[2]*xi + 15/64*Zeta[2]*xi^2 - 21/
         4*Zeta[2]*Log[2] + 29/16*Zeta[3] )

       + M1^(-1)*eM1^2*nh*cf*tf*as0pi^2*dmf * ( 25/9 - 15/8*Zeta[2] )

       + M1^(-1)*eM1^2*nl*cf*tf*as0pi^2*dmf * (  - 11/8 - 5/16*Zeta[2]
          )

       + eM1^2*cf*ca*as0pi^2 * ( 139/128 - 5/8*Zeta[2] + 3/4*Zeta[2]*
         Log[2] - 1/16*Zeta[3] )

       + eM1^2*cf^2*as0pi^2 * ( 71/128 + 13/8*xi + 3/2*Zeta[2] + 9/32*
         Zeta[2]*xi - 3/2*Zeta[2]*Log[2] - 1/8*Zeta[3] )

       + eM1^2*nh*cf*tf*as0pi^2 * ( 29/32 - 1/2*Zeta[2] )

       + eM1^2*nl*cf*tf*as0pi^2 * (  - 11/32 )

       + eM1^3*cf*ca^2*as0pi^3 * ( 303575/20736 + 13/1536*xi - 120845/
         10368*Zeta[2] + 3/256*Zeta[2]*xi + 54571/11520*Zeta[2]^2 - 17/
         1536*Zeta[2]^2*xi + 11/24*Zeta[2]*Zeta[3] - 1/48*Zeta[2]*Zeta[
         3]*xi + 827/48*Zeta[2]*Log[2] - 179/48*Zeta[2]*Log[2]^2 - 
         19409/3456*Zeta[3] + 13/512*Zeta[3]*xi - 187/384*Zeta[5] - 7/
         768*Zeta[5]*xi - 161/576*Log[2]^4 - 161/24*PolyLog[4,1/2] )

       + eM1^3*cf^2*ca*as0pi^3 * ( 1044841/41472 + 18029/2048*xi - 179/
         1024*xi^2 + 1012271/41472*Zeta[2] - 25/256*Zeta[2]*xi - 231/
         2048*Zeta[2]*xi^2 - 105277/11520*Zeta[2]^2 + 519/640*Zeta[2]^2
         *xi - 55/48*Zeta[2]*Zeta[3] - 1247/96*Zeta[2]*Log[2] + 3*Zeta[
         2]*Log[2]*xi + 217/48*Zeta[2]*Log[2]^2 - 9/8*Zeta[2]*Log[2]^2*
         xi + 29105/13824*Zeta[3] - 709/512*Zeta[3]*xi - 63/512*Zeta[3]
         *xi^2 + 31/24*Zeta[5] + 37/576*Log[2]^4 - 3/32*Log[2]^4*xi + 
         37/24*PolyLog[4,1/2] - 9/4*PolyLog[4,1/2]*xi )

       + eM1^3*cf^3*as0pi^3 * ( 19967/1536 + 957/2048*xi - 1/64*xi^2 + 
         24893/768*Zeta[2] + 4503/1024*Zeta[2]*xi + 33/256*Zeta[2]*xi^2
          - 637/96*Zeta[2]^2 - 189/80*Zeta[2]^2*xi + 71/24*Zeta[2]*
         Zeta[3] - 687/16*Zeta[2]*Log[2] - 6*Zeta[2]*Log[2]*xi + 47/8*
         Zeta[2]*Log[2]^2 + 9/4*Zeta[2]*Log[2]^2*xi + 9311/384*Zeta[3]
          + 915/256*Zeta[3]*xi + 63/256*Zeta[3]*xi^2 - 35/12*Zeta[5] + 
         95/96*Log[2]^4 + 3/16*Log[2]^4*xi + 95/4*PolyLog[4,1/2] + 9/2*
         PolyLog[4,1/2]*xi )

       + eM1^3*nh*cf*ca*tf*as0pi^3 * ( 233/324 - 407/1728*xi - 2881/270
         *Zeta[2] - 3/128*Zeta[2]*xi - 479/240*Zeta[2]^2 + 31/32*Zeta[2
         ]*Zeta[3] + 101/6*Zeta[2]*Log[2] + 4/3*Zeta[2]*Log[2]^2 - 157/
         64*Zeta[3] + 7/192*Zeta[3]*xi - 55/64*Zeta[5] + 1/36*Log[2]^4
          + 2/3*PolyLog[4,1/2] )

       + eM1^3*nh*cf^2*tf*as0pi^3 * (  - 1721/3456 - 1385/512*xi - 3953/
         1728*Zeta[2] + 365/256*Zeta[2]*xi + 47/12*Zeta[2]^2 + 149/24*
         Zeta[2]*Log[2] - 9/4*Zeta[2]*Log[2]*xi - 8/3*Zeta[2]*Log[2]^2
          - 263/24*Zeta[3] + 87/64*Zeta[3]*xi - 1/18*Log[2]^4 - 4/3*
         PolyLog[4,1/2] )

       + eM1^3*nh^2*cf*tf^2*as0pi^3 * (  - 104/81 + 343/360*Zeta[2] - 2
         *Zeta[2]*Log[2] + 11/6*Zeta[3] )

       + eM1^3*nl*cf*ca*tf*as0pi^3 * (  - 4045/432 + 2501/1296*Zeta[2]
          - 1217/720*Zeta[2]^2 - 16/3*Zeta[2]*Log[2] + 4/3*Zeta[2]*Log[
         2]^2 + 451/216*Zeta[3] + 1/9*Log[2]^4 + 8/3*PolyLog[4,1/2] )

       + eM1^3*nl*cf^2*tf*as0pi^3 * (  - 78143/10368 - 929/512*xi - 
         1813/162*Zeta[2] - 187/256*Zeta[2]*xi + 463/90*Zeta[2]^2 + 32/
         3*Zeta[2]*Log[2] - 8/3*Zeta[2]*Log[2]^2 - 9383/1728*Zeta[3] - 
         21/64*Zeta[3]*xi - 2/9*Log[2]^4 - 16/3*PolyLog[4,1/2] )

       + eM1^3*nl*nh*cf*tf^2*as0pi^3 * ( 62/81 + 35/36*Zeta[2] - 2*
         Zeta[2]*Log[2] + 7/6*Zeta[3] )

       + eM1^3*nl^2*cf*tf^2*as0pi^3 * ( 457/324 + 3/8*Zeta[2] )

       + ep*M1^(-2)*eM1*cf*as0pi*dmf^2 * ( 1 + 2*xi + 3/8*Zeta[2] + 3/
         16*Zeta[2]*xi )

       + ep*M1^(-1)*eM1*cf*as0pi*dmf * (  - 1/2 + 3/2*xi + 3/16*Zeta[2]
         *xi )

       + ep*M1^(-1)*eM1^2*cf*ca*as0pi^2*dmf * ( 143/8 + 165/128*xi + 75/
         256*xi^2 - 313/64*Zeta[2] - 11/128*Zeta[2]*xi + 5/256*Zeta[2]*
         xi^2 + 993/160*Zeta[2]^2 - 37/80*Zeta[2]^2*xi + 105/8*Zeta[2]*
         Log[2] - 21/4*Zeta[2]*Log[2]^2 - 1765/192*Zeta[3] + 1/64*Zeta[
         3]*xi - 11/64*Zeta[3]*xi^2 - 7/16*Log[2]^4 - 21/2*PolyLog[4,1/
         2] )

       + ep*M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * ( 241/16 - 55/8*xi - 15/16
         *xi^2 + 409/16*Zeta[2] - 13/16*Zeta[2]*xi - 5/32*Zeta[2]*xi^2
          - 367/40*Zeta[2]^2 - 105/4*Zeta[2]*Log[2] + 21/2*Zeta[2]*Log[
         2]^2 + 175/8*Zeta[3] + 13/16*Zeta[3]*xi + 11/32*Zeta[3]*xi^2
          + 7/8*Log[2]^4 + 21*PolyLog[4,1/2] )

       + ep*M1^(-1)*eM1^2*nh*cf*tf*as0pi^2*dmf * (  - 155/108 - 8/3*
         Zeta[2] + 21/2*Zeta[2]*Log[2] - 145/24*Zeta[3] )

       + ep*M1^(-1)*eM1^2*nl*cf*tf*as0pi^2*dmf * (  - 37/8 - 15/8*Zeta[
         2] - 11/24*Zeta[3] )

       + ep*eM1^2*cf*ca*as0pi^2 * ( 1129/256 - 171/64*Zeta[2] + 163/80*
         Zeta[2]^2 + 9/2*Zeta[2]*Log[2] - 3/2*Zeta[2]*Log[2]^2 - 19/8*
         Zeta[3] - 1/8*Log[2]^4 - 3*PolyLog[4,1/2] )

       + ep*eM1^2*cf^2*as0pi^2 * ( 933/256 + 7/4*xi + 497/64*Zeta[2] + 
         3/16*Zeta[2]*xi - 5*Zeta[2]^2 - 9*Zeta[2]*Log[2] + 3*Zeta[2]*
         Log[2]^2 + 21/4*Zeta[3] - 3/16*Zeta[3]*xi + 1/4*Log[2]^4 + 6*
         PolyLog[4,1/2] )

       + ep*eM1^2*nh*cf*tf*as0pi^2 * ( 63/64 - 25/16*Zeta[2] + 3*Zeta[2
         ]*Log[2] - 7/4*Zeta[3] )

       + ep*eM1^2*nl*cf*tf*as0pi^2 * (  - 89/64 - 5/16*Zeta[2] )

       + ep^2*M1^(-2)*eM1*cf*as0pi*dmf^2 * ( 2 - 4*xi + 1/8*Zeta[2] - 1/
         2*Zeta[2]*xi - 1/4*Zeta[3] - 1/8*Zeta[3]*xi )

       + ep^2*M1^(-1)*eM1*cf*as0pi*dmf * ( 3 - xi + 3/8*Zeta[2] - 1/8*
         Zeta[2]*xi - 1/8*Zeta[3]*xi )

       + 1;
 );
