(
   Gam["scalar"] =
       + ep^(-3)*eM1^3*cf*ca^2*as0pi^3 * ( 121/576 )

       + ep^(-3)*eM1^3*cf^2*ca*as0pi^3 * ( 453/512 + 93/512*xi - 9/512*
         xi^2 )

       + ep^(-3)*eM1^3*cf^3*as0pi^3 * ( 81/64 + 171/256*xi + 9/256*xi^2
          )

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

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

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

       + ep^(-3)*eM1^3*nl*cf*ca*tf*as0pi^3 * (  - 11/72 )

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

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

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

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

       + ep^(-2)*M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * ( 39/16 + 21/16*xi + 
         3/64*xi^2 )

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

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

       + ep^(-2)*eM1^2*cf*ca*as0pi^2 * ( 11/32 )

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

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

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

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

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

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

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

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

       + ep^(-2)*eM1^3*nh^2*cf*tf^2*as0pi^3 * ( 5/27 )

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

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

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

       + ep^(-2)*eM1^3*nl^2*cf*tf^2*as0pi^3 * ( 5/27 )

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

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

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

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

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

       + ep^(-1)*eM1*cf*as0pi * ( 3/4 )

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

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

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

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

       + ep^(-1)*eM1^3*cf*ca^2*as0pi^3 * ( 211999/20736 + 1/1536*xi - 
         5779/3456*Zeta[2] - 1/120*Zeta[2]^2 - 1/240*Zeta[2]^2*xi + 33/
         8*Zeta[2]*Log[2] - 1865/2304*Zeta[3] + 3/512*Zeta[3]*xi )

       + ep^(-1)*eM1^3*cf^2*ca*as0pi^3 * ( 222203/13824 + 3117/1024*xi
          - 281/1024*xi^2 + 81005/27648*Zeta[2] + 47/1024*Zeta[2]*xi - 
         99/1024*Zeta[2]*xi^2 - 1/60*Zeta[2]^2 - 57/32*Zeta[2]*Log[2]
          + 9/16*Zeta[2]*Log[2]*xi - 203/1152*Zeta[3] - 15/64*Zeta[3]*
         xi )

       + ep^(-1)*eM1^3*cf^3*as0pi^3 * ( 12499/1536 + 5263/1024*xi + 31/
         64*xi^2 + 5677/384*Zeta[2] + 1161/512*Zeta[2]*xi + 99/512*
         Zeta[2]*xi^2 - 1/12*Zeta[2]^2 - 207/16*Zeta[2]*Log[2] - 9/8*
         Zeta[2]*Log[2]*xi + 225/64*Zeta[3] + 9/32*Zeta[3]*xi )

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

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

       + ep^(-1)*eM1^3*nh^2*cf*tf^2*as0pi^3 * ( 193/324 - 7/24*Zeta[2]
          )

       + ep^(-1)*eM1^3*nl*cf*ca*tf*as0pi^3 * (  - 16451/2592 - 185/432*
         Zeta[2] - 3/2*Zeta[2]*Log[2] + 7/144*Zeta[3] )

       + ep^(-1)*eM1^3*nl*cf^2*tf*as0pi^3 * (  - 19105/3456 - 155/256*
         xi - 13565/3456*Zeta[2] - 33/128*Zeta[2]*xi + 3*Zeta[2]*Log[2]
          - 7/36*Zeta[3] )

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

       + ep^(-1)*eM1^3*nl^2*cf*tf^2*as0pi^3 * ( 301/324 + 3/8*Zeta[2] )

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

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

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

       + M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * ( 81/16 + 57/8*xi + 9/16*xi^2
          + 221/16*Zeta[2] + 33/16*Zeta[2]*xi + 15/64*Zeta[2]*xi^2 - 39/
         4*Zeta[2]*Log[2] + 55/16*Zeta[3] )

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

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

       + eM1*cf*as0pi * ( 1 )

       + eM1^2*cf*ca*as0pi^2 * ( 93/16 - 25/32*Zeta[2] + 9/4*Zeta[2]*
         Log[2] - 7/16*Zeta[3] )

       + eM1^2*cf^2*as0pi^2 * ( 4 + 13/8*xi + 135/32*Zeta[2] + 9/32*
         Zeta[2]*xi - 9/2*Zeta[2]*Log[2] + 5/8*Zeta[3] )

       + eM1^2*nh*cf*tf*as0pi^2 * (  - 5/4 + 3/8*Zeta[2] )

       + eM1^2*nl*cf*tf*as0pi^2 * (  - 7/4 - 5/8*Zeta[2] )

       + eM1^3*cf*ca^2*as0pi^3 * ( 774857/15552 + 13/1536*xi - 106751/
         10368*Zeta[2] + 3/256*Zeta[2]*xi + 142579/11520*Zeta[2]^2 - 17/
         1536*Zeta[2]^2*xi - 415/96*Zeta[2]*Zeta[3] - 1/48*Zeta[2]*
         Zeta[3]*xi + 693/16*Zeta[2]*Log[2] - 177/16*Zeta[2]*Log[2]^2
          - 56843/3456*Zeta[3] + 13/512*Zeta[3]*xi + 593/384*Zeta[5] - 
         7/768*Zeta[5]*xi - 57/64*Log[2]^4 - 171/8*PolyLog[4,1/2] )

       + eM1^3*cf^2*ca*as0pi^3 * ( 4893479/82944 + 18029/2048*xi - 179/
         1024*xi^2 + 700225/20736*Zeta[2] - 25/256*Zeta[2]*xi - 231/
         2048*Zeta[2]*xi^2 - 191869/11520*Zeta[2]^2 + 519/640*Zeta[2]^2
         *xi + 287/48*Zeta[2]*Zeta[3] - 519/32*Zeta[2]*Log[2] + 3*Zeta[
         2]*Log[2]*xi + 187/16*Zeta[2]*Log[2]^2 - 9/8*Zeta[2]*Log[2]^2*
         xi + 185237/13824*Zeta[3] - 709/512*Zeta[3]*xi - 63/512*Zeta[3
         ]*xi^2 - 73/48*Zeta[5] + 55/192*Log[2]^4 - 3/32*Log[2]^4*xi + 
         55/8*PolyLog[4,1/2] - 9/4*PolyLog[4,1/2]*xi )

       + eM1^3*cf^3*as0pi^3 * ( 30571/1024 + 10045/2048*xi - 1/64*xi^2
          + 127063/1536*Zeta[2] + 5151/1024*Zeta[2]*xi + 33/256*Zeta[2]
         *xi^2 - 11897/480*Zeta[2]^2 - 189/80*Zeta[2]^2*xi + 10/3*Zeta[
         2]*Zeta[3] - 2253/16*Zeta[2]*Log[2] - 6*Zeta[2]*Log[2]*xi + 
         167/8*Zeta[2]*Log[2]^2 + 9/4*Zeta[2]*Log[2]^2*xi + 5273/96*
         Zeta[3] + 861/256*Zeta[3]*xi + 63/256*Zeta[3]*xi^2 - 85/24*
         Zeta[5] + 287/96*Log[2]^4 + 3/16*Log[2]^4*xi + 287/4*PolyLog[4
         ,1/2] + 9/2*PolyLog[4,1/2]*xi )

       + eM1^3*nh*cf*ca*tf*as0pi^3 * (  - 462365/15552 - 407/1728*xi + 
         149173/8640*Zeta[2] - 3/128*Zeta[2]*xi - 213/80*Zeta[2]^2 + 7/
         32*Zeta[2]*Zeta[3] - 43/2*Zeta[2]*Log[2] + 3*Zeta[2]*Log[2]^2
          + 4055/576*Zeta[3] + 7/192*Zeta[3]*xi - 15/64*Zeta[5] + 1/4*
         Log[2]^4 + 6*PolyLog[4,1/2] )

       + eM1^3*nh*cf^2*tf*as0pi^3 * (  - 167911/6912 - 1385/512*xi + 
         32039/3456*Zeta[2] + 365/256*Zeta[2]*xi + 33/5*Zeta[2]^2 - 41/
         8*Zeta[2]*Log[2] - 9/4*Zeta[2]*Log[2]*xi - 6*Zeta[2]*Log[2]^2
          - 27/8*Zeta[3] + 87/64*Zeta[3]*xi - 1/2*Log[2]^4 - 12*
         PolyLog[4,1/2] )

       + eM1^3*nh^2*cf*tf^2*as0pi^3 * ( 4193/972 - 209/90*Zeta[2] + 2*
         Zeta[2]*Log[2] - 41/36*Zeta[3] )

       + eM1^3*nl*cf*ca*tf*as0pi^3 * (  - 463001/15552 - 18913/5184*
         Zeta[2] - 2957/720*Zeta[2]^2 - 15*Zeta[2]*Log[2] + 4*Zeta[2]*
         Log[2]^2 + 845/432*Zeta[3] + 1/3*Log[2]^4 + 8*PolyLog[4,1/2] )

       + eM1^3*nl*cf^2*tf*as0pi^3 * (  - 451333/20736 - 929/512*xi - 
         289021/10368*Zeta[2] - 187/256*Zeta[2]*xi + 2039/180*Zeta[2]^2
          + 30*Zeta[2]*Log[2] - 8*Zeta[2]*Log[2]^2 - 26249/1728*Zeta[3]
          - 21/64*Zeta[3]*xi - 2/3*Log[2]^4 - 16*PolyLog[4,1/2] )

       + eM1^3*nl*nh*cf*tf^2*as0pi^3 * ( 2218/243 - 13/9*Zeta[2] + 2*
         Zeta[2]*Log[2] - 7/9*Zeta[3] )

       + eM1^3*nl^2*cf*tf^2*as0pi^3 * ( 2029/486 + 5/2*Zeta[2] + 37/36*
         Zeta[3] )

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

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

       + ep*M1^(-1)*eM1^2*cf*ca*as0pi^2*dmf * ( 89/8 + 165/128*xi + 75/
         256*xi^2 - 329/64*Zeta[2] - 11/128*Zeta[2]*xi + 5/256*Zeta[2]*
         xi^2 + 335/32*Zeta[2]^2 - 37/80*Zeta[2]^2*xi + 159/8*Zeta[2]*
         Log[2] - 39/4*Zeta[2]*Log[2]^2 - 871/64*Zeta[3] + 1/64*Zeta[3]
         *xi - 11/64*Zeta[3]*xi^2 - 13/16*Log[2]^4 - 39/2*PolyLog[4,1/2
         ] )

       + ep*M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * ( 197/16 - 41/8*xi - 15/16
         *xi^2 + 401/16*Zeta[2] - 5/8*Zeta[2]*xi - 5/32*Zeta[2]*xi^2 - 
         671/40*Zeta[2]^2 - 159/4*Zeta[2]*Log[2] + 39/2*Zeta[2]*Log[2]^
         2 + 497/16*Zeta[3] + 5/8*Zeta[3]*xi + 11/32*Zeta[3]*xi^2 + 13/
         8*Log[2]^4 + 39*PolyLog[4,1/2] )

       + ep*M1^(-1)*eM1^2*nh*cf*tf*as0pi^2*dmf * (  - 199/36 + 9/2*
         Zeta[2] - 9/2*Zeta[2]*Log[2] + 23/8*Zeta[3] )

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

       + ep*eM1*cf*as0pi * ( 2 + 3/8*Zeta[2] )

       + ep*eM1^2*cf*ca*as0pi^2 * ( 77/4 - 29/8*Zeta[2] + 83/16*Zeta[2]
         ^2 + 27/2*Zeta[2]*Log[2] - 9/2*Zeta[2]*Log[2]^2 - 281/48*Zeta[
         3] - 3/8*Log[2]^4 - 9*PolyLog[4,1/2] )

       + ep*eM1^2*cf^2*as0pi^2 * ( 13 + 7/4*xi + 77/4*Zeta[2] + 3/16*
         Zeta[2]*xi - 113/10*Zeta[2]^2 - 27*Zeta[2]*Log[2] + 9*Zeta[2]*
         Log[2]^2 + 207/16*Zeta[3] - 3/16*Zeta[3]*xi + 3/4*Log[2]^4 + 
         18*PolyLog[4,1/2] )

       + ep*eM1^2*nh*cf*tf*as0pi^2 * (  - 25/4 + 3*Zeta[2] - 3*Zeta[2]*
         Log[2] + 11/6*Zeta[3] )

       + ep*eM1^2*nl*cf*tf*as0pi^2 * (  - 23/4 - 5/2*Zeta[2] - 11/12*
         Zeta[3] )

       + ep^2*M1^(-2)*eM1*cf*as0pi*dmf^2 * (  - 4 - 4*xi - 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 + 1/8*Zeta[2] - 1/8*
         Zeta[2]*xi - 1/4*Zeta[3] - 1/8*Zeta[3]*xi )

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

       + 1;
 );
