(* ::Package:: *)

(*
Results for \zeta_m, C_ 1 and C_ 2

apinf -> alpha_s^{nf}/Pi
apinl -> alpha_s^{nl}/Pi 
lmMS  -> Log[(\mu/m_h)^2]
lmOS  -> Log[(\mu/M_h)^2]

l2      -> Log[2]
zn      -> Zeta[n] 
an      -> PolyLog[n,1/2]
*)


Print["... zetamMSnf"];

(zetamMSnf =

       + apinf^2 * (
          - 89/1152*nc^-1
          + 5/96*nc^-1*lmMS
          - 1/32*nc^-1*lmMS^2
          + 89/1152*nc
          - 5/96*nc*lmMS
          + 1/32*nc*lmMS^2
          )

       + apinf^3 * (
          + 907/31104
          - 683/4608*nc^-2
          - 13/512*nc^-2*lmMS
          - 1/32*nc^-2*lmMS^2
          + 1/2*nc^-2*a4
          - 1/8*nc^-2*l2^2*z2
          + 1/48*nc^-2*l2^4
          - 11/32*nc^-2*z4
          + 57/256*nc^-2*z3
          + 3/32*nc^-2*z3*lmMS
          + 1685/62208*nc^-1
          + 31/864*nc^-1*lmMS
          - 5/576*nc^-1*lmMS^2
          + 1/144*nc^-1*lmMS^3
          - 1327/31104*nc^-1*nl
          + 53/1152*nc^-1*nl*lmMS
          + 1/288*nc^-1*nl*lmMS^3
          - 7/144*nc^-1*z3
          + 1/36*nc^-1*z3*nl
          + 1/32*lmMS
          - 5/384*lmMS^2
          - 11/576*lmMS^3
          - 1/2*a4
          + 1/8*l2^2*z2
          - 1/48*l2^4
          + 13/64*z4
          - 397/2304*z3
          - 1685/62208*nc
          - 31/864*nc*lmMS
          + 5/576*nc*lmMS^2
          - 1/144*nc*lmMS^3
          + 1327/31104*nc*nl
          - 53/1152*nc*nl*lmMS
          - 1/288*nc*nl*lmMS^3
          + 7/144*nc*z3
          - 1/36*nc*z3*nl
          + 14813/124416*nc^2
          - 3/512*nc^2*lmMS
          + 17/384*nc^2*lmMS^2
          + 11/576*nc^2*lmMS^3
          + 9/64*nc^2*z4
          - 29/576*nc^2*z3
          - 3/32*nc^2*z3*lmMS
          )

       + apinf^4 * (
          - 2112311/81285120
          + 38579/163840*nc^-3
          - 793/12288*nc^-3*lmMS
          - 15/256*nc^-3*lmMS^2
          + 655/1024*nc^-3*z6
          + 7/5*nc^-3*a5
          - 32027/3840*nc^-3*a4
          + 32027/15360*nc^-3*l2^2*z2
          + 7/60*nc^-3*l2^3*z2
          - 32027/92160*nc^-3*l2^4
          - 7/600*nc^-3*l2^5
          - 5953/5120*nc^-3*z5
          + 15/128*nc^-3*z5*lmMS
          + 150947/30720*nc^-3*z4
          + 253/160*nc^-3*z4*l2
          - 46639/15360*nc^-3*z3
          + 1271/4096*nc^-3*z3*lmMS
          - 447/512*nc^-3*z3^2
          + 257128337/3251404800*nc^-2
          + 8951/110592*nc^-2*lmMS
          + 193/36864*nc^-2*lmMS^2
          + 31/3072*nc^-2*lmMS^3
          + 1/2048*nc^-2*lmMS^4
          + 5041/26880*nc^-2*a4
          - 1/4*nc^-2*a4*lmMS
          - 5041/107520*nc^-2*l2^2*z2
          + 1/16*nc^-2*l2^2*z2*lmMS
          + 5041/645120*nc^-2*l2^4
          - 1/96*nc^-2*l2^4*lmMS
          - 15557/82944*nc^-2*nl
          + 319/3072*nc^-2*nl*lmMS
          - 9/1024*nc^-2*nl*lmMS^2
          + 1/128*nc^-2*nl*lmMS^3
          - 1/2*nc^-2*nl*a5
          - 49/96*nc^-2*nl*a4
          - 1/4*nc^-2*nl*a4*lmMS
          + 49/384*nc^-2*nl*l2^2*z2
          + 1/16*nc^-2*nl*l2^2*z2*lmMS
          - 1/24*nc^-2*nl*l2^3*z2
          - 49/2304*nc^-2*nl*l2^4
          - 1/96*nc^-2*nl*l2^4*lmMS
          + 1/240*nc^-2*nl*l2^5
          - 63/256*nc^-2*z5
          + 157/256*nc^-2*z5*nl
          - 90269/430080*nc^-2*z4
          + 25/128*nc^-2*z4*lmMS
          + 1639/3072*nc^-2*z4*nl
          + 7/32*nc^-2*z4*nl*lmMS
          - 49/64*nc^-2*z4*nl*l2
          + 7671973/30105600*nc^-2*z3
          - 595/4096*nc^-2*z3*lmMS
          - 3/128*nc^-2*z3*lmMS^2
          - 5/128*nc^-2*z3*nl
          - 97/512*nc^-2*z3*nl*lmMS
          + 272586767/2508226560*nc^-1
          - 468247/995328*nc^-1*lmMS
          - 1301/27648*nc^-1*lmMS^2
          - 1247/27648*nc^-1*lmMS^3
          - 1/768*nc^-1*lmMS^4
          - 595/512*nc^-1*z6
          + 1/12*nc^-1*a5
          + 8675/512*nc^-1*a4
          + 11/8*nc^-1*a4*lmMS
          - 8675/2048*nc^-1*l2^2*z2
          - 11/32*nc^-1*l2^2*z2*lmMS
          + 1/144*nc^-1*l2^3*z2
          + 8675/12288*nc^-1*l2^4
          + 11/192*nc^-1*l2^4*lmMS
          - 1/1440*nc^-1*l2^5
          + 57/4096*nc^-1*nl
          + 143/10368*nc^-1*nl*lmMS
          - 283/13824*nc^-1*nl*lmMS^2
          - 1/768*nc^-1*nl*lmMS^4
          - 17671/331776*nc^-1*nl^2
          + 3401/124416*nc^-1*nl^2*lmMS
          - 31/3456*nc^-1*nl^2*lmMS^2
          - 1/2304*nc^-1*nl^2*lmMS^4
          - 415/1536*nc^-1*z5
          - 5/16*nc^-1*z5*lmMS
          - 297695/24576*nc^-1*z4
          - 275/256*nc^-1*z4*lmMS
          + 73/96*nc^-1*z4*l2
          - 9/256*nc^-1*z4*nl
          + 7/256*nc^-1*z4*nl^2
          + 2055427/344064*nc^-1*z3
          + 18773/24576*nc^-1*z3*lmMS
          + 33/256*nc^-1*z3*lmMS^2
          + 5/2304*nc^-1*z3*nl
          + 5/2304*nc^-1*z3*nl^2
          - 7/288*nc^-1*z3*nl^2*lmMS
          + 719/512*nc^-1*z3^2
          + 6203/82944*lmMS
          + 185/6912*lmMS^2
          + 7/512*lmMS^3
          + 19/3072*lmMS^4
          - 12913/8960*a4
          + 1/4*a4*lmMS
          + 12913/35840*l2^2*z2
          - 1/16*l2^2*z2*lmMS
          - 12913/215040*l2^4
          + 1/96*l2^4*lmMS
          + 386477/663552*nl
          - 46523/248832*nl*lmMS
          + 3067/27648*nl*lmMS^2
          + 1/768*nl*lmMS^3
          + 11/2304*nl*lmMS^4
          + 1/2*nl*a5
          + 49/96*nl*a4
          + 1/4*nl*a4*lmMS
          - 49/384*nl*l2^2*z2
          - 1/16*nl*l2^2*z2*lmMS
          + 1/24*nl*l2^3*z2
          + 49/2304*nl*l2^4
          + 1/96*nl*l2^4*lmMS
          - 1/240*nl*l2^5
          - 5/384*z5
          - 689/768*z5*nl
          + 127895/86016*z4
          - 13/128*z4*lmMS
          - 1649/3072*z4*nl
          - 13/128*z4*nl*lmMS
          + 49/64*z4*nl*l2
          - 3038363/3612672*z3
          - 773/12288*z3*lmMS
          + 185/1536*z3*nl
          + 749/4608*z3*nl*lmMS
          - 5207600543/2508226560*nc
          + 425053/497664*nc*lmMS
          - 1813/18432*nc*lmMS^2
          + 43/13824*nc*lmMS^3
          - 109/9216*nc*lmMS^4
          + 975/1024*nc*z6
          - 89/60*nc*a5
          - 66071/7680*nc*a4
          - 11/8*nc*a4*lmMS
          + 66071/30720*nc*l2^2*z2
          + 11/32*nc*l2^2*z2*lmMS
          - 89/720*nc*l2^3*z2
          - 66071/184320*nc*l2^4
          - 11/192*nc*l2^4*lmMS
          + 89/7200*nc*l2^5
          - 57/4096*nc*nl
          - 143/10368*nc*nl*lmMS
          + 283/13824*nc*nl*lmMS^2
          + 1/768*nc*nl*lmMS^4
          + 17671/331776*nc*nl^2
          - 3401/124416*nc*nl^2*lmMS
          + 31/3456*nc*nl^2*lmMS^2
          + 1/2304*nc*nl^2*lmMS^4
          + 47999/15360*nc*z5
          - 15/128*nc*z5*lmMS
          + 260479/40960*nc*z4
          + 143/256*nc*z4*lmMS
          - 281/120*nc*z4*l2
          + 9/256*nc*z4*nl
          - 7/256*nc*z4*nl^2
          - 13823281/5160960*nc*z3
          - 41293/73728*nc*z3*lmMS
          - 5/2304*nc*z3*nl
          - 5/2304*nc*z3*nl^2
          + 7/288*nc*z3*nl^2*lmMS
          - 21/32*nc*z3^2
          - 8220757/154828800*nc^2
          - 51665/331776*nc^2*lmMS
          - 3539/110592*nc^2*lmMS^2
          - 73/3072*nc^2*lmMS^3
          - 41/6144*nc^2*lmMS^4
          + 2407/1920*nc^2*a4
          - 2407/7680*nc^2*l2^2*z2
          + 2407/46080*nc^2*l2^4
          - 262021/663552*nc^2*nl
          + 5171/62208*nc^2*nl*lmMS
          - 353/3456*nc^2*nl*lmMS^2
          - 7/768*nc^2*nl*lmMS^3
          - 11/2304*nc^2*nl*lmMS^4
          + 199/768*nc^2*z5
          + 109/384*nc^2*z5*nl
          - 39229/30720*nc^2*z4
          - 3/32*nc^2*z4*lmMS
          + 5/1536*nc^2*z4*nl
          - 15/128*nc^2*z4*nl*lmMS
          + 1890827/3225600*nc^2*z3
          + 1279/6144*nc^2*z3*lmMS
          + 3/128*nc^2*z3*lmMS^2
          - 125/1536*nc^2*z3*nl
          + 31/1152*nc^2*z3*nl*lmMS
          + 1532419/884736*nc^3
          - 158813/497664*nc^3*lmMS
          + 11281/55296*nc^3*lmMS^2
          + 43/1024*nc^3*lmMS^3
          + 121/9216*nc^3*lmMS^4
          - 55/128*nc^3*z6
          - 2599/1536*nc^3*z5
          + 5/16*nc^3*z5*lmMS
          + 10325/12288*nc^3*z4
          + 33/64*nc^3*z4*lmMS
          - 9553/36864*nc^3*z3
          - 2369/4608*nc^3*z3*lmMS
          - 33/256*nc^3*z3*lmMS^2
          + 1/8*nc^3*z3^2
          )

       + 1);


