Zetaa3lmubmuc = 1 + api*((-11*CA*lcb)/6 - (2*lb1*nb*TF)/3 + 
       ((2*l21)/3 - (2*lb1)/3)*nc*TF + (2*lcb*nl*TF)/3) + 
     api^2*((-17*CA^2*lcb)/12 + (121*CA^2*lcb^2)/36 + 
       (((5*CA)/6 + CF/2)*lcb - (22*CA*lcb^2)/9)*nl*TF + 
       (4*lb1^2*nb^2*TF^2)/9 + ((-8*l21*lb1)/9 + (8*lb1^2)/9)*nb*nc*TF^2 + 
       ((4*l21^2)/9 - (8*l21*lb1)/9 + (4*lb1^2)/9)*nc^2*TF^2 + 
       (4*lcb^2*nl^2*TF^2)/9 + 
       nc*(((2*CA)/9 - (13*CF)/48 + ((5*CA)/6 - CF/2)*l21 + 
           ((-5*CA)/6 + CF/2)*lb1 + ((-22*CA*l21)/9 + (22*CA*lb1)/9)*lcb)*
          TF + ((8*l21)/9 - (8*lb1)/9)*lcb*nl*TF^2) + 
       nb*(((2*CA)/9 - (13*CF)/48 + ((-5*CA)/6 + CF/2)*lb1 + 
           (22*CA*lb1*lcb)/9)*TF - (8*lb1*lcb*nl*TF^2)/9)) + 
     api^3*((-2857*CA^3*lcb)/1728 + (935*CA^3*lcb^2)/144 - 
       (1331*CA^3*lcb^3)/216 + (((1415*CA^2)/864 + (205*CA*CF)/288 - CF^2/16)*
          lcb + ((-445*CA^2)/72 - (55*CA*CF)/24)*lcb^2 + (121*CA^2*lcb^3)/18)*
        nl*TF + (((-79*CA)/432 - (11*CF)/72)*lcb + ((25*CA)/18 + (5*CF)/6)*
          lcb^2 - (22*CA*lcb^3)/9)*nl^2*TF^2 - (8*lb1^3*nb^3*TF^3)/27 + 
       ((8*l21^3)/27 - (8*l21^2*lb1)/9 + (8*l21*lb1^2)/9 - (8*lb1^3)/27)*nc^3*
        TF^3 + (8*lcb^3*nl^3*TF^3)/27 + 
       nb*nc*(((-16*l21*lb1)/9 + (16*lb1^2)/9)*lcb*nl*TF^3 + 
         TF^2*((83*CA)/2592 + (7*CF)/2592 + ((-113*CA)/216 + (5*CF)/6 + 
             ((-25*CA)/9 + (5*CF)/3)*l21)*lb1 + ((25*CA)/9 - (5*CF)/3)*
            lb1^2 + ((44*CA*l21*lb1)/9 - (44*CA*lb1^2)/9)*lcb + 
           (CA/32 + (5*CF)/64)/x^2 + (CA/32 + (5*CF)/64)*x^2 + 
           l21^2*(CA/48 + (7*CF)/192 + (CA/64 + (5*CF)/128)/x^2 + 
             (CA/64 + (5*CF)/128)*x^2) + l21*((113*CA)/432 - (5*CF)/12 + 
             (-CA/32 - (5*CF)/64)/x^2 + (CA/32 + (5*CF)/64)*x^2) + 
           Lmx*((-CA/64 - (5*CF)/128)/x^3 + (-CA/64 - (3*CF)/128)/x + 
             (-CA/64 - (3*CF)/128)*x + (-CA/64 - (5*CF)/128)*x^3))) + 
       nb^2*(((8*l21*lb1^2)/9 - (8*lb1^3)/9)*nc*TF^3 + (8*lb1^2*lcb*nl*TF^3)/
          9 + TF^2*((245*CA)/5184 + (103*CF)/1296 + 
           ((-113*CA)/432 + (5*CF)/12)*lb1 + ((25*CA)/18 - (5*CF)/6)*lb1^2 - 
           (22*CA*lb1^2*lcb)/9 - (7*CA*z3)/128 - (7*CF*z3)/64)) + 
       nc^2*(((8*l21^2)/9 - (16*l21*lb1)/9 + (8*lb1^2)/9)*lcb*nl*TF^3 + 
         TF^2*((245*CA)/5184 + (103*CF)/1296 + ((113*CA)/432 - (5*CF)/12)*
            l21 + ((25*CA)/18 - (5*CF)/6)*l21^2 + 
           ((-113*CA)/432 + (5*CF)/12 + ((-25*CA)/9 + (5*CF)/3)*l21)*lb1 + 
           ((25*CA)/18 - (5*CF)/6)*lb1^2 + ((-22*CA*l21^2)/9 + 
             (44*CA*l21*lb1)/9 - (22*CA*lb1^2)/9)*lcb - (7*CA*z3)/128 - 
           (7*CF*z3)/64)) + nb*((-CA/2592 - (41*CF)/162 + 
           ((47*CA)/216 + (5*CF)/24)*lb1 - (CF*lb1^2)/3 + 
           ((8*CA)/27 - (13*CF)/36 + ((-25*CA)/9 - CF/3)*lb1)*lcb + 
           (44*CA*lb1*lcb^2)/9)*nl*TF^2 + ((-8*l21^2*lb1)/9 + 
           (16*l21*lb1^2)/9 - (8*lb1^3)/9)*nc^2*TF^3 - 
         (8*lb1*lcb^2*nl^2*TF^3)/9 + TF*((11347*CA^2)/20736 + 
           (2999*CA*CF)/2592 - (97*CF^2)/288 + ((-1063*CA^2)/864 + 
             (25*CA*CF)/18 - (9*CF^2)/16)*lb1 + 
           ((-7*CA^2)/24 + (11*CA*CF)/24)*lb1^2 + 
           ((-22*CA^2)/27 + (143*CA*CF)/144 + ((53*CA^2)/9 - (11*CA*CF)/6)*
              lb1)*lcb - (121*CA^2*lb1*lcb^2)/18 - (5*CA^2*z3)/1536 - 
           (1273*CA*CF*z3)/768 + (95*CF^2*z3)/192)) + 
       nc*((-CA/2592 - (41*CF)/162 + ((-47*CA)/216 - (5*CF)/24)*l21 - 
           (CF*l21^2)/3 + ((47*CA)/216 + (5*CF)/24 + (2*CF*l21)/3)*lb1 - 
           (CF*lb1^2)/3 + ((8*CA)/27 - (13*CF)/36 + ((25*CA)/9 + CF/3)*l21 + 
             ((-25*CA)/9 - CF/3)*lb1)*lcb + ((-44*CA*l21)/9 + (44*CA*lb1)/9)*
            lcb^2)*nl*TF^2 + ((8*l21)/9 - (8*lb1)/9)*lcb^2*nl^2*TF^3 + 
         TF*((11347*CA^2)/20736 + (2999*CA*CF)/2592 - (97*CF^2)/288 + 
           ((1063*CA^2)/864 - (25*CA*CF)/18 + (9*CF^2)/16)*l21 + 
           ((-7*CA^2)/24 + (11*CA*CF)/24)*l21^2 + 
           ((-1063*CA^2)/864 + (25*CA*CF)/18 - (9*CF^2)/16 + 
             ((7*CA^2)/12 - (11*CA*CF)/12)*l21)*lb1 + 
           ((-7*CA^2)/24 + (11*CA*CF)/24)*lb1^2 + 
           ((-22*CA^2)/27 + (143*CA*CF)/144 + ((-53*CA^2)/9 + (11*CA*CF)/6)*
              l21 + ((53*CA^2)/9 - (11*CA*CF)/6)*lb1)*lcb + 
           ((121*CA^2*l21)/18 - (121*CA^2*lb1)/18)*lcb^2 - (5*CA^2*z3)/1536 - 
           (1273*CA*CF*z3)/768 + (95*CF^2*z3)/192)))
