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

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

       + cf*ca^2*asmpi^3 * ( 2621401/248832 + 24277/13824*lmM1 - 1901/
         2304*lmM1^2 + 121/1152*lmM1^3 - 73919/10368*Zeta[2] + 65/288*
         Zeta[2]*lmM1 - 11/96*Zeta[2]*lmM1^2 + 281/1440*Zeta[2]^2 - 3/
         40*Zeta[2]^2*lmM1 + 31/32*Zeta[2]*Zeta[3] + 283/24*Zeta[2]*
         Log[2] - 1/6*Zeta[2]*Log[2]^2 - 1453/1728*Zeta[3] - 11/192*
         Zeta[3]*lmM1 - 197/128*Zeta[5] - 1/18*Log[2]^4 - 4/3*PolyLog[4
         ,1/2] )

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

       + cf^2*ca*asmpi^3 * ( 117305/82944 + 7211/18432*lmM1 + 61/1536*
         lmM1^2 - 33/512*lmM1^3 + 119443/5184*Zeta[2] - 1193/1152*Zeta[
         2]*lmM1 + 97/192*Zeta[2]*lmM1^2 - 1787/576*Zeta[2]^2 - 1/20*
         Zeta[2]^2*lmM1 - 59/24*Zeta[2]*Zeta[3] - 515/24*Zeta[2]*Log[2]
          + 9/32*Zeta[2]*Log[2]*lmM1 + 7/12*Zeta[2]*Log[2]^2 + 467/216*
         Zeta[3] + 133/384*Zeta[3]*lmM1 + 289/96*Zeta[5] - 1/72*Log[2]^
         4 - 1/3*PolyLog[4,1/2] )

       + cf^3*asmpi^3 * ( 2249/768 + 3051/2048*lmM1 + 51/512*lmM1^2 + 9/
         1024*lmM1^3 - 575/192*Zeta[2] + 77/128*Zeta[2]*lmM1 - 3/16*
         Zeta[2]*lmM1^2 - 251/480*Zeta[2]^2 - 1/4*Zeta[2]^2*lmM1 + 89/
         24*Zeta[2]*Zeta[3] - 17/4*Zeta[2]*Log[2] - 9/16*Zeta[2]*Log[2]
         *lmM1 - 1/2*Zeta[2]*Log[2]^2 + 47/96*Zeta[3] - 21/64*Zeta[3]*
         lmM1 - 355/96*Zeta[5] + 1/4*Log[2]^4 + 6*PolyLog[4,1/2] )

       + nh*cf*tf*asmpi^2 * ( 149/144 - 5/24*lmM1 + 1/8*lmM1^2 - 1/2*
         Zeta[2] )

       + nh*cf*ca*tf*asmpi^3 * ( 6481/3888 - 211/432*lmM1 + 197/288*
         lmM1^2 - 11/72*lmM1^3 - 12839/1080*Zeta[2] - 1/30*Zeta[2]^2 + 
         29/32*Zeta[2]*Zeta[3] + 59/3*Zeta[2]*Log[2] + Zeta[2]*Log[2]^2
          - 983/384*Zeta[3] - 3/4*Zeta[3]*lmM1 - 45/64*Zeta[5] - 1/6*
         Log[2]^4 - 4*PolyLog[4,1/2] )

       + nh*cf^2*tf*asmpi^3 * (  - 29/16 + 103/384*lmM1 + 7/64*lmM1^2
          + 3/64*lmM1^3 + 95/18*Zeta[2] - 3/16*Zeta[2]*lmM1 - 29/120*
         Zeta[2]^2 - 22/3*Zeta[2]*Log[2] - 3/2*Zeta[2]*Log[2]^2 - 163/
         576*Zeta[3] + 3/4*Zeta[3]*lmM1 + 1/4*Log[2]^4 + 6*PolyLog[4,1/
         2] )

       + nh^2*cf*tf^2*asmpi^3 * (  - 4045/3888 - 35/432*lmM1 - 5/72*
         lmM1^2 + 1/36*lmM1^3 - 2/15*Zeta[2] + 11/9*Zeta[3] )

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

       + nl*cf*ca*tf*asmpi^3 * (  - 150227/31104 - 203/432*lmM1 + 245/
         576*lmM1^2 - 11/144*lmM1^3 - 581/648*Zeta[2] - 7/72*Zeta[2]*
         lmM1 + 1/24*Zeta[2]*lmM1^2 - 239/720*Zeta[2]^2 - 25/12*Zeta[2]
         *Log[2] + 1/3*Zeta[2]*Log[2]^2 - 77/216*Zeta[3] - 17/48*Zeta[3
         ]*lmM1 + 1/36*Log[2]^4 + 2/3*PolyLog[4,1/2] )

       + nl*cf^2*tf*asmpi^3 * (  - 9709/5184 - 2527/4608*lmM1 + 37/384*
         lmM1^2 + 3/128*lmM1^3 - 2581/648*Zeta[2] + 85/288*Zeta[2]*lmM1
          - 1/6*Zeta[2]*lmM1^2 + 679/360*Zeta[2]^2 + 25/6*Zeta[2]*Log[2
         ] - 2/3*Zeta[2]*Log[2]^2 - 449/432*Zeta[3] + 7/24*Zeta[3]*lmM1
          - 1/18*Log[2]^4 - 4/3*PolyLog[4,1/2] )

       + nl*nh*cf*tf^2*asmpi^3 * ( 353/1944 - 35/216*lmM1 - 5/36*lmM1^2
          + 1/18*lmM1^3 + 1/12*Zeta[2] - 2/9*Zeta[3] )

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

       + 1;
 );
