matadSymbols = {D3 -> -Pi^4/24 - 6*Im[PolyLog[2, E^((I/3)*Pi)]]^2 + 
       6*Zeta[3], B4 -> (-13*Pi^4)/180 - (2*Pi^2*Log[2]^2)/3 + 
       (2*Log[2]^4)/3 + 16*PolyLog[4, 1/2], 
     S2 -> (4*Im[PolyLog[2, E^((I/3)*Pi)]])/(9*Sqrt[3]), 
     OepS2 -> -763/32 + (65*Pi^2)/32 - (35*Pi^3)/(16*Sqrt[3]) + 
       (19*Pi^4)/480 + (45*Sqrt[3]*Im[PolyLog[2, E^((I/3)*Pi)]])/2 - 
       27*Sqrt[3]*Im[PolyLog[3, 1/(Sqrt[3]*E^((I/6)*Pi))]] - 
       (9*Sqrt[3]*Pi*Log[3]^2)/16 - (15*Zeta[3])/4, 
     T1ep -> -45/2 - (3*Pi^2)/4 - (35*Pi^3)/(72*Sqrt[3]) + 
       6*Sqrt[3]*Im[PolyLog[2, E^((I/3)*Pi)]] - 
       6*Sqrt[3]*Im[PolyLog[3, 1/(Sqrt[3]*E^((I/6)*Pi))]] - 
       (Sqrt[3]*Pi*Log[3]^2)/8 + Zeta[3]}

hplsToKnownFunctions = {HPL[{0}, x] -> Log[x], HPL[{-1, 0}, x] -> 
      Log[x]*Log[1 + x] + PolyLog[2, -x], HPL[{0, 0}, x] -> Log[x]^2/2, 
     HPL[{1, 0}, x] -> -(Log[1 - x]*Log[x]) - PolyLog[2, x], 
     HPL[{-1, 0, 0}, x] -> (Log[x]^2*Log[1 + x])/2 + Log[x]*PolyLog[2, -x] - 
       PolyLog[3, -x], HPL[{0, 0, 0}, x] -> Log[x]^3/6, 
     HPL[{1, 0, 0}, x] -> -(Log[1 - x]*Log[x]^2)/2 - Log[x]*PolyLog[2, x] + 
       PolyLog[3, x],HPL[{0}, X] -> Log[X], HPL[{-1, 0}, X] -> 
      Log[X]*Log[1 + X] + PolyLog[2, -X], HPL[{0, 0}, X] -> Log[X]^2/2, 
     HPL[{1, 0}, X] -> -(Log[1 - X]*Log[X]) - PolyLog[2, X], 
     HPL[{-1, 0, 0}, X] -> (Log[X]^2*Log[1 + X])/2 + Log[X]*PolyLog[2, -X] - 
       PolyLog[3, -X], HPL[{0, 0, 0}, X] -> Log[X]^3/6, 
     HPL[{1, 0, 0}, X] -> -(Log[1 - X]*Log[X]^2)/2 - Log[X]*PolyLog[2, X] + 
       PolyLog[3, X]}


rhoMsBar0 = 2/3*(nc*(1 + x^2 - (4*x^2*HPL[{0}, x])/(-1 + x^2)))
 
rhoMsBar1 = 2/3*(cf*nc*((3*lm*(1 + x^2))/2 + 
      (Pi^2*(-1 + x^2) + 9*(1 + x^2))/6 + ((-6*lm*x^2)/(-1 + x^2) + 
        (x^2 + 3*x^4)/(1 - x^2))*HPL[{0}, x] + (2 - 2*x^2)*HPL[{-1, 0}, x] + 
      (4*x^2*(-3 + x^4)*HPL[{0, 0}, x])/(-1 + x^2)^2 + 
      2*(-1 + x^2)*HPL[{1, 0}, x]))
 