ZetaaInv3lmubmuc = 1 + anlpi*((11*CA*lcb)/6 + (2*lb1*nb*TF)/3 + 
       ((-2*l21)/3 + (2*lb1)/3)*nc*TF - (2*lcb*nl*TF)/3) + 
     anlpi^2*((17*CA^2*lcb)/12 + (121*CA^2*lcb^2)/36 + 
       (((-5*CA)/6 - CF/2)*lcb - (22*CA*lcb^2)/9)*nl*TF + 
       (4*lb1^2*nb^2*TF^2)/9 + ((-8*l21*lb1)/9 + (8*lb1^2)/9)*nb*nc*TF^2 + 
       ((4*l21^2)/9 - (8*l21*lb1)/9 + (4*lb1^2)/9)*nc^2*TF^2 + 
       (4*lcb^2*nl^2*TF^2)/9 + 
       nc*(((-2*CA)/9 + (13*CF)/48 + ((-5*CA)/6 + CF/2)*l21 + 
           ((5*CA)/6 - CF/2)*lb1 + ((-22*CA*l21)/9 + (22*CA*lb1)/9)*lcb)*TF + 
         ((8*l21)/9 - (8*lb1)/9)*lcb*nl*TF^2) + 
       nb*(((-2*CA)/9 + (13*CF)/48 + ((5*CA)/6 - CF/2)*lb1 + 
           (22*CA*lb1*lcb)/9)*TF - (8*lb1*lcb*nl*TF^2)/9)) + 
     anlpi^3*((2857*CA^3*lcb)/1728 + (935*CA^3*lcb^2)/144 + 
       (1331*CA^3*lcb^3)/216 + (((-1415*CA^2)/864 - (205*CA*CF)/288 + 
           CF^2/16)*lcb + ((-445*CA^2)/72 - (55*CA*CF)/24)*lcb^2 - 
         (121*CA^2*lcb^3)/18)*nl*TF + (((79*CA)/432 + (11*CF)/72)*lcb + 
         ((25*CA)/18 + (5*CF)/6)*lcb^2 + (22*CA*lcb^3)/9)*nl^2*TF^2 + 
       (8*lb1^3*nb^3*TF^3)/27 + ((-8*l21^3)/27 + (8*l21^2*lb1)/9 - 
         (8*l21*lb1^2)/9 + (8*lb1^3)/27)*nc^3*TF^3 - (8*lcb^3*nl^3*TF^3)/27 + 
       nb*nc*(((16*l21*lb1)/9 - (16*lb1^2)/9)*lcb*nl*TF^3 + 
         TF^2*((-83*CA)/2592 - (7*CF)/2592 + ((-23*CA)/24 + (35*CF)/36 + 
             ((-25*CA)/9 + (5*CF)/3)*l21)*lb1 + ((25*CA)/9 - (5*CF)/3)*
            lb1^2 + ((-44*CA*l21*lb1)/9 + (44*CA*lb1^2)/9)*lcb + 
           (-CA/32 - (5*CF)/64)/x^2 + (-CA/32 - (5*CF)/64)*x^2 + 
           l21*((23*CA)/48 - (35*CF)/72 + (CA/32 + (5*CF)/64)/x^2 + 
             (-CA/32 - (5*CF)/64)*x^2) + l21^2*(-CA/48 - (7*CF)/192 + 
             (-CA/64 - (5*CF)/128)/x^2 + (-CA/64 - (5*CF)/128)*x^2) + 
           Lmx*((CA/64 + (5*CF)/128)/x^3 + (CA/64 + (3*CF)/128)/x + 
             (CA/64 + (3*CF)/128)*x + (CA/64 + (5*CF)/128)*x^3))) + 
       nb^2*(((-8*l21*lb1^2)/9 + (8*lb1^3)/9)*nc*TF^3 - (8*lb1^2*lcb*nl*TF^3)/
          9 + TF^2*((-245*CA)/5184 - (103*CF)/1296 + 
           ((-23*CA)/48 + (35*CF)/72)*lb1 + ((25*CA)/18 - (5*CF)/6)*lb1^2 + 
           (22*CA*lb1^2*lcb)/9 + (7*CA*z3)/128 + (7*CF*z3)/64)) + 
       nc^2*(((-8*l21^2)/9 + (16*l21*lb1)/9 - (8*lb1^2)/9)*lcb*nl*TF^3 + 
         TF^2*((-245*CA)/5184 - (103*CF)/1296 + ((23*CA)/48 - (35*CF)/72)*
            l21 + ((25*CA)/18 - (5*CF)/6)*l21^2 + ((-23*CA)/48 + (35*CF)/72 + 
             ((-25*CA)/9 + (5*CF)/3)*l21)*lb1 + ((25*CA)/18 - (5*CF)/6)*
            lb1^2 + ((22*CA*l21^2)/9 - (44*CA*l21*lb1)/9 + (22*CA*lb1^2)/9)*
            lcb + (7*CA*z3)/128 + (7*CF*z3)/64)) + 
       nb*((CA/2592 + (41*CF)/162 + ((-47*CA)/216 - (5*CF)/24)*lb1 + 
           (CF*lb1^2)/3 + ((4*CA)/9 - (13*CF)/24 + ((-25*CA)/9 + CF/3)*lb1)*
            lcb - (44*CA*lb1*lcb^2)/9)*nl*TF^2 + 
         ((8*l21^2*lb1)/9 - (16*l21*lb1^2)/9 + (8*lb1^3)/9)*nc^2*TF^3 + 
         (8*lb1*lcb^2*nl^2*TF^3)/9 + TF*((-11347*CA^2)/20736 - 
           (2999*CA*CF)/2592 + (97*CF^2)/288 + ((1063*CA^2)/864 - 
             (25*CA*CF)/18 + (9*CF^2)/16)*lb1 + ((7*CA^2)/24 - (11*CA*CF)/24)*
            lb1^2 + ((-11*CA^2)/9 + (143*CA*CF)/96 + 
             ((233*CA^2)/36 - (11*CA*CF)/4)*lb1)*lcb + (121*CA^2*lb1*lcb^2)/
            18 + (5*CA^2*z3)/1536 + (1273*CA*CF*z3)/768 - 
           (95*CF^2*z3)/192)) + 
       nc*((CA/2592 + (41*CF)/162 + ((47*CA)/216 + (5*CF)/24)*l21 + 
           (CF*l21^2)/3 + ((-47*CA)/216 - (5*CF)/24 - (2*CF*l21)/3)*lb1 + 
           (CF*lb1^2)/3 + ((4*CA)/9 - (13*CF)/24 + ((25*CA)/9 - CF/3)*l21 + 
             ((-25*CA)/9 + CF/3)*lb1)*lcb + ((44*CA*l21)/9 - (44*CA*lb1)/9)*
            lcb^2)*nl*TF^2 + ((-8*l21)/9 + (8*lb1)/9)*lcb^2*nl^2*TF^3 + 
         TF*((-11347*CA^2)/20736 - (2999*CA*CF)/2592 + (97*CF^2)/288 + 
           ((-1063*CA^2)/864 + (25*CA*CF)/18 - (9*CF^2)/16)*l21 + 
           ((7*CA^2)/24 - (11*CA*CF)/24)*l21^2 + ((1063*CA^2)/864 - 
             (25*CA*CF)/18 + (9*CF^2)/16 + ((-7*CA^2)/12 + (11*CA*CF)/12)*
              l21)*lb1 + ((7*CA^2)/24 - (11*CA*CF)/24)*lb1^2 + 
           ((-11*CA^2)/9 + (143*CA*CF)/96 + ((-233*CA^2)/36 + (11*CA*CF)/4)*
              l21 + ((233*CA^2)/36 - (11*CA*CF)/4)*lb1)*lcb + 
           ((-121*CA^2*l21)/18 + (121*CA^2*lb1)/18)*lcb^2 + 
           (5*CA^2*z3)/1536 + (1273*CA*CF*z3)/768 - (95*CF^2*z3)/192)))