Print["... zetamOSnf"];

(zetamOSnf =

       + apinf^2 * (
          - 89/1152*nc^-1
          + 5/96*nc^-1*lmOS
          - 1/32*nc^-1*lmOS^2
          + 89/1152*nc
          - 5/96*nc*lmOS
          + 1/32*nc*lmOS^2
          )

       + apinf^3 * (
          + 4147/31104
          - 923/4608*nc^-2
          - 1/512*nc^-2*lmOS
          + 1/64*nc^-2*lmOS^2
          + 1/2*nc^-2*a4
          - 1/8*nc^-2*l2^2*z2
          + 1/48*nc^-2*l2^4
          - 11/32*nc^-2*z4
          + 57/256*nc^-2*z3
          + 3/32*nc^-2*z3*lmOS
          + 1685/62208*nc^-1
          + 31/864*nc^-1*lmOS
          - 5/576*nc^-1*lmOS^2
          + 1/144*nc^-1*lmOS^3
          - 1327/31104*nc^-1*nl
          + 53/1152*nc^-1*nl*lmOS
          + 1/288*nc^-1*nl*lmOS^3
          - 7/144*nc^-1*z3
          + 1/36*nc^-1*z3*nl
          - 1/64*lmOS
          - 41/384*lmOS^2
          - 11/576*lmOS^3
          - 1/2*a4
          + 1/8*l2^2*z2
          - 1/48*l2^4
          + 13/64*z4
          - 397/2304*z3
          - 1685/62208*nc
          - 31/864*nc*lmOS
          + 5/576*nc*lmOS^2
          - 1/144*nc*lmOS^3
          + 1327/31104*nc*nl
          - 53/1152*nc*nl*lmOS
          - 1/288*nc*nl*lmOS^3
          + 7/144*nc*z3
          - 1/36*nc*z3*nl
          + 8333/124416*nc^2
          + 9/512*nc^2*lmOS
          + 35/384*nc^2*lmOS^2
          + 11/576*nc^2*lmOS^3
          + 9/64*nc^2*z4
          - 29/576*nc^2*z3
          - 3/32*nc^2*z3*lmOS
          )

       + apinf^4 * (
          - 2585651/81285120
          + 39519/163840*nc^-3
          - 127/3072*nc^-3*lmOS
          - 33/1024*nc^-3*lmOS^2
          + 655/1024*nc^-3*z6
          + 25/512*nc^-3*z2
          - 15/256*nc^-3*z2*lmOS
          + 7/5*nc^-3*a5
          - 32027/3840*nc^-3*a4
          - 5/64*nc^-3*l2*z2
          + 3/32*nc^-3*l2*z2*lmOS
          + 32027/15360*nc^-3*l2^2*z2
          + 7/60*nc^-3*l2^3*z2
          - 32027/92160*nc^-3*l2^4
          - 7/600*nc^-3*l2^5
          - 5953/5120*nc^-3*z5
          + 15/128*nc^-3*z5*lmOS
          + 150947/30720*nc^-3*z4
          + 253/160*nc^-3*z4*l2
          - 47779/15360*nc^-3*z3
          + 887/4096*nc^-3*z3*lmOS
          - 447/512*nc^-3*z3^2
          + 266595137/3251404800*nc^-2
          + 4307/110592*nc^-2*lmOS
          - 599/36864*nc^-2*lmOS^2
          - 29/3072*nc^-2*lmOS^3
          + 1/2048*nc^-2*lmOS^4
          - 5/192*nc^-2*z2
          + 1/32*nc^-2*z2*lmOS
          + 5041/26880*nc^-2*a4
          - 1/4*nc^-2*a4*lmOS
          - 5041/107520*nc^-2*l2^2*z2
          + 1/16*nc^-2*l2^2*z2*lmOS
          + 5041/645120*nc^-2*l2^4
          - 1/96*nc^-2*l2^4*lmOS
          - 35551/165888*nc^-2*nl
          + 139/2304*nc^-2*nl*lmOS
          - 101/3072*nc^-2*nl*lmOS^2
          - 1/256*nc^-2*nl*lmOS^3
          + 5/384*nc^-2*nl*z2
          - 1/64*nc^-2*nl*z2*lmOS
          - 1/2*nc^-2*nl*a5
          - 49/96*nc^-2*nl*a4
          - 1/4*nc^-2*nl*a4*lmOS
          + 49/384*nc^-2*nl*l2^2*z2
          + 1/16*nc^-2*nl*l2^2*z2*lmOS
          - 1/24*nc^-2*nl*l2^3*z2
          - 49/2304*nc^-2*nl*l2^4
          - 1/96*nc^-2*nl*l2^4*lmOS
          + 1/240*nc^-2*nl*l2^5
          - 63/256*nc^-2*z5
          + 157/256*nc^-2*z5*nl
          - 90269/430080*nc^-2*z4
          + 25/128*nc^-2*z4*lmOS
          + 1639/3072*nc^-2*z4*nl
          + 7/32*nc^-2*z4*nl*lmOS
          - 49/64*nc^-2*z4*nl*l2
          + 7671973/30105600*nc^-2*z3
          - 595/4096*nc^-2*z3*lmOS
          - 3/128*nc^-2*z3*lmOS^2
          - 5/128*nc^-2*z3*nl
          - 97/512*nc^-2*z3*nl*lmOS
          - 99558013/2508226560*nc^-1
          - 36563/124416*nc^-1*lmOS
          + 4063/27648*nc^-1*lmOS^2
          + 535/27648*nc^-1*lmOS^3
          - 1/768*nc^-1*lmOS^4
          - 595/512*nc^-1*z6
          - 185/1536*nc^-1*z2
          + 37/256*nc^-1*z2*lmOS
          + 1/12*nc^-1*a5
          + 8675/512*nc^-1*a4
          + 11/8*nc^-1*a4*lmOS
          + 5/32*nc^-1*l2*z2
          - 3/16*nc^-1*l2*z2*lmOS
          - 8675/2048*nc^-1*l2^2*z2
          - 11/32*nc^-1*l2^2*z2*lmOS
          + 1/144*nc^-1*l2^3*z2
          + 8675/12288*nc^-1*l2^4
          + 11/192*nc^-1*l2^4*lmOS
          - 1/1440*nc^-1*l2^5
          + 57/4096*nc^-1*nl
          + 143/10368*nc^-1*nl*lmOS
          - 283/13824*nc^-1*nl*lmOS^2
          - 1/768*nc^-1*nl*lmOS^4
          - 17671/331776*nc^-1*nl^2
          + 3401/124416*nc^-1*nl^2*lmOS
          - 31/3456*nc^-1*nl^2*lmOS^2
          - 1/2304*nc^-1*nl^2*lmOS^4
          - 415/1536*nc^-1*z5
          - 5/16*nc^-1*z5*lmOS
          - 297695/24576*nc^-1*z4
          - 275/256*nc^-1*z4*lmOS
          + 73/96*nc^-1*z4*l2
          - 9/256*nc^-1*z4*nl
          + 7/256*nc^-1*z4*nl^2
          + 2074243/344064*nc^-1*z3
          + 21653/24576*nc^-1*z3*lmOS
          + 33/256*nc^-1*z3*lmOS^2
          + 5/2304*nc^-1*z3*nl
          + 5/2304*nc^-1*z3*nl^2
          - 7/288*nc^-1*z3*nl^2*lmOS
          + 719/512*nc^-1*z3^2
          + 13169/82944*lmOS
          + 241/3456*lmOS^2
          + 27/512*lmOS^3
          + 19/3072*lmOS^4
          + 5/96*z2
          - 1/16*z2*lmOS
          - 12913/8960*a4
          + 1/4*a4*lmOS
          + 12913/35840*l2^2*z2
          - 1/16*l2^2*z2*lmOS
          - 12913/215040*l2^4
          + 1/96*l2^4*lmOS
          + 421973/663552*nl
          - 24869/248832*nl*lmOS
          + 4399/27648*nl*lmOS^2
          + 19/768*nl*lmOS^3
          + 11/2304*nl*lmOS^4
          - 5/192*nl*z2
          + 1/32*nl*z2*lmOS
          + 1/2*nl*a5
          + 49/96*nl*a4
          + 1/4*nl*a4*lmOS
          - 49/384*nl*l2^2*z2
          - 1/16*nl*l2^2*z2*lmOS
          + 1/24*nl*l2^3*z2
          + 49/2304*nl*l2^4
          + 1/96*nl*l2^4*lmOS
          - 1/240*nl*l2^5
          - 5/384*z5
          - 689/768*z5*nl
          + 127895/86016*z4
          - 13/128*z4*lmOS
          - 1649/3072*z4*nl
          - 13/128*z4*nl*lmOS
          + 49/64*z4*nl*l2
          - 3038363/3612672*z3
          - 773/12288*z3*lmOS
          + 185/1536*z3*nl
          + 749/4608*z3*nl*lmOS
          - 4506482363/2508226560*nc
          + 429365/995328*nc*lmOS
          - 10423/18432*nc*lmOS^2
          - 1739/13824*nc*lmOS^3
          - 109/9216*nc*lmOS^4
          + 975/1024*nc*z6
          + 145/1536*nc*z2
          - 29/256*nc*z2*lmOS
          - 89/60*nc*a5
          - 66071/7680*nc*a4
          - 11/8*nc*a4*lmOS
          - 5/64*nc*l2*z2
          + 3/32*nc*l2*z2*lmOS
          + 66071/30720*nc*l2^2*z2
          + 11/32*nc*l2^2*z2*lmOS
          - 89/720*nc*l2^3*z2
          - 66071/184320*nc*l2^4
          - 11/192*nc*l2^4*lmOS
          + 89/7200*nc*l2^5
          - 57/4096*nc*nl
          - 143/10368*nc*nl*lmOS
          + 283/13824*nc*nl*lmOS^2
          + 1/768*nc*nl*lmOS^4
          + 17671/331776*nc*nl^2
          - 3401/124416*nc*nl^2*lmOS
          + 31/3456*nc*nl^2*lmOS^2
          + 1/2304*nc*nl^2*lmOS^4
          + 47999/15360*nc*z5
          - 15/128*nc*z5*lmOS
          + 260479/40960*nc*z4
          + 143/256*nc*z4*lmOS
          - 281/120*nc*z4*l2
          + 9/256*nc*z4*nl
          - 7/256*nc*z4*nl^2
          - 13238641/5160960*nc*z3
          - 37837/73728*nc*z3*lmOS
          - 5/2304*nc*z3*nl
          - 5/2304*nc*z3*nl^2
          + 7/288*nc*z3*nl^2*lmOS
          - 21/32*nc*z3^2
          - 7769957/154828800*nc^2
          - 65597/331776*nc^2*lmOS
          - 5915/110592*nc^2*lmOS^2
          - 133/3072*nc^2*lmOS^3
          - 41/6144*nc^2*lmOS^4
          - 5/192*nc^2*z2
          + 1/32*nc^2*z2*lmOS
          + 2407/1920*nc^2*a4
          - 2407/7680*nc^2*l2^2*z2
          + 2407/46080*nc^2*l2^4
          - 279769/663552*nc^2*nl
          + 9857/248832*nc^2*nl*lmOS
          - 1745/13824*nc^2*nl*lmOS^2
          - 1/48*nc^2*nl*lmOS^3
          - 11/2304*nc^2*nl*lmOS^4
          + 5/384*nc^2*nl*z2
          - 1/64*nc^2*nl*z2*lmOS
          + 199/768*nc^2*z5
          + 109/384*nc^2*z5*nl
          - 39229/30720*nc^2*z4
          - 3/32*nc^2*z4*lmOS
          + 5/1536*nc^2*z4*nl
          - 15/128*nc^2*z4*nl*lmOS
          + 1890827/3225600*nc^2*z3
          + 1279/6144*nc^2*z3*lmOS
          + 3/128*nc^2*z3*lmOS^2
          - 125/1536*nc^2*z3*nl
          + 31/1152*nc^2*z3*nl*lmOS
          + 1411303/884736*nc^3
          - 95713/995328*nc^3*lmOS
          + 24925/55296*nc^3*lmOS^2
          + 109/1024*nc^3*lmOS^3
          + 121/9216*nc^3*lmOS^4
          - 55/128*nc^3*z6
          - 35/1536*nc^3*z2
          + 7/256*nc^3*z2*lmOS
          - 2599/1536*nc^3*z5
          + 5/16*nc^3*z5*lmOS
          + 10325/12288*nc^3*z4
          + 33/64*nc^3*z4*lmOS
          - 13009/36864*nc^3*z3
          - 2693/4608*nc^3*z3*lmOS
          - 33/256*nc^3*z3*lmOS^2
          + 1/8*nc^3*z3^2
          )

       + 1); 