rhoMsBar2Exp0 = 2/3*(nc*(cf*(77/48 - lm^2/4 + (7*Pi^2)/72 + 
        (lm*(-33 + 2*Pi^2))/36 + (243*S2)/8 + 
        nl*(-25/32 - lm^2/8 + (7*Pi^2)/72 + (lm*(-33 + 2*Pi^2))/72 - 
          Zeta[3]/3) - (34*Zeta[3])/3) + ca*cf*(1013/192 + (3*B4)/4 - D3/8 + 
        (11*lm^2)/16 - (73*Pi^2)/144 + (67*Pi^4)/720 + (lm*(489 - 22*Pi^2))/
         144 + (1053*S2)/32 + (2*Pi^2*Log[2]^2)/3 - (2*Log[2]^4)/3 - 
        16*PolyLog[4, 1/2] - (113*Zeta[3])/12) + 
      cf^2*(-1591/576 - (3*B4)/2 + D3/4 + (9*lm^2)/8 - (259*Pi^2)/432 - 
        (17*Pi^4)/120 + (lm*(3 - 4*Pi^2))/16 - (1053*S2)/16 - 
        (4*Pi^2*Log[2]^2)/3 + (4*Log[2]^4)/3 + 32*PolyLog[4, 1/2] + 
        (271*Zeta[3])/12) + 
      x^30*(ca*cf*(-8899970915383290187877291494307757785379543857650573/
           3121998464581776894187075146770545920000000 + 
          (6593691705196*OepS2)/1412147682405 + (413483134619419040907632123*
            Pi^2)/506619843103325952 - (31320035599681*Pi^4)/
           169457721888600 - (3348766661833717351*S2)/18847464401165400 - 
          (3296845852598*T1ep)/470715894135 + (11*lm^2*Log[x])/4 + 
          (2294203017781615387123*Log[x]^2)/234247513500 + 
          (347694541078*Log[x]^3)/105 + lm*(-319/529200 + (3704*Log[x])/315 - 
            (176*Log[x]^2)/3) + Log[x]*
           (29698798127367637331192848128939527765983/
             19436340473748897676230052800000 + (2086167213029*Pi^2)/2520 + 
            (3296845852598*S2)/17433922005 - 6*Zeta[3]) - 
          (5455511847590575223147*Zeta[3])/3295011258945) + 
        cf*(901924402408539819525293391294800564608159570056821974861/
           1518990301798016212650897723418457600087838000000 - 
          (22181171683297*OepS2)/847288609443 - (2001691411612598071015673*
            Pi^2)/12336522153490080 + (421442261982643*Pi^4)/
           406698532532640 + (2827584827341465903*S2)/1130847864069924 + 
          (22181171683297*T1ep)/564859072962 - lm^2*Log[x] + 
          (-15027507493917431110138734984452913552292477/
             84152516824416068843715026266140000 - (3071706077*Pi^2)/18 - 
            (22181171683297*S2)/20920706406)*Log[x] - 
          (26359368501267207855601*Log[x]^2)/13537833083775 - 
          (6143411876*Log[x]^3)/9 + lm*(29/132300 - (946*Log[x])/315 + 
            (64*Log[x]^2)/3) + (385573119530223220973*Zeta[3])/
           1129718145924 + nl*(-1553833941721187149/140388324043968000 + 
            (209*Pi^2)/3780 - (lm^2*Log[x])/2 + (719538439547/55653998400 + 
              (29*Pi^2)/18)*Log[x] + (49*Log[x]^2)/3 - (64*Log[x]^3)/9 + 
            lm*(29/264600 - (473*Log[x])/315 + (32*Log[x]^2)/3) + 
            (29*Zeta[3])/3)) + 
        cf^2*(202272065472312835132853932766505398599381519158317/
           35477255279338373797580399395119840000000 - (13187383410392*OepS2)/
           1412147682405 - (8786516610057672557722283153*Pi^2)/
           5382835832972838240 + (31320035599681*Pi^4)/84728860944300 + 
          (3348766661833717351*S2)/9423732200582700 + (6593691705196*T1ep)/
           470715894135 + (9*lm^2*Log[x])/2 - (9176812182289817897647*
            Log[x]^2)/468495027000 - (139077806027*Log[x]^3)/21 + 
          lm*(-29/29400 - (313*Log[x])/140 - 96*Log[x]^2) + 
          (10911023758884701452564*Zeta[3])/3295011258945 + 
          Log[x]*(-2121342685929056528063721300232099122197/
             694155016919603488436787600000 - (2086167208969*Pi^2)/1260 - 
            (6593691705196*S2)/17433922005 + 12*Zeta[3]))) + 
      x^26*(ca*cf*((11*lm^2*Log[x])/4 + lm*(-275/292032 + (11009*Log[x])/
             936 - (154*Log[x]^2)/3) + 
          (3*(-110978356169892890645527755360666118362937/
              2888088818741710305352386852096000 + (3236352234550*OepS2)/
              1223861324751 + (4207088299364469421996517*Pi^2)/
              392027259544240320 - (6149069245645*Pi^4)/58745343588048 - 
             (1578034643561982247*S2)/16334469147676680 - 
             (1618176117275*T1ep)/407953774917 + (438441725754534083*Log[x]^
                2)/3404591190 + (10404265529*Log[x]^3)/234 + 
             Log[x]*(127637794895084006147588729717539/
                8270850863222634220934400 + (10404258235*Pi^2)/936 + 
               (1618176117275*S2)/15109399071 - 4*Zeta[3]) - 
             (18138756229428738805*Zeta[3])/815907549834))/2) + 
        cf*(31203413708499950191374997282323795141677579/
           2244362825354969586918474771948000000 - (225140300071*OepS2)/
           10460353203 - (6194807462902180040437*Pi^2)/1675330168992480 + 
          (4277665701349*Pi^4)/5020969537440 + (11335546977284381*S2)/
           5369647977540 + (225140300071*T1ep)/6973568802 - lm^2*Log[x] + 
          (-26471229115627911114258739896959/13972524495161100813945000 - 
            (72027481*Pi^2)/18 - (225140300071*S2)/258280326)*Log[x] - 
          (38736184790982661*Log[x]^2)/872972100 - (144054724*Log[x]^3)/9 + 
          lm*(25/73008 - (703*Log[x])/234 + (56*Log[x]^2)/3) + 
          (111615012295820339*Zeta[3])/13947137604 + 
          nl*(-817625476519723/86392814796288 + (155*Pi^2)/2808 - 
            (lm^2*Log[x])/2 + (64318606241/5993507520 + (25*Pi^2)/18)*
             Log[x] + (43*Log[x]^2)/3 - (56*Log[x]^3)/9 + 
            lm*(25/146016 - (703*Log[x])/468 + (28*Log[x]^2)/3) + 
            (25*Zeta[3])/3)) + 
        cf^2*(1220761716977221333767274420556694302311427/
           10589659002052937786292085124352000 - (3236352234550*OepS2)/
           407953774917 - (4207088284622530936403237*Pi^2)/
           130675753181413440 + (6149069245645*Pi^4)/19581781196016 + 
          (1578034643561982247*S2)/5444823049225560 + (1618176117275*T1ep)/
           135984591639 + (9*lm^2*Log[x])/2 - (1753767845559363617*Log[x]^2)/
           4539454920 - (2400977227*Log[x]^3)/18 + 
          lm*(-25/16224 - (29*Log[x])/13 - 84*Log[x]^2) + 
          (18138760762248460105*Zeta[3])/271969183278 + 
          Log[x]*(-1404014238663981623071664146859849/
             30326453165149658810092800 - (31212772105*Pi^2)/936 - 
            (1618176117275*S2)/5036466357 + 12*Zeta[3]))) + 
      x^24*(ca*cf*(-1687072701766452148409556505474395410333/
           199121793200995411366175959603200 + (201786670873*OepS2)/
           55788550416 + (32743248382726845607073*Pi^2)/14040862368698880 - 
          (3833946746587*Pi^4)/26778504199680 - (191333864404573*S2)/
           1487694677760 - (201786670873*T1ep)/37192366944 + 
          (11*lm^2*Log[x])/4 + (61422302388488201*Log[x]^2)/2194901280 + 
          (969989888*Log[x]^3)/99 + lm*(-23/19008 + (847*Log[x])/72 - 
            (143*Log[x]^2)/3) + Log[x]*(478434629021524093398149297101/
             178200312402054619076400 + (7759884575*Pi^2)/3168 + 
            (201786670873*S2)/1377495072 - 6*Zeta[3]) - 
          (4008525293737465447*Zeta[3])/818232072768) + 
        cf*(38171294945949999940668098554044217301195701/
           17206781661054766833041639918268000000 - (67181973811*OepS2)/
           3486784401 - (24960809970064015741*Pi^2)/42957183820320 + 
          (1276457502409*Pi^4)/1673656512480 + (312389104158763*S2)/
           162716605380 + (67181973811*T1ep)/2324522934 - lm^2*Log[x] + 
          (-3678761218474568157770096289173/53561343898117553120122500 - 
            (11541479*Pi^2)/18 - (67181973811*S2)/86093442)*Log[x] - 
          (140015739251848667*Log[x]^2)/20078358300 - (23082740*Log[x]^3)/9 + 
          lm*(23/52272 - (595*Log[x])/198 + (52*Log[x]^2)/3) + 
          (5960453171867951*Zeta[3])/4649045868 + 
          nl*(-50316679953851/5809091904000 + (131*Pi^2)/2376 - 
            (lm^2*Log[x])/2 + (403698833/41912640 + (23*Pi^2)/18)*Log[x] + 
            (40*Log[x]^2)/3 - (52*Log[x]^3)/9 + lm*(23/104544 - 
              (595*Log[x])/396 + (26*Log[x]^2)/3) + (23*Zeta[3])/3)) + 
        cf^2*(8435354859457052965431815787530262899497/
           497804483002488528415439899008000 - (201786670873*OepS2)/
           27894275208 - (229202733120776839153319*Pi^2)/49143018290446080 + 
          (3833946746587*Pi^4)/13389252099840 + (191333864404573*S2)/
           743847338880 + (201786670873*T1ep)/18596183472 + 
          (9*lm^2*Log[x])/2 - (20474171336258273*Log[x]^2)/365816880 - 
          (646648372*Log[x]^3)/33 + lm*(-23/11616 - (49*Log[x])/22 - 
            78*Log[x]^2) + (4008531566850023335*Zeta[3])/409116036384 + 
          Log[x]*(-1913721593089873659204941652991/356400624804109238152800 - 
            (7759880527*Pi^2)/1584 - (201786670873*S2)/688747536 + 
            12*Zeta[3]))) + 
      x^22*(ca*cf*(-103252976687230269785307215319329773/
           80640923199444125061223680000 + (13930593748*OepS2)/4261625379 + 
          (1820728422400692493*Pi^2)/5250322466928 - (66170320303*Pi^4)/
           511395045480 - (491679138747457*S2)/4375268722440 - 
          (6965296874*T1ep)/1420541793 + (11*lm^2*Log[x])/4 + 
          (19474547100968*Log[x]^2)/4679675 + (97873681*Log[x]^3)/66 + 
          lm*(-7/4400 + (353*Log[x])/30 - 44*Log[x]^2) + 
          Log[x]*(604303801072088090300814941/2193909663306304944000 + 
            (24467753*Pi^2)/66 + (6965296874*S2)/52612659 - 6*Zeta[3]) - 
          (4213497840187199*Zeta[3])/5682167172) + 
        cf*(84895205792624436237300097986043541327/
           231741167152252751960156766576000 - (2203827571*OepS2)/129140163 - 
          (150222751927236661*Pi^2)/1591006808160 + (41872723849*Pi^4)/
           61987278240 + (372614130181*S2)/215233605 + (2203827571*T1ep)/
           86093442 - lm^2*Log[x] + (22924832893190425913076256757/
             634801112866578407349600 - (641429*Pi^2)/6 - (2203827571*S2)/
             3188646)*Log[x] - (42152335920852*Log[x]^2)/37182145 - 
          (1282792*Log[x]^3)/3 + lm*(7/12100 - (496*Log[x])/165 + 
            16*Log[x]^2) + (36767651927735*Zeta[3])/172186884 + 
          nl*(-2174001437191/276623424000 + (109*Pi^2)/1980 - 
            (lm^2*Log[x])/2 + (3405575/399168 + (7*Pi^2)/6)*Log[x] + 
            (37*Log[x]^2)/3 - (16*Log[x]^3)/3 + 
            lm*(7/24200 - (248*Log[x])/165 + 8*Log[x]^2) + 7*Zeta[3])) + 
        cf^2*(18773152790283817003023353588235463/
           7330993018131284096474880000 - (27861187496*OepS2)/4261625379 - 
          (946778624745064602721*Pi^2)/1365083841401280 + 
          (66170320303*Pi^4)/255697522740 + (491679138747457*S2)/
           2187634361220 + (13930593748*T1ep)/1420541793 + 
          (9*lm^2*Log[x])/2 - (5608789542587981*Log[x]^2)/673873200 - 
          (44483432*Log[x]^3)/15 + lm*(-63/24200 - (489*Log[x])/220 - 
            72*Log[x]^2) + (4213537615357403*Zeta[3])/2841083586 + 
          Log[x]*(-604254058361228332395829763/1096954831653152472000 - 
            (8155892*Pi^2)/11 - (13930593748*S2)/52612659 + 12*Zeta[3]))) + 
      x^18*(ca*cf*(-239115413391483975052159567/7368702352419792384000 + 
          (36990160*OepS2)/14348907 + (437011363411423*Pi^2)/51426482688 - 
          (8785163*Pi^4)/86093442 - (567978619*S2)/7085880 - 
          (18495080*T1ep)/4782969 + (11*lm^2*Log[x])/4 + 
          (1765274666509*Log[x]^2)/17297280 + (4155323*Log[x]^3)/108 + 
          lm*(-187/62208 + (5087*Log[x])/432 - (110*Log[x]^2)/3) + 
          Log[x]*(-102492934161788597929/40896338952268800 + 
            (461293*Pi^2)/48 + (18495080*S2)/177147 - 6*Zeta[3]) - 
          (1476065031239*Zeta[3])/76527504) + 
        cf*(695376268664998050425665068152107/62078124828842564960149632000 - 
          (184255631*OepS2)/14348907 - (45048776518267*Pi^2)/16070775840 + 
          (3500856989*Pi^4)/6887475360 + (202772737451*S2)/148803480 + 
          (184255631*T1ep)/9565938 - lm^2*Log[x] + 
          (218416509411378713237927/38095304879921889600 - (63575*Pi^2)/18 - 
            (184255631*S2)/354294)*Log[x] - (2382247160213*Log[x]^2)/
           69837768 - (126992*Log[x]^3)/9 + lm*(17/15552 - (325*Log[x])/108 + 
            (40*Log[x]^2)/3) + (130964617843*Zeta[3])/19131876 + 
          nl*(-8337740303/1333584000 + (71*Pi^2)/1296 - (lm^2*Log[x])/2 + 
            (3348547/529200 + (17*Pi^2)/18)*Log[x] + (31*Log[x]^2)/3 - 
            (40*Log[x]^3)/9 + lm*(17/31104 - (325*Log[x])/216 + 
              (20*Log[x]^2)/3) + (17*Zeta[3])/3)) + 
        cf^2*(119534477349041168133479531/1842175588104948096000 - 
          (73980320*OepS2)/14348907 - (24035462944890073*Pi^2)/
           1414228273920 + (8785163*Pi^4)/43046721 + (567978619*S2)/3542940 + 
          (36990160*T1ep)/4782969 + (9*lm^2*Log[x])/2 - 
          (5299678078327*Log[x]^2)/25945920 - (230165*Log[x]^3)/3 + 
          lm*(-17/3456 - (53*Log[x])/24 - 60*Log[x]^2) + 
          (1476498687095*Zeta[3])/38263752 + 
          Log[x]*(103334087302961851123/20448169476134400 - 
            (1383743*Pi^2)/72 - (36990160*S2)/177147 + 12*Zeta[3]))) + 
      x^16*(ca*cf*(-1053266389635314097149/193233633094656000 + 
          (38023127*OepS2)/17006112 + (6758790850741*Pi^2)/4761711360 - 
          (722439413*Pi^4)/8162933760 - (19410594103*S2)/302330880 - 
          (38023127*T1ep)/11337408 + (11*lm^2*Log[x])/4 + 
          (26546935249*Log[x]^2)/1552320 + (572213*Log[x]^3)/84 + 
          lm*(-55/12544 + (3959*Log[x])/336 - 33*Log[x]^2) + 
          Log[x]*(-173017545026353/98028425880 + (759481*Pi^2)/448 + 
            (38023127*S2)/419904 - 6*Zeta[3]) - (547635759565*Zeta[3])/
           158723712) + cf*(74433472231187106510494431/
           37245303145512798336000 - (639595*OepS2)/59049 - 
          (32860090919*Pi^2)/66134880 + (2430461*Pi^4)/5668704 + 
          (310580299*S2)/262440 + (639595*T1ep)/39366 - lm^2*Log[x] + 
          (232027348870899763/121595081864256 - (1473*Pi^2)/2 - 
            (639595*S2)/1458)*Log[x] - (38963249749*Log[x]^2)/6126120 - 
          (8792*Log[x]^3)/3 + lm*(5/3136 - (253*Log[x])/84 + 12*Log[x]^2) + 
          (100997471*Zeta[3])/78732 + nl*(-7687021/1411200 + (55*Pi^2)/1008 - 
            (lm^2*Log[x])/2 + (26321/5040 + (5*Pi^2)/6)*Log[x] + 
            (28*Log[x]^2)/3 - 4*Log[x]^3 + lm*(5/6272 - (253*Log[x])/168 + 
              6*Log[x]^2) + 5*Zeta[3])) + 
        cf^2*(1052200087927301109659/96616816547328000 - 
          (38023127*OepS2)/8503056 - (6758516076061*Pi^2)/2380855680 + 
          (722439413*Pi^4)/4081466880 + (19410594103*S2)/151165440 + 
          (38023127*T1ep)/5668704 + (9*lm^2*Log[x])/2 - 
          (26650826509*Log[x]^2)/776160 - (282040*Log[x]^3)/21 + 
          lm*(-45/6272 - (123*Log[x])/56 - 54*Log[x]^2) + 
          (548429378125*Zeta[3])/79361856 + Log[x]*
           (11195734117734871/3136909628160 - (2277323*Pi^2)/672 - 
            (38023127*S2)/209952 + 12*Zeta[3]))) + 
      x^28*(ca*cf*((11*lm^2*Log[x])/4 + lm*(-99/132496 + (6421*Log[x])/546 - 
            55*Log[x]^2) + 
          (3*(-120665361104369465170829838016645242311764966369/
              450955333772923329160355298977967744000 + (7587098764289*OepS2)/
              2636009007156 + (415537493983381987564697531*Pi^2)/
              5488381633619364480 - (144154876521491*Pi^4)/1265284323434880 - 
             (96987482178757741*S2)/902100860226720 - (7587098764289*T1ep)/
              1757339338104 + (83001664842879085993*Log[x]^2)/91356530265 + 
             (84663923971*Log[x]^3)/273 + Log[x]*
              (3050093838227336307514918844804485691/
                24064040586546254265808636800 + (507983489317*Pi^2)/6552 + 
               (7587098764289*S2)/65086642152 - 4*Zeta[3]) - 
             (7084918861283308273633*Zeta[3])/45690822790704))/2) + 
        cf*(508330969704585008245733066933357077672661322670261/
           5665358672074773561929209505549616402000000 - 
          (2243258641711*OepS2)/94143178827 - (364937600950108250332381*Pi^2)/
           15077971520932320 + (42621914192509*Pi^4)/45188725836960 + 
          (111375157711161629*S2)/48326831797860 + (2243258641711*T1ep)/
           62762119218 - lm^2*Log[x] + 
          (-28030677785990493777198295816406033941/
             1389962783879491747084138980000 - (154717081*Pi^2)/6 - 
            (2243258641711*S2)/2324522934)*Log[x] - 
          (42278977219885005868*Log[x]^2)/145568097675 - 103144692*Log[x]^3 + 
          lm*(9/33124 - (820*Log[x])/273 + 20*Log[x]^2) + 
          (6473535426560925155*Zeta[3])/125524238436 + 
          nl*(-(lm^2*Log[x])/2 + lm*(9/66248 - (410*Log[x])/273 + 
              10*Log[x]^2) + (3*(-4106121836154133/599950102752000 + 
               (181*Pi^2)/4914 + (1125475189/142702560 + Pi^2)*Log[x] + 
               (92*Log[x]^2)/9 - (40*Log[x]^3)/9 + 6*Zeta[3]))/2)) + 
        cf^2*((9*lm^2*Log[x])/2 + lm*(-81/66248 - (813*Log[x])/364 - 
            90*Log[x]^2) + 
          (3*(965322864101212371341511622053863118866889515417/
              1803821335091693316641421195911870976000 - 
             (7587098764289*OepS2)/1318004503578 - 
             (103884373444355350552891253*Pi^2)/686047704202420560 + 
             (144154876521491*Pi^4)/632642161717440 + (96987482178757741*S2)/
              451050430113360 + (7587098764289*T1ep)/878669669052 - 
             (13551293325220121831*Log[x]^2)/7457675940 - 
             (72569043032*Log[x]^3)/117 + (7084919135428245017857*Zeta[3])/
              22845411395352 + Log[x]*
              (-42701307919075849317640868790938680033/
                168448284105823779860660457600 - (507983482765*Pi^2)/3276 - 
               (7587098764289*S2)/32543321076 + 8*Zeta[3])))/2)) + 
      x^20*(ca*cf*((11*lm^2*Log[x])/4 + lm*(-209/97200 + (1589*Log[x])/135 - 
            (121*Log[x]^2)/3) + (3*(-22106251774616673741252427537/
              165795802929445328640000 + (1677317401*OepS2)/860934420 + 
             (75290696502318989*Pi^2)/2121342410880 - (31869030619*Pi^4)/
              413248521600 - (5154463799657*S2)/80353879200 - 
             (1677317401*T1ep)/573956280 + (41445408955207*Log[x]^2)/
              97297200 + (62918603*Log[x]^3)/405 + Log[x]*
              (64879476896879392469/5477188252536000 + (125817073*Pi^2)/
                3240 + (1677317401*S2)/21257640 - 4*Zeta[3]) - 
             (89210062236917*Zeta[3])/1147912560))/2) + 
        cf*(162680414770750104877157055584957/2586588534535106873339568000 - 
          (214104749*OepS2)/14348907 - (256374291157057*Pi^2)/16070775840 + 
          (4067990231*Pi^4)/6887475360 + (86235152342*S2)/55801305 + 
          (214104749*T1ep)/9565938 - lm^2*Log[x] + 
          (2240897395271650570797499/133333567079726613600 - 
            (337705*Pi^2)/18 - (214104749*S2)/354294)*Log[x] - 
          (1675980103913*Log[x]^2)/8729721 - (675232*Log[x]^3)/9 + 
          lm*(19/24300 - (406*Log[x])/135 + (44*Log[x]^2)/3) + 
          (713151237097*Zeta[3])/19131876 + nl*(-12546276727/1778112000 + 
            (89*Pi^2)/1620 - (lm^2*Log[x])/2 + (1048577/141120 + 
              (19*Pi^2)/18)*Log[x] + (34*Log[x]^2)/3 - (44*Log[x]^3)/9 + 
            lm*(19/48600 - (203*Log[x])/135 + (22*Log[x]^2)/3) + 
            (19*Zeta[3])/3)) + cf^2*((9*lm^2*Log[x])/2 + 
          lm*(-19/5400 - (133*Log[x])/60 - 66*Log[x]^2) + 
          (3*(1105273369952377886669846201/4144895073236133216000 - 
             (1677317401*OepS2)/430467210 - (75290615924171687*Pi^2)/
              1060671205440 + (31869030619*Pi^4)/206624260800 + 
             (5154463799657*S2)/40176939600 + (1677317401*T1ep)/286978140 - 
             (20725352206759*Log[x]^2)/24324300 - (25153048*Log[x]^3)/81 + 
             (89214908978837*Zeta[3])/573956280 + Log[x]*
              (-453603886594730245709/19170158883876000 - (125814793*Pi^2)/
                1620 - (1677317401*S2)/10628820 + 8*Zeta[3])))/2)) + 
      x^12*(ca*cf*(-8347346586907/80015040000 + (27319*OepS2)/17496 + 
          (4003489507*Pi^2)/97977600 - (519061*Pi^4)/8398080 - 
          (127285*S2)/3888 - (27319*T1ep)/11664 + (11*lm^2*Log[x])/4 + 
          (2686829*Log[x]^2)/5040 + (9487*Log[x]^3)/30 + 
          lm*(-121/10800 + (1063*Log[x])/90 - (77*Log[x]^2)/3) + 
          Log[x]*(-28321662139/95256000 + (52477*Pi^2)/720 + (27319*S2)/432 - 
            6*Zeta[3]) - (22140823*Zeta[3])/116640) + 
        cf*(-8879402758922099/438713512637400 - (46331*OepS2)/6561 + 
          (5469587*Pi^2)/1049760 + (880289*Pi^4)/3149280 + 
          (4882415*S2)/5832 + (46331*T1ep)/4374 - lm^2*Log[x] + 
          (4788765199784/18261468225 - (941*Pi^2)/18 - (46331*S2)/162)*
           Log[x] - (52318379*Log[x]^2)/270270 - (1784*Log[x]^3)/9 + 
          lm*(11/2700 - (136*Log[x])/45 + (28*Log[x]^2)/3) - 
          (265025*Zeta[3])/8748 + nl*(-496307/129600 + (29*Pi^2)/540 - 
            (lm^2*Log[x])/2 + ((2159 + 440*Pi^2)*Log[x])/720 + 
            (22*Log[x]^2)/3 - (28*Log[x]^3)/9 + 
            lm*(11/5400 - (68*Log[x])/45 + (14*Log[x]^2)/3) + 
            (11*Zeta[3])/3)) + cf^2*((9*lm^2*Log[x])/2 + 
          lm*(-11/600 - (43*Log[x])/20 - 42*Log[x]^2) + 
          (3*(6423179858993/48009024000 - (27319*OepS2)/13122 - 
             (3997719877*Pi^2)/73483200 + (519061*Pi^4)/6298560 + 
             (127285*S2)/2916 + (27319*T1ep)/8748 - (7376777*Log[x]^2)/9450 - 
             (15032*Log[x]^3)/45 + (22568503*Zeta[3])/87480 + 
             Log[x]*(12000383443/28576800 - (52037*Pi^2)/540 - (27319*S2)/
                324 + 8*Zeta[3])))/2)) + 
      x^2*(ca*cf*(1667/192 + (3*B4)/4 + (17*Pi^2)/32 + (2*Pi^4)/45 + 
          (1377*S2)/32 + (2*Pi^2*Log[2]^2)/3 - (2*Log[2]^4)/3 - 
          (43*Log[x]^2)/6 + (11*Log[x]^3)/3 + (11*lm^2*(1 + 4*Log[x]))/16 + 
          (lm*(621 + 22*Pi^2 + 1032*Log[x] - 792*Log[x]^2))/144 - 
          16*PolyLog[4, 1/2] + Log[x]*(85/48 - (9*Zeta[3])/2) - 
          (365*Zeta[3])/24) + cf^2*(-2827/576 - (3*B4)/2 + (17*Pi^2)/108 - 
          (2*Pi^4)/15 - (1377*S2)/16 - (4*Pi^2*Log[2]^2)/3 + (4*Log[2]^4)/3 + 
          (39*Log[x]^2)/4 + 6*Log[x]^3 + (9*lm^2*(1 + 4*Log[x]))/8 + 
          lm*(27/16 + Pi^2/4 - (39*Log[x])/4 - 9*Log[x]^2) + 
          32*PolyLog[4, 1/2] + (191*Zeta[3])/6 + 
          Log[x]*(-101/16 + 9*Zeta[3])) + cf*(lm^2*(-1/4 - Log[x]) + 
          lm*(-5/4 - Pi^2/18 - (4*Log[x])/3 + 2*Log[x]^2) + 
          (3221 - 210*Pi^2 + 39366*S2 - 2148*Log[x] + 1584*Log[x]^2 - 
            288*Log[x]^3 - 8640*Zeta[3])/432 + nl*((lm^2*(-1 - 4*Log[x]))/8 + 
            lm*(-5/8 - Pi^2/36 - (2*Log[x])/3 + Log[x]^2) + 
            (-393 + 4*Pi^2 - 60*Log[x] + 192*Log[x]^2 - 192*Log[x]^3 + 
              96*Zeta[3])/288))) + 
      x^4*(cf*(-(lm^2*Log[x]) + lm*(1/4 - (10*Log[x])/3 + 4*Log[x]^2) + 
          (3*(-1725443/162000 - (2*OepS2)/3 + (773*Pi^2)/432 + 
             (19*Pi^4)/720 + (585*S2)/4 + T1ep + (5491/675 - (5*Pi^2)/9 - 27*
                S2)*Log[x] + (109*Log[x]^2)/45 - (8*Log[x]^3)/9 - 
             (121*Zeta[3])/6))/2 + nl*(-7/48 + Pi^2/36 - (lm^2*Log[x])/2 + 
            ((-45 + 4*Pi^2)*Log[x])/24 + (10*Log[x]^2)/3 - (4*Log[x]^3)/3 + 
            lm*(1/8 - (5*Log[x])/3 + 2*Log[x]^2) + Zeta[3])) + 
        ca*cf*((11*lm^2*Log[x])/4 + lm*(-11/16 + (38*Log[x])/3 - 
            11*Log[x]^2) + (-29445 + 3360*OepS2 - 10985*Pi^2 - 133*Pi^4 + 
            311040*S2 - 5040*T1ep - 291840*Log[x]^2 + 84480*Log[x]^3 - 
            102360*Zeta[3] - 240*Log[x]*(-687 + 61*Pi^2 - 567*S2 + 
              288*Zeta[3]))/11520) + cf^2*((9*lm^2*Log[x])/2 - 
          (3*lm*(3 + 2*Log[x] + 48*Log[x]^2))/8 + 
          (3*(-71/72 - (7*OepS2)/18 + (899*Pi^2)/576 + (133*Pi^4)/8640 - 
             36*S2 + (7*T1ep)/12 + Log[x]^2 + 16*Log[x]^3 + 
             (937*Zeta[3])/72 + (Log[x]*(165 + 71*Pi^2 - 567*S2 + 
                288*Zeta[3]))/36))/2)) + 
      x^6*(ca*cf*(166849/6912 + (16*OepS2)/27 - (5425*Pi^2)/1728 - 
          (19*Pi^4)/810 + (101*S2)/8 - (8*T1ep)/9 + (11*lm^2*Log[x])/4 - 
          (104*Log[x]^2)/3 + (160*Log[x]^3)/9 + 
          (lm*(-55 + 5208*Log[x] - 6336*Log[x]^2))/432 + 
          Log[x]*(-3227/288 + (43*Pi^2)/72 + 24*S2 - 6*Zeta[3]) - 
          (383*Zeta[3])/18) + cf^2*((9*lm^2*Log[x])/2 + 
          lm*(-5/24 - (7*Log[x])/4 - 24*Log[x]^2) + 
          (-700315 - 15360*OepS2 + 84480*Pi^2 + 608*Pi^4 - 327240*S2 + 
            23040*T1ep + 95760*Log[x]^2 + 282240*Log[x]^3 - 
            90*Log[x]*(-8011 + 112*Pi^2 + 6912*S2 - 1728*Zeta[3]) + 
            600120*Zeta[3])/12960) + cf*(-3233892733/111132000 - 
          (61*OepS2)/27 + (6409*Pi^2)/864 + (1159*Pi^4)/12960 + (2885*S2)/8 + 
          (61*T1ep)/18 - lm^2*Log[x] + (1453253/66150 - (41*Pi^2)/18 - 
            (183*S2)/2)*Log[x] + (8321*Log[x]^2)/630 - (44*Log[x]^3)/9 + 
          (lm*(5 - 336*Log[x] + 576*Log[x]^2))/108 - (207*Zeta[3])/4 + 
          nl*(-(lm^2*Log[x])/2 + lm*(5/216 - (14*Log[x])/9 + 
              (8*Log[x]^2)/3) + (-553 + 20*Pi^2 + 30*(-7 + 4*Pi^2)*Log[x] + 
              1872*Log[x]^2 - 768*Log[x]^3 + 720*Zeta[3])/432))) + 
      x^8*(ca*cf*((11*lm^2*Log[x])/4 + lm*(-77/1728 + (857*Log[x])/72 - 
            (55*Log[x]^2)/3) + (76673505 + 1692960*OepS2 - 6097407*Pi^2 - 
            67013*Pi^4 - 4026510*S2 - 2539440*T1ep - 41938560*Log[x]^2 + 
            71539200*Log[x]^3 + 3240*Log[x]*(-30725 + 2862*Pi^2 + 21162*S2 - 
              3456*Zeta[3]) - 71165880*Zeta[3])/1866240) + 
        cf^2*((9*lm^2*Log[x])/2 + lm*(-7/96 - 2*Log[x] - 30*Log[x]^2) + 
          (-81731640 - 1692960*OepS2 + 6197631*Pi^2 + 67013*Pi^4 + 
            4026510*S2 + 2539440*T1ep - 29535840*Log[x]^2 + 
            2903040*Log[x]^3 - 1080*Log[x]*(-120727 + 7914*Pi^2 + 63486*S2 - 
              10368*Zeta[3]) + 75520440*Zeta[3])/933120) + 
        cf*(-27367344769/571536000 - (301*OepS2)/81 + (33545*Pi^2)/2592 + 
          (5719*Pi^4)/38880 + (4097*S2)/8 + (301*T1ep)/54 - lm^2*Log[x] + 
          (10349461/226800 - (121*Pi^2)/18 - (301*S2)/2)*Log[x] + 
          (4471*Log[x]^2)/180 - (184*Log[x]^3)/9 + 
          lm*(7/432 - (55*Log[x])/18 + (20*Log[x]^2)/3) - 
          (7375*Zeta[3])/108 + nl*(-(lm^2*Log[x])/2 + 
            lm*(7/864 - (55*Log[x])/36 + (10*Log[x]^2)/3) + 
            (-3757 + 88*Pi^2 + 96*(13 + 7*Pi^2)*Log[x] + 9216*Log[x]^2 - 
              3840*Log[x]^3 + 4032*Zeta[3])/1728))) + 
      x^14*(cf*(21882602296492111/74780712381375 - (175207*OepS2)/19683 - 
          (322591579*Pi^2)/4408992 + (3328933*Pi^4)/9447840 + 
          (88186271*S2)/87480 + (175207*T1ep)/13122 - lm^2*Log[x] + 
          (3374718939341/4980400425 - (3211*Pi^2)/18 - (175207*S2)/486)*
           Log[x] - (29563757*Log[x]^2)/24570 - (6304*Log[x]^3)/9 + 
          lm*(13/5292 - (190*Log[x])/63 + (32*Log[x]^2)/3) + 
          (5104187*Zeta[3])/26244 + nl*(-21047903/4536000 + (41*Pi^2)/756 - 
            (lm^2*Log[x])/2 + ((14809 + 2600*Pi^2)*Log[x])/3600 + 
            (25*Log[x]^2)/3 - (32*Log[x]^3)/9 + lm*(13/10584 - 
              (95*Log[x])/63 + (16*Log[x]^2)/3) + (13*Zeta[3])/3)) + 
        cf^2*((9*lm^2*Log[x])/2 - (325055279*Log[x]^2)/52920 - 
          (22738*Log[x]^3)/9 + lm*(-13/1176 - (61*Log[x])/28 - 48*Log[x]^2) + 
          Log[x]*(336265978891/225042300 - (27575*Pi^2)/42 - 
            (784016*S2)/5103 + 12*Zeta[3]) + (-5163215769600*OepS2 - 
            670365909190740*Pi^2 + 204377290880*Pi^4 + 
            9*(271415005289977 + 14634589257120*S2 + 860535961600*T1ep + 
              216020894555200*Zeta[3]))/1361055830400) + 
        ca*cf*((11*lm^2*Log[x])/4 + (159372823*Log[x]^2)/52920 + 
          (169321*Log[x]^3)/126 + lm*(-143/21168 + (743*Log[x])/63 - 
            (88*Log[x]^2)/3) + Log[x]*(-749520229441/1028764800 + 
            (82907*Pi^2)/252 + (392008*S2)/5103 - 6*Zeta[3]) + 
          (206528630784000*OepS2 + 26820978659756400*Pi^2 - 
            8175091635200*Pi^4 - 9*(10913948420987887 + 585383570284800*S2 + 
              34421438464000*T1ep + 8588409928000000*Zeta[3]))/
           108884466432000)) + x^10*(cf*(-323092231612649/5325000912000 - 
          (3883*OepS2)/729 + (1869547*Pi^2)/116640 + (73777*Pi^4)/349920 + 
          (145045*S2)/216 + (3883*T1ep)/486 - lm^2*Log[x] + 
          (6893875081/64033200 - (109*Pi^2)/6 - (3883*S2)/18)*Log[x] + 
          (19631*Log[x]^2)/13860 - 64*Log[x]^3 + lm*(3/400 - (91*Log[x])/30 + 
            8*Log[x]^2) - (69625*Zeta[3])/972 + 
          nl*(-43361/14400 + (19*Pi^2)/360 - (lm^2*Log[x])/2 + 
            ((337 + 90*Pi^2)*Log[x])/180 + (19*Log[x]^2)/3 - (8*Log[x]^3)/3 + 
            lm*(3/800 - (91*Log[x])/60 + 4*Log[x]^2) + 3*Zeta[3])) + 
        cf^2*((9*lm^2*Log[x])/2 - (270319*Log[x]^2)/1200 - (265*Log[x]^3)/3 - 
          (3*lm*(9 + 560*Log[x] + 9600*Log[x]^2))/800 + 
          Log[x]*(188917963/648000 - (1091*Pi^2)/30 - (4487*S2)/45 + 
            12*Zeta[3]) + (-6030528000*OepS2 - 15061354500*Pi^2 + 
            238708400*Pi^4 + 21*(-8020055831 + 4037306400*S2 + 
              430752000*T1ep + 17852664000*Zeta[3]))/2449440000) + 
        ca*cf*((11*lm^2*Log[x])/4 + lm*(-33/1600 + (1421*Log[x])/120 - 
            22*Log[x]^2) + (2359583460*Log[x]^2 + 3369729600*Log[x]^3 + 
            27*Log[x]*(-167089039 + 24213600*Pi^2 + 64612800*S2 - 
              7776000*Zeta[3]) - 5*(-8615040*OepS2 - 21930015*Pi^2 + 
              341012*Pi^4 + 24*(-9075694 + 5046633*S2 + 538440*T1ep + 
                21441030*Zeta[3])))/34992000))))
 