ZetaaOS3lmubmuc = 1 + api*((-11*CA*lcb)/6 - (2*lb1*nb*TF)/3 + 
       ((2*l21)/3 - (2*lb1)/3)*nc*TF + (2*lcb*nl*TF)/3) + 
     api^2*((-17*CA^2*lcb)/12 + (121*CA^2*lcb^2)/36 + 
       (((5*CA)/6 + CF/2)*lcb - (22*CA*lcb^2)/9)*nl*TF + 
       (4*lb1^2*nb^2*TF^2)/9 + ((-8*l21*lb1)/9 + (8*lb1^2)/9)*nb*nc*TF^2 + 
       ((4*l21^2)/9 - (8*l21*lb1)/9 + (4*lb1^2)/9)*nc^2*TF^2 + 
       (4*lcb^2*nl^2*TF^2)/9 + 
       nc*(((2*CA)/9 - (15*CF)/16 + ((5*CA)/6 + CF/2)*l21 + 
           ((-5*CA)/6 - CF/2)*lb1 + ((-22*CA*l21)/9 + (22*CA*lb1)/9)*lcb)*
          TF + ((8*l21)/9 - (8*lb1)/9)*lcb*nl*TF^2) + 
       nb*(((2*CA)/9 - (15*CF)/16 + ((-5*CA)/6 - CF/2)*lb1 + 
           (22*CA*lb1*lcb)/9)*TF - (8*lb1*lcb*nl*TF^2)/9)) + 
     api^3*((-2857*CA^3*lcb)/1728 + (935*CA^3*lcb^2)/144 - 
       (1331*CA^3*lcb^3)/216 + (((1415*CA^2)/864 + (205*CA*CF)/288 - CF^2/16)*
          lcb + ((-445*CA^2)/72 - (55*CA*CF)/24)*lcb^2 + (121*CA^2*lcb^3)/18)*
        nl*TF + (((-79*CA)/432 - (11*CF)/72)*lcb + ((25*CA)/18 + (5*CF)/6)*
          lcb^2 - (22*CA*lcb^3)/9)*nl^2*TF^2 - (8*lb1^3*nb^3*TF^3)/27 + 
       ((8*l21^3)/27 - (8*l21^2*lb1)/9 + (8*l21*lb1^2)/9 - (8*lb1^3)/27)*nc^3*
        TF^3 + (8*lcb^3*nl^3*TF^3)/27 + 
       nb*nc*(((-16*l21*lb1)/9 + (16*lb1^2)/9)*lcb*nl*TF^3 + 
         TF^2*((83*CA)/2592 + (2563*CF)/2592 - (CF*Hm10xp1)/3 + 
           (CF*Hp10xp1)/3 + ((-113*CA)/216 + (73*CF)/18 + 
             ((-25*CA)/9 - (5*CF)/3)*l21)*lb1 + ((25*CA)/9 + (5*CF)/3)*
            lb1^2 + ((44*CA*l21*lb1)/9 - (44*CA*lb1^2)/9)*lcb + 
           Hm10xm1*(-CF/3 - CF/(3*x^4) - CF/(3*x^3) - CF/(3*x)) + 
           Hp10xm1*(CF/3 + CF/(3*x^4) - CF/(3*x^3) - CF/(3*x)) + 
           (CA/32 + (37*CF)/64)/x^2 + (CA/32 + (37*CF)/64)*x^2 + 
           l21*((113*CA)/432 - (73*CF)/36 + (-CA/32 - (79*CF)/192)/x^2 + 
             (CA/32 + (79*CF)/192)*x^2) + Lmx*((-CA/64 - (5*CF)/128)/x^3 + 
             (-CA/64 - (3*CF)/128)/x + (-CA/64 - (3*CF)/128)*x + 
             (-CA/64 - (5*CF)/128)*x^3) + l21^2*(CA/48 + (71*CF)/192 + 
             CF/(3*x^4) + (CA/64 + (5*CF)/128)/x^2 + (CA/64 + (5*CF)/128)*
              x^2 + (CF*x^4)/3) + (2*CF*z2)/3 + (CF*z2)/(3*x^4) - 
           (CF*z2)/x^3 - (CF*z2)/x + x*(-(CF*Hm10xp1)/3 - (CF*Hp10xp1)/3 - 
             CF*z2) + x^3*(-(CF*Hm10xp1)/3 - (CF*Hp10xp1)/3 - CF*z2) + 
           x^4*(-(CF*Hm10xp1)/3 + (CF*Hp10xp1)/3 + (CF*z2)/3))) + 
       nb^2*(((8*l21*lb1^2)/9 - (8*lb1^3)/9)*nc*TF^3 + (8*lb1^2*lcb*nl*TF^3)/
          9 + TF^2*((245*CA)/5184 + (695*CF)/648 + 
           ((-113*CA)/432 + (73*CF)/36)*lb1 + ((25*CA)/18 + (5*CF)/6)*lb1^2 - 
           (22*CA*lb1^2*lcb)/9 - (2*CF*z2)/3 - (7*CA*z3)/128 - 
           (7*CF*z3)/64)) + nc^2*(((8*l21^2)/9 - (16*l21*lb1)/9 + 
           (8*lb1^2)/9)*lcb*nl*TF^3 + TF^2*((245*CA)/5184 + (695*CF)/648 + 
           ((113*CA)/432 - (73*CF)/36)*l21 + ((25*CA)/18 + (5*CF)/6)*l21^2 + 
           ((-113*CA)/432 + (73*CF)/36 + ((-25*CA)/9 - (5*CF)/3)*l21)*lb1 + 
           ((25*CA)/18 + (5*CF)/6)*lb1^2 + ((-22*CA*l21^2)/9 + 
             (44*CA*l21*lb1)/9 - (22*CA*lb1^2)/9)*lcb - (2*CF*z2)/3 - 
           (7*CA*z3)/128 - (7*CF*z3)/64)) + 
       nb*(((-8*l21^2*lb1)/9 + (16*l21*lb1^2)/9 - (8*lb1^3)/9)*nc^2*TF^3 - 
         (8*lb1*lcb^2*nl^2*TF^3)/9 + nl*TF^2*(-CA/2592 + (311*CF)/1296 + 
           ((47*CA)/216 + (67*CF)/72)*lb1 + ((8*CA)/27 - (5*CF)/4 + 
             ((-25*CA)/9 - (5*CF)/3)*lb1)*lcb + (44*CA*lb1*lcb^2)/9 + 
           (CF*z2)/3) + TF*((11347*CA^2)/20736 - (8321*CA*CF)/5184 - 
           (77*CF^2)/576 + ((-1063*CA^2)/864 - (175*CA*CF)/72 + CF^2/16)*
            lb1 + ((-7*CA^2)/24 - (11*CA*CF)/24)*lb1^2 + 
           ((-22*CA^2)/27 + (55*CA*CF)/16 + ((53*CA^2)/9 + (11*CA*CF)/6)*lb1)*
            lcb - (121*CA^2*lb1*lcb^2)/18 + (CA*CF*z2)/3 - (5*CF^2*z2)/4 - 
           CA*CF*log2*z2 + 2*CF^2*log2*z2 - (5*CA^2*z3)/1536 - 
           (1081*CA*CF*z3)/768 - (CF^2*z3)/192)) + 
       nc*(((8*l21)/9 - (8*lb1)/9)*lcb^2*nl^2*TF^3 + 
         nl*TF^2*(-CA/2592 + (311*CF)/1296 + ((-47*CA)/216 - (67*CF)/72)*
            l21 + ((47*CA)/216 + (67*CF)/72)*lb1 + 
           ((8*CA)/27 - (5*CF)/4 + ((25*CA)/9 + (5*CF)/3)*l21 + 
             ((-25*CA)/9 - (5*CF)/3)*lb1)*lcb + ((-44*CA*l21)/9 + 
             (44*CA*lb1)/9)*lcb^2 + (CF*z2)/3) + 
         TF*((11347*CA^2)/20736 - (8321*CA*CF)/5184 - (77*CF^2)/576 + 
           ((1063*CA^2)/864 + (175*CA*CF)/72 - CF^2/16)*l21 + 
           ((-7*CA^2)/24 - (11*CA*CF)/24)*l21^2 + 
           ((-1063*CA^2)/864 - (175*CA*CF)/72 + CF^2/16 + 
             ((7*CA^2)/12 + (11*CA*CF)/12)*l21)*lb1 + 
           ((-7*CA^2)/24 - (11*CA*CF)/24)*lb1^2 + 
           ((-22*CA^2)/27 + (55*CA*CF)/16 + ((-53*CA^2)/9 - (11*CA*CF)/6)*
              l21 + ((53*CA^2)/9 + (11*CA*CF)/6)*lb1)*lcb + 
           ((121*CA^2*l21)/18 - (121*CA^2*lb1)/18)*lcb^2 + (CA*CF*z2)/3 - 
           (5*CF^2*z2)/4 - CA*CF*log2*z2 + 2*CF^2*log2*z2 - 
           (5*CA^2*z3)/1536 - (1081*CA*CF*z3)/768 - (CF^2*z3)/192)))
