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

       + ep^(-3)*eM1^3*cf^3*as0pi^3 * (  - 75/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 * ( 5/64 - 3/64*xi )

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

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

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

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

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

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

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

       + ep^(-2)*eM1^3*cf^3*as0pi^3 * (  - 91/64 - 51/512*xi + 3/128*
         xi^2 - 7/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 * ( 7/12 - 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 * ( 5/36 - 1/36*Zeta[2] )

       + ep^(-2)*eM1^3*nl*cf^2*tf*as0pi^3 * ( 151/288 - 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 * (  - 1/2 + 3/8*xi )

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

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

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

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

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

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

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

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

       + ep^(-1)*eM1^3*cf*ca^2*as0pi^3 * (  - 14041/6912 + 1/1536*xi + 
         659/864*Zeta[2] - 1/120*Zeta[2]^2 - 1/240*Zeta[2]^2*xi - 11/24
         *Zeta[2]*Log[2] + 775/2304*Zeta[3] + 3/512*Zeta[3]*xi )

       + ep^(-1)*eM1^3*cf^2*ca*as0pi^3 * (  - 3575/432 + 3117/1024*xi
          - 281/1024*xi^2 - 55135/27648*Zeta[2] + 47/1024*Zeta[2]*xi - 
         99/1024*Zeta[2]*xi^2 - 1/60*Zeta[2]^2 - 47/96*Zeta[2]*Log[2]
          + 9/16*Zeta[2]*Log[2]*xi - 359/1152*Zeta[3] - 15/64*Zeta[3]*
         xi )

       + ep^(-1)*eM1^3*cf^3*as0pi^3 * (  - 1499/384 + 2863/1024*xi + 31/
         64*xi^2 - 3427/768*Zeta[2] + 999/512*Zeta[2]*xi + 99/512*Zeta[
         2]*xi^2 - 1/12*Zeta[2]^2 + 45/16*Zeta[2]*Log[2] - 9/8*Zeta[2]*
         Log[2]*xi - 107/64*Zeta[3] + 9/32*Zeta[3]*xi )

       + ep^(-1)*eM1^3*nh*cf*ca*tf*as0pi^3 * ( 5/9 + 35/576*xi + 55/
         1152*Zeta[2] + 1/128*Zeta[2]*xi + 1/6*Zeta[2]*Log[2] - 1/8*
         Zeta[3] )

       + ep^(-1)*eM1^3*nh*cf^2*tf*as0pi^3 * ( 7/6 - 227/256*xi + 493/
         576*Zeta[2] + 39/128*Zeta[2]*xi - 1/3*Zeta[2]*Log[2] + 5/12*
         Zeta[3] )

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

       + ep^(-1)*eM1^3*nl*cf*ca*tf*as0pi^3 * ( 1189/864 - 61/216*Zeta[2
         ] + 1/6*Zeta[2]*Log[2] - 17/144*Zeta[3] )

       + ep^(-1)*eM1^3*nl*cf^2*tf*as0pi^3 * ( 1613/864 - 155/256*xi + 
         4621/3456*Zeta[2] - 33/128*Zeta[2]*xi - 1/3*Zeta[2]*Log[2] + 7/
         18*Zeta[3] )

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

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

       + 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 * (  - 4525/864 + 143/128*xi
          - 79/256*xi^2 + 197/384*Zeta[2] + 13/128*Zeta[2]*xi - 15/128*
         Zeta[2]*xi^2 - 7/8*Zeta[2]*Log[2] + 21/32*Zeta[3] - 1/8*Zeta[3
         ]*xi )

       + M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * (  - 835/144 + 21/4*xi + 9/16
         *xi^2 - 1345/288*Zeta[2] + 57/32*Zeta[2]*xi + 15/64*Zeta[2]*
         xi^2 + 7/4*Zeta[2]*Log[2] - 77/48*Zeta[3] )

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

       + M1^(-1)*eM1^2*nl*cf*tf*as0pi^2*dmf * ( 191/216 + 5/48*Zeta[2]
          )

       + eM1*cf*as0pi * (  - 1/2 )

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

       + eM1^2*cf^2*as0pi^2 * (  - 1075/384 + 13/8*xi - 37/24*Zeta[2]
          + 9/32*Zeta[2]*xi + 1/2*Zeta[2]*Log[2] - 5/8*Zeta[3] )

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

       + eM1^2*nl*cf*tf*as0pi^2 * ( 149/288 )

       + eM1^3*cf*ca^2*as0pi^3 * (  - 286945/20736 + 13/1536*xi + 88873/
         10368*Zeta[2] + 3/256*Zeta[2]*xi + 2681/3840*Zeta[2]^2 - 17/
         1536*Zeta[2]^2*xi - 1/48*Zeta[2]*Zeta[3]*xi - 1447/144*Zeta[2]
         *Log[2] + 119/144*Zeta[2]*Log[2]^2 + 13639/3456*Zeta[3] + 13/
         512*Zeta[3]*xi + 133/384*Zeta[5] - 7/768*Zeta[5]*xi + 173/1728
         *Log[2]^4 + 173/72*PolyLog[4,1/2] )

       + eM1^3*cf^2*ca*as0pi^3 * (  - 382945/20736 + 18029/2048*xi - 
         179/1024*xi^2 - 588139/41472*Zeta[2] - 25/256*Zeta[2]*xi - 231/
         2048*Zeta[2]*xi^2 - 12709/2304*Zeta[2]^2 + 519/640*Zeta[2]^2*
         xi - 31/48*Zeta[2]*Zeta[3] + 1567/288*Zeta[2]*Log[2] + 3*Zeta[
         2]*Log[2]*xi + 407/144*Zeta[2]*Log[2]^2 - 9/8*Zeta[2]*Log[2]^2
         *xi - 106783/13824*Zeta[3] - 709/512*Zeta[3]*xi - 63/512*Zeta[
         3]*xi^2 - 53/48*Zeta[5] - 133/1728*Log[2]^4 - 3/32*Log[2]^4*xi
          - 133/72*PolyLog[4,1/2] - 9/4*PolyLog[4,1/2]*xi )

       + eM1^3*cf^3*as0pi^3 * (  - 3085/256 - 2243/2048*xi - 1/64*xi^2
          - 13967/384*Zeta[2] + 4287/1024*Zeta[2]*xi + 33/256*Zeta[2]*
         xi^2 + 187/160*Zeta[2]^2 - 189/80*Zeta[2]^2*xi + 55/24*Zeta[2]
         *Zeta[3] + 469/16*Zeta[2]*Log[2] - 6*Zeta[2]*Log[2]*xi - 215/
         24*Zeta[2]*Log[2]^2 + 9/4*Zeta[2]*Log[2]^2*xi - 2593/384*Zeta[
         3] + 915/256*Zeta[3]*xi + 63/256*Zeta[3]*xi^2 - 5/6*Zeta[5] - 
         71/288*Log[2]^4 + 3/16*Log[2]^4*xi - 71/12*PolyLog[4,1/2] + 9/
         2*PolyLog[4,1/2]*xi )

       + eM1^3*nh*cf*ca*tf*as0pi^3 * ( 1481/324 - 407/1728*xi + 673/360
         *Zeta[2] - 3/128*Zeta[2]*xi + 83/720*Zeta[2]^2 - 31/96*Zeta[2]
         *Zeta[3] - 9/2*Zeta[2]*Log[2] - 4/9*Zeta[2]*Log[2]^2 + 235/576
         *Zeta[3] + 7/192*Zeta[3]*xi + 55/192*Zeta[5] - 1/108*Log[2]^4
          - 2/9*PolyLog[4,1/2] )

       + eM1^3*nh*cf^2*tf*as0pi^3 * ( 53/32 - 1385/512*xi + 319/120*
         Zeta[2] + 365/256*Zeta[2]*xi + 161/180*Zeta[2]^2 - 31/24*Zeta[
         2]*Log[2] - 9/4*Zeta[2]*Log[2]*xi + 8/9*Zeta[2]*Log[2]^2 + 625/
         144*Zeta[3] + 87/64*Zeta[3]*xi + 1/54*Log[2]^4 + 4/9*PolyLog[4
         ,1/2] )

       + eM1^3*nh^2*cf*tf^2*as0pi^3 * ( 11/162 - 71/360*Zeta[2] + 2/3*
         Zeta[2]*Log[2] - 11/18*Zeta[3] )

       + eM1^3*nl*cf*ca*tf*as0pi^3 * ( 1219/144 - 2341/1296*Zeta[2] - 7/
         48*Zeta[2]^2 + 26/9*Zeta[2]*Log[2] - 4/9*Zeta[2]*Log[2]^2 - 
         269/216*Zeta[3] - 1/27*Log[2]^4 - 8/9*PolyLog[4,1/2] )

       + eM1^3*nl*cf^2*tf*as0pi^3 * ( 6011/1296 - 929/512*xi + 58501/
         5184*Zeta[2] - 187/256*Zeta[2]*xi + 101/90*Zeta[2]^2 - 52/9*
         Zeta[2]*Log[2] + 8/9*Zeta[2]*Log[2]^2 + 7489/1728*Zeta[3] - 21/
         64*Zeta[3]*xi + 2/27*Log[2]^4 + 16/9*PolyLog[4,1/2] )

       + eM1^3*nl*nh*cf*tf^2*as0pi^3 * (  - 112/81 - 49/108*Zeta[2] + 2/
         3*Zeta[2]*Log[2] - 7/18*Zeta[3] )

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

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

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

       + ep*M1^(-1)*eM1^2*cf*ca*as0pi^2*dmf * (  - 3379/324 + 165/128*
         xi + 75/256*xi^2 + 3443/576*Zeta[2] - 11/128*Zeta[2]*xi + 5/
         256*Zeta[2]*xi^2 - 7/32*Zeta[2]^2 - 37/80*Zeta[2]^2*xi - 197/
         24*Zeta[2]*Log[2] + 7/4*Zeta[2]*Log[2]^2 + 2089/576*Zeta[3] + 
         1/64*Zeta[3]*xi - 11/64*Zeta[3]*xi^2 + 7/48*Log[2]^4 + 7/2*
         PolyLog[4,1/2] )

       + ep*M1^(-1)*eM1^2*cf^2*as0pi^2*dmf * (  - 3533/432 - 65/8*xi - 
         15/16*xi^2 - 2377/108*Zeta[2] - Zeta[2]*xi - 5/32*Zeta[2]*xi^2
          - 77/120*Zeta[2]^2 + 197/12*Zeta[2]*Log[2] - 7/2*Zeta[2]*Log[
         2]^2 - 337/36*Zeta[3] + 13/16*Zeta[3]*xi + 11/32*Zeta[3]*xi^2
          - 7/24*Log[2]^4 - 7*PolyLog[4,1/2] )

       + ep*M1^(-1)*eM1^2*nh*cf*tf*as0pi^2*dmf * ( 25/108 + 19/18*Zeta[
         2] - 7/2*Zeta[2]*Log[2] + 145/72*Zeta[3] )

       + ep*M1^(-1)*eM1^2*nl*cf*tf*as0pi^2*dmf * ( 1867/648 + 95/72*
         Zeta[2] + 11/72*Zeta[3] )

       + ep*eM1*cf*as0pi * (  - 1 )

       + ep*eM1^2*cf*ca*as0pi^2 * (  - 46885/6912 + 551/192*Zeta[2] - 1/
         16*Zeta[2]^2 - 19/6*Zeta[2]*Log[2] + 1/2*Zeta[2]*Log[2]^2 + 37/
         24*Zeta[3] + 1/24*Log[2]^4 + PolyLog[4,1/2] )

       + ep*eM1^2*cf^2*as0pi^2 * (  - 6331/2304 + 7/4*xi - 5551/576*
         Zeta[2] + 3/16*Zeta[2]*xi - 4/5*Zeta[2]^2 + 19/3*Zeta[2]*Log[2
         ] - Zeta[2]*Log[2]^2 - 11/3*Zeta[3] - 3/16*Zeta[3]*xi - 1/12*
         Log[2]^4 - 2*PolyLog[4,1/2] )

       + ep*eM1^2*nh*cf*tf*as0pi^2 * ( 989/1728 + 55/144*Zeta[2] - 
         Zeta[2]*Log[2] + 7/12*Zeta[3] )

       + ep*eM1^2*nl*cf*tf*as0pi^2 * ( 3413/1728 + 25/48*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 * (  - 1 - xi - 3/8*Zeta[2] - 1/
         8*Zeta[2]*xi + 1/6*Zeta[3] - 1/8*Zeta[3]*xi )

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

       + 1;
 );