rhoMsBar2Exp1 = 2/3*(nc*((ca*cf*(-13/36 + (11*lm)/6) + cf^2*(-31/4 + 3*lm) + 
        cf*(5/9 - (2*lm)/3 + (5/18 - lm/3)*nl))*(1 - x)^3 + 
      ca*cf*((3*B4)/2 + (13*Pi^4)/120 + Pi^2*Log[2]^2 - Log[2]^4 - 
        24*PolyLog[4, 1/2]) + cf^2*(-3*B4 - (13*Pi^4)/60 - 2*Pi^2*Log[2]^2 + 
        2*Log[2]^4 + 48*PolyLog[4, 1/2]) + 
      (1 - x)*(ca*cf*((-3*B4)/2 - (13*Pi^4)/120 - Pi^2*Log[2]^2 + Log[2]^4 + 
          24*PolyLog[4, 1/2]) + cf^2*(3*B4 + (13*Pi^4)/60 + 2*Pi^2*Log[2]^2 - 
          2*(Log[2]^4 + 24*PolyLog[4, 1/2]))) + 
      (1 - x)^2*(ca*cf*(-1213/288 + (3*B4)/4 - (13*lm)/36 + (11*lm^2)/12 + 
          (13*Pi^4)/240 + (Pi^2*Log[2]^2)/2 - Log[2]^4/2 - 
          12*PolyLog[4, 1/2] + (91*Zeta[3])/64) + 
        cf^2*(-21/4 - (3*B4)/2 - (31*lm)/4 + (3*lm^2)/2 - (13*Pi^4)/120 - 
          Pi^2*Log[2]^2 + Log[2]^4 + 24*PolyLog[4, 1/2] + (293*Zeta[3])/32) + 
        cf*((5*lm)/9 - lm^2/3 + (-5/72 + (5*lm)/18 - lm^2/6)*nl - 
          (7*(26 + 45*Zeta[3]))/144)) + (1 - x)^4*
       (cf*(-7627/12960 - (131*lm)/540 + lm^2/60 + 
          (-101/864 - (131*lm)/1080 + lm^2/120)*nl + (91*Zeta[3])/192) + 
        cf^2*((65*lm)/48 - (3*lm^2)/40 + (79*(-703 + 561*Zeta[3]))/2880) + 
        ca*cf*((263*lm)/432 - (11*lm^2)/240 + (-8278 + 12321*Zeta[3])/
           5760)) + (1 - x)^5*(cf*((-53*lm)/540 + lm^2/60 + 
          (-509/4320 - (53*lm)/1080 + lm^2/120)*nl + 
          (7*(-2186 + 1755*Zeta[3]))/25920) + 
        ca*cf*((457*lm)/2160 - (11*lm^2)/240 + (-29114 + 36963*Zeta[3])/
           17280) + cf^2*((169*lm)/240 - (3*lm^2)/40 + 
          (-52237 + 44319*Zeta[3])/2880)) + 
      (1 - x)^7*(cf*((-26*lm)/4725 + lm^2/105 + 
          (-33101/567000 - (13*lm)/4725 + lm^2/210)*nl + 
          (-756074 + 635775*Zeta[3])/1612800) + 
        ca*cf*((-86*lm)/4725 - (11*lm^2)/420 + (-90348254 + 84241485*Zeta[3])/
           29030400) + cf^2*((367*lm)/2100 - (3*lm^2)/70 + 
          (-733535842 + 616358835*Zeta[3])/14515200)) + 
      (1 - x)^6*(cf*(-510019/967680 - (433*lm)/12600 + (11*lm^2)/840 + 
          (-131303/1512000 - (433*lm)/25200 + (11*lm^2)/1680)*nl + 
          (889*Zeta[3])/2048) + ca*cf*((2453*lm)/50400 - (121*lm^2)/3360 + 
          (-47145514 + 48779775*Zeta[3])/19353600) + 
        cf^2*((2021*lm)/5600 - (33*lm^2)/560 + 
          (-993148114 + 839726595*Zeta[3])/29030400)) + 
      (1 - x)^9*(cf*(-262399397/677376000 + (991*lm)/88200 + lm^2/210 + 
          (-282467/12348000 + (991*lm)/176400 + lm^2/420)*nl + 
          (2737*Zeta[3])/8192) + cf^2*((479*lm)/19600 - (3*lm^2)/140 + 
          (-1279260563662 + 1067386092885*Zeta[3])/13005619200) + 
        ca*cf*((-16781*lm)/352800 - (11*lm^2)/840 + 
          (-1292345766898 + 1109132937675*Zeta[3])/303464448000)) + 
      (1 - x)^8*(cf*(-122710351/290304000 + (601*lm)/88200 + (17*lm^2)/2520 + 
          (-8287577/222264000 + (601*lm)/176400 + (17*lm^2)/5040)*nl + 
          (8869*Zeta[3])/24576) + cf^2*((2963*lm)/39200 - (17*lm^2)/560 + 
          (-935626646042 + 782653708935*Zeta[3])/13005619200) + 
        ca*cf*((-14941*lm)/352800 - (187*lm^2)/10080 + 
          (-3380681634418 + 2986079035275*Zeta[3])/910393344000)) + 
      (1 - x)^11*(cf*(-357954044629/1072963584000 + (21041*lm)/1871100 + 
          (17*lm^2)/6930 + (-1102054/147349125 + (21041*lm)/3742200 + 
            (17*lm^2)/13860)*nl + (38059*Zeta[3])/131072) + 
        cf^2*((-622*lm)/51975 - (17*lm^2)/1540 + 
          (-3609653661671174 + 3005497866947625*Zeta[3])/21459271680000) + 
        ca*cf*((-295711*lm)/7484400 - (17*lm^2)/2520 + 
          (-14176793751969158 + 11881516779995625*Zeta[3])/
           2703868231680000)) + (1 - x)^10*
       (cf*(-220050891107/613122048000 + (89927*lm)/7484400 + 
          (47*lm^2)/13860 + (-126901261/9430344000 + (89927*lm)/14968800 + 
            (47*lm^2)/27720)*nl + (163275*Zeta[3])/524288) + 
        cf^2*((-1097*lm)/1663200 - (47*lm^2)/3080 + 
          (-11183318144847134 + 9318330609202125*Zeta[3])/85837086720000) + 
        ca*cf*((-1344517*lm)/29937600 - (47*lm^2)/5040 + 
          (-51530312707196318 + 43547329208746125*Zeta[3])/
           10815472926720000)) + (1 - x)^13*
       (cf*(-1113628911766537/3965673406464000 + (10660483*lm)/1248647400 + 
          (491*lm^2)/360360 + (-10349611/6837038208 + (10660483*lm)/
             2497294800 + (491*lm^2)/720720)*nl + (1540819*Zeta[3])/
           6291456) + ca*cf*((-12825793*lm)/454053600 - (491*lm^2)/131040 + 
          (-16731023084911613246 + 13927839611505535485*Zeta[3])/
           2722038226196889600) + cf^2*((-9411761*lm)/554954400 - 
          (491*lm^2)/80080 + (-163373772891394952606 + 135948152227010525325*
             Zeta[3])/618645051408384000)) + 
      (1 - x)^12*(cf*((70667*lm)/7114800 + (163*lm^2)/90090 + 
          (-9665553277/2563889328000 + (70667*lm)/14229600 + 
            (163*lm^2)/180180)*nl + (-748777737854 + 652572892125*Zeta[3])/
           2427715584000) + ca*cf*((-3394973*lm)/100900800 - 
          (163*lm^2)/32760 + (-388050392650421524678 + 323785976238457511625*
             Zeta[3])/68050955654922240000) + cf^2*((-998873*lm)/61661600 - 
          (163*lm^2)/20020 + (-657742848253666457582 + 547442048279497573125*
             Zeta[3])/3093225257041920000)) + 
      (1 - x)^14*(cf*(-74545579155169961/299949115834368000 + 
          (436837*lm)/60231600 + (151*lm^2)/144144 + 
          (-40375393397/236343979872000 + (436837*lm)/120463200 + 
            (151*lm^2)/288288)*nl + (87375925*Zeta[3])/402653184) + 
        ca*cf*(-847226943938028730504133/128806848863636815872000 - 
          (5688557*lm)/240926400 - (151*lm^2)/52416 + (2204331044549*Zeta[3])/
           403055837184) + cf^2*((-431999*lm)/26769600 - (151*lm^2)/32032 + 
          (-640164402839332182006518 + 532638645898900080164625*Zeta[3])/
           1981643828671335628800)) + (1 - x)^15*
       (cf*(-498838977901586797/2382929086906368000 + (377563*lm)/61486425 + 
          (37*lm^2)/45045 + (1279369009/2110214106000 + (377563*lm)/
             122972850 + (37*lm^2)/90090)*nl + (12304341*Zeta[3])/67108864) + 
        ca*cf*((-4860263*lm)/245945700 - (37*lm^2)/16380 + 
          (-4505176686361864059922534 + 3743878904526521749613625*Zeta[3])/
           644034244318184079360000) + cf^2*((-401591*lm)/27327300 - 
          (37*lm^2)/10010 + (-9661494228831316850595202 + 
            8038179841327467312355875*Zeta[3])/24770547858391695360000)) + 
      (1 - x)^17*(cf*(-125018145037487663293/1250016515302883328000 + 
          (102682121*lm)/23156733600 + (163*lm^2)/306306 + 
          (136211853551/108373513248000 + (102682121*lm)/46313467200 + 
            (163*lm^2)/612612)*nl + (145762939*Zeta[3])/1610612736) + 
        ca*cf*((-1302022531*lm)/92626934400 - (163*lm^2)/111384 + 
          (-6232010015073126041460974554 + 5177495152044716636686926375*
             Zeta[3])/800810580363639175249920000) + 
        cf^2*((-59552321*lm)/5145940800 - (163*lm^2)/68068 + 
          (-73427308876932911424544950626 + 61086635943007537208693221875*
             Zeta[3])/133468430060606529208320000)) + 
      (1 - x)^16*(cf*(-108235840658362742917/673085815932321792000 + 
          (120577199*lm)/23156733600 + (1607*lm^2)/2450448 + 
          (4325529748559/4172380260048000 + (120577199*lm)/46313467200 + 
            (1607*lm^2)/4900896)*nl + (228676273*Zeta[3])/1610612736) + 
        ca*cf*((-1538955289*lm)/92626934400 - (1607*lm^2)/891072 + 
          (-41470481093262514649240858074 + 34455001941803377527904686375*
             Zeta[3])/5605674062545474226749440000) + 
        cf^2*((-33712837*lm)/2572970400 - (1607*lm^2)/544544 + 
          (-144990584103703819535139249782 + 120625045552201408471248620625*
             Zeta[3])/311426336808081901486080000)) + 
      (1 - x)^19*(cf*(185874336617187398132527/2512255414310950404096000 + 
          (17002792*lm)/5194183995 + (5314*lm^2)/14549535 + 
          (7849522227690809/5727635001980892000 + (8501396*lm)/5194183995 + 
            (2657*lm^2)/14549535)*nl - (718615037*Zeta[3])/12884901888) + 
        ca*cf*((-53397521*lm)/5194183995 - (2657*lm^2)/2645370 + 
          (-35475322809822928859089221338 + 29475829883039267899184148135*
             Zeta[3])/4187857996930231153616486400) + 
        cf^2*((-10362949*lm)/1154263110 - (2657*lm^2)/1616615 + 
          (-9680488858421347267246628318518 + 8053363739483860411665588762225*
             Zeta[3])/12933090872872772680286208000)) + 
      (1 - x)^18*(cf*(-122496326626501010177/5254390409016366858240 + 
          (541064597*lm)/142469046720 + (5099*lm^2)/11639628 + 
          (167074541339851/123673630272192000 + (541064597*lm)/284938093440 + 
            (5099*lm^2)/23279256)*nl + (889419041*Zeta[3])/34359738368) + 
        cf^2*((-322623437*lm)/31659788160 - (5099*lm^2)/2586584 + 
          (-13332193579047943940995016718802 + 
            11091361248965454541520278117275*Zeta[3])/
           20692945396596436288457932800) + 
        ca*cf*((-6825475207*lm)/569876186880 - (5099*lm^2)/4232592 + 
          (-28643323027511346870819731901766 + 
            23797462282664945413752972383625*Zeta[3])/
           3517800717421394169037848576000)) + 
      (1 - x)^20*(cf*(30900501384055440263078591/
           156055395147786095689728000 + (7254102841*lm)/2554316364600 + 
          (17921*lm^2)/58198140 + (14416858809063620089/
             10703265220172292768000 + (7254102841*lm)/5108632729200 + 
            (17921*lm^2)/116396280)*nl - (65793708253*Zeta[3])/
           412316860416) + ca*cf*((-750469991*lm)/84440210400 - 
          (17921*lm^2)/21162960 + (-481421417309581371948660393930646 + 
            400039758769646868538293085082625*Zeta[3])/
           54975154068793216234747330560000) + 
        cf^2*((-2250584213*lm)/283812929400 - (17921*lm^2)/12932920 + 
          (-288515930313138147003626063348261506 + 
            240020238330689676893503479458323875*Zeta[3])/
           334191068155032446058595614720000))))
 