Print["... C2MSnf"];

   (C2MSnf =

        + apinf^2 * (
          + 5/18
          - 1/3*lmMS
          )

       + apinf^3 * (
          + 311/1296
          - 175/108*lmMS
          - 29/36*lmMS^2
          + 53/216*nl
          + 1/18*nl*lmMS^2
          + 5/3*z3
          )

       + apinf^4 * (
          + 2800175/186624
          - 50885/2592*lmMS
          - 1219/216*lmMS^2
          - 301/144*lmMS^3
          + 62/9*a4
          - 31/18*l2^2*z2
          + 31/108*l2^4
          - 16669/15552*nl
          + 7825/2592*nl*lmMS
          + 23/48*nl*lmMS^2
          + 5/18*nl*lmMS^3
          - 4/9*nl*a4
          + 1/9*nl*l2^2*z2
          - 1/54*nl*l2^4
          + 3401/23328*nl^2
          - 31/324*nl^2*lmMS
          - 1/108*nl^2*lmMS^3
          - 575/36*z5
          - 4123/144*z4
          + 163/72*z4*nl
          + 373261/13824*z3
          + 155/12*z3*lmMS
          - 221/288*z3*nl
          - 7/54*z3*nl^2
          )

       + apinf^5 * (
          - 6555736362659/12345177600
          + 26299405/1119744*lmMS
          - 359921/3888*lmMS^2
          - 315317/15552*lmMS^3
          - 27985/5184*lmMS^4
          + 829565/2592*z6
          + 56875/432*z7
          + 852973/5670*a5
          + 4299029/12960*a4
          + 1922/27*a4*lmMS
          - 4299029/51840*l2^2*z2
          - 961/54*l2^2*z2*lmMS
          + 852973/68040*l2^3*z2
          + 4299029/311040*l2^4
          + 961/324*l2^4*lmMS
          - 852973/680400*l2^5
          + 516545017993/4938071040*nl
          - 2081491/559872*nl*lmMS
          + 701407/31104*nl*lmMS^2
          + 527/162*nl*lmMS^3
          + 307/288*nl*lmMS^4
          - 16495/648*nl*z6
          - 1216/81*nl*a5
          - 6661933/272160*nl*a4
          - 248/27*nl*a4*lmMS
          + 6661933/1088640*nl*l2^2*z2
          + 62/27*nl*l2^2*z2*lmMS
          - 304/243*nl*l2^3*z2
          - 6661933/6531840*nl*l2^4
          - 31/81*nl*l2^4*lmMS
          + 152/1215*nl*l2^5
          - 8035069/1119744*nl^2
          + 256807/139968*nl^2*lmMS
          - 7867/5184*nl^2*lmMS^2
          - 163/1296*nl^2*lmMS^3
          - 91/1296*nl^2*lmMS^4
          + 16/27*nl^2*a5
          + 49/81*nl^2*a4
          + 8/27*nl^2*a4*lmMS
          - 49/324*nl^2*l2^2*z2
          - 2/27*nl^2*l2^2*z2*lmMS
          + 4/81*nl^2*l2^3*z2
          + 49/1944*nl^2*l2^4
          + 1/81*nl^2*l2^4*lmMS
          - 2/405*nl^2*l2^5
          + 18451/93312*nl^3
          - 197/2187*nl^3*lmMS
          + 31/972*nl^3*lmMS^2
          + 1/648*nl^3*lmMS^4
          + 639619703/2177280*z5
          - 12875/108*z5*lmMS
          - 10847/162*z5*nl
          + 737/216*z5*nl^2
          - 1122122327/1451520*z4
          - 127813/432*z4*lmMS
          + 34490467/181440*z4*l2
          + 111452531/1741824*z4*nl
          + 4093/108*z4*nl*lmMS
          - 7457/324*z4*nl*l2
          + 1853/2592*z4*nl^2
          - 133/108*z4*nl^2*lmMS
          + 49/54*z4*nl^2*l2
          - 25/216*z4*nl^3
          + 188941403323/609638400*z3
          + 18568279/82944*z3*lmMS
          + 4805/72*z3*lmMS^2
          - 285285281/13547520*z3*nl
          - 236143/41472*z3*nl*lmMS
          - 55/12*z3*nl*lmMS^2
          - 2141/1296*z3*nl^2
          - 1673/1296*z3*nl^2*lmMS
          + 5/1944*z3*nl^3
          + 2/27*z3*nl^3*lmMS
          - 21601/864*z3^2
          - 613/144*z3^2*nl
          )

       + 1);



