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

       + cf*ca*asmpi^2 * (  - 6821/2304 - 49/192*lmM1 + 11/64*lmM1^2 + 
         21/16*Zeta[2] + 1/8*Zeta[2]*lmM1 - 7/4*Zeta[2]*Log[2] + 9/16*
         Zeta[3] )

       + cf*ca^2*asmpi^3 * (  - 2333689/248832 + 1451/13824*lmM1 + 845/
         2304*lmM1^2 - 121/1152*lmM1^3 - 6905/10368*Zeta[2] + 65/288*
         Zeta[2]*lmM1 - 11/96*Zeta[2]*lmM1^2 - 223/120*Zeta[2]^2 - 3/40
         *Zeta[2]^2*lmM1 + 127/24*Zeta[2]*Zeta[3] - 793/72*Zeta[2]*Log[
         2] + 23/9*Zeta[2]*Log[2]^2 + 10721/1728*Zeta[3] - 11/192*Zeta[
         3]*lmM1 - 1051/384*Zeta[5] + 37/216*Log[2]^4 + 37/9*PolyLog[4,
         1/2] )

       + cf^2*asmpi^2 * ( 1453/768 + 29/64*lmM1 + 9/128*lmM1^2 - 173/48
         *Zeta[2] - 1/2*Zeta[2]*lmM1 + 7/2*Zeta[2]*Log[2] - 11/8*Zeta[3
         ] )

       + cf^2*ca*asmpi^3 * ( 660965/82944 + 22547/18432*lmM1 - 227/1536
         *lmM1^2 - 33/512*lmM1^3 + 38611/5184*Zeta[2] - 1895/1152*Zeta[
         2]*lmM1 + 79/192*Zeta[2]*lmM1^2 + 13201/2880*Zeta[2]^2 - 1/20*
         Zeta[2]^2*lmM1 - 109/12*Zeta[2]*Zeta[3] - 139/18*Zeta[2]*Log[2
         ] + 21/32*Zeta[2]*Log[2]*lmM1 - 205/36*Zeta[2]*Log[2]^2 - 5255/
         432*Zeta[3] + 61/384*Zeta[3]*lmM1 + 329/96*Zeta[5] - 35/216*
         Log[2]^4 - 35/9*PolyLog[4,1/2] )

       + cf^3*asmpi^3 * (  - 2395/512 - 1909/2048*lmM1 - 51/512*lmM1^2
          - 9/1024*lmM1^3 - 6331/192*Zeta[2] + 301/128*Zeta[2]*lmM1 + 3/
         16*Zeta[2]*lmM1^2 - 433/160*Zeta[2]^2 - 1/4*Zeta[2]^2*lmM1 + 8/
         3*Zeta[2]*Zeta[3] + 119/2*Zeta[2]*Log[2] - 21/16*Zeta[2]*Log[2
         ]*lmM1 + 7/6*Zeta[2]*Log[2]^2 - 265/96*Zeta[3] + 15/64*Zeta[3]
         *lmM1 - 95/96*Zeta[5] - 13/36*Log[2]^4 - 26/3*PolyLog[4,1/2] )

       + nh*cf*tf*asmpi^2 * ( 133/96 - 5/6*Zeta[2] )

       + nh*cf*ca*tf*asmpi^3 * ( 49145/5184 - 4799/270*Zeta[2] - 34/45*
         Zeta[2]^2 + 35/96*Zeta[2]*Zeta[3] + 197/9*Zeta[2]*Log[2] + 8/9
         *Zeta[2]*Log[2]^2 - 815/384*Zeta[3] - 35/192*Zeta[5] - 4/27*
         Log[2]^4 - 32/9*PolyLog[4,1/2] )

       + nh*cf^2*tf*asmpi^3 * (  - 2987/576 - 133/256*lmM1 - 316/45*
         Zeta[2] + 5/16*Zeta[2]*lmM1 + 379/360*Zeta[2]^2 + 74/9*Zeta[2]
         *Log[2] - 23/18*Zeta[2]*Log[2]^2 + 1445/576*Zeta[3] + 23/108*
         Log[2]^4 + 46/9*PolyLog[4,1/2] )

       + nh^2*cf*tf^2*asmpi^3 * (  - 3641/2592 + 22/135*Zeta[2] + 17/18
         *Zeta[3] )

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

       + nl*cf*ca*tf*asmpi^3 * ( 153347/31104 - 4/27*lmM1 - 149/576*
         lmM1^2 + 11/144*lmM1^3 + 19/162*Zeta[2] - 7/72*Zeta[2]*lmM1 + 
         1/24*Zeta[2]*lmM1^2 + 31/240*Zeta[2]^2 + 169/36*Zeta[2]*Log[2]
          - 7/9*Zeta[2]*Log[2]^2 - 199/432*Zeta[3] + 19/48*Zeta[3]*lmM1
          - 7/108*Log[2]^4 - 14/9*PolyLog[4,1/2] )

       + nl*cf^2*tf*asmpi^3 * (  - 1145/10368 + 65/4608*lmM1 + 1/384*
         lmM1^2 + 3/128*lmM1^3 + 2383/324*Zeta[2] + 85/288*Zeta[2]*lmM1
          - 1/6*Zeta[2]*lmM1^2 - 95/72*Zeta[2]^2 - 169/18*Zeta[2]*Log[2
         ] + 14/9*Zeta[2]*Log[2]^2 + 1741/432*Zeta[3] - 11/24*Zeta[3]*
         lmM1 + 7/54*Log[2]^4 + 28/9*PolyLog[4,1/2] )

       + nl*nh*cf*tf^2*asmpi^3 * (  - 2545/1296 + 127/108*Zeta[2] )

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

       + 1;
 );