rhoOnShell0 = nc*((2*(1 + X^2))/3 + (8*X^2*HPL[{0}, X])/(3 - 3*X^2))
 
rhoOnShell1 = cf*nc*((Pi^2*(-1 + X^2) - 3*(1 + X^2))/9 + 
      (4*X^2*HPL[{0}, X])/(3 - 3*X^2) - (4*(-1 + X^2)*HPL[{-1, 0}, X])/3 + 
      (8*X^2*(3 + X^4)*HPL[{0, 0}, X])/(3*(-1 + X^2)^2) + 
      (4*(-1 + X^2)*HPL[{1, 0}, X])/3)
 
 
rhoOnShell2Exp0 = nc*(cf*(23/9 + (3473*X^2)/648 - (2121653*X^4)/162000 - 
        (8026884581*X^6)/333396000 - (31170910639*X^8)/857304000 - 
        (358269233365619*X^10)/7987501368000 - (94171548247865317*X^12)/
         5264562151648800 + (5475853419488548141*X^14)/28715793554448000 + 
        (74188250820529333980733651*X^16)/55867954718269197504000 + 
        (694967637685143695400792659308567*X^18)/
         93117187243263847440224448000 + (325326781326746571173736104040709*
          X^20)/7759765603605320620018704000 + 
        (84893680669019796983048383657759158437*X^22)/
         347611750728379127940235149864000 + 
        (38171181711603557366796976820978658097711801*X^24)/
         25810172491582150249562459877402000000 + 
        (31203398939350390844210413466207004324047879*X^26)/
         3366544238032454380377712157922000000 + 
        (508330932424376156856495560238081937168899844988011*X^28)/
         8498038008112160342893814258324424603000000 + 
        (7215395139305612173143122949088609537471991725612106150813*X^30)/
         18227883621576194551810772681021491201054056000000 - 
        (2*OepS2*X^4*(847288609443 + 1914244636149*X^2 + 3148566314103*X^4 + 
           4513061276361*X^6 + 5983192891953*X^8 + 7542086846247*X^10 + 
           9177489172665*X^12 + 10880110754919*X^14 + 12642671323701*X^16 + 
           14459312693331*X^18 + 16325219636073*X^20 + 18236364305751*X^22 + 
           20189327775399*X^24 + 22181171683297*X^26))/2541865828329 + 
        (19*Pi^4*X^4*(847288609443 + 1914244636149*X^2 + 3148566314103*X^4 + 
           4513061276361*X^6 + 5983192891953*X^8 + 7542086846247*X^10 + 
           9177489172665*X^12 + 10880110754919*X^14 + 12642671323701*X^16 + 
           14459312693331*X^18 + 16325219636073*X^20 + 18236364305751*X^22 + 
           20189327775399*X^24 + 22181171683297*X^26))/610047798798960 + 
        (T1ep*X^4*(847288609443 + 1914244636149*X^2 + 3148566314103*X^4 + 
           4513061276361*X^6 + 5983192891953*X^8 + 7542086846247*X^10 + 
           9177489172665*X^12 + 10880110754919*X^14 + 12642671323701*X^16 + 
           14459312693331*X^18 + 16325219636073*X^20 + 18236364305751*X^22 + 
           20189327775399*X^24 + 22181171683297*X^26))/847288609443 + 
        Pi^2*(1/108 - X/6 - (59*X^2)/108 + X^3/6 + (749*X^4)/432 + 
          (2*X^5)/3 + (6121*X^6)/1296 + (2*X^7)/3 + (32681*X^8)/3888 + 
          (2*X^9)/3 + (1830667*X^10)/174960 + (2*X^11)/3 + 
          (5119667*X^12)/1574640 + (2*X^13)/3 - (324061243*X^14)/6613488 + 
          (2*X^15)/3 - (32882135879*X^16)/99202320 + (2*X^17)/3 - 
          (45054133443547*X^18)/24106163760 + (2*X^19)/3 - 
          (256379648082337*X^20)/24106163760 + (2*X^21)/3 - 
          (150223282262839381*X^22)/2386510212240 + (2*X^23)/3 - 
          (24960824289125289181*X^24)/64435775730480 + (2*X^25)/3 - 
          (6194808021345569704597*X^26)/2512995253488720 + (2*X^27)/3 - 
          (364937605976098757309821*X^28)/22616957281398480 + (2*X^29)/3 - 
          (2001691415724772122179033*X^30)/18504783230235120) - 
        (X^2*(2356270471083649927624020735451920000 - 
           6453937023653883039853984414437830400*X^2 - 
           8857441036644921391062257534020294800*X^4 - 
           15251968800499859148754465807466005200*X^6 - 
           48011990573810312623011244019647215600*X^8 - 
           143160813285002430485944094355666207600*X^10 - 
           413908397381618540620572192266775541500*X^12 - 
           1233496643938168392990243997814573513100*X^14 - 
           3799847359412774052610574080142545918500*X^16 - 
           11245735179410156093800975581930459023100*X^18 - 
           24234536659689127237688479513944013078800*X^20 + 
           46325472662224571648553778805659375681464*X^22 + 
           1275524568660245504448489040987341287069448*X^24 + 
           13576594925243562245511698554947669482161344*X^26 + 
           120220173135794183117054801510596445259787085*X^28 + 
           37401118588629363930540011673840000*Pi^2*(12 + 33*X^2 - 18*X^3 + 
             65*X^4 - 54*X^5 + 157*X^6 - 90*X^7 + 375*X^8 - 126*X^9 + 
             1001*X^10 - 162*X^11 + 3283*X^12 - 198*X^13 + 13341*X^14 - 
             234*X^15 + 63671*X^16 - 270*X^17 + 337813*X^18 - 306*X^19 + 
             1924407*X^20 - 342*X^21 + 11541611*X^22 - 378*X^23 + 
             72027625*X^24 - 414*X^25 + 464151399*X^26 - 450*X^27 + 
             3071706245*X^28))*Log[X])/
         1009830201892992826124580315193680000 + 
        ((5*X^2)/9 - (122*X^4)/45 - (622*X^6)/945 + (4171*X^8)/630 - 
          (204613*X^10)/20790 - (57073364*X^12)/405405 - 
          (18877411*X^14)/23166 - (58958191*X^16)/13860 - 
          (11918951656369*X^18)/523783260 - (33523728150643*X^20)/261891630 - 
          (758758896597934*X^22)/1003917915 - (140016275873962409*X^24)/
           30117537450 - (296977606008255923*X^26)/10039179150 - 
          (84557963126865574466*X^28)/436704293025 - 
          (2317307158359742099337*X^30)/1785208758300)*Log[X]^2 - 
        (4*X^2*(-6 - 12*X^2 - 5*X^4 + 65*X^6 + 261*X^8 + 865*X^10 + 
           3125*X^12 + 13161*X^14 + 63469*X^16 + 337589*X^18 + 1924161*X^20 + 
           11541343*X^22 + 72027335*X^24 + 464151087*X^26 + 3071705911*X^28)*
          Log[X]^3)/27 + S2*(81/4 + (243*X^2)/4 + (585*X^4)/4 + 
          (2885*X^6)/12 + (4097*X^8)/12 + (145045*X^10)/324 + 
          (4882415*X^12)/8748 + (88186271*X^14)/131220 + 
          (310580299*X^16)/393660 + (202772737451*X^18)/223205220 + 
          (172470304684*X^20)/167403915 + (745228260362*X^22)/645700815 + 
          (312389104158763*X^24)/244074908070 + (11335546977284381*X^26)/
           8054471966310 + (111375157711161629*X^28)/72490247696790 + 
          (2827584827341465903*X^30)/1696271796104886 - 
          (X^4*(847288609443 + 1914244636149*X^2 + 3148566314103*X^4 + 
             4513061276361*X^6 + 5983192891953*X^8 + 7542086846247*X^10 + 
             9177489172665*X^12 + 10880110754919*X^14 + 12642671323701*X^16 + 
             14459312693331*X^18 + 16325219636073*X^20 + 18236364305751*
              X^22 + 20189327775399*X^24 + 22181171683297*X^26)*Log[X])/
           31381059609) + lM*(1/9 - X^2/9 + X^4/6 + (5*X^6)/162 + 
          (7*X^8)/648 + X^10/200 + (11*X^12)/4050 + (13*X^14)/7938 + 
          (5*X^16)/4704 + (17*X^18)/23328 + (19*X^20)/36450 + 
          (7*X^22)/18150 + (23*X^24)/78408 + (25*X^26)/109512 + 
          (3*X^28)/16562 + (29*X^30)/198450 - (Pi^2*(-1 + X^2))/27 - 
          (X^4*(540540 + 420420*X^2 + 390390*X^4 + 378378*X^6 + 372372*X^8 + 
             368940*X^10 + 366795*X^12 + 365365*X^14 + 364364*X^16 + 
             363636*X^18 + 363090*X^20 + 362670*X^22 + 362340*X^24 + 
             362076*X^26)*Log[X])/810810 - (4*X^2*(3 + 6*X^2 + 10*X^4 + 
             14*X^6 + 18*X^8 + 22*X^10 + 26*X^12 + 30*X^14 + 34*X^16 + 
             38*X^18 + 42*X^20 + 46*X^22 + 50*X^24 + 54*X^26 + 58*X^28)*
            Log[X]^2)/9) - (68*Zeta[3])/9 - (40*X^2*Zeta[3])/3 - 
        (121*X^4*Zeta[3])/6 - (69*X^6*Zeta[3])/2 - (7375*X^8*Zeta[3])/162 - 
        (69625*X^10*Zeta[3])/1458 - (265025*X^12*Zeta[3])/13122 + 
        (5104187*X^14*Zeta[3])/39366 + (100997471*X^16*Zeta[3])/118098 + 
        (130964617843*X^18*Zeta[3])/28697814 + (713151237097*X^20*Zeta[3])/
         28697814 + (36767651927735*X^22*Zeta[3])/258280326 + 
        (5960453171867951*X^24*Zeta[3])/6973568802 + 
        (111615012295820339*X^26*Zeta[3])/20920706406 + 
        (6473535426560925155*X^28*Zeta[3])/188286357654 + 
        (385573119530223220973*X^30*Zeta[3])/1694577218886 + 
        nl*(-1/36 - (5*X^2)/12 - (7*X^4)/72 - (553*X^6)/648 - 
          (3757*X^8)/2592 - (43361*X^10)/21600 - (496307*X^12)/194400 - 
          (21047903*X^14)/6804000 - (7687021*X^16)/2116800 - 
          (8337740303*X^18)/2000376000 - (12546276727*X^20)/2667168000 - 
          (2174001437191*X^22)/414935136000 - (50316679953851*X^24)/
           8713637856000 - (817625476519723*X^26)/129589222194432 - 
          (4106121836154133*X^28)/599950102752000 - 
          (1553833941721187149*X^30)/210582486065952000 + 
          (Pi^2*(1171170 + 630630*X^2 + 180180*X^4 + 300300*X^6 + 
             330330*X^8 + 342342*X^10 + 348348*X^12 + 351780*X^14 + 
             353925*X^16 + 355355*X^18 + 356356*X^20 + 357084*X^22 + 
             357630*X^24 + 358050*X^26 + 358380*X^28 + 358644*X^30))/
           9729720 + (X^2*(-194788994400 - 422042821200*X^2 + 
             119037718800*X^4 + 589777788600*X^6 + 1037792475720*X^8 + 
             1476608793660*X^10 + 1910988251172*X^12 + 2342956035420*X^14 + 
             2773495661508*X^16 + 3203133302655*X^18 + 3632171968425*X^20 + 
             4060796560535*X^22 + 4489125313465*X^24 + 4917236806755*X^26 + 
             5345184984629*X^28 + 21643221600*Pi^2*(6 + 9*X^2 + 11*X^4 + 13*
                X^6 + 15*X^8 + 17*X^10 + 19*X^12 + 21*X^14 + 23*X^16 + 25*
                X^18 + 27*X^20 + 29*X^22 + 31*X^24 + 33*X^26 + 35*X^28))*
            Log[X])/584366983200 - (X^2*(13 + 26*X^2 + 46*X^4 + 66*X^6 + 
             86*X^8 + 106*X^10 + 126*X^12 + 146*X^14 + 166*X^16 + 186*X^18 + 
             206*X^20 + 226*X^22 + 246*X^24 + 266*X^26 + 286*X^28)*Log[X]^2)/
           9 + (4*X^2*(6 + 12*X^2 + 19*X^4 + 26*X^6 + 33*X^8 + 40*X^10 + 
             47*X^12 + 54*X^14 + 61*X^16 + 68*X^18 + 75*X^20 + 82*X^22 + 
             89*X^24 + 96*X^26 + 103*X^28)*Log[X]^3)/27 + 
          lM*(1/18 - X^2/18 + X^4/12 + (5*X^6)/324 + (7*X^8)/1296 + 
            X^10/400 + (11*X^12)/8100 + (13*X^14)/15876 + (5*X^16)/9408 + 
            (17*X^18)/46656 + (19*X^20)/72900 + (7*X^22)/36300 + 
            (23*X^24)/156816 + (25*X^26)/219024 + (3*X^28)/33124 + 
            (29*X^30)/396900 - (Pi^2*(-1 + X^2))/54 - 
            (X^4*(540540 + 420420*X^2 + 390390*X^4 + 378378*X^6 + 372372*
                X^8 + 368940*X^10 + 366795*X^12 + 365365*X^14 + 364364*X^16 + 
               363636*X^18 + 363090*X^20 + 362670*X^22 + 362340*X^24 + 362076*
                X^26)*Log[X])/1621620 - (2*X^2*(3 + 6*X^2 + 10*X^4 + 14*X^6 + 
               18*X^8 + 22*X^10 + 26*X^12 + 30*X^14 + 34*X^16 + 38*X^18 + 42*
                X^20 + 46*X^22 + 50*X^24 + 54*X^26 + 58*X^28)*Log[X]^2)/9) - 
          (2*Zeta[3])/9 + (2*X^2*Zeta[3])/9 + (2*X^4*Zeta[3])/3 + 
          (10*X^6*Zeta[3])/9 + (14*X^8*Zeta[3])/9 + 2*X^10*Zeta[3] + 
          (22*X^12*Zeta[3])/9 + (26*X^14*Zeta[3])/9 + (10*X^16*Zeta[3])/3 + 
          (34*X^18*Zeta[3])/9 + (38*X^20*Zeta[3])/9 + (14*X^22*Zeta[3])/3 + 
          (46*X^24*Zeta[3])/9 + (50*X^26*Zeta[3])/9 + 6*X^28*Zeta[3] + 
          (58*X^30*Zeta[3])/9)) + cf^2*(-119/108 + D3/6 - (835*X^2)/216 + 
        X^4/72 - (139343*X^6)/3888 - (680341*X^8)/11664 - 
        (8014807031*X^10)/174960000 + (6423962228273*X^12)/48009024000 + 
        (271417234280377*X^14)/226842638400 + (1052201012195316549659*X^16)/
         144925224820992000 + (119534489431211306167475531*X^18)/
         2763263382157422144000 + (1105273382915835728889851321*X^20)/
         4144895073236133216000 + (18773152815730239049429463675173063*X^22)/
         10996489527196926144712320000 + 
        (8435354860771275085387880543447976131497*X^24)/
         746706724503732792623159848512000 + 
        (1220761716998978522709690877218208283911427*X^26)/
         15884488503079406679438127686528000 + 
        (965322864103172809937795312411280434500654379417*X^28)/
         1803821335091693316641421195911870976000 + 
        (202272065472359494561565987770592919306068978358317*X^30)/
         53215882919007560696370599092679760000000 + B4*(-1 - X^2) + 
        OepS2*((-7*X^4)/18 - (64*X^6)/81 - (3527*X^8)/2916 - 
          (17948*X^10)/10935 - (27319*X^12)/13122 - (3136064*X^14)/1240029 - 
          (38023127*X^16)/12754584 - (147960640*X^18)/43046721 - 
          (1677317401*X^20)/430467210 - (55722374992*X^22)/12784876137 - 
          (201786670873*X^24)/41841412812 - (6472704469100*X^26)/
           1223861324751 - (7587098764289*X^28)/1318004503578 - 
          (26374766820784*X^30)/4236443047215) + 
        Pi^4*(-17/180 - (4*X^2)/45 + (133*X^4)/8640 + (38*X^6)/1215 + 
          (67013*X^8)/1399680 + (85253*X^10)/1312200 + (519061*X^12)/
           6298560 + (1862038*X^14)/18600435 + (722439413*X^16)/6122200320 + 
          (17570326*X^18)/129140163 + (31869030619*X^20)/206624260800 + 
          (66170320303*X^22)/383546284110 + (3833946746587*X^24)/
           20083878149760 + (6149069245645*X^26)/29372671794024 + 
          (144154876521491*X^28)/632642161717440 + (31320035599681*X^30)/
           127093291416450) + T1ep*((7*X^4)/12 + (32*X^6)/27 + 
          (3527*X^8)/1944 + (8974*X^10)/3645 + (27319*X^12)/8748 + 
          (1568032*X^14)/413343 + (38023127*X^16)/8503056 + 
          (73980320*X^18)/14348907 + (1677317401*X^20)/286978140 + 
          (27861187496*X^22)/4261625379 + (201786670873*X^24)/27894275208 + 
          (3236352234550*X^26)/407953774917 + (7587098764289*X^28)/
           878669669052 + (13187383410392*X^30)/1412147682405) + 
        (8*Log[2]^4)/9 + (8*X^2*Log[2]^4)/9 + 
        Pi^2*(-385/648 + (899*X^4)/576 + (352*X^6)/81 + 
          (2065877*X^8)/466560 - (10040903*X^10)/2449440 - 
          (3997719877*X^12)/73483200 - (228015615371*X^14)/694416240 - 
          (6758516076061*X^16)/3571283520 - (24035462944890073*X^18)/
           2121342410880 - (75290615924171687*X^20)/1060671205440 - 
          (946778624745064602721*X^22)/2047625762101920 - 
          (229202733120776839153319*X^24)/73714527435669120 - 
          (4207088284622530936403237*X^26)/196013629772120160 - 
          (103884373444355350552891253*X^28)/686047704202420560 - 
          (8786516610057672557722283153*X^30)/8074253749459257360 + 
          (2*Log[2])/3 - (8*Log[2]^2)/9 + X^2*(-173/324 + (2*Log[2])/3 - 
            (8*Log[2]^2)/9)) + (X^2/2 + 7*X^4 + (304*X^6)/27 - 
          (1415*X^8)/108 - (252319*X^10)/1800 - (7262117*X^12)/9450 - 
          (323917499*X^14)/79380 - (26631533389*X^16)/1164240 - 
          (5298945106087*X^18)/38918880 - (20724838693759*X^20)/24324300 - 
          (5608765889638661*X^22)/1010809800 - (20474157235680353*X^24)/
           548725320 - (1753767654902256977*X^26)/6809182380 - 
          (13551293099195174111*X^28)/7457675940 - 
          (9176812159367025505147*X^30)/702742540500)*Log[X]^2 + 
        (4*X^2 + 16*X^4 + (716*X^6)/27 + (920*X^8)/27 + (10*X^10)/9 - 
          (10712*X^12)/45 - (41696*X^14)/27 - (551984*X^16)/63 - 
          (458062*X^18)/9 - (25127128*X^20)/81 - (88949044*X^22)/45 - 
          (1293249224*X^24)/99 - (2400961783*X^26)/27 - (72568964408*X^28)/
           117 - (278155562914*X^30)/63)*Log[X]^3 + 
        S2*(-351/8 - (459*X^2)/8 - 36*X^4 - (101*X^6)/6 + (1657*X^8)/576 + 
          (560737*X^10)/24300 + (127285*X^12)/2916 + (622218931*X^14)/
           9644670 + (19410594103*X^16)/226748160 + (567978619*X^18)/
           5314410 + (5154463799657*X^20)/40176939600 + 
          (491679138747457*X^22)/3281451541830 + (191333864404573*X^24)/
           1115771008320 + (1578034643561982247*X^26)/8167234573838340 + 
          (96987482178757741*X^28)/451050430113360 + 
          (3348766661833717351*X^30)/14135598300874050 + 
          ((-63*X^4)/4 - 32*X^6 - (3527*X^8)/72 - (8974*X^10)/135 - 
            (27319*X^12)/324 - (1568032*X^14)/15309 - (38023127*X^16)/
             314928 - (73980320*X^18)/531441 - (1677317401*X^20)/10628820 - 
            (27861187496*X^22)/157837977 - (201786670873*X^24)/1033121304 - 
            (3236352234550*X^26)/15109399071 - (7587098764289*X^28)/
             32543321076 - (13187383410392*X^30)/52301766015)*Log[X]) + 
        (64*PolyLog[4, 1/2])/3 + (64*X^2*PolyLog[4, 1/2])/3 + 
        (253*Zeta[3])/18 + (182*X^2*Zeta[3])/9 + (937*X^4*Zeta[3])/72 + 
        (1667*X^6*Zeta[3])/54 + (209779*X^8*Zeta[3])/3888 + 
        (743861*X^10*Zeta[3])/7290 + (22568503*X^12*Zeta[3])/87480 + 
        (393624079*X^14*Zeta[3])/413343 + (548429378125*X^16*Zeta[3])/
         119042784 + (1476498687095*X^18*Zeta[3])/57395628 + 
        (89214908978837*X^20*Zeta[3])/573956280 + 
        (4213537615357403*X^22*Zeta[3])/4261625379 + 
        (4008531566850023335*X^24*Zeta[3])/613674054576 + 
        (18138760762248460105*X^26*Zeta[3])/407953774917 + 
        (7084919135428245017857*X^28*Zeta[3])/22845411395352 + 
        (21822047517769402905128*X^30*Zeta[3])/9885033776835 + 
        Log[X]*(X^30*(-303048954820040604806044041982567598171/
             148747503625629318950740200000 + Pi^2*(-2086167212119/1890 + 
              (8*Log[2])/3) + 4*Zeta[3]) + 
          X^28*(-42701307565584823137410993196475073233/
             168448284105823779860660457600 + Pi^2*(-507983488225/3276 + 
              (8*Log[2])/3) + 4*Zeta[3]) + 
          X^26*(-1404014143407814629973376858747849/
             45489679747724488215139200 + Pi^2*(-31212774445/1404 + 
              (8*Log[2])/3) + 4*Zeta[3]) + 
          X^24*(-1913720476699280156939209952941/534600937206163857229200 + 
            Pi^2*(-7759884487/2376 + (8*Log[2])/3) + 4*Zeta[3]) + 
          X^20*(-453564191349066071609/19170158883876000 + 
            Pi^2*(-125817493/1620 + (8*Log[2])/3) + 4*Zeta[3]) + 
          X^18*(103397171325807494923/30672254214201600 + 
            Pi^2*(-1383923/108 + (8*Log[2])/3) + 4*Zeta[3]) + 
          X^16*(11205316896509671/4705364442240 + Pi^2*(-2279003/1008 + 
              (8*Log[2])/3) + 4*Zeta[3]) + X^12*(12056282839/28576800 + 
            Pi^2*(-52937/540 + (8*Log[2])/3) + 4*Zeta[3]) + 
          X^14*(1347772460389/1350253800 + (8*Pi^2*(-3460 + 21*Log[2]))/63 + 
            4*Zeta[3]) + X^10*(190732363/972000 + (2*Pi^2*(-583 + 60*Log[2]))/
             45 + 4*Zeta[3]) + X^22*(-604250634615613993136080223/
             1645432247479728708000 + (Pi^2*(-16311839 + 88*Log[2]))/33 + 
            4*Zeta[3]) + X^4*(27/8 + (Pi^2*(11 + 96*Log[2]))/36 + 
            4*Zeta[3]) + (X^2*(5 + 8*Pi^2*(-1 + Log[4]) + 12*Zeta[3]))/6 + 
          (X^6*(4133 + 4*Pi^2*(-59 + 72*Log[2]) + 432*Zeta[3]))/108 + 
          (X^8*(122905 + 6*Pi^2*(-1679 + 576*Log[2]) + 5184*Zeta[3]))/
           1296)) + ca*cf*(-49/144 - D3/12 + (139*X^2)/72 - (1963*X^4)/1152 + 
        (166849*X^6)/10368 + (5111567*X^8)/186624 + (4537847*X^10)/218700 - 
        (8347346586907*X^12)/120022560000 - (10913948420987887*X^14)/
         18147411072000 - (1053266389635314097149*X^16)/289850449641984000 - 
        (239115413391483975052159567*X^18)/11053053528629688576000 - 
        (22106251774616673741252427537*X^20)/165795802929445328640000 - 
        (103252976687230269785307215319329773*X^22)/
         120961384799166187591835520000 - 
        (1687072701766452148409556505474395410333*X^24)/
         298682689801493117049263939404800 - 
        (110978356169892890645527755360666118362937*X^26)/
         2888088818741710305352386852096000 - 
        (120665361104369465170829838016645242311764966369*X^28)/
         450955333772923329160355298977967744000 - 
        (8899970915383290187877291494307757785379543857650573*X^30)/
         4682997696872665341280612720155818880000000 + (B4*(1 + X^2))/2 + 
        T1ep*((-7*X^4)/24 - (16*X^6)/27 - (3527*X^8)/3888 - 
          (4487*X^10)/3645 - (27319*X^12)/17496 - (784016*X^14)/413343 - 
          (38023127*X^16)/17006112 - (36990160*X^18)/14348907 - 
          (1677317401*X^20)/573956280 - (13930593748*X^22)/4261625379 - 
          (201786670873*X^24)/55788550416 - (1618176117275*X^26)/
           407953774917 - (7587098764289*X^28)/1757339338104 - 
          (6593691705196*X^30)/1412147682405) + 
        Pi^4*(67/1080 + (4*X^2)/135 - (133*X^4)/17280 - (19*X^6)/1215 - 
          (67013*X^8)/2799360 - (85253*X^10)/2624400 - (519061*X^12)/
           12597120 - (931019*X^14)/18600435 - (722439413*X^16)/12244400640 - 
          (8785163*X^18)/129140163 - (31869030619*X^20)/413248521600 - 
          (66170320303*X^22)/767092568220 - (3833946746587*X^24)/
           40167756299520 - (6149069245645*X^26)/58745343588048 - 
          (144154876521491*X^28)/1265284323434880 - (31320035599681*X^30)/
           254186582832900) + OepS2*((7*X^4)/36 + (32*X^6)/81 + 
          (3527*X^8)/5832 + (8974*X^10)/10935 + (27319*X^12)/26244 + 
          (1568032*X^14)/1240029 + (38023127*X^16)/25509168 + 
          (73980320*X^18)/43046721 + (1677317401*X^20)/860934420 + 
          (27861187496*X^22)/12784876137 + (201786670873*X^24)/83682825624 + 
          (3236352234550*X^26)/1223861324751 + (7587098764289*X^28)/
           2636009007156 + (13187383410392*X^30)/4236443047215) - 
        (4*Log[2]^4)/9 - (4*X^2*Log[2]^4)/9 + 
        Pi^2*((-2197*X^4)/3456 - (5425*X^6)/2592 - (2032469*X^8)/933120 + 
          (1462001*X^10)/699840 + (4003489507*X^12)/146966400 + 
          (456139092853*X^14)/2777664960 + (6758790850741*X^16)/7142567040 + 
          (437011363411423*X^18)/77139724032 + (75290696502318989*X^20)/
           2121342410880 + (1820728422400692493*X^22)/7875483700392 + 
          (32743248382726845607073*X^24)/21061293553048320 + 
          (4207088299364469421996517*X^26)/392027259544240320 + 
          (415537493983381987564697531*X^28)/5488381633619364480 + 
          (413483134619419040907632123*X^30)/759929764654988928 + 
          (X^2*(67 - 48*Log[2] + 64*Log[2]^2))/144 + 
          (-49 - 72*Log[2] + 96*Log[2]^2)/216) + 
        ((185*X^2)/18 + (185*X^4)/9 + (314*X^6)/9 + (3433*X^8)/54 + 
          (172879*X^10)/1200 + (3591509*X^12)/7560 + (170503663*X^14)/79380 + 
          (26921303089*X^16)/2328480 + (1769979526669*X^18)/25945920 + 
          (41465052178807*X^20)/97297200 + (116856651315158*X^22)/42117075 + 
          (61423102430004761*X^24)/3292351920 + (438442623053456603*X^26)/
           3404591190 + (83001690798284406838*X^28)/91356530265 + 
          (2294203124832729056623*X^30)/351371270250)*Log[X]^2 + 
        ((-44*X^2)/9 - (88*X^4)/9 - (274*X^6)/27 - (34*X^8)/9 + 
          (413*X^10)/15 + (7507*X^12)/45 + (159619*X^14)/189 + 
          (564821*X^16)/126 + (4144631*X^18)/162 + (62888903*X^20)/405 + 
          (97865695*X^22)/99 + (1939953640*X^24)/297 + (10404243221*X^26)/
           234 + (84663895943*X^28)/273 + (695389047506*X^30)/315)*Log[X]^3 + 
        S2*(351/16 + (459*X^2)/16 + 18*X^4 + (101*X^6)/12 - (1657*X^8)/1152 - 
          (560737*X^10)/48600 - (127285*X^12)/5832 - (622218931*X^14)/
           19289340 - (19410594103*X^16)/453496320 - (567978619*X^18)/
           10628820 - (5154463799657*X^20)/80353879200 - 
          (491679138747457*X^22)/6562903083660 - (191333864404573*X^24)/
           2231542016640 - (1578034643561982247*X^26)/16334469147676680 - 
          (96987482178757741*X^28)/902100860226720 - 
          (3348766661833717351*X^30)/28271196601748100 + 
          ((63*X^4)/8 + 16*X^6 + (3527*X^8)/144 + (4487*X^10)/135 + 
            (27319*X^12)/648 + (784016*X^14)/15309 + (38023127*X^16)/629856 + 
            (36990160*X^18)/531441 + (1677317401*X^20)/21257640 + 
            (13930593748*X^22)/157837977 + (201786670873*X^24)/2066242608 + 
            (1618176117275*X^26)/15109399071 + (7587098764289*X^28)/
             65086642152 + (6593691705196*X^30)/52301766015)*Log[X]) + 
        lM*(-11/36 + (11*X^2)/36 - (11*X^4)/24 - (55*X^6)/648 - 
          (77*X^8)/2592 - (11*X^10)/800 - (121*X^12)/16200 - 
          (143*X^14)/31752 - (55*X^16)/18816 - (187*X^18)/93312 - 
          (209*X^20)/145800 - (7*X^22)/6600 - (23*X^24)/28512 - 
          (275*X^26)/438048 - (33*X^28)/66248 - (319*X^30)/793800 + 
          (11*Pi^2*(-1 + X^2))/108 + (X^4*(540540 + 420420*X^2 + 390390*X^4 + 
             378378*X^6 + 372372*X^8 + 368940*X^10 + 366795*X^12 + 
             365365*X^14 + 364364*X^16 + 363636*X^18 + 363090*X^20 + 
             362670*X^22 + 362340*X^24 + 362076*X^26)*Log[X])/294840 + 
          (11*X^2*(3 + 6*X^2 + 10*X^4 + 14*X^6 + 18*X^8 + 22*X^10 + 26*X^12 + 
             30*X^14 + 34*X^16 + 38*X^18 + 42*X^20 + 46*X^22 + 50*X^24 + 
             54*X^26 + 58*X^28)*Log[X]^2)/9) - (32*PolyLog[4, 1/2])/3 - 
        (32*X^2*PolyLog[4, 1/2])/3 - (52*Zeta[3])/9 - (347*X^2*Zeta[3])/36 - 
        (853*X^4*Zeta[3])/144 - (383*X^6*Zeta[3])/27 - (197683*X^8*Zeta[3])/
         7776 - (714701*X^10*Zeta[3])/14580 - (22140823*X^12*Zeta[3])/
         174960 - (391235875*X^14*Zeta[3])/826686 - 
        (547635759565*X^16*Zeta[3])/238085568 - (1476065031239*X^18*Zeta[3])/
         114791256 - (89210062236917*X^20*Zeta[3])/1147912560 - 
        (4213497840187199*X^22*Zeta[3])/8523250758 - 
        (4008525293737465447*X^24*Zeta[3])/1227348109152 - 
        (18138756229428738805*X^26*Zeta[3])/815907549834 - 
        (7084918861283308273633*X^28*Zeta[3])/45690822790704 - 
        (10911023695181150446294*X^30*Zeta[3])/9885033776835 + 
        Log[X]*(X^10*(-177106039/1944000 + Pi^2*(129/10 - (4*Log[2])/3) - 
            2*Zeta[3]) + X^12*(-29057911639/142884000 + 
            Pi^2*(52957/1080 - (4*Log[2])/3) - 2*Zeta[3]) + 
          X^14*(-757471724041/1543147200 + Pi^2*(83075/378 - (4*Log[2])/3) - 
            2*Zeta[3]) + X^16*(-347550446136101/294085277640 + 
            Pi^2*(2279339/2016 - (4*Log[2])/3) - 2*Zeta[3]) + 
          X^18*(-102809028781607175529/61344508428403200 + 
            Pi^2*(153775/24 - (4*Log[2])/3) - 2*Zeta[3]) + 
          X^20*(64851254162967019469/5477188252536000 + 
            Pi^2*(125818513/3240 - (4*Log[2])/3) - 2*Zeta[3]) + 
          X^22*(604286843978648785318851941/3290864494959457416000 + 
            Pi^2*(24467797/99 - (4*Log[2])/3) - 2*Zeta[3]) + 
          X^24*(239216625840804742925494676213/133650234301540964307300 + 
            Pi^2*(7759886687/4752 - (4*Log[2])/3) - 2*Zeta[3]) + 
          X^26*(127637752277227474819848507958339/8270850863222634220934400 + 
            Pi^2*(1156028739/104 - (4*Log[2])/3) - 2*Zeta[3]) + 
          X^28*(3050093714230682729616858669596093291/
             24064040586546254265808636800 + Pi^2*(169327830743/2184 - 
              (4*Log[2])/3) - 2*Zeta[3]) + 
          X^30*(29698797977140922419508659839744744665983/
             29154510710623346514345079200000 + Pi^2*(2086167214709/3780 - 
              (4*Log[2])/3) - 2*Zeta[3]) - 
          (X^6*(5453 + 4*Pi^2*(-91 + 144*Log[2]) + 864*Zeta[3]))/432 - 
          (X^8*(35177 + 6*Pi^2*(-541 + 192*Log[2]) + 1728*Zeta[3]))/864 - 
          (X^2*(8*Pi^2*(-1 + Log[8]) + 3*(-7 + 6*Zeta[3])))/18 - 
          (X^4*(Pi^2*(29 + 96*Log[2]) + 4*(-79 + 36*Zeta[3])))/72)))
 