Print["... C2OSnl"];

   (C2OSnl =

       + apinl^2 * (
          + 5/18
          - 1/3*lmOS
          )

       + apinl^3 * (
          - 841/1296
          - 79/36*lmOS
          - 11/12*lmOS^2
          + 53/216*nl
          + 1/18*nl*lmOS^2
          + 5/3*z3
          )

       + apinl^4 * (
          + 609215/186624
          - 4645/144*lmOS
          - 91/8*lmOS^2
          - 121/48*lmOS^3
          - 4/3*z2
          + 62/9*a4
          - 4/9*l2*z2
          - 31/18*l2^2*z2
          + 31/108*l2^4
          - 11557/15552*nl
          + 9535/2592*nl*lmOS
          + 109/144*nl*lmOS^2
          + 11/36*nl*lmOS^3
          + 2/9*nl*z2
          - 4/9*nl*a4
          + 1/9*nl*l2^2*z2
          - 1/54*nl*l2^4
          + 3401/23328*nl^2
          - 31/324*nl^2*lmOS
          - 1/108*nl^2*lmOS^3
          - 575/36*z5
          - 4123/144*z4
          + 163/72*z4*nl
          + 374797/13824*z3
          + 55/4*z3*lmOS
          - 221/288*z3*nl
          - 7/54*z3*nl^2
          )

       + apinl^5 * (
          - 4344482572567/6172588800
          - 63682459/373248*lmOS
          - 360983/1728*lmOS^2
          - 27599/576*lmOS^3
          - 1331/192*lmOS^4
          + 829565/2592*z6
          - 720881/9720*z2
          - 44/3*z2*lmOS
          + 56875/432*z7
          + 852973/5670*a5
          + 4369429/12960*a4
          + 682/9*a4*lmOS
          + 961/81*l2*z2
          - 44/9*l2*z2*lmOS
          - 4242709/51840*l2^2*z2
          - 341/18*l2^2*z2*lmOS
          + 852973/68040*l2^3*z2
          + 4369429/311040*l2^4
          + 341/108*l2^4*lmOS
          - 852973/680400*l2^5
          + 603951764833/4938071040*nl
          + 9284417/559872*nl*lmOS
          + 356177/10368*nl*lmOS^2
          + 2723/432*nl*lmOS^3
          + 121/96*nl*lmOS^4
          - 16495/648*nl*z6
          + 590/81*nl*z2
          + 10/3*nl*z2*lmOS
          - 1216/81*nl*a5
          - 6715693/272160*nl*a4
          - 256/27*nl*a4*lmOS
          + 44/81*nl*l2*z2
          + 8/27*nl*l2*z2*lmOS
          + 6554413/1088640*nl*l2^2*z2
          + 64/27*nl*l2^2*z2*lmOS
          - 304/243*nl*l2^3*z2
          - 6715693/6531840*nl*l2^4
          - 32/81*nl*l2^4*lmOS
          + 152/1215*nl*l2^5
          - 2798687/373248*nl^2
          + 71825/46656*nl^2*lmOS
          - 27617/15552*nl^2*lmOS^2
          - 89/432*nl^2*lmOS^3
          - 11/144*nl^2*lmOS^4
          - 13/81*nl^2*z2
          - 4/27*nl^2*z2*lmOS
          + 16/27*nl^2*a5
          + 49/81*nl^2*a4
          + 8/27*nl^2*a4*lmOS
          - 49/324*nl^2*l2^2*z2
          - 2/27*nl^2*l2^2*z2*lmOS
          + 4/81*nl^2*l2^3*z2
          + 49/1944*nl^2*l2^4
          + 1/81*nl^2*l2^4*lmOS
          - 2/405*nl^2*l2^5
          + 18451/93312*nl^3
          - 197/2187*nl^3*lmOS
          + 31/972*nl^3*lmOS^2
          + 1/648*nl^3*lmOS^4
          + 626347703/2177280*z5
          - 4675/36*z5*lmOS
          - 10847/162*z5*nl
          + 737/216*z5*nl^2
          - 371446109/483840*z4
          - 45353/144*z4*lmOS
          + 34490467/181440*z4*l2
          + 108173171/1741824*z4*nl
          + 1064/27*z4*nl*lmOS
          - 7457/324*z4*nl*l2
          + 1853/2592*z4*nl^2
          - 133/108*z4*nl^2*lmOS
          + 49/54*z4*nl^2*l2
          - 25/216*z4*nl^3
          + 70438900691/203212800*z3
          + 22515733/82944*z3*lmOS
          + 605/8*z3*lmOS^2
          + 1439/108*z3*z2
          - 255054241/13547520*z3*nl
          - 237391/41472*z3*nl*lmOS
          - 55/12*z3*nl*lmOS^2
          - 751/432*z3*nl^2
          - 595/432*z3*nl^2*lmOS
          + 5/1944*z3*nl^3
          + 2/27*z3*nl^3*lmOS
          - 21601/864*z3^2
          - 613/144*z3^2*nl
          )

       + 1);