ZetaaInvOS3lmubmuc = 1 + anlpi*((11*CA*lcb)/6 + (2*lb1*nb*TF)/3 + 
       ((-2*l21)/3 + (2*lb1)/3)*nc*TF - (2*lcb*nl*TF)/3) + 
     anlpi^2*((17*CA^2*lcb)/12 + (121*CA^2*lcb^2)/36 + 
       (((-5*CA)/6 - CF/2)*lcb - (22*CA*lcb^2)/9)*nl*TF + 
       (4*lb1^2*nb^2*TF^2)/9 + ((-8*l21*lb1)/9 + (8*lb1^2)/9)*nb*nc*TF^2 + 
       ((4*l21^2)/9 - (8*l21*lb1)/9 + (4*lb1^2)/9)*nc^2*TF^2 + 
       (4*lcb^2*nl^2*TF^2)/9 + 
       nc*(((-2*CA)/9 + (15*CF)/16 + ((-5*CA)/6 - CF/2)*l21 + 
           ((5*CA)/6 + CF/2)*lb1 + ((-22*CA*l21)/9 + (22*CA*lb1)/9)*lcb)*TF + 
         ((8*l21)/9 - (8*lb1)/9)*lcb*nl*TF^2) + 
       nb*(((-2*CA)/9 + (15*CF)/16 + ((5*CA)/6 + CF/2)*lb1 + 
           (22*CA*lb1*lcb)/9)*TF - (8*lb1*lcb*nl*TF^2)/9)) + 
     anlpi^3*((2857*CA^3*lcb)/1728 + (935*CA^3*lcb^2)/144 + 
       (1331*CA^3*lcb^3)/216 + (((-1415*CA^2)/864 - (205*CA*CF)/288 + 
           CF^2/16)*lcb + ((-445*CA^2)/72 - (55*CA*CF)/24)*lcb^2 - 
         (121*CA^2*lcb^3)/18)*nl*TF + (((79*CA)/432 + (11*CF)/72)*lcb + 
         ((25*CA)/18 + (5*CF)/6)*lcb^2 + (22*CA*lcb^3)/9)*nl^2*TF^2 + 
       (8*lb1^3*nb^3*TF^3)/27 + ((-8*l21^3)/27 + (8*l21^2*lb1)/9 - 
         (8*l21*lb1^2)/9 + (8*lb1^3)/27)*nc^3*TF^3 - (8*lcb^3*nl^3*TF^3)/27 + 
       nb*nc*(((16*l21*lb1)/9 - (16*lb1^2)/9)*lcb*nl*TF^3 + 
         TF^2*((-83*CA)/2592 - (2563*CF)/2592 + (CF*Hm10xp1)/3 - 
           (CF*Hp10xp1)/3 + ((-23*CA)/24 + (79*CF)/36 + 
             ((-25*CA)/9 - (5*CF)/3)*l21)*lb1 + ((25*CA)/9 + (5*CF)/3)*
            lb1^2 + ((-44*CA*l21*lb1)/9 + (44*CA*lb1^2)/9)*lcb + 
           Hp10xm1*(-CF/3 - CF/(3*x^4) + CF/(3*x^3) + CF/(3*x)) + 
           Hm10xm1*(CF/3 + CF/(3*x^4) + CF/(3*x^3) + CF/(3*x)) + 
           (-CA/32 - (37*CF)/64)/x^2 + (-CA/32 - (37*CF)/64)*x^2 + 
           l21*((23*CA)/48 - (79*CF)/72 + (CA/32 + (79*CF)/192)/x^2 + 
             (-CA/32 - (79*CF)/192)*x^2) + Lmx*((CA/64 + (5*CF)/128)/x^3 + 
             (CA/64 + (3*CF)/128)/x + (CA/64 + (3*CF)/128)*x + 
             (CA/64 + (5*CF)/128)*x^3) + l21^2*(-CA/48 - (71*CF)/192 - 
             CF/(3*x^4) + (-CA/64 - (5*CF)/128)/x^2 + (-CA/64 - (5*CF)/128)*
              x^2 - (CF*x^4)/3) - (2*CF*z2)/3 - (CF*z2)/(3*x^4) + 
           (CF*z2)/x^3 + (CF*z2)/x + x^4*((CF*Hm10xp1)/3 - (CF*Hp10xp1)/3 - 
             (CF*z2)/3) + x*((CF*Hm10xp1)/3 + (CF*Hp10xp1)/3 + CF*z2) + 
           x^3*((CF*Hm10xp1)/3 + (CF*Hp10xp1)/3 + CF*z2))) + 
       nb^2*(((-8*l21*lb1^2)/9 + (8*lb1^3)/9)*nc*TF^3 - (8*lb1^2*lcb*nl*TF^3)/
          9 + TF^2*((-245*CA)/5184 - (695*CF)/648 + 
           ((-23*CA)/48 + (79*CF)/72)*lb1 + ((25*CA)/18 + (5*CF)/6)*lb1^2 + 
           (22*CA*lb1^2*lcb)/9 + (2*CF*z2)/3 + (7*CA*z3)/128 + 
           (7*CF*z3)/64)) + nc^2*(((-8*l21^2)/9 + (16*l21*lb1)/9 - 
           (8*lb1^2)/9)*lcb*nl*TF^3 + TF^2*((-245*CA)/5184 - (695*CF)/648 + 
           ((23*CA)/48 - (79*CF)/72)*l21 + ((25*CA)/18 + (5*CF)/6)*l21^2 + 
           ((-23*CA)/48 + (79*CF)/72 + ((-25*CA)/9 - (5*CF)/3)*l21)*lb1 + 
           ((25*CA)/18 + (5*CF)/6)*lb1^2 + ((22*CA*l21^2)/9 - 
             (44*CA*l21*lb1)/9 + (22*CA*lb1^2)/9)*lcb + (2*CF*z2)/3 + 
           (7*CA*z3)/128 + (7*CF*z3)/64)) + 
       nb*(((8*l21^2*lb1)/9 - (16*l21*lb1^2)/9 + (8*lb1^3)/9)*nc^2*TF^3 + 
         (8*lb1*lcb^2*nl^2*TF^3)/9 + nl*TF^2*(CA/2592 - (311*CF)/1296 + 
           ((-47*CA)/216 - (67*CF)/72)*lb1 + ((4*CA)/9 - (15*CF)/8 + 
             ((-25*CA)/9 - (5*CF)/3)*lb1)*lcb - (44*CA*lb1*lcb^2)/9 - 
           (CF*z2)/3) + TF*((-11347*CA^2)/20736 + (8321*CA*CF)/5184 + 
           (77*CF^2)/576 + ((1063*CA^2)/864 + (175*CA*CF)/72 - CF^2/16)*lb1 + 
           ((7*CA^2)/24 + (11*CA*CF)/24)*lb1^2 + 
           ((-11*CA^2)/9 + (165*CA*CF)/32 + ((233*CA^2)/36 + (11*CA*CF)/4)*
              lb1)*lcb + (121*CA^2*lb1*lcb^2)/18 - (CA*CF*z2)/3 + 
           (5*CF^2*z2)/4 + CA*CF*log2*z2 - 2*CF^2*log2*z2 + 
           (5*CA^2*z3)/1536 + (1081*CA*CF*z3)/768 + (CF^2*z3)/192)) + 
       nc*(((-8*l21)/9 + (8*lb1)/9)*lcb^2*nl^2*TF^3 + 
         nl*TF^2*(CA/2592 - (311*CF)/1296 + ((47*CA)/216 + (67*CF)/72)*l21 + 
           ((-47*CA)/216 - (67*CF)/72)*lb1 + ((4*CA)/9 - (15*CF)/8 + 
             ((25*CA)/9 + (5*CF)/3)*l21 + ((-25*CA)/9 - (5*CF)/3)*lb1)*lcb + 
           ((44*CA*l21)/9 - (44*CA*lb1)/9)*lcb^2 - (CF*z2)/3) + 
         TF*((-11347*CA^2)/20736 + (8321*CA*CF)/5184 + (77*CF^2)/576 + 
           ((-1063*CA^2)/864 - (175*CA*CF)/72 + CF^2/16)*l21 + 
           ((7*CA^2)/24 + (11*CA*CF)/24)*l21^2 + ((1063*CA^2)/864 + 
             (175*CA*CF)/72 - CF^2/16 + ((-7*CA^2)/12 - (11*CA*CF)/12)*l21)*
            lb1 + ((7*CA^2)/24 + (11*CA*CF)/24)*lb1^2 + 
           ((-11*CA^2)/9 + (165*CA*CF)/32 + ((-233*CA^2)/36 - (11*CA*CF)/4)*
              l21 + ((233*CA^2)/36 + (11*CA*CF)/4)*lb1)*lcb + 
           ((-121*CA^2*l21)/18 + (121*CA^2*lb1)/18)*lcb^2 - (CA*CF*z2)/3 + 
           (5*CF^2*z2)/4 + CA*CF*log2*z2 - 2*CF^2*log2*z2 + 
           (5*CA^2*z3)/1536 + (1081*CA*CF*z3)/768 + (CF^2*z3)/192)))