rhoOnShell2Exp1 = nc*(cf^2*((-13*Pi^4)/90 - (13*Pi^4*(-1 + X))/90 - 
        (365*(-1 + X)^2)/72 - (5*Pi^2*(-1 + X)^2)/9 - (13*Pi^4*(-1 + X)^2)/
         180 - (26347*(-1 + X)^4)/2160 + (Pi^2*(-1 + X)^4)/36 + 
        (26347*(-1 + X)^5)/2160 - (Pi^2*(-1 + X)^5)/36 - 
        (502124321*(-1 + X)^6)/21772800 + (11*Pi^2*(-1 + X)^6)/504 + 
        (369335441*(-1 + X)^7)/10886400 - (Pi^2*(-1 + X)^7)/63 - 
        (469560837901*(-1 + X)^8)/9754214400 + (17*Pi^2*(-1 + X)^8)/1512 + 
        (640860544103*(-1 + X)^9)/9754214400 - (Pi^2*(-1 + X)^9)/126 + 
        B4*(-1 - X^2) + (8*Pi^2*(-1 + X)^2*Log[2])/9 - 
        (2*Pi^2*(-1 + X)^4*Log[2])/45 + (2*Pi^2*(-1 + X)^5*Log[2])/45 - 
        (11*Pi^2*(-1 + X)^6*Log[2])/315 + (8*Pi^2*(-1 + X)^7*Log[2])/315 - 
        (17*Pi^2*(-1 + X)^8*Log[2])/945 + (4*Pi^2*(-1 + X)^9*Log[2])/315 - 
        (4*Pi^2*Log[2]^2)/3 - (4*Pi^2*(-1 + X)*Log[2]^2)/3 - 
        (2*Pi^2*(-1 + X)^2*Log[2]^2)/3 + (4*Log[2]^4)/3 + 
        (4*(-1 + X)*Log[2]^4)/3 + (2*(-1 + X)^2*Log[2]^4)/3 + 
        32*PolyLog[4, 1/2] + 32*(-1 + X)*PolyLog[4, 1/2] + 
        16*(-1 + X)^2*PolyLog[4, 1/2] - 
        ((-1 + X)^20*(5888169433143304996682348083092514 + 
           1050078817023126348496896000*Pi^2*(-5 + Log[256]) - 
           4898384811776205846771355012921875*Zeta[3]))/
         10230338821072421818120273920000 + 
        ((-1 + X)^19*(9680681715220846533929619729718 + 
           2361808982158052268441600*Pi^2*(-5 + Log[256]) - 
           8053392081191646308292810061425*Zeta[3]))/
         19399636309309159020429312000 - 
        ((-1 + X)^18*(9523248801163625344205941865366 + 
           3237503091570489950208000*Pi^2*(-5 + Log[256]) - 
           7922439742155280661251029722625*Zeta[3]))/
         22171012924924753166204928000 - 
        ((-1 + X)^16*(144997935385059933320077418102 + 
           102116454470894223360000*Pi^2*(-5 + Log[256]) - 
           120626270949655059201928940625*Zeta[3]))/
         467139505212122852229120000 + 
        ((-1 + X)^17*(73429970943259673297782557026 + 35512451763724615680000*
            Pi^2*(-5 + Log[256]) - 61087062092428701904081381875*Zeta[3]))/
         200202645090909793812480000 + 
        ((-1 + X)^15*(9662203069589382378116482 + 10173274178715648000*Pi^2*
            (-5 + Log[256]) - 8038301920617611900131875*Zeta[3]))/
         37155821787587543040000 - ((-1 + X)^14*(640235393637704822397686 + 
           1037948919585177600*Pi^2*(-5 + Log[256]) - 
           532651101285935102295825*Zeta[3]))/2972465743007003443200 - 
        ((-1 + X)^12*(657938937809299404782 + 2798289027072000*Pi^2*
            (-5 + Log[256]) - 547475627747822437125*Zeta[3]))/
         4639837885562880000 + ((-1 + X)^13*(163402632771817033118 + 
           421460095795200*Pi^2*(-5 + Log[256]) - 135953209748160067725*
            Zeta[3]))/927967577112576000 - 
        ((-1 + X)^10*(11194461929540894 + 145539072000*Pi^2*(-5 + Log[256]) - 
           9320077078066125*Zeta[3]))/128755630080000 + 
        ((-1 + X)^11*(3611575184645894 + 26320896000*Pi^2*(-5 + Log[256]) - 
           3005813717699625*Zeta[3]))/32188907520000 + 
        (229*(-1 + X)^2*Zeta[3])/48 + (14869*(-1 + X)^4*Zeta[3])/1440 - 
        (14869*(-1 + X)^5*Zeta[3])/1440 + (2079031*(-1 + X)^6*Zeta[3])/
         107520 - (914353*(-1 + X)^7*Zeta[3])/32256 + 
        (13154953*(-1 + X)^8*Zeta[3])/327680 - (125544701*(-1 + X)^9*Zeta[3])/
         2293760) + ca*cf*((13*Pi^4)/180 + lM*((-11*(-1 + X)^2)/9 - 
          (121*(-1 + X)^4)/1620 + (121*(-1 + X)^5)/1620 - 
          (209*(-1 + X)^6)/4200 + (352*(-1 + X)^7)/14175 - 
          (11209*(-1 + X)^8)/1587600 - (209*(-1 + X)^9)/58800 + 
          (36559*(-1 + X)^10)/4082400 - (5657*(-1 + X)^11)/510300 + 
          (21521*(-1 + X)^12)/1891890 - (7315409*(-1 + X)^13)/681080400 + 
          (15645143*(-1 + X)^14)/1609826400 - (3845*(-1 + X)^15)/447174 + 
          (7310077*(-1 + X)^16)/971611200 - (2365553*(-1 + X)^17)/360884160 + 
          (575157263*(-1 + X)^18)/100735689600 - (56486071*(-1 + X)^19)/
           11332765080 + (36505713779*(-1 + X)^20)/8359580829600) + 
        (13*Pi^4*(-1 + X))/180 - (475*(-1 + X)^2)/432 + 
        (4*Pi^2*(-1 + X)^2)/27 + (13*Pi^4*(-1 + X)^2)/360 + 
        (11*(-1 + X)^3)/9 - (8191*(-1 + X)^4)/3240 - (Pi^2*(-1 + X)^4)/135 + 
        (2591*(-1 + X)^5)/1080 + (Pi^2*(-1 + X)^5)/135 - 
        (36852517*(-1 + X)^6)/14515200 - (11*Pi^2*(-1 + X)^6)/1890 + 
        (58396207*(-1 + X)^7)/21772800 + (4*Pi^2*(-1 + X)^7)/945 - 
        (1949049527609*(-1 + X)^8)/682795008000 - (17*Pi^2*(-1 + X)^8)/5670 + 
        (77416253041*(-1 + X)^9)/25288704000 + (2*Pi^2*(-1 + X)^9)/945 - 
        (2430966149861069*(-1 + X)^10)/737418608640000 - 
        (47*Pi^2*(-1 + X)^10)/31185 + (654828276680489*(-1 + X)^11)/
         184354652160000 + (34*Pi^2*(-1 + X)^11)/31185 - 
        (14996426883951776303*(-1 + X)^12)/3926016672399360000 - 
        (326*Pi^2*(-1 + X)^12)/405405 + (41786027378128344347*(-1 + X)^13)/
         10207643348238336000 + (491*Pi^2*(-1 + X)^13)/810810 - 
        (844068292023640928736197*(-1 + X)^14)/193210273295455223808000 - 
        (151*Pi^2*(-1 + X)^14)/324324 + (2241672813672368462713907*
          (-1 + X)^15)/483025683238638059520000 + (148*Pi^2*(-1 + X)^15)/
         405405 - (20628924146767834082596697837*(-1 + X)^16)/
         4204255546909105670062080000 - (1607*Pi^2*(-1 + X)^16)/5513508 + 
        (3100431420968415483477643757*(-1 + X)^17)/
         600607935272729381437440000 + (326*Pi^2*(-1 + X)^17)/1378377 - 
        (14254554221409886440141133104803*(-1 + X)^18)/
         2638350538066045626778386432000 - (5099*Pi^2*(-1 + X)^18)/26189163 + 
        (618141811785520774803098615783*(-1 + X)^19)/
         109931272419418567782432768000 + (21256*Pi^2*(-1 + X)^19)/
         130945815 - (1678343513966954775746582656878221*(-1 + X)^20)/
         288619558861164385232423485440000 - (17921*Pi^2*(-1 + X)^20)/
         130945815 + (B4*(1 + X^2))/2 - (4*Pi^2*(-1 + X)^2*Log[2])/9 + 
        (Pi^2*(-1 + X)^4*Log[2])/45 - (Pi^2*(-1 + X)^5*Log[2])/45 + 
        (11*Pi^2*(-1 + X)^6*Log[2])/630 - (4*Pi^2*(-1 + X)^7*Log[2])/315 + 
        (17*Pi^2*(-1 + X)^8*Log[2])/1890 - (2*Pi^2*(-1 + X)^9*Log[2])/315 + 
        (47*Pi^2*(-1 + X)^10*Log[2])/10395 - (34*Pi^2*(-1 + X)^11*Log[2])/
         10395 + (326*Pi^2*(-1 + X)^12*Log[2])/135135 - 
        (491*Pi^2*(-1 + X)^13*Log[2])/270270 + (151*Pi^2*(-1 + X)^14*Log[2])/
         108108 - (148*Pi^2*(-1 + X)^15*Log[2])/135135 + 
        (1607*Pi^2*(-1 + X)^16*Log[2])/1837836 - 
        (326*Pi^2*(-1 + X)^17*Log[2])/459459 + (5099*Pi^2*(-1 + X)^18*Log[2])/
         8729721 - (21256*Pi^2*(-1 + X)^19*Log[2])/43648605 + 
        (17921*Pi^2*(-1 + X)^20*Log[2])/43648605 + (2*Pi^2*Log[2]^2)/3 + 
        (2*Pi^2*(-1 + X)*Log[2]^2)/3 + (Pi^2*(-1 + X)^2*Log[2]^2)/3 - 
        (2*Log[2]^4)/3 - (2*(-1 + X)*Log[2]^4)/3 - ((-1 + X)^2*Log[2]^4)/3 - 
        16*PolyLog[4, 1/2] - 16*(-1 + X)*PolyLog[4, 1/2] - 
        8*(-1 + X)^2*PolyLog[4, 1/2] + (155*(-1 + X)^2*Zeta[3])/96 + 
        (1337*(-1 + X)^4*Zeta[3])/960 - (1337*(-1 + X)^5*Zeta[3])/960 + 
        (1067099*(-1 + X)^6*Zeta[3])/645120 - (617867*(-1 + X)^7*Zeta[3])/
         322560 + (89725001*(-1 + X)^8*Zeta[3])/41287680 - 
        (100208411*(-1 + X)^9*Zeta[3])/41287680 + 
        (29184380059*(-1 + X)^10*Zeta[3])/10899947520 - 
        (7969509151*(-1 + X)^11*Zeta[3])/2724986880 + 
        (598609105021*(-1 + X)^12*Zeta[3])/188932423680 - 
        (643958218201*(-1 + X)^13*Zeta[3])/188932423680 + 
        (734354788247*(-1 + X)^14*Zeta[3])/201527918592 - 
        (390339643549*(-1 + X)^15*Zeta[3])/100763959296 + 
        (7698594913313131*(-1 + X)^16*Zeta[3])/1879391789383680 - 
        (56690038719832699*(-1 + X)^17*Zeta[3])/13155742525685760 + 
        (2704979535900269951*(-1 + X)^18*Zeta[3])/599901859171270656 - 
        (65149589276001959*(-1 + X)^19*Zeta[3])/13886617110446080 + 
        (1939899600092708999*(-1 + X)^20*Zeta[3])/399934572780847104) + 
      cf*(lM*((4*(-1 + X)^2)/9 + (11*(-1 + X)^4)/405 - (11*(-1 + X)^5)/405 + 
          (19*(-1 + X)^6)/1050 - (128*(-1 + X)^7)/14175 + 
          (1019*(-1 + X)^8)/396900 + (19*(-1 + X)^9)/14700 - 
          (36559*(-1 + X)^10)/11226600 + (5657*(-1 + X)^11)/1403325 - 
          (43042*(-1 + X)^12)/10405395 + (7315409*(-1 + X)^13)/1872971100 - 
          (15645143*(-1 + X)^14)/4427022600 + (7690*(-1 + X)^15)/2459457 - 
          (7310077*(-1 + X)^16)/2671930800 + (2365553*(-1 + X)^17)/
           992431440 - (575157263*(-1 + X)^18)/277023146400 + 
          (56486071*(-1 + X)^19)/31165103970 - (36505713779*(-1 + X)^20)/
           22988847281400) + nl*(lM*((2*(-1 + X)^2)/9 + (11*(-1 + X)^4)/810 - 
            (11*(-1 + X)^5)/810 + (19*(-1 + X)^6)/2100 - (64*(-1 + X)^7)/
             14175 + (1019*(-1 + X)^8)/793800 + (19*(-1 + X)^9)/29400 - 
            (36559*(-1 + X)^10)/22453200 + (5657*(-1 + X)^11)/2806650 - 
            (21521*(-1 + X)^12)/10405395 + (7315409*(-1 + X)^13)/3745942200 - 
            (15645143*(-1 + X)^14)/8854045200 + (3845*(-1 + X)^15)/2459457 - 
            (7310077*(-1 + X)^16)/5343861600 + (2365553*(-1 + X)^17)/
             1984862880 - (575157263*(-1 + X)^18)/554046292800 + 
            (56486071*(-1 + X)^19)/62330207940 - (36505713779*(-1 + X)^20)/
             45977694562800) - (19*(-1 + X)^2)/54 + (2*Pi^2*(-1 + X)^2)/27 - 
          (2*(-1 + X)^3)/9 + (473*(-1 + X)^4)/3240 - (Pi^2*(-1 + X)^4)/270 - 
          (397*(-1 + X)^5)/3240 + (Pi^2*(-1 + X)^5)/270 + 
          (1697*(-1 + X)^6)/18000 - (11*Pi^2*(-1 + X)^6)/3780 - 
          (28337*(-1 + X)^7)/425250 + (2*Pi^2*(-1 + X)^7)/945 + 
          (459269*(-1 + X)^8)/10418625 - (17*Pi^2*(-1 + X)^8)/11340 - 
          (1530079*(-1 + X)^9)/55566000 + (Pi^2*(-1 + X)^9)/945 + 
          (228978179*(-1 + X)^10)/14145516000 - (47*Pi^2*(-1 + X)^10)/62370 - 
          (30965651*(-1 + X)^11)/3536379000 + (17*Pi^2*(-1 + X)^11)/31185 + 
          (46860660169*(-1 + X)^12)/11537501976000 - (163*Pi^2*(-1 + X)^12)/
           405405 - (13705010779*(-1 + X)^13)/11537501976000 + 
          (491*Pi^2*(-1 + X)^13)/1621620 - (181990589119*(-1 + X)^14)/
           354515969808000 - (151*Pi^2*(-1 + X)^14)/648648 + 
          (935479733*(-1 + X)^15)/633064231800 + (74*Pi^2*(-1 + X)^15)/
           405405 - (3110398765219*(-1 + X)^16)/1564642597518000 - 
          (1607*Pi^2*(-1 + X)^16)/11027016 + (361413161401*(-1 + X)^17)/
           162560269872000 + (163*Pi^2*(-1 + X)^17)/1378377 - 
          (8090581442816351*(-1 + X)^18)/3524698462757472000 - 
          (5099*Pi^2*(-1 + X)^18)/52378326 + (1562350718792519*(-1 + X)^19)/
           687316200237707040 + (10628*Pi^2*(-1 + X)^19)/130945815 - 
          (3529383873768458831*(-1 + X)^20)/1605489783025843915200 - 
          (17921*Pi^2*(-1 + X)^20)/261891630) - 
        ((-1 + X)^2*(-5708138177915703969985553248706 + 11396161962875289600*
            Pi^2*(51310044002 - 712212463163*X + 4845876627983*X^2 - 
             21467714647963*X^3 + 68739680256418*X^4 - 167810084748115*X^5 + 
             322170220231012*X^6 - 495700941042245*X^7 + 618155049669620*
              X^8 - 628302464336413*X^9 + 521032080603162*X^10 - 
             351307041698935*X^11 + 190998316632298*X^12 - 82540068673935*
              X^13 + 27713469350228*X^14 - 6972648526153*X^15 + 
             1237432685694*X^16 - 138186680422*X^17 + 7305936767*X^18) - 
           29329280717360601515821408125*Zeta[3] + 2320500*X^12*
            (-9284876576624485287398562074 + 21975185659658255647325703*
              Zeta[3]) + 46512*X^3*(51855869507367859137241511954 + 
             61244332492010070579301125*Zeta[3]) + 189618*X^8*
            (-367372057803075655609274351806 + 553841560570242003119416125*
              Zeta[3]) + 3*X^18*(-274948226954499408678787226974 + 
             846663768239383320795694125*Zeta[3]) - 
           1292*X^4*(5983237588631646978860484090094 + 
             1791605814510686347589290875*Zeta[3]) - 
           18*X^17*(-866647514728684828080725604538 + 
             2589311142438826358418240375*Zeta[3]) + 
           153*X^16*(-912908074687227007882198369282 + 
             2635937932104078943478011875*Zeta[3]) - 
           425*X^2*(1279233063257340665183887435322 + 
             3287048074302674112371966025*Zeta[3]) - 6664*X^13*
            (-1397450616051460924167517527638 + 3520250362582573547405102625*
              Zeta[3]) + 6188*X^6*(-5862340904979013913044976010686 + 
             4746853410560273391822568125*Zeta[3]) + 24310*X^10*
            (-2416676319186126528976116374594 + 4829795635014096641044625475*
              Zeta[3]) - 10608*X^7*(-5263987672686752149029186460502 + 
             6278527723148215732995422625*Zeta[3]) - 
           952*X^5*(-19837114849984755206526625113466 + 
             6536064515079789213865716375*Zeta[3]) - 
           272*X^15*(-2893134932498202412977097828982 + 
             8034318800827139495238188625*Zeta[3]) - 
           9724*X^9*(-7283614993606111395585532268354 + 
             12908167842588263492732359875*Zeta[3]) - 
           3536*X^11*(-11205062346405922777695320895602 + 
             24592131981307896018006868875*Zeta[3]) + 
           68*X^14*(-45989504410980675492116012783426 + 
             122099047976067723113894215875*Zeta[3]) + 
           X*(79726106207833950481284967047572 + 
             413770010728119647887900865250*Zeta[3])))/
         23876475457611272640528384000))