Print["... C1OSnl"];

 (C1OSnl =
       + apinl * (
          - 1/12
          )

       + apinl^2 * (
          - 11/48
          )

       + apinl^3 * (
          - 2777/3456
          - 19/192*lmOS
          + 67/1152*nl
          - 1/36*nl*lmOS
          )

       + apinl^4 * (
          + 2761331/497664
          - 2417/3456*lmOS
          - 209/768*lmOS^2
          - 58723/248832*nl
          - 91/648*nl*lmOS
          - 23/384*nl*lmOS^2
          + 6865/373248*nl^2
          - 77/20736*nl^2*lmOS
          + 1/216*nl^2*lmOS^2
          - 897943/110592*z3
          + 110779/165888*z3*nl
          )

       + apinl^5 * (
          + 330267654131/1567641600
          + 46543033/663552*lmOS
          - 14747/4608*lmOS^2
          - 2299/3072*lmOS^3
          + 1/2*b4nf
          - 1/2*b4nl
          - 19/48*z2
          + 10566443/90720*a5
          + 93970579/435456*a4
          - 19/144*l2*z2
          - 93970579/1741824*l2^2*z2
          + 10566443/1088640*l2^3*z2
          + 93970579/10450944*l2^4
          - 10566443/10886400*l2^5
          - 40412995727/1567641600*nl
          - 5581849/497664*nl*lmOS
          - 12547/27648*nl*lmOS^2
          - 275/2304*nl*lmOS^3
          - 13/288*nl*z2
          - 340853/45360*nl*a5
          - 11679301/870912*nl*a4
          - 1/27*nl*l2*z2
          + 11679301/3483648*nl*l2^2*z2
          - 340853/544320*nl*l2^3*z2
          - 11679301/20901888*nl*l2^4
          + 340853/5443200*nl*l2^5
          + 1272247/3981312*nl^2
          + 22687/41472*nl^2*lmOS
          + 5107/165888*nl^2*lmOS^2
          + 157/6912*nl^2*lmOS^3
          + 1/54*nl^2*z2
          - 685/20736*nl^2*a4
          + 685/82944*nl^2*l2^2*z2
          - 685/497664*nl^2*l2^4
          + 255947/17915904*nl^3
          - 481/62208*nl^3*lmOS
          + 77/82944*nl^3*lmOS^2
          - 1/1296*nl^3*lmOS^3
          - 130269691/1290240*z5
          + 5376961/645120*z5*nl
          - 115/2304*z5*nl^2
          - 2375772389/17418240*z4
          + 288872477/2903040*z4*l2
          - 461541691/139345920*z4*nl
          - 9318467/1451520*z4*nl*l2
          + 518509/663552*z4*nl^2
          - 66998691323/696729600*z3
          - 9364157/147456*z3*lmOS
          + 8402310869/464486400*z3*nl
          + 32257/3456*z3*nl*lmOS
          - 2900089/3981312*z3*nl^2
          - 28297/110592*z3*nl^2*lmOS
          - 5/768*z3*nl^3
		   )
  );