Zetam3lmubmuc = 1 - (3*api*CF*lcb)/2 + 
     api^2*(((-97*CA*CF)/48 - (3*CF^2)/16)*lcb + ((11*CA*CF)/8 + (9*CF^2)/8)*
        lcb^2 + ((89*CF)/288 - (5*CF*lb1)/12 + (CF*lb1^2)/2 + CF*lb1*lcb)*nb*
        TF + ((89*CF)/288 + (5*CF*l21)/12 + (CF*l21^2)/2 + 
         ((-5*CF)/12 - CF*l21)*lb1 + (CF*lb1^2)/2 + (-(CF*l21) + CF*lb1)*lcb)*
        nc*TF + ((5*CF*lcb)/12 - (CF*lcb^2)/2)*nl*TF) + 
     api^3*(((-11413*CA^2*CF)/3456 + (129*CA*CF^2)/128 - (129*CF^3)/64)*lcb + 
       ((1373*CA^2*CF)/288 + (27*CA*CF^2)/8 + (9*CF^3)/32)*lcb^2 + 
       ((-121*CA^2*CF)/72 - (33*CA*CF^2)/16 - (9*CF^3)/16)*lcb^3 + 
       ((35*CF*lcb)/216 + (5*CF*lcb^2)/18 - (2*CF*lcb^3)/9)*nl^2*TF^2 + 
       nb*nc*TF^2*((-875*CF)/3888 + (CF*l21^3)/9 + 
         ((-31*CF)/27 - (5*CF*l21)/9 - (2*CF*l21^2)/3)*lb1 + 
         ((5*CF)/9 + (4*CF*l21)/3)*lb1^2 - (8*CF*lb1^3)/9 + 
         ((4*CF*l21*lb1)/3 - (4*CF*lb1^2)/3)*lcb - (2*CF*Lpx)/3 - 
         (5*CF)/(48*x^2) - (5*CF*x^2)/48 + l21*((31*CF)/54 + 
           (5*CF)/(48*x^2) - (5*CF*x^2)/48) + l21^2*((-43*CF)/144 - 
           (5*CF)/(96*x^2) - (5*CF*x^2)/96) + Lmx*((5*CF)/(96*x^3) + 
           (9*CF)/(32*x) + (9*CF*x)/32 + (5*CF*x^3)/96) - (2*CF*z3)/9) + 
       nb^2*TF^2*((-1685*CF)/7776 - (31*CF*lb1)/54 + (5*CF*lb1^2)/18 - 
         (4*CF*lb1^3)/9 - (2*CF*lb1^2*lcb)/3 + (7*CF*z3)/18) + 
       nc^2*TF^2*((-1685*CF)/7776 + (31*CF*l21)/54 + (5*CF*l21^2)/18 + 
         (4*CF*l21^3)/9 + ((-31*CF)/54 - (5*CF*l21)/9 - (4*CF*l21^2)/3)*lb1 + 
         ((5*CF)/18 + (4*CF*l21)/3)*lb1^2 - (4*CF*lb1^3)/9 + 
         ((-2*CF*l21^2)/3 + (4*CF*l21*lb1)/3 - (2*CF*lb1^2)/3)*lcb + 
         (7*CF*z3)/18) + nl*TF*(((-197*CA*CF)/72 - (9*CF^2)/8)*lcb^2 + 
         ((11*CA*CF)/9 + (3*CF^2)/4)*lcb^3 + lcb*((139*CA*CF)/216 + 
           (23*CF^2)/16 + (3*CA*CF*z3)/2 - (3*CF^2*z3)/2)) + 
       nb*(nl*TF^2*((1327*CF)/3888 - (53*CF*lb1)/72 - (2*CF*lb1^3)/9 - 
           (5*CF*lb1*lcb)/9 + (2*CF*lb1*lcb^2)/3 - (2*CF*z3)/9) + 
         TF*((16627*CA*CF)/15552 - (B4*CA*CF)/8 - (683*CF^2)/576 + 
           (B4*CF^2)/4 + ((29*CA*CF)/24 - CF^2)*lb1^2 + (11*CA*CF*lb1^3)/18 + 
           (-(CA*CF)/3 - (11*CF^2)/192 + ((71*CA*CF)/18 + CF^2/8)*lb1 - 
             (3*CF^2*lb1^2)/4)*lcb + ((-11*CA*CF)/6 - (3*CF^2)/2)*lb1*lcb^2 - 
           (629*CA*CF*z3)/576 + (57*CF^2*z3)/32 + 
           lb1*((5*CA*CF)/32 - (13*CF^2)/32 - (3*CA*CF*z3)/2 + 
             (3*CF^2*z3)/2) + (9*CA*CF*z4)/8 - (9*CF^2*z4)/8)) + 
       nc*(nl*TF^2*((1327*CF)/3888 + (53*CF*l21)/72 + (2*CF*l21^3)/9 + 
           ((-53*CF)/72 - (2*CF*l21^2)/3)*lb1 + (2*CF*l21*lb1^2)/3 - 
           (2*CF*lb1^3)/9 + ((5*CF*l21)/9 - (5*CF*lb1)/9)*lcb + 
           ((-2*CF*l21)/3 + (2*CF*lb1)/3)*lcb^2 - (2*CF*z3)/9) + 
         TF*((16627*CA*CF)/15552 - (B4*CA*CF)/8 - (683*CF^2)/576 + 
           (B4*CF^2)/4 + ((29*CA*CF)/24 - CF^2)*l21^2 - (11*CA*CF*l21^3)/18 + 
           ((29*CA*CF)/24 - CF^2 - (11*CA*CF*l21)/6)*lb1^2 + 
           (11*CA*CF*lb1^3)/18 + (-(CA*CF)/3 - (11*CF^2)/192 + 
             ((-71*CA*CF)/18 - CF^2/8)*l21 - (3*CF^2*l21^2)/4 + 
             ((71*CA*CF)/18 + CF^2/8 + (3*CF^2*l21)/2)*lb1 - 
             (3*CF^2*lb1^2)/4)*lcb + (((11*CA*CF)/6 + (3*CF^2)/2)*l21 + 
             ((-11*CA*CF)/6 - (3*CF^2)/2)*lb1)*lcb^2 - (629*CA*CF*z3)/576 + 
           (57*CF^2*z3)/32 + l21*((-5*CA*CF)/32 + (13*CF^2)/32 + 
             (3*CA*CF*z3)/2 - (3*CF^2*z3)/2) + lb1*((5*CA*CF)/32 - 
             (13*CF^2)/32 + ((-29*CA*CF)/12 + 2*CF^2)*l21 + (11*CA*CF*l21^2)/
              6 - (3*CA*CF*z3)/2 + (3*CF^2*z3)/2) + (9*CA*CF*z4)/8 - 
           (9*CF^2*z4)/8)))