CbarDbMsBar = cf*nc*((lm^2*(-1 - x^2))/4 + (1 + x^2)/(9*ep^3) + 
      (lm^3*(1 + x^2))/2 + (lm*(-1 + Pi^2)*(1 + x^2))/12 + 
      (lm*x^2 - 3*lm^2*x^2 + (2 - (1 + Pi^2)*x^2)/6)*HPL[{0}, x] + 
      ((-1 - x^2)/18 + (lm*(1 + x^2))/3 - (2*x^2*HPL[{0}, x])/3)/ep^2 + 
      (12*lm*x^2 + (-1 - 7*x^2)/3)*HPL[{0, 0}, x] + 
      ((lm*(-1 - x^2))/6 + (lm^2*(1 + x^2))/2 + ((-1 + Pi^2)*(1 + x^2))/36 + 
        (x^2/3 - 2*lm*x^2)*HPL[{0}, x] + 4*x^2*HPL[{0, 0}, x])/ep + 
      ((1 - 8*x - 6*x^2 - 8*x^3 + x^4)*HPL[{-1, 0, 0}, x])/(6*x) - 
      (64*x^2*HPL[{0, 0, 0}, x])/3 + ((1 + 8*x - 6*x^2 + 8*x^3 + x^4)*
        HPL[{1, 0, 0}, x])/(6*x) - 
      ((1 + x^2)*(49 + Pi^2 - 4*PolyGamma[2, 1] - 46*Zeta[3]))/72)
 
