(
   Co["gamma0"] =
       + cf*asmpi * (  - 1/2 - 3/8*lmM1 )

       + cf*ca*asmpi^2 * (  - 871/768 - 49/192*lmM1 + 11/64*lmM1^2 + 5/
         16*Zeta[2] + 1/8*Zeta[2]*lmM1 - 3/4*Zeta[2]*Log[2] + 5/16*
         Zeta[3] )

       + cf*ca^2*asmpi^3 * (  - 23689/248832 + 1451/13824*lmM1 + 845/
         2304*lmM1^2 - 121/1152*lmM1^3 - 109109/10368*Zeta[2] + 65/288*
         Zeta[2]*lmM1 - 11/96*Zeta[2]*lmM1^2 - 601/360*Zeta[2]^2 - 3/40
         *Zeta[2]^2*lmM1 + 23/4*Zeta[2]*Zeta[3] + 53/24*Zeta[2]*Log[2]
          + 5/3*Zeta[2]*Log[2]^2 + 6605/1728*Zeta[3] - 11/192*Zeta[3]*
         lmM1 - 457/128*Zeta[5] + 7/72*Log[2]^4 + 7/3*PolyLog[4,1/2] )

       + cf^2*asmpi^2 * ( 255/256 + 17/64*lmM1 + 9/128*lmM1^2 - 15/16*
         Zeta[2] - 1/2*Zeta[2]*lmM1 + 3/2*Zeta[2]*Log[2] - 7/8*Zeta[3]
          )

       + cf^2*ca*asmpi^3 * ( 88811/82944 + 7571/18432*lmM1 - 95/1536*
         lmM1^2 - 33/512*lmM1^3 + 149197/5184*Zeta[2] - 1391/1152*Zeta[
         2]*lmM1 + 79/192*Zeta[2]*lmM1^2 + 1333/576*Zeta[2]^2 - 1/20*
         Zeta[2]^2*lmM1 - 115/12*Zeta[2]*Zeta[3] - 613/24*Zeta[2]*Log[2
         ] + 9/32*Zeta[2]*Log[2]*lmM1 - 55/12*Zeta[2]*Log[2]^2 - 4873/
         864*Zeta[3] + 97/384*Zeta[3]*lmM1 + 559/96*Zeta[5] - 5/72*Log[
         2]^4 - 5/3*PolyLog[4,1/2] )

       + cf^3*asmpi^3 * (  - 301/128 - 1141/2048*lmM1 - 33/512*lmM1^2
          - 9/1024*lmM1^3 - 4613/192*Zeta[2] + 141/128*Zeta[2]*lmM1 + 3/
         16*Zeta[2]*lmM1^2 - 611/480*Zeta[2]^2 - 1/4*Zeta[2]^2*lmM1 + 
         10/3*Zeta[2]*Zeta[3] + 169/4*Zeta[2]*Log[2] - 9/16*Zeta[2]*
         Log[2]*lmM1 + 5/2*Zeta[2]*Log[2]^2 - 463/96*Zeta[3] + 3/64*
         Zeta[3]*lmM1 - 295/96*Zeta[5] - 1/4*Log[2]^4 - 6*PolyLog[4,1/2
         ] )

       + nh*cf*tf*asmpi^2 * ( 727/288 - 3/2*Zeta[2] )

       + nh*cf*ca*tf*asmpi^3 * ( 6457/576 - 8261/270*Zeta[2] - 22/15*
         Zeta[2]^2 + 53/32*Zeta[2]*Zeta[3] + 41*Zeta[2]*Log[2] + 4/3*
         Zeta[2]*Log[2]^2 - 3821/1152*Zeta[3] - 85/64*Zeta[5] - 2/9*
         Log[2]^4 - 16/3*PolyLog[4,1/2] )

       + nh*cf^2*tf*asmpi^3 * (  - 19/16 - 727/768*lmM1 + 5/9*Zeta[2]
          + 9/16*Zeta[2]*lmM1 + 51/40*Zeta[2]^2 - 2*Zeta[2]*Log[2] - 13/
         6*Zeta[2]*Log[2]^2 - 1435/576*Zeta[3] + 13/36*Log[2]^4 + 26/3*
         PolyLog[4,1/2] )

       + nh^2*cf*tf^2*asmpi^3 * (  - 5857/2592 + 2/45*Zeta[2] + 11/6*
         Zeta[3] )

       + nl*cf*tf*asmpi^2 * ( 47/192 + 5/48*lmM1 - 1/16*lmM1^2 + 1/4*
         Zeta[2] )

       + nl*cf*ca*tf*asmpi^3 * (  - 8701/31104 - 4/27*lmM1 - 149/576*
         lmM1^2 + 11/144*lmM1^3 + 497/324*Zeta[2] - 7/72*Zeta[2]*lmM1
          + 1/24*Zeta[2]*lmM1^2 - 11/720*Zeta[2]^2 + 19/12*Zeta[2]*Log[
         2] - 1/3*Zeta[2]*Log[2]^2 + 113/432*Zeta[3] + 19/48*Zeta[3]*
         lmM1 - 1/36*Log[2]^4 - 2/3*PolyLog[4,1/2] )

       + nl*cf^2*tf*asmpi^3 * ( 395/324 + 1217/4608*lmM1 - 11/384*
         lmM1^2 + 3/128*lmM1^3 + 551/648*Zeta[2] + 85/288*Zeta[2]*lmM1
          - 1/6*Zeta[2]*lmM1^2 - 7/72*Zeta[2]^2 - 19/6*Zeta[2]*Log[2]
          + 2/3*Zeta[2]*Log[2]^2 + 541/432*Zeta[3] - 11/24*Zeta[3]*lmM1
          + 1/18*Log[2]^4 + 4/3*PolyLog[4,1/2] )

       + nl*nh*cf*tf^2*asmpi^3 * (  - 193/144 + 29/36*Zeta[2] )

       + nl^2*cf*tf^2*asmpi^3 * ( 1751/15552 + 35/864*lmM1 + 5/144*
         lmM1^2 - 1/72*lmM1^3 - 13/36*Zeta[2] - 7/36*Zeta[3] )

       + 1;
 );