ZetamInv3lmubmuc = 1 + (3*anlpi*CF*lcb)/2 + 
     anlpi^2*(((97*CA*CF)/48 + (3*CF^2)/16)*lcb + ((11*CA*CF)/8 + (9*CF^2)/8)*
        lcb^2 + ((-89*CF)/288 + (5*CF*lb1)/12 - (CF*lb1^2)/2)*nb*TF + 
       ((-89*CF)/288 - (5*CF*l21)/12 - (CF*l21^2)/2 + ((5*CF)/12 + CF*l21)*
          lb1 - (CF*lb1^2)/2)*nc*TF + ((-5*CF*lcb)/12 - (CF*lcb^2)/2)*nl*
        TF) + anlpi^3*(((11413*CA^2*CF)/3456 - (129*CA*CF^2)/128 + 
         (129*CF^3)/64)*lcb + ((1373*CA^2*CF)/288 + (27*CA*CF^2)/8 + 
         (9*CF^3)/32)*lcb^2 + ((121*CA^2*CF)/72 + (33*CA*CF^2)/16 + 
         (9*CF^3)/16)*lcb^3 + ((-35*CF*lcb)/216 + (5*CF*lcb^2)/18 + 
         (2*CF*lcb^3)/9)*nl^2*TF^2 + nb^2*TF^2*((1685*CF)/7776 + 
         (35*CF*lb1)/216 + (5*CF*lb1^2)/18 - (2*CF*lb1^3)/9 - (7*CF*z3)/18) + 
       nc^2*TF^2*((1685*CF)/7776 - (35*CF*l21)/216 + (5*CF*l21^2)/18 + 
         (2*CF*l21^3)/9 + ((35*CF)/216 - (5*CF*l21)/9 - (2*CF*l21^2)/3)*lb1 + 
         ((5*CF)/18 + (2*CF*l21)/3)*lb1^2 - (2*CF*lb1^3)/9 - (7*CF*z3)/18) + 
       nb*nc*TF^2*((875*CF)/3888 - (CF*l21^3)/9 + 
         ((35*CF)/108 - (5*CF*l21)/9)*lb1 + ((5*CF)/9 + (2*CF*l21)/3)*lb1^2 - 
         (4*CF*lb1^3)/9 + (2*CF*Lpx)/3 + (5*CF)/(48*x^2) + (5*CF*x^2)/48 + 
         l21^2*((43*CF)/144 + (5*CF)/(96*x^2) + (5*CF*x^2)/96) + 
         l21*((-35*CF)/216 - (5*CF)/(48*x^2) + (5*CF*x^2)/48) + 
         Lmx*((-5*CF)/(96*x^3) - (9*CF)/(32*x) - (9*CF*x)/32 - 
           (5*CF*x^3)/96) + (2*CF*z3)/9) + 
       nl*TF*(((-197*CA*CF)/72 - (9*CF^2)/8)*lcb^2 + 
         ((-11*CA*CF)/9 - (3*CF^2)/4)*lcb^3 + lcb*((-139*CA*CF)/216 - 
           (23*CF^2)/16 - (3*CA*CF*z3)/2 + (3*CF^2*z3)/2)) + 
       nb*(nl*TF^2*((-1327*CF)/3888 + (53*CF*lb1)/72 + (2*CF*lb1^3)/9 + 
           ((89*CF)/216 - (5*CF*lb1)/9 + (2*CF*lb1^2)/3)*lcb + (2*CF*z3)/9) + 
         TF*((-16627*CA*CF)/15552 + (B4*CA*CF)/8 + (683*CF^2)/576 - 
           (B4*CF^2)/4 + ((-29*CA*CF)/24 + CF^2)*lb1^2 - (11*CA*CF*lb1^3)/
            18 + ((-979*CA*CF)/864 - (89*CF^2)/192 + 
             ((55*CA*CF)/36 + (5*CF^2)/8)*lb1 + ((-11*CA*CF)/6 - (3*CF^2)/4)*
              lb1^2)*lcb + (629*CA*CF*z3)/576 - (57*CF^2*z3)/32 + 
           lb1*((-5*CA*CF)/32 + (13*CF^2)/32 + (3*CA*CF*z3)/2 - 
             (3*CF^2*z3)/2) - (9*CA*CF*z4)/8 + (9*CF^2*z4)/8)) + 
       nc*(nl*TF^2*((-1327*CF)/3888 - (53*CF*l21)/72 - (2*CF*l21^3)/9 + 
           ((53*CF)/72 + (2*CF*l21^2)/3)*lb1 - (2*CF*l21*lb1^2)/3 + 
           (2*CF*lb1^3)/9 + ((89*CF)/216 + (5*CF*l21)/9 + (2*CF*l21^2)/3 + 
             ((-5*CF)/9 - (4*CF*l21)/3)*lb1 + (2*CF*lb1^2)/3)*lcb + 
           (2*CF*z3)/9) + TF*((-16627*CA*CF)/15552 + (B4*CA*CF)/8 + 
           (683*CF^2)/576 - (B4*CF^2)/4 + ((-29*CA*CF)/24 + CF^2)*l21^2 + 
           (11*CA*CF*l21^3)/18 + ((-29*CA*CF)/24 + CF^2 + (11*CA*CF*l21)/6)*
            lb1^2 - (11*CA*CF*lb1^3)/18 + ((-979*CA*CF)/864 - (89*CF^2)/192 + 
             ((-55*CA*CF)/36 - (5*CF^2)/8)*l21 + ((-11*CA*CF)/6 - (3*CF^2)/4)*
              l21^2 + ((55*CA*CF)/36 + (5*CF^2)/8 + ((11*CA*CF)/3 + 
                 (3*CF^2)/2)*l21)*lb1 + ((-11*CA*CF)/6 - (3*CF^2)/4)*lb1^2)*
            lcb + (629*CA*CF*z3)/576 - (57*CF^2*z3)/32 + 
           lb1*((-5*CA*CF)/32 + (13*CF^2)/32 + ((29*CA*CF)/12 - 2*CF^2)*l21 - 
             (11*CA*CF*l21^2)/6 + (3*CA*CF*z3)/2 - (3*CF^2*z3)/2) + 
           l21*((5*CA*CF)/32 - (13*CF^2)/32 - (3*CA*CF*z3)/2 + 
             (3*CF^2*z3)/2) - (9*CA*CF*z4)/8 + (9*CF^2*z4)/8)))