CbarSingMsBar = cf*nc*(2*x*HPL[{-1, 0, 0}, x] + 2*x*HPL[{1, 0, 0}, x] - 
      (7*(1 + x^2)*Zeta[3])/4)
 
CbarDiagMsBar = nc*((-4*(1 + x^2))/(3*ep) - (4*lm*(1 + x^2))/3 + 
       (8*x^2*HPL[{0}, x])/3) + api*cf*nc*((-1 - x^2)/12 + lm^2*(-1 - x^2) + 
       (1 + x^2)/ep^2 - (5*(1 + x^2))/(6*ep) - (2*x^2*HPL[{0}, x])/3 + 
       lm*((1 + x^2)/3 + 4*x^2*HPL[{0}, x]) - 8*x^2*HPL[{0, 0}, x]) + 
     api^2*nc*(cf*((-7*lm^3*(1 + x^2))/18 + (nl*(1 + x^2))/(18*ep^3) + 
         ((-8 + 3*Pi^2)*x^2*HPL[{0}, x])/18 + ((lm*(-1 - x^2))/3 - 
           (2*(1 + x^2))/9 - (5*nl*(1 + x^2))/36 + (2*x^2*HPL[{0}, x])/3)/
          ep^2 + lm^2*((17*(1 + x^2))/36 + (7*x^2*HPL[{0}, x])/3) + 
         (34*x^2*HPL[{0, 0}, x])/9 + lm*(-((-8 + 3*Pi^2)*(1 + x^2))/36 - 
           (17*x^2*HPL[{0}, x])/9 - (28*x^2*HPL[{0, 0}, x])/3) + 
         ((lm^2*(-1 - x^2))/2 + (nl*(1 + x^2))/9 - ((-9 + Pi^2)*(1 + x^2))/
            36 - (x^2*HPL[{0}, x])/3 + lm*((1 + x^2)/6 + 2*x^2*HPL[{0}, x]) - 
           4*x^2*HPL[{0, 0}, x])/ep + (56*x^2*HPL[{0, 0, 0}, x])/3 + 
         ((1 + x^2)*(-1 + Pi^2 + 8*Zeta[3]))/72 + 
         nl*((lm^3*(1 + x^2))/18 - (5*x^2*HPL[{0}, x])/36 + 
           lm^2*((1 + x^2)/9 - (x^2*HPL[{0}, x])/3) + (8*x^2*HPL[{0, 0}, x])/
            9 + lm*((5*(1 + x^2))/72 - (4*x^2*HPL[{0}, x])/9 + 
             (4*x^2*HPL[{0, 0}, x])/3) - (8*x^2*HPL[{0, 0, 0}, x])/3 + 
           ((1 + x^2)*(-1 + 32*Zeta[3]))/72)) + 
       cf^2*((-1 - x^2)/(2*ep^3) + (lm^3*(-1 - x^2))/2 + (1 + x^2)/(2*ep^2) + 
         lm^2*((13*(1 + x^2))/8 + 3*x^2*HPL[{0}, x]) + 
         13*x^2*HPL[{0, 0}, x] + 24*x^2*HPL[{0, 0, 0}, x] + 
         ((1 + x^2)*(85 - 160*B4 + 8*Pi^4 - 960*Zeta[3]))/160 - 
         ((1 + x^2)*(-119 + 144*Zeta[3]))/(144*ep) + 
         (x^2*HPL[{0}, x]*(-101 + 144*Zeta[3]))/24 + 
         lm*((-13*x^2*HPL[{0}, x])/2 - 12*x^2*HPL[{0, 0}, x] - 
           ((1 + x^2)*(-101 + 144*Zeta[3]))/48)) + 
       ca*cf*((-11*(1 + x^2))/(36*ep^3) + (83*(1 + x^2))/(72*ep^2) - 
         (11*lm^3*(1 + x^2))/36 + lm^2*((-43*(1 + x^2))/36 + 
           (11*x^2*HPL[{0}, x])/6) - (86*x^2*HPL[{0, 0}, x])/9 + 
         (44*x^2*HPL[{0, 0, 0}, x])/3 + x^2*HPL[{0}, x]*(85/72 - 3*Zeta[3]) + 
         ((1 + x^2)*(-77 + 36*Zeta[3]))/(72*ep) - 
         ((1 + x^2)*(-180*B4 + 9*Pi^4 - 20*(9 + 28*Zeta[3])))/360 + 
         lm*((43*x^2*HPL[{0}, x])/9 - (22*x^2*HPL[{0, 0}, x])/3 + 
           ((1 + x^2)*(-85 + 216*Zeta[3]))/144)))