ZetamOS3lmubmuc = 1 - (3*api*CF*lcb)/2 + 
     api^2*(((-97*CA*CF)/48 - (3*CF^2)/16)*lcb + ((11*CA*CF)/8 + (9*CF^2)/8)*
        lcb^2 + ((89*CF)/288 - (5*CF*lb1)/12 + (CF*lb1^2)/2 + CF*lb1*lcb)*nb*
        TF + ((89*CF)/288 + (5*CF*l21)/12 + (CF*l21^2)/2 + 
         ((-5*CF)/12 - CF*l21)*lb1 + (CF*lb1^2)/2 + (-(CF*l21) + CF*lb1)*lcb)*
        nc*TF + ((5*CF*lcb)/12 - (CF*lcb^2)/2)*nl*TF) + 
     api^3*(((-11413*CA^2*CF)/3456 + (129*CA*CF^2)/128 - (129*CF^3)/64)*lcb + 
       ((1373*CA^2*CF)/288 + (27*CA*CF^2)/8 + (9*CF^3)/32)*lcb^2 + 
       ((-121*CA^2*CF)/72 - (33*CA*CF^2)/16 - (9*CF^3)/16)*lcb^3 + 
       ((35*CF*lcb)/216 + (5*CF*lcb^2)/18 - (2*CF*lcb^3)/9)*nl^2*TF^2 + 
       nb*nc*TF^2*((-875*CF)/3888 + (CF*l21^3)/9 + 
         ((-31*CF)/27 - (5*CF*l21)/9 - (2*CF*l21^2)/3)*lb1 + 
         ((5*CF)/9 + (4*CF*l21)/3)*lb1^2 - (8*CF*lb1^3)/9 + 
         ((4*CF*l21*lb1)/3 - (4*CF*lb1^2)/3)*lcb - (2*CF*Lpx)/3 - 
         (5*CF)/(48*x^2) - (5*CF*x^2)/48 + l21*((31*CF)/54 + 
           (5*CF)/(48*x^2) - (5*CF*x^2)/48) + l21^2*((-43*CF)/144 - 
           (5*CF)/(96*x^2) - (5*CF*x^2)/96) + Lmx*((5*CF)/(96*x^3) + 
           (9*CF)/(32*x) + (9*CF*x)/32 + (5*CF*x^3)/96) - (2*CF*z3)/9) + 
       nb^2*TF^2*((-1685*CF)/7776 - (31*CF*lb1)/54 + (5*CF*lb1^2)/18 - 
         (4*CF*lb1^3)/9 - (2*CF*lb1^2*lcb)/3 + (7*CF*z3)/18) + 
       nc^2*TF^2*((-1685*CF)/7776 + (31*CF*l21)/54 + (5*CF*l21^2)/18 + 
         (4*CF*l21^3)/9 + ((-31*CF)/54 - (5*CF*l21)/9 - (4*CF*l21^2)/3)*lb1 + 
         ((5*CF)/18 + (4*CF*l21)/3)*lb1^2 - (4*CF*lb1^3)/9 + 
         ((-2*CF*l21^2)/3 + (4*CF*l21*lb1)/3 - (2*CF*lb1^2)/3)*lcb + 
         (7*CF*z3)/18) + nl*TF*(((-197*CA*CF)/72 - (9*CF^2)/8)*lcb^2 + 
         ((11*CA*CF)/9 + (3*CF^2)/4)*lcb^3 + lcb*((139*CA*CF)/216 + 
           (23*CF^2)/16 + (3*CA*CF*z3)/2 - (3*CF^2*z3)/2)) + 
       nb*(nl*TF^2*((1327*CF)/3888 - (53*CF*lb1)/72 - (2*CF*lb1^3)/9 - 
           (5*CF*lb1*lcb)/9 + (2*CF*lb1*lcb^2)/3 - (2*CF*z3)/9) + 
         TF*((16627*CA*CF)/15552 - (B4*CA*CF)/8 - (923*CF^2)/576 + 
           (B4*CF^2)/4 + ((29*CA*CF)/24 + CF^2/2)*lb1^2 + (11*CA*CF*lb1^3)/
            18 + (-(CA*CF)/3 + (181*CF^2)/192 + ((71*CA*CF)/18 + (13*CF^2)/8)*
              lb1 - (3*CF^2*lb1^2)/4)*lcb + ((-11*CA*CF)/6 - (3*CF^2)/2)*lb1*
            lcb^2 - (629*CA*CF*z3)/576 + (57*CF^2*z3)/32 + 
           lb1*((5*CA*CF)/32 - CF^2/32 - (3*CA*CF*z3)/2 + (3*CF^2*z3)/2) + 
           (9*CA*CF*z4)/8 - (9*CF^2*z4)/8)) + 
       nc*(nl*TF^2*((1327*CF)/3888 + (53*CF*l21)/72 + (2*CF*l21^3)/9 + 
           ((-53*CF)/72 - (2*CF*l21^2)/3)*lb1 + (2*CF*l21*lb1^2)/3 - 
           (2*CF*lb1^3)/9 + ((5*CF*l21)/9 - (5*CF*lb1)/9)*lcb + 
           ((-2*CF*l21)/3 + (2*CF*lb1)/3)*lcb^2 - (2*CF*z3)/9) + 
         TF*((16627*CA*CF)/15552 - (B4*CA*CF)/8 - (923*CF^2)/576 + 
           (B4*CF^2)/4 + ((29*CA*CF)/24 + CF^2/2)*l21^2 - (11*CA*CF*l21^3)/
            18 + ((29*CA*CF)/24 + CF^2/2 - (11*CA*CF*l21)/6)*lb1^2 + 
           (11*CA*CF*lb1^3)/18 + (-(CA*CF)/3 + (181*CF^2)/192 + 
             ((-71*CA*CF)/18 - (13*CF^2)/8)*l21 - (3*CF^2*l21^2)/4 + 
             ((71*CA*CF)/18 + (13*CF^2)/8 + (3*CF^2*l21)/2)*lb1 - 
             (3*CF^2*lb1^2)/4)*lcb + (((11*CA*CF)/6 + (3*CF^2)/2)*l21 + 
             ((-11*CA*CF)/6 - (3*CF^2)/2)*lb1)*lcb^2 - (629*CA*CF*z3)/576 + 
           (57*CF^2*z3)/32 + l21*((-5*CA*CF)/32 + CF^2/32 + (3*CA*CF*z3)/2 - 
             (3*CF^2*z3)/2) + lb1*((5*CA*CF)/32 - CF^2/32 + 
             ((-29*CA*CF)/12 - CF^2)*l21 + (11*CA*CF*l21^2)/6 - 
             (3*CA*CF*z3)/2 + (3*CF^2*z3)/2) + (9*CA*CF*z4)/8 - 
           (9*CF^2*z4)/8)))
ZetamInvOS3lmubmuc = 1 + (3*anlpi*CF*lcb)/2 + 
     anlpi^2*(((97*CA*CF)/48 + (3*CF^2)/16)*lcb + ((11*CA*CF)/8 + (9*CF^2)/8)*
        lcb^2 + ((-89*CF)/288 + (5*CF*lb1)/12 - (CF*lb1^2)/2)*nb*TF + 
       ((-89*CF)/288 - (5*CF*l21)/12 - (CF*l21^2)/2 + ((5*CF)/12 + CF*l21)*
          lb1 - (CF*lb1^2)/2)*nc*TF + ((-5*CF*lcb)/12 - (CF*lcb^2)/2)*nl*
        TF) + anlpi^3*(((11413*CA^2*CF)/3456 - (129*CA*CF^2)/128 + 
         (129*CF^3)/64)*lcb + ((1373*CA^2*CF)/288 + (27*CA*CF^2)/8 + 
         (9*CF^3)/32)*lcb^2 + ((121*CA^2*CF)/72 + (33*CA*CF^2)/16 + 
         (9*CF^3)/16)*lcb^3 + ((-35*CF*lcb)/216 + (5*CF*lcb^2)/18 + 
         (2*CF*lcb^3)/9)*nl^2*TF^2 + nb^2*TF^2*((1685*CF)/7776 + 
         (35*CF*lb1)/216 + (5*CF*lb1^2)/18 - (2*CF*lb1^3)/9 - (7*CF*z3)/18) + 
       nc^2*TF^2*((1685*CF)/7776 - (35*CF*l21)/216 + (5*CF*l21^2)/18 + 
         (2*CF*l21^3)/9 + ((35*CF)/216 - (5*CF*l21)/9 - (2*CF*l21^2)/3)*lb1 + 
         ((5*CF)/18 + (2*CF*l21)/3)*lb1^2 - (2*CF*lb1^3)/9 - (7*CF*z3)/18) + 
       nb*nc*TF^2*((875*CF)/3888 - (CF*l21^3)/9 + 
         ((35*CF)/108 - (5*CF*l21)/9)*lb1 + ((5*CF)/9 + (2*CF*l21)/3)*lb1^2 - 
         (4*CF*lb1^3)/9 + (2*CF*Lpx)/3 + (5*CF)/(48*x^2) + (5*CF*x^2)/48 + 
         l21^2*((43*CF)/144 + (5*CF)/(96*x^2) + (5*CF*x^2)/96) + 
         l21*((-35*CF)/216 - (5*CF)/(48*x^2) + (5*CF*x^2)/48) + 
         Lmx*((-5*CF)/(96*x^3) - (9*CF)/(32*x) - (9*CF*x)/32 - 
           (5*CF*x^3)/96) + (2*CF*z3)/9) + 
       nl*TF*(((-197*CA*CF)/72 - (9*CF^2)/8)*lcb^2 + 
         ((-11*CA*CF)/9 - (3*CF^2)/4)*lcb^3 + lcb*((-139*CA*CF)/216 - 
           (23*CF^2)/16 - (3*CA*CF*z3)/2 + (3*CF^2*z3)/2)) + 
       nb*(nl*TF^2*((-1327*CF)/3888 + (53*CF*lb1)/72 + (2*CF*lb1^3)/9 + 
           ((89*CF)/216 - (5*CF*lb1)/9 + (2*CF*lb1^2)/3)*lcb + (2*CF*z3)/9) + 
         TF*((-16627*CA*CF)/15552 + (B4*CA*CF)/8 + (923*CF^2)/576 - 
           (B4*CF^2)/4 + ((-29*CA*CF)/24 - CF^2/2)*lb1^2 - 
           (11*CA*CF*lb1^3)/18 + ((-979*CA*CF)/864 - (89*CF^2)/192 + 
             ((55*CA*CF)/36 + (5*CF^2)/8)*lb1 + ((-11*CA*CF)/6 - (3*CF^2)/4)*
              lb1^2)*lcb + (629*CA*CF*z3)/576 - (57*CF^2*z3)/32 + 
           lb1*((-5*CA*CF)/32 + CF^2/32 + (3*CA*CF*z3)/2 - (3*CF^2*z3)/2) - 
           (9*CA*CF*z4)/8 + (9*CF^2*z4)/8)) + 
       nc*(nl*TF^2*((-1327*CF)/3888 - (53*CF*l21)/72 - (2*CF*l21^3)/9 + 
           ((53*CF)/72 + (2*CF*l21^2)/3)*lb1 - (2*CF*l21*lb1^2)/3 + 
           (2*CF*lb1^3)/9 + ((89*CF)/216 + (5*CF*l21)/9 + (2*CF*l21^2)/3 + 
             ((-5*CF)/9 - (4*CF*l21)/3)*lb1 + (2*CF*lb1^2)/3)*lcb + 
           (2*CF*z3)/9) + TF*((-16627*CA*CF)/15552 + (B4*CA*CF)/8 + 
           (923*CF^2)/576 - (B4*CF^2)/4 + ((-29*CA*CF)/24 - CF^2/2)*l21^2 + 
           (11*CA*CF*l21^3)/18 + ((-29*CA*CF)/24 - CF^2/2 + (11*CA*CF*l21)/6)*
            lb1^2 - (11*CA*CF*lb1^3)/18 + ((-979*CA*CF)/864 - (89*CF^2)/192 + 
             ((-55*CA*CF)/36 - (5*CF^2)/8)*l21 + ((-11*CA*CF)/6 - (3*CF^2)/4)*
              l21^2 + ((55*CA*CF)/36 + (5*CF^2)/8 + ((11*CA*CF)/3 + 
                 (3*CF^2)/2)*l21)*lb1 + ((-11*CA*CF)/6 - (3*CF^2)/4)*lb1^2)*
            lcb + (629*CA*CF*z3)/576 - (57*CF^2*z3)/32 + 
           lb1*((-5*CA*CF)/32 + CF^2/32 + ((29*CA*CF)/12 + CF^2)*l21 - 
             (11*CA*CF*l21^2)/6 + (3*CA*CF*z3)/2 - (3*CF^2*z3)/2) + 
           l21*((5*CA*CF)/32 - CF^2/32 - (3*CA*CF*z3)/2 + (3*CF^2*z3)/2) - 
           (9*CA*CF*z4)/8 + (9*CF^2*z4)/8)))
funcLabels = {l21 -> Log[M2/M1], lcb -> Log[muc/mub], lb1 -> Log[mub/M1], 
     Lpx -> Log[M2/M1]^3/6 - (Log[M2/M1]^2*Log[1 - M2/M1])/2 - 
       (Log[M2/M1]^2*Log[1 + M2/M1])/2 - Log[M2/M1]*PolyLog[2, -(M2/M1)] - 
       Log[M2/M1]*PolyLog[2, M2/M1] + PolyLog[3, -(M2/M1)] + 
       PolyLog[3, M2/M1], Lmx -> -(Log[M2/M1]^2*Log[1 - M2/M1])/2 + 
       (Log[M2/M1]^2*Log[1 + M2/M1])/2 + Log[M2/M1]*PolyLog[2, -(M2/M1)] - 
       Log[M2/M1]*PolyLog[2, M2/M1] - PolyLog[3, -(M2/M1)] + 
       PolyLog[3, M2/M1], Hp10xm1 -> -(Log[1 - M1/M2]*Log[M1/M2]) - 
       PolyLog[2, M1/M2], Hm10xm1 -> Log[1 + M1/M2]*Log[M1/M2] + 
       PolyLog[2, -(M1/M2)], Hp10xp1 -> -(Log[M2/M1]*Log[1 - M2/M1]) - 
       PolyLog[2, M2/M1], Hm10xp1 -> Log[M2/M1]*Log[1 + M2/M1] + 
       PolyLog[2, -(M2/M1)], x -> M2/M1}
constantLabels = {B4 -> (-13*Pi^4)/180 - (2*Pi^2*Log[2]^2)/3 + 
       (2*Log[2]^4)/3 + 16*PolyLog[4, 1/2], z2 -> Pi^2/6, z3 -> Zeta[3], 
     z4 -> Pi^4/90, log2 -> Log[2], TF -> 1/2, CA -> 3, CF -> 4/3, nb -> 1, 
     nc -> 1}
