
Meql     = {"T[6,11]",  "T[7, 5]",  "T[7, 7]",  "T[7, 8]",  "T[7, 9]",  "T[7,10]",  "T[7,12]",  "T[7,13]",  "T[7,16]",  "T[8, 1]",  "T[8, 2]",  "T[8, 3]",  "T[8, 4]",  "T[8, 5]",  "T[8, 6]",  "T[8, 7]",  "T[8, 8]",  "T[8, 9]",  "T[8,10]",  "T[9, 1]",  "T[9, 3]",  "T[9, 5]",  "T[9, 6]",  "T[9, 7]"};
Mstd     = {"T[5, 2]",  "T[5, 4]",  "T[5, 8]",  "T[5, 9]",  "T[5,10]",  "T[6, 1]",  "T[6, 2]",  "T[6, 4]",  "T[6, 9]",  "T[6,10]",  "T[6,12]",  "T[6,13]",  "T[6,14]",  "T[6,15]",  "T[6,16]",  "T[6,17]",  "T[6,18]",  "T[6,19]",  "T[7, 1]",  "T[7, 2]",  "T[7, 6]",  "T[7,11]",  "T[7,14]",  "T[7,15]",  "T[9, 4]"};
Mepf     = {"Tf[5, 2]", "Tf[5, 4]", "Tf[5, 8]", "Tf[5, 9]", "Tf[5,10]", "Tf[6, 1]", "Tf[6, 2]", "Tf[6, 4]", "Tf[6, 9]", "Tf[6,10]", "Tf[6,12]", "Tf[6,13]", "Tf[6,14]", "Tf[6,15]", "Tf[6,16]", "Tf[6,17]", "Tf[6,18]", "Tf[6,19]", "Tf[7, 1]", "Tf[7, 2]", "Tf[7, 6]", "Tf[7,11]", "Tf[7,14]", "Tf[7,15]", "Tf[9, 4]"};
Mana     = {"T[4, 1]",  "T[5, 1]",  "T[5, 3]",  "T[5, 5]",  "T[5, 6]",  "T[5, 7]",  "T[6, 3]",  "T[6, 5]",  "T[6, 6]",  "T[6, 7]",  "T[6, 8]",  "T[7, 3]",  "T[7, 4]",  "T[9, 2]"};
Mfullana = {"Ta[4, 1]", "Ta[5, 1]", "Ta[5, 3]", "Ta[5, 5]", "Ta[5, 6]", "Ta[5, 7]", "Ta[6, 3]", "Ta[6, 5]", "Ta[6, 6]", "Ta[6, 7]", "Ta[6, 8]", "Ta[7, 3]"};
StdBasis = {"T[4, 1]",  "T[5, 1]",  "T[5, 2]",  "T[5, 3]",  "T[5, 4]",  "T[5, 5]",  "T[5, 6]",  "T[5, 7]",  "T[5, 8]",  "T[5, 9]",  "T[5,10]",  "T[6, 1]",  "T[6, 2]",  "T[6, 3]",  "T[6, 4]",  "T[6, 5]",  "T[6, 6]",  "T[6, 7]",  "T[6, 8]",  "T[6, 9]",  "T[6,10]",  "T[6,11]",  "T[6,12]",  "T[6,13]",  "T[6,14]",  "T[6,15]",  "T[6,16]",  "T[6,17]",  "T[6,18]",  "T[6,19]",  "T[7, 1]",  "T[7, 2]",  "T[7, 3]",  "T[7, 4]",  "T[7, 5]",  "T[7, 6]",  "T[7, 7]",  "T[7, 8]",  "T[7, 9]",  "T[7,10]",  "T[7,11]",  "T[7,12]",  "T[7,13]",  "T[7,14]",  "T[7,15]",  "T[7,16]",  "T[8, 1]",  "T[8, 2]",  "T[8, 3]",  "T[8, 4]",  "T[8, 5]",  "T[8, 6]",  "T[8, 7]",  "T[8, 8]",  "T[8, 9]",  "T[8,10]",  "T[9, 1]",  "T[9, 2]",  "T[9, 3]",  "T[9, 4]",  "T[9, 5]",  "T[9, 6]",  "T[9, 7]"};
EpfBasis = {"T[4, 1]",  "T[5, 1]",  "Tf[5, 2]", "T[5, 3]",  "Tf[5, 4]", "T[5, 5]",  "T[5, 6]",  "T[5, 7]",  "Tf[5, 8]", "Tf[5, 9]", "Tf[5,10]", "Tf[6, 1]", "Tf[6, 2]", "T[6, 3]",  "Tf[6, 4]", "T[6, 5]",  "T[6, 6]",  "T[6, 7]",  "T[6, 8]",  "Tf[6, 9]", "Tf[6,10]", "T[6,11]",  "Tf[6,12]", "Tf[6,13]", "Tf[6,14]", "Tf[6,15]", "Tf[6,16]", "Tf[6,17]", "Tf[6,18]", "Tf[6,19]", "Tf[7, 1]", "Tf[7, 2]", "T[7, 3]",  "T[7, 4]",  "T[7, 5]",  "Tf[7, 6]", "T[7, 7]",  "T[7, 8]",  "T[7, 9]",  "T[7,10]",  "Tf[7,11]", "T[7,12]",  "T[7,13]",  "Tf[7,14]", "Tf[7,15]", "T[7,16]",  "T[8, 1]",  "T[8, 2]",  "T[8, 3]",  "T[8, 4]",  "T[8, 5]",  "T[8, 6]",  "T[8, 7]",  "T[8, 8]",  "T[8, 9]",  "T[8,10]",  "T[9, 1]",  "T[9, 2]",  "T[9, 3]",  "Tf[9, 4]", "T[9, 5]",  "T[9, 6]",  "T[9, 7]"};

(* s2 = Im[PolyLog[2, E^(I/3*Pi)]]; *)

T[4, 1] = (
  +(1)*ep^-4
  +(4)*ep^-3
  +(10)*ep^-2
  +(20)*ep^-1
  +(35)
  +(56)*ep
  +(84)*ep^2
  +(120)*ep^3
  +(165)*ep^4
  +(220)*ep^5
  +O[ep]^6
);

T[5, 1] = (
  +(-1/3)*ep^-4
  +(-3/2)*ep^-3
  +(-43/12)*ep^-2
  +(-27/8-(8*Zeta[3])/3)*ep^-1
  +(797/48+(2*Pi^4)/15-12*Zeta[3])
  +(4029/32+(3*Pi^4)/5-(86*Zeta[3])/3-48*Zeta[5])*ep
  +(109349/192+(43*Pi^4)/30+(4*Pi^6)/21-(Zeta[3]*(81+32*Zeta[3]))/3-216*Zeta[5])*ep^2
  +(276261/128+(6*Pi^6)/7+(797*Zeta[3])/6-48*Zeta[3]^2+(Pi^4*(81+64*Zeta[3]))/60-516*Zeta[5]-664*Zeta[7])*ep^3
  +(5745485/768+(43*Pi^6)/21+(52*Pi^8)/225-(344*Zeta[3]^2)/3+(Pi^4*(-797+576*Zeta[3]))/120+Zeta[3]*(4029/4-384*Zeta[5])-486*Zeta[5]-2988*Zeta[7])*ep^4
  +(12607821/512+(26*Pi^8)/25+(Pi^6*(81+64*Zeta[3]))/42+2391*Zeta[5]+Pi^4*(-4029/80+(172*Zeta[3])/15+(96*Zeta[5])/5)-(Zeta[3]*(-328047+32*Zeta[3]*(243+64*Zeta[3])+124416*Zeta[5]))/72-7138*Zeta[7]-(80480*Zeta[9])/9)*ep^5
  +O[ep]^6
);

T[5, 2] = (
  +(-2)*ep^-4
  +(-29/3)*ep^-3
  +(-163/6)*ep^-2
  +(-601/12)*ep^-1
  +((-635-896*Zeta[3])/24)
  +(4549/16+(8*(17*Pi^4+60*Pi^2*Log[2]^2-60*Log[2]^4))/45-256*PolyLog[4,1/2]-(952*Zeta[3])/3)*ep
  +774.005466874027463354204247639953443434699271881776`51*ep^2
  +4459.06024625572626904843082548365792864849946480554`51*ep^3
  +14510.0356501941068565889763674855766391082040565046`51*ep^4
  +56830.0101966073686659496400867473486043301135134363`51*ep^5
  +O[ep]^6
);

T[5, 3] = (
  +(1/12)*ep^-3
  +(41/72)*ep^-2
  +(805/432)*ep^-1
  +(2093/2592+(8*Zeta[3])/3)
  +(-539483/15552-Pi^4/5+(164*Zeta[3])/9)*ep
  +(-28264243/93312-(41*Pi^4)/30+(1610*Zeta[3])/27+104*Zeta[5])*ep^2
  +(-1025549195/559872-(161*Pi^4)/36-(38*Pi^6)/63+(Zeta[3]*(2093+3456*Zeta[3]))/81+(2132*Zeta[5])/3)*ep^3
  +(-31881942499/3359232-(779*Pi^6)/189-(539483*Zeta[3])/486+(2624*Zeta[3]^2)/9-(Pi^4*(2093+6912*Zeta[3]))/1080+(20930*Zeta[5])/9+3084*Zeta[7])*ep^4
  +(-909890808827/20155392-(2185*Pi^6)/162-(343*Pi^8)/225+Pi^4*(539483/6480-(656*Zeta[3])/15)+(7*Zeta[3]*(-4037749+397440*Zeta[3]))/2916+(27209*Zeta[5])/27+3328*Zeta[3]*Zeta[5]+21074*Zeta[7])*ep^5
  +O[ep]^6
);

T[5, 4] = (
  +(-1)*ep^-4
  +(-9/2)*ep^-3
  +(-415/36)*ep^-2
  +(-4991/216)*ep^-1
  +(-89383/1296)
  +(-3367679/7776+(1792*Zeta[3])/9)*ep
  +(-137735095/46656-(256*(17*Pi^4+60*Pi^2*Log[2]^2-60*Log[2]^4))/135+(8192*PolyLog[4,1/2])/3+(68992*Zeta[3])/27)*ep^2
  +(-4908181487/279936+(128*(12*(385-72*Log[2])*Log[2]^4+60*Pi^2*Log[2]^2*(-77+24*Log[2])+17*Pi^4*(-77+72*Log[2])))/405+(315392*PolyLog[4,1/2])/9+32768*PolyLog[5,1/2]+(1602496*Zeta[3])/81-(87296*Zeta[5])/3)*ep^3
  -63850.1424496280567654758079403655066039675528545348`51*ep^4
  -235075.471229695213915120640095263103657996339489490`51*ep^5
  +O[ep]^6
);

T[5, 5] = (
  +(-1/144)*ep^-2
  +(-77/864)*ep^-1
  +((-3577-72*Pi^2)/5184)
  +((-7*(18695+792*Pi^2-864*Zeta[3]))/31104)*ep
  +((-20858285-72*Pi^2*(17885+264*Pi^2)+2328480*Zeta[3])/933120)*ep^2
  +((7*(-87082435-72*Pi^2*(93475+2904*Pi^2)+4320*(3577+72*Pi^2)*Zeta[3]))/5598720+(7*Zeta[5])/4)*ep^3
  +(-3363405529/6718464-(39347*Pi^4)/19440-(319*Pi^6)/11340+(49*(18695-432*Zeta[3])*Zeta[3])/7776+(7*Pi^2*(-595951+66528*Zeta[3]))/93312+(539*Zeta[5])/24)*ep^4
  +(-89200656545/40310784-(3509*Pi^6)/9720+(77*Pi^4*(-18695+864*Zeta[3]))/116640-(49*Zeta[3]*(-595951+33264*Zeta[3]))/46656+(25039*Zeta[5])/144+Pi^2*(-121915409/559872+(25039*Zeta[3])/648+(7*Zeta[5])/2)+(221*Zeta[7])/12)*ep^5
  +O[ep]^6
);

T[5, 6] = (
  +(1/12)*ep^-3
  +(17/24)*ep^-2
  +((179+4*Pi^2)/48)*ep^-1
  +((1513+68*Pi^2-80*Zeta[3])/96)
  +(3769/64+(179*Pi^2)/48+(23*Pi^4)/360-(85*Zeta[3])/12)*ep
  +(26123/128+(1513*Pi^2)/96+(391*Pi^4)/720-(5*(179+4*Pi^2)*Zeta[3])/24-(9*Zeta[5])/2)*ep^2
  +(517603/768+(3769*Pi^2)/64+(4117*Pi^4)/1440+(361*Pi^6)/7560-(85*(89+4*Pi^2)*Zeta[3])/48+(25*Zeta[3]^2)/6-(153*Zeta[5])/4)*ep^3
  +((6137*Pi^6)/15120+(23*Pi^4*(1513-80*Zeta[3]))/2880+Pi^2*(26123/128-(895*Zeta[3])/24-(9*Zeta[5])/2)+(3312041+80*Zeta[3]*(-11307+680*Zeta[3])-309312*Zeta[5]-42240*Zeta[7])/1536)*ep^4
  +(20750843/3072+(517603*Pi^2)/768+(86687*Pi^4)/1920+(64619*Pi^6)/30240+(5549*Pi^8)/151200-((1175535+90780*Pi^2+3128*Pi^4)*Zeta[3])/576+(25*(179+4*Pi^2)*Zeta[3]^2)/24-(153*(89+4*Pi^2)*Zeta[5])/16+45*Zeta[3]*Zeta[5]-(935*Zeta[7])/4)*ep^5
  +O[ep]^6
);

T[5, 7] = (
  +(-1/2)*ep^-4
  +(-5/2)*ep^-3
  +((-48-Pi^2)/6)*ep^-2
  +(-21-(5*Pi^2)/6+Zeta[3])*ep^-1
  +(-99/2-(8*Pi^2)/3-Pi^4/20+5*Zeta[3])
  +(-219/2-Pi^4/4+(Pi^2*(-21+Zeta[3]))/3+16*Zeta[3]+3*Zeta[5])*ep
  +(-233-(4*Pi^4)/5-(61*Pi^6)/3780-(-42+Zeta[3])*Zeta[3]+(Pi^2*(-99+10*Zeta[3]))/6+15*Zeta[5])*ep^2
  +(-484-(61*Pi^6)/756+(Pi^4*(-21+Zeta[3]))/10+(99-5*Zeta[3])*Zeta[3]+48*Zeta[5]+Pi^2*(-73/2+(16*Zeta[3])/3+Zeta[5])+9*Zeta[7])*ep^3
  +(-1981/2-(244*Pi^6)/945-(1261*Pi^8)/226800+(Pi^4*(-99+10*Zeta[3]))/20-(Pi^2*(233+(-42+Zeta[3])*Zeta[3]-15*Zeta[5]))/3+Zeta[3]*(219-16*Zeta[3]-6*Zeta[5])+126*Zeta[5]+45*Zeta[7])*ep^4
  +(-4017/2-(1261*Pi^8)/45360+(61*Pi^6*(-21+Zeta[3]))/1890+(2*Zeta[3]*(699+(-63+Zeta[3])*Zeta[3]-45*Zeta[5]))/3+297*Zeta[5]+(Pi^4*(-219+32*Zeta[3]+6*Zeta[5]))/20+144*Zeta[7]+(Pi^2*(-484+(99-5*Zeta[3])*Zeta[3]+48*Zeta[5]+9*Zeta[7]))/3+(85*Zeta[9])/3)*ep^5
  +O[ep]^6
);

T[5, 8] = (
  +(-1/4)*ep^-4
  +(-1)*ep^-3
  +(-97/48)*ep^-2
  +(-833/288)*ep^-1
  +(-26509/1728-Pi^4/15+(11*Sqrt[3]*s2)/2-4*s2^2)
  -69.6505441108505972304230794567366852451294478408788`51*ep
  -1026.45331711805694242567716850884259174534898355292`51*ep^2
  -4327.78556120107883777346998423861816859446980360435`51*ep^3
  -32631.2716737648748204746041955661927417100593304520`51*ep^4
  -121932.824454653195548790807105637775405715708473525`51*ep^5
  +O[ep]^6
);

T[5, 9] = (
  +(1/4)*ep^-3
  +(43/48)*ep^-2
  +(51/32)*ep^-1
  +(-161/192+(3*Sqrt[3]*s2)/2)
  -12.8201652257708830630523252025004250241866940303171`51*ep
  -109.565386271893281088571418372818911763660869123789`51*ep^2
  -376.660003044759726976545173478756733537037725480716`51*ep^3
  -2568.87433557725167824587651296555800886720384266300`51*ep^4
  -5759.42807235680031654293317419589161488863596136549`51*ep^5
  +O[ep]^6
);

T[5,10] = (
  +(-1)*ep^-4
  +(-19/4)*ep^-3
  +(-103/8)*ep^-2
  +(-341/16-3*Sqrt[3]*s2)*ep^-1
  -15.0423370505166027480948058347299032544022499112429`51
  +90.9453247668764379063683744128614170639466869823317`51*ep
  +786.412048166605103186941550664331466717690455280466`51*ep^2
  +3063.98510867182615803842655166426250231694118784335`51*ep^3
  +12328.7828669327819255448978864758705237188193073124`51*ep^4
  +39105.7901592708895392871664914561600411166146411504`51*ep^5
  +O[ep]^6
);

T[6, 1] = (
  +(3/2)*ep^-4
  +(19/2)*ep^-3
  +(67/2)*ep^-2
  +(127/2-3*Zeta[3])*ep^-1
  +(-237/2-Pi^4/20+97*Zeta[3])
  +(-3969/2-(559*Pi^4)/60-32*Pi^2*Log[2]^2+32*Log[2]^4+768*PolyLog[4,1/2]+1261*Zeta[3]+189*Zeta[5])*ep
  -7607.58935727331640672412576730534525017830137117901`51*ep^2
  -36043.9315389277686247446140343608059543249254100018`51*ep^3
  -199692.431290421928115148395110471119663343842784435`51*ep^4
  -784516.116603849522465732718633141346131423652952883`51*ep^5
  +O[ep]^6
);

T[6, 2] = (
  +(2/3)*ep^-4
  +(4)*ep^-3
  +(38/3)*ep^-2
  +((4*(11+4*Zeta[3]))/3)*ep^-1
  +(-118-(4*Pi^4)/15+88*Zeta[3])
  +((-374*Pi^4-960*Pi^2*Log[2]^2+60*(-867+16*Log[2]^4+384*PolyLog[4,1/2]+538*Zeta[3]+72*Zeta[5]))/45)*ep
  +((-2*(60*Pi^6+Pi^4*(12425-8456*Log[2])+672*Log[2]^4*(-55+Log[256])-1120*Pi^2*Log[2]^2*(-33+Log[256])-105*(-10469+8448*PolyLog[4,1/2]+6144*PolyLog[5,1/2]+4*Zeta[3]*(1619+8*Zeta[3])-3600*Zeta[5])))/315)*ep^2
  -28107.5715670960878830893417996758562446024390712307`51*ep^3
  -116811.754347202793190833719877891648542359421317239`51*ep^4
  -600956.150548688478765386244911222880217775654083441`51*ep^5
  +O[ep]^6
);

T[6, 3] = (
  +(1/6)*ep^-4
  +(5/6)*ep^-3
  +(3/2)*ep^-2
  +((-39+38*Zeta[3])/6)*ep^-1
  +(-157/2-(23*Pi^4)/60+(95*Zeta[3])/3)
  +(-989/2-(23*Pi^4)/12+57*Zeta[3]+215*Zeta[5])*ep
  +(-15191/6-(69*Pi^4)/20-(25*Pi^6)/21+(19*Zeta[3]*(-39+19*Zeta[3]))/3+1075*Zeta[5])*ep^2
  +(-70231/6-(125*Pi^6)/21-(23*Pi^4*(-39+38*Zeta[3]))/60+(19*Zeta[3]*(-471+95*Zeta[3]))/3+1935*Zeta[5]+6217*Zeta[7])*ep^3
  +(-102173/2-(75*Pi^6)/7-(11071*Pi^8)/3600+(23*Pi^4*(471-190*Zeta[3]))/60-8385*Zeta[5]+19*Zeta[3]*(-989+57*Zeta[3]+430*Zeta[5])+31085*Zeta[7])*ep^4
  +(-430877/2-(11071*Pi^8)/720+(25*Pi^6*(39-38*Zeta[3]))/21-101265*Zeta[5]-(23*Pi^4*(-989+114*Zeta[3]+430*Zeta[5]))/20+(19*Zeta[3]*(-45573+19*Zeta[3]*(-117+38*Zeta[3])+19350*Zeta[5]))/9+55953*Zeta[7]+(1632475*Zeta[9])/9)*ep^5
  +O[ep]^6
);

T[6, 4] = (
  +(1/4)*ep^-4
  +(3/2)*ep^-3
  +(9/2)*ep^-2
  +(2+(13*Zeta[3])/2)*ep^-1
  +(-147/2-(5*Pi^4)/8+30*Zeta[3])
  +(-624-(27*Pi^4)/10+45*Zeta[3]+(693*Zeta[5])/2)*ep
  -4609.64010449680326204100165437574570967564954728942`51*ep^2
  -16580.4910776068099036452615032781775337054142273815`51*ep^3
  -109190.181936620344099843280795304690763732312776615`51*ep^4
  -362516.971463352268504307890461417628856149005805110`51*ep^5
  +O[ep]^6
);

T[6, 5] = (
  +(-1/12)*ep^-3
  +(-11/12)*ep^-2
  +((-75-2*Pi^2)/12)*ep^-1
  +((-411-22*Pi^2+22*Zeta[3])/12)
  +(-665/4-(25*Pi^2)/2-(91*Pi^4)/360+(121*Zeta[3])/6)*ep
  +((-1001*Pi^4)/360+(Pi^2*(-411+22*Zeta[3]))/6+(-3001+550*Zeta[3]+70*Zeta[5])/4)*ep^2
  +((-455*Pi^4)/24-(1369*Pi^6)/3780+(Pi^2*(-1995+242*Zeta[3]))/6+(-38827+22*(411-11*Zeta[3])*Zeta[3]+2310*Zeta[5])/12)*ep^3
  +((-15059*Pi^6)/3780+(91*Pi^4*(-411+22*Zeta[3]))/360+Pi^2*(-3001/2+275*Zeta[3]+35*Zeta[5])+(-162731+22*(1995-121*Zeta[3])*Zeta[3]+15750*Zeta[5]+2358*Zeta[7])/12)*ep^4
  +((-6845*Pi^6)/252-(77989*Pi^8)/151200+(91*Pi^4*(-1995+242*Zeta[3]))/360+(Pi^2*(-38827+22*(411-11*Zeta[3])*Zeta[3]+2310*Zeta[5]))/6+(-223289+28770*Zeta[5]-22*Zeta[3]*(-3001+275*Zeta[3]+70*Zeta[5])+8646*Zeta[7])/4)*ep^5
  +O[ep]^6
);

T[6, 6] = (
  +(-1/16)*ep^-3
  +(-19/32)*ep^-2
  +((-221-8*Pi^2)/64)*ep^-1
  +((-2039-152*Pi^2+96*Zeta[3])/128)
  +(-16437/256-(221*Pi^2)/32-Pi^4/5+(57*Zeta[3])/8)*ep
  +((-607515-81560*Pi^2-4864*Pi^4+480*(221+8*Pi^2)*Zeta[3]+9600*Zeta[5])/2560)*ep^2
  +(-846877/1024-(221*Pi^4)/20-(1193*Pi^6)/3780+(3*(2039-48*Zeta[3])*Zeta[3])/32+(3*Pi^2*(-5479+608*Zeta[3]))/128+(285*Zeta[5])/8)*ep^3
  +(-5664343/2048-(22667*Pi^6)/7560+(9*(5479-304*Zeta[3])*Zeta[3])/64+(Pi^4*(-2039+96*Zeta[3]))/40+(3315*Zeta[5])/16+(3*Pi^2*(-40501+7072*Zeta[3]+640*Zeta[5]))/256+(87*Zeta[7])/4)*ep^4
  +((-263653*Pi^6)/15120-(19057*Pi^8)/37800+(3*Pi^4*(-5479+608*Zeta[3]))/80+Pi^2*(-846877/512+(6117*Zeta[3])/16-9*Zeta[3]^2+(285*Zeta[5])/4)+(-36784133+3914880*Zeta[5]-288*Zeta[3]*(-40501+3536*Zeta[3]+640*Zeta[5])+846336*Zeta[7])/4096)*ep^5
  +O[ep]^6
);

T[6, 7] = (
  +(1/3)*ep^-4
  +(2)*ep^-3
  +((23+Pi^2)/3)*ep^-2
  +(24+2*Pi^2-(4*Zeta[3])/3)*ep^-1
  +(67+(23*Pi^2)/3+(13*Pi^4)/45-8*Zeta[3])
  +(174+24*Pi^2+(26*Pi^4)/15-(4*(23+Pi^2)*Zeta[3])/3-4*Zeta[5])*ep
  +(1291/3+(299*Pi^4)/45+(242*Pi^6)/945+Pi^2*(67-8*Zeta[3])+(8*(-36+Zeta[3])*Zeta[3])/3-24*Zeta[5])*ep^2
  +((484*Pi^6)/315-(52*Pi^4*(-18+Zeta[3]))/45+Pi^2*(174-(92*Zeta[3])/3-4*Zeta[5])+4*(257+Zeta[3]*(-67+4*Zeta[3])-23*Zeta[5]-3*Zeta[7]))*ep^3
  +(7181/3+(5566*Pi^6)/945+(1093*Pi^8)/4725-(13*Pi^4*(-67+8*Zeta[3]))/15+(Pi^2*(1291+8*(-36+Zeta[3])*Zeta[3]-72*Zeta[5]))/3-288*Zeta[5]+(8*Zeta[3]*(-261+23*Zeta[3]+6*Zeta[5]))/3-72*Zeta[7])*ep^4
  +(5466+(2186*Pi^8)/1575-(968*Pi^6*(-18+Zeta[3]))/945-804*Zeta[5]-(26*Pi^4*(-261+46*Zeta[3]+6*Zeta[5]))/45+(4*Zeta[3]*(-3873-8*(-54+Zeta[3])*Zeta[3]+216*Zeta[5]))/9+4*Pi^2*(257+Zeta[3]*(-67+4*Zeta[3])-23*Zeta[5]-3*Zeta[7])-276*Zeta[7]-(340*Zeta[9])/9)*ep^5
  +O[ep]^6
);

T[6, 8] = (
  +(1/4)*ep^-4
  +(3/2)*ep^-3
  +(23/4+Pi^2/6)*ep^-2
  +(18+Pi^2-Zeta[3])*ep^-1
  +(201/4+(23*Pi^2)/6+(7*Pi^4)/90-6*Zeta[3])
  +(261/2+(7*Pi^4)/15-(2*Pi^2*(-18+Zeta[3]))/3-23*Zeta[3]-3*Zeta[5])*ep
  +(1291/4+(161*Pi^4)/90+(31*Pi^6)/945+Pi^2*(67/2-4*Zeta[3])+2*(-36+Zeta[3])*Zeta[3]-18*Zeta[5])*ep^2
  +((62*Pi^6)/315-(14*Pi^4*(-18+Zeta[3]))/45+Pi^2*(87-(46*Zeta[3])/3-2*Zeta[5])+3*(257+Zeta[3]*(-67+4*Zeta[3])-23*Zeta[5]-3*Zeta[7]))*ep^3
  +(7181/4+(713*Pi^6)/945+(127*Pi^8)/9450-(7*Pi^4*(-67+8*Zeta[3]))/30+(Pi^2*(1291+8*(-36+Zeta[3])*Zeta[3]-72*Zeta[5]))/6-216*Zeta[5]+2*Zeta[3]*(-261+23*Zeta[3]+6*Zeta[5])-54*Zeta[7])*ep^4
  +(8199/2+(127*Pi^8)/1575-(124*Pi^6*(-18+Zeta[3]))/945-603*Zeta[5]-(7*Pi^4*(-261+46*Zeta[3]+6*Zeta[5]))/45+(Zeta[3]*(-3873-8*(-54+Zeta[3])*Zeta[3]+216*Zeta[5]))/3+2*Pi^2*(257+Zeta[3]*(-67+4*Zeta[3])-23*Zeta[5]-3*Zeta[7])-207*Zeta[7]-(85*Zeta[9])/3)*ep^5
  +O[ep]^6
);

T[6, 9] = (
  +(1)*ep^-4
  +(25/4)*ep^-3
  +(257/12)*ep^-2
  +(433/12+6*Sqrt[3]*s2-2*Zeta[3])*ep^-1
  -53.3431573078311113878444050309141441704781215668838`51
  -740.732476883867657524902616106453155974798891291779`51*ep
  -6741.86165372823991163645454267615958732687467604314`51*ep^2
  -29171.1801952361837741528205756035238340511898969484`51*ep^3
  -170737.762505082563952818195086853029518923443551139`51*ep^4
  -640207.942056859484083975040325762869424889351199004`51*ep^5
  +O[ep]^6
);

T[6,10] = (
  +(-5/12)*ep^-4
  +(-7/4)*ep^-3
  +(-43/12)*ep^-2
  +(13/4-2*Sqrt[3]*s2+Zeta[3]/2)*ep^-1
  +49.1471677206145743293211843554330513425580705902593`51
  +298.047852285403108719760320923150121416370488257760`51*ep
  +1878.41451599274340307816701182114565581014315742770`51*ep^2
  +7332.54241105293400661103386319347693434012791991672`51*ep^3
  +40241.5145930986868154666035802713583681613265821298`51*ep^4
  +139863.910164380149573650856671369573977634414383713`51*ep^5
  +O[ep]^6
);

T[6,11] = (
  +(7/12)*ep^-4
  +(43/12)*ep^-3
  +(47/4)*ep^-2
  +(187/12+6*Sqrt[3]*s2+(7*Zeta[3])/6)*ep^-1
  -87.9261424095009363612687893654620107473927960195902`51
  -569.406457016489862087645486512240183239542065549736`51*ep
  -5855.66085223584587057683649755963886313717162935362`51*ep^2
  -22160.5520141217692429128547995974118424901910011306`51*ep^3
  -142927.585167354942005841699030740423215261989147644`51*ep^4
  -489366.566721874883382713601349244423443225591005301`51*ep^5
  +O[ep]^6
);

T[6,12] = (
  +(-1/3)*ep^-4
  +(-4/3)*ep^-3
  +(-7/3)*ep^-2
  +(16/3-4*Sqrt[3]*s2)*ep^-1
  +65.2891346352764051182112101341904605154404461333644`51
  +235.378897804749678044597587462774656737157101747014`51*ep
  +2150.74945309205496063232566551798107820292227024787`51*ep^2
  +6222.87086902451243107690707730700055230770221858007`51*ep^3
  +44770.6542262683613422489020019241678350088606452694`51*ep^4
  +121556.575316814885572650674390186850862381110081173`51*ep^5
  +O[ep]^6
);

T[6,13] = (
  +(1/3)*ep^-4
  +(23/12)*ep^-3
  +(65/12)*ep^-2
  +((13-Pi^2+36*Sqrt[3]*s2+32*Zeta[3])/12)*ep^-1
  -76.1116203247541030218866146841167132872551134655119`51
  -579.527720797210645153228527443070648586673034785390`51*ep
  -4056.41938134389585102475465812870949463404507745492`51*ep^2
  -19103.3405264140566016820578331943647873538263031573`51*ep^3
  -97977.6398166848758327755615165768204653977860963905`51*ep^4
  -410630.835681254547240771070277894707075844502423836`51*ep^5
  +O[ep]^6
);

T[6,14] = (
  +(-1/4)*ep^-4
  +(-11/12)*ep^-3
  +(-13/12)*ep^-2
  +(89/12-2*Sqrt[3]*s2-(19*Zeta[3])/6)*ep^-1
  +66.6447233294059950135392875885872818918490991587217`51
  +261.073616883964081237073348445916210693974952247406`51*ep
  +2031.08779600320272948973995231519991824419025879602`51*ep^2
  +6891.05242086483155421426678030349172605988636207033`51*ep^3
  +41668.2229555379844251011754342412682628891095870607`51*ep^4
  +135551.585073576832754485113344699652902365571106400`51*ep^5
  +O[ep]^6
);

T[6,15] = (
  +(1/12)*ep^-4
  +(5/12)*ep^-3
  +(7/12)*ep^-2
  +((-67-2*Pi^2+26*Zeta[3])/12)*ep^-1
  +(-235/4-2*Pi^2-(13*Pi^4)/72+(101*Zeta[3])/6)
  +((-501*Pi^4+60*Pi^2*(-89+16*Zeta[3])+30*(-4491+1046*Zeta[3]+726*Zeta[5]))/360)*ep
  -2848.65364899561089848554539142820121003334361464887`51*ep^2
  -13086.0433346924542908439247754691849713271184772043`51*ep^3
  -66103.2927374676695704447815055258476452325092458040`51*ep^4
  -274798.301615962817291679989927277620778418503167198`51*ep^5
  +O[ep]^6
);

T[6,16] = (
  +(7/8)*ep^-4
  +(85/16)*ep^-3
  +((601+8*Pi^2)/32)*ep^-2
  +(2747/64+(35*Pi^2)/24-(3*Zeta[3])/2)*ep^-1
  +(2329/128+(245*Pi^2)/48+(13*Pi^4)/120+(17*Sqrt[3]*s2)/2+2*s2^2+(343*Zeta[3])/12)
  -171.318541252613251780835728791802769743004355180802`51*ep
  -1484.23423452676423430283578012278823316465768368487`51*ep^2
  -15756.8227036156414833491686451093387252075407345472`51*ep^3
  -61784.9262105249516867263782178284342470787649723502`51*ep^4
  -379234.671466128300538566659065130195616335047078208`51*ep^5
  +O[ep]^6
);

T[6,17] = (
  +(5/8)*ep^-4
  +(33/8)*ep^-3
  +(775/48+Pi^2/4)*ep^-2
  +(1389/32+(3*Pi^2)/2+3*Sqrt[3]*s2-(3*Zeta[3])/2)*ep^-1
  +89.1253843662396089100751397605580246619255586179741`51
  +271.853771423526536244060281201618507436482160881522`51*ep
  -2285.75335829698084263732602170670403858153054650404`51*ep^2
  -5143.23119756650402939123588964287773512593212038633`51*ep^3
  -76626.3302154978090854065032397049297354652037994107`51*ep^4
  -176028.154555601635614737265459017339129994011889373`51*ep^5
  +O[ep]^6
);

T[6,18] = (
  +(1/2)*ep^-4
  +(47/16)*ep^-3
  +(317/32+Pi^2/6)*ep^-2
  +(1315/64+(23*Pi^2)/24+3*Sqrt[3]*s2-Zeta[3])*ep^-1
  +46.5951023350029333262768262200087593931078424251702`51
  -82.4458756013779331454227175918532458320048394468667`51*ep
  -1716.05600327076619038320639527224996002910110174859`51*ep^2
  -9542.34374186199562911508547710953924345386356556808`51*ep^3
  -54334.1548724053092390424686707923952435094025222036`51*ep^4
  -234597.610608431240251139948227857914331916668095494`51*ep^5
  +O[ep]^6
);

T[6,19] = (
  +(5/24)*ep^-4
  +(55/48)*ep^-3
  +((331+8*Pi^2)/96)*ep^-2
  +((977+88*Pi^2+416*Zeta[3])/192)*ep^-1
  +(-1831/128+(55*Pi^2)/48-(2*Pi^4)/15+(143*Zeta[3])/12)
  +((-8*Pi^2*(535+352*Pi^2-1120*Zeta[3])+25*(-25491+6368*Zeta[3]+8832*Zeta[5]))/3840)*ep
  -1306.48849778486874986108799144773529718560066405195`51*ep^2
  -5997.65155434525980602082748734519239317681020886212`51*ep^3
  -35213.4371437164874838456913460567631842709841245385`51*ep^4
  -138627.971084958086694167619293369246707690280964378`51*ep^5
  +O[ep]^6
);

T[7, 1] = (
  +(-1/6)*ep^-4
  +(-3/2)*ep^-3
  +(-26/3-Zeta[3])*ep^-2
  +(-41-Pi^4/60-(10*Zeta[3])/3)*ep^-1
  +(-1039/6-(7*Pi^4)/30+3*Zeta[3]+53*Zeta[5])
  -929.858212294016976382457232976107077779091868024097`51*ep
  -1698.90639250653423244336742929239614618805875415910`51*ep^2
  -16655.8998994692523769449673612842926626244934637372`51*ep^3
  -21751.2175089760729706512346573100940966374203326102`51*ep^4
  -283551.819410167943346104338229079644742630341423269`51*ep^5
  +O[ep]^6
);

T[7, 2] = (
  +(-1/6)*ep^-4
  +(-3/2)*ep^-3
  +(-26/3-Zeta[3]/2)*ep^-2
  +(-41-Pi^4/120+Zeta[3]/6)*ep^-1
  +(-1039/6-(29*Pi^4)/120+(17*Zeta[3])/2+(43*Zeta[5])/2)
  -758.669486731967257046185950411227644240941810349848`51*ep
  -2529.44744210570103969358709741446931936388651107443`51*ep^2
  -12777.5114312429869652348444531738491309289357459403`51*ep^3
  -38865.8166684274098583605018809462774767419496261688`51*ep^4
  -210003.011997666017499226365581434745569104737812603`51*ep^5
  +O[ep]^6
);

T[7, 3] = (
  +(-1/4)*ep^-4
  +(-7/4)*ep^-3
  +((-31-2*Pi^2)/4)*ep^-2
  +((-111-14*Pi^2+6*Zeta[3])/4)*ep^-1
  +((-3510-620*Pi^2-33*Pi^4+420*Zeta[3])/40)
  +((-3*(3410+740*Pi^2+77*Pi^4-20*(31+2*Pi^2)*Zeta[3]-60*Zeta[5]))/40)*ep
  +((-1023*Pi^4)/40-(5051*Pi^6)/3780+Pi^2*(-351/2+21*Zeta[3])+(-2815-18*(-37+Zeta[3])*Zeta[3]+126*Zeta[5])/4)*ep^2
  +((-5051*Pi^6)/540+(99*Pi^4*(-37+2*Zeta[3]))/40+(3*Pi^2*(-341+62*Zeta[3]+6*Zeta[5]))/2+(-7423+18*(117-7*Zeta[3])*Zeta[3]+558*Zeta[5]+54*Zeta[7])/4)*ep^3
  +((-156581*Pi^6)/3780-(65459*Pi^8)/30240+(99*Pi^4*(-117+14*Zeta[3]))/40+Pi^2*(-2815/2-9*(-37+Zeta[3])*Zeta[3]+63*Zeta[5])+(-18943+1998*Zeta[5]-18*Zeta[3]*(-341+31*Zeta[3]+6*Zeta[5])+378*Zeta[7])/4)*ep^4
  +((-65459*Pi^8)/4320+(5051*Pi^6*(-37+2*Zeta[3]))/1260+(99*Pi^4*(-341+62*Zeta[3]+6*Zeta[5]))/40+Pi^2*(-7423/2+9*(117-7*Zeta[3])*Zeta[3]+279*Zeta[5]+27*Zeta[7])+(3*(-15701+2*Zeta[3]*(2815-333*Zeta[3]+6*Zeta[3]^2-126*Zeta[5])+2106*Zeta[5]+558*Zeta[7]))/4+(85*Zeta[9])/2)*ep^5
  +O[ep]^6
);

T[7, 4] = (
  +(-1/12)*ep^-4
  +(-2/3)*ep^-3
  +((-39-2*Pi^2)/12)*ep^-2
  +((-71-8*Pi^2-4*Zeta[3])/6)*ep^-1
  +((-5655-1170*Pi^2-53*Pi^4-960*Zeta[3])/180)
  +(-33-(106*Pi^4)/45-32*Zeta[3]-(Pi^2*(71+4*Zeta[3]))/3-30*Zeta[5])*ep
  +(3557/12-(139*Pi^4)/12-(1067*Pi^6)/1890-(Pi^2*(377+64*Zeta[3]))/6+(Zeta[3]*(-554+73*Zeta[3]))/3-240*Zeta[5])*ep^2
  +O[ep]^3
);

T[7, 5] = (
  +(-1/8)*ep^-4
  +(-5/4)*ep^-3
  +(-65/8-Pi^2/12-Zeta[3])*ep^-2
  +(-(Pi^2*(20+Pi^2))/60+6*Sqrt[3]*s2-(7*(25+2*Zeta[3]))/4)*ep^-1
  -156.049584217013908703316644017514331570234624101191`51
  -835.421750060231317129023644997856214596933401364486`51*ep
  -2120.36177575634886782298722795385479202254270566468`51*ep^2
  -14808.2235938865963602830098820547190080295069712873`51*ep^3
  -29612.2181184316730879620573005673680545959915664177`51*ep^4
  -250634.846174409858661438600401104662221407239068637`51*ep^5
  +O[ep]^6
);

T[7, 6] = (
  +(-1/12)*ep^-4
  +(-5/6)*ep^-3
  +(-65/12-Pi^2/6-Zeta[3])*ep^-2
  +(-175/6-Pi^2-Pi^4/60+(4*Zeta[3])/3)*ep^-1
  +(-567/4-(7*Pi^4)/15-4*s2*(2*Sqrt[3]+s2)+(119*Zeta[3])/3+(Pi^2*(-21+4*Zeta[3]))/6+(89*Zeta[5])/3)
  -795.845747921992172541378362705390925979643448630938`51*ep
  -2275.19417932292514912539410182373663629703908558971`51*ep^2
  -13954.9926720189379646826119246076008930914330603312`51*ep^3
  -33402.6350498761995592224304839342514651987201430382`51*ep^4
  -233861.011428556045093870531508261144138372439654607`51*ep^5
  +O[ep]^6
);

T[7, 7] = (
  +(-1/8)*ep^-4
  +(-13/12)*ep^-3
  +((-143-2*Pi^2-24*Zeta[3])/24)*ep^-2
  +((-1585-Pi^2*(50+Pi^2)-390*Zeta[3])/60)*ep^-1
  -132.927909227859234562180140709595508846351935604449`51
  -902.030459371319167831276699395406228702636845534783`51*ep
  -1807.19566262899838014518254691893272656274869334805`51*ep^2
  -16006.6445345049084243422219691236555956831640164622`51*ep^3
  -24699.2696603309058067048927288212817342739539069764`51*ep^4
  -270263.055586806460208406029078469699507232709731666`51*ep^5
  +O[ep]^6
);

T[7, 8] = (
  +(-1/6)*ep^-4
  +(-3/2)*ep^-3
  +(-26/3-Zeta[3]/2)*ep^-2
  +(-41-Pi^4/120+Zeta[3]/6)*ep^-1
  -162.976245715464622145364548222119432508271661334372`51
  -764.817094193383828292354323113322264422954387069122`51*ep
  -2493.08813759412388791047787495152565072349971695351`51*ep^2
  -12948.2074113991637243025468651240872264533079113639`51*ep^3
  -38080.9422548261554852729826583236440069036961982664`51*ep^4
  -213427.317295268612465071650676274936043883153799733`51*ep^5
  +O[ep]^6
);

T[7, 9] = (
  +(-1/8)*ep^-4
  +(-5/4)*ep^-3
  +((-2*Pi^2-3*(65+4*Zeta[3]))/24)*ep^-2
  +(-175/4-(Pi^2*(40+Pi^2))/120+6*Sqrt[3]*s2)*ep^-1
  -180.629629484096466921122165628045571803328138558462`51
  -676.134169798872450171140948157526349869316884831680`51*ep
  -2881.63688000058492263949888583848283707409747138069`51*ep^2
  -11256.7214766441906930294998522061124206744264207114`51*ep^3
  -45229.4474118318170962413127328228664592904872525285`51*ep^4
  -183615.889345881491682532270842342286637710809154608`51*ep^5
  +O[ep]^6
);

T[7,10] = (
  +(-1/12)*ep^-4
  +(-5/6)*ep^-3
  +((-65-2*Pi^2-6*Zeta[3])/12)*ep^-2
  +(-175/6-Pi^2-Pi^4/120+(29*Zeta[3])/6)*ep^-1
  -169.68281024338906010112803809`29
  -665.57628415369448105619377231`29*ep
  -2850.6848464132727472947198342`29*ep^2
  -11181.718497095853379575724011`29*ep^3
  -45469.633236247619900460707251`29*ep^4
  +O[ep]^5
);

T[7,11] = (
  +(Zeta[3]/2)*ep^-2
  +(Pi^4/120+(3*Zeta[3])/2)*ep^-1
  -6.572315129519925781565283243`28
  +71.30962362014960763888308696`28*ep
  -318.6084779218926893927388846`28*ep^2
  +1602.676470257724982973999546`28*ep^3
  -7082.780107764081057021038489`28*ep^4
  +O[ep]^5
);

T[7,12] = (
  +(-1/12)*ep^-4
  +(-2/3)*ep^-3
  +(-13/4-Pi^2/6-Zeta[3])*ep^-2
  +((-710-80*Pi^2-Pi^4-460*Zeta[3])/60)*ep^-1
  -135.07366864912844484857884400`29
  -782.87577036928136782322700285`29*ep
  -2183.7114630194528880793176923`29*ep^2
  +O[ep]^3
);

T[7,13] = (
  +(1/12)*ep^-4
  +(1/3)*ep^-3
  +((7+2*Pi^2)/12)*ep^-2
  +((-7+4*Pi^2+4*Zeta[3])/6)*ep^-1
  +2.6657462440150049239997240951`29
  +102.86025722858080565176524103`29*ep
  -350.37032650986523973158094743`29*ep^2
  +O[ep]^3
);

T[7,14] = (
  +(-1/6)*ep^-4
  +(-3/2)*ep^-3
  +(-26/3)*ep^-2
  +(-41+(11*Zeta[3])/3)*ep^-1
  +(-1039/6-Pi^2/6-(107*Pi^4)/360-2*s2*(Sqrt[3]+5*s2)+(88*Zeta[3])/3)
  -620.899055909975738516131952589545183485195239852635`51*ep
  -3152.26135128742202236124626320684984343505741826017`51*ep^2
  -9799.32345191744160458636896689769344726956162669185`51*ep^3
  -51868.5239490901307210153664262839525771395359693890`51*ep^4
  -153962.632662747997028853728610167591277850604211675`51*ep^5
  +O[ep]^6
);

T[7,15] = (
  +(-1/8)*ep^-4
  +(-13/12)*ep^-3
  +((-143-2*Pi^2-12*Zeta[3])/24)*ep^-2
  +(-317/12-(Pi^2*(100+Pi^2))/120-3*Zeta[3])*ep^-1
  +((3*Pi^4+10*Pi^2*(-67+4*Zeta[3])-5*(2455+48*(Sqrt[3]-4*s2)*s2+260*Zeta[3]-236*Zeta[5]))/120)
  -772.326304945612408521125582223746870415454230297447`51*ep
  -2380.54278631017463578393120978523351441128237277803`51*ep^2
  -13248.1221180512548321610604536200317660965438346575`51*ep^3
  -36698.4781126567203746779252235629935102098568501546`51*ep^4
  -218523.209479111011385577827885740307876954923621396`51*ep^5
  +O[ep]^6
);

T[7,16] = (
  +(-2*Zeta[3])*ep^-2
  +(Pi^4/24+6*s2^2-8*Zeta[3])*ep^-1
  -22.1816263809632850054323152554315622943194170350060`51
  +30.9206490544864176839986405789579378893338357050830`51*ep
  -170.922690887707549260516228781859603075752767258103`51*ep^2
  +406.817552793235719075535960219234386388596335675407`51*ep^3
  -1435.41175550174209817970147238332294930617652441141`51*ep^4
  +4083.62192628794169653335489322322732268110319991762`51*ep^5
  +O[ep]^6
);

T[8, 1] = (
  +(5*Zeta[5])*ep^-1
  -22.5376034664572420845745563240437167867299150237220`51
  +139.691001197759679495462990254105927707407615819754`51*ep
  -625.606683387386379323122696026107572459578150052571`51*ep^2
  +2893.63704617442450174181954152246120985053321265976`51*ep^3
  -12453.6923484070619614497904070996213676331975501436`51*ep^4
  +53077.1123829600559498689276663763067610864397505338`51*ep^5
  +O[ep]^6
);

T[8, 2] = (
  +(5*Zeta[5])*ep^-1
  -20.8585976739540823192137169923590025170233079442547`51
  +131.757549396226936476090999955288889156274264880900`51*ep
  -580.487243803379244598846616689356519101143295986134`51*ep^2
  +2681.90756289415256901070042568239562056161835574962`51*ep^3
  -11507.8282524391926755587904461935426353147554901259`51*ep^4
  +O[ep]^5
);

T[8, 3] = (
  +(5*Zeta[5])*ep^-1
  -19.1118657059490732311103013898`30.
  +124.087397612137003608129119969`29.9*ep
  -534.246031506076674961968630976`29.9*ep^2
  +2473.04910756305895256706669346`29.8*ep^3
  +O[ep]^4
);

T[8, 4] = (
  +(5*Zeta[5])*ep^-1
  -19.0793750927079960314257405875623941525789104197843`51
  +124.195470293149892759222880468915672769830502954819`51*ep
  -533.95135196587595436445960426058731048205531150251`50*ep^2
  +2473.65262820126583714929043301059337082260959851`48*ep^3
  -10553.010275894406552790616481496857780968660270`47*ep^4
  +44849.59018769062572562708741799166955140730510`46*ep^5
  +O[ep]^6
);

T[8, 5] = (
  +(5*Zeta[5])*ep^-1
  -18.7699726985410818615637417243`30.
  +125.135163171638095080201914343`29.9*ep
  -530.195134968427619086820803930`29.9*ep^2
  +2481.75523870799122066812875155`29.8*ep^3
  +O[ep]^4
);

T[8, 6] = (
  +(5*Zeta[5])*ep^-1
  -7.85876279518922242360822990783004542937882100674500`51
  +111.474608783473449039353499987996944146578608976659`51*ep
  -315.715241817850688287306967430566714060665077044280`51*ep^2
  +1919.83292772839679713249994717844876805487305298533`51*ep^3
  -7046.39837794066951989835293569918154711673923387054`51*ep^4
  +32505.2820150069394250113591520370747074441854612974`51*ep^5
  +O[ep]^6
);

T[8, 7] = (
  +(5*Zeta[5])*ep^-1
  -13.7206302625920162070487342658557388782474967767706`51.
  +118.872128215995017633768624758569462618371336937765`51*ep
  -444.315939169987245888368415779776502167496370107304`51*ep^2
  +2248.61825955351812510018837480736599668191219456295`51*ep^3
  -9231.94617741275567691090539582690062244307382318077`51*ep^4
  +40129.0862267827692308107936676652484481609245345321`51*ep^5
  +O[ep]^6
);

T[8, 8] = (
  +(5*Zeta[5])*ep^-1
  +0.37484419192685127491080516587109756699794227132128`51
  +124.626355435305786046350805293568847947516422168178`51*ep
  -139.607561640075995777546540864167505584519986790149`51*ep^2
  +2005.70645205831280590173139698951241558666286126275`51*ep^3
  -4733.46716382038114637049761573464137229252675381718`51*ep^4
  +30617.9126908797386004266727157562061757703730503775`51*ep^5
  +O[ep]^6
);

T[8, 9] = (
  +(5*Zeta[5])*ep^-1
  +7.30445968508659934705130294982328536505673135057504`51
  +124.386639963921089710170469017755092329624975818697`51*ep
  +24.7563113389735947896296575370552341503182665984131`51*ep^2
  +1813.21342231067481693321346202497779047941437636689`51*ep^3
  -2086.67145447038092994514592376266295688002312310178`51*ep^4
  +24386.3391667760896788727645914721248654230367448070`51*ep^5
  +O[ep]^6
);

T[8,10] = (
  +((3*Zeta[3])/2)*ep^-2
  +(-Pi^4/20-6*s2^2+(15*Zeta[3])/2)*ep^-1
  +34.5454770641910364260696584117360443212020763884804`51
  -91.0420216826153216549385757172524696346537853406570`51*ep
  +550.434498485040019731708881873691740291433975212660`51*ep^2
  -2008.49820920768977583338914584426701518133569091105`51*ep^3
  +9045.06764562757780161516231952869079163872372129316`51*ep^4
  -36073.9750971761986186794428165377704378284936790190`51*ep^5
  +O[ep]^6
);

T[9, 1] = (
  +((-53*Pi^4*Log[2]-80*Pi^2*Log[2]^3+48*Log[2]^5)/15-384*PolyLog[5,1/2]+(873*Zeta[5])/2)
  -5.5459654251191013`17*ep
  +30.77339207827551`16*ep^2
  -119.8656758489377`16*ep^3
  +506.1570434525610`16*ep^4
  -2031.403298527746`16*ep^5
  +O[ep]^6
);

T[9, 2] = (
  +(-5*Zeta[5])*ep^-1
  +((-5*Pi^6)/378-17*Zeta[3]^2-15*Zeta[5])
  +((-5*Pi^6)/126-(17*Pi^4*Zeta[3])/30-51*Zeta[3]^2-5*(7+2*Pi^2)*Zeta[5]-(225*Zeta[7])/2)*ep
  +O[ep]^2
);

T[9, 3] = (
  +2.42695639537700735`18
  -4.01554669961524192`18*ep
  +35.7118868533108389`17.8*ep^2
  -111.835752497273400`17.8*ep^3
  +O[ep]^4
);

T[9, 4] = (
  +0.473611472272364450`18
  +1.09585342206826990`18*ep
  +3.81952404214692972`17.7*ep^2
  +5.98282857632256558`17.6*ep^3
  +O[ep]^4
);

T[9, 5] = (
  +3.71140264536682392682628965373`30
  +2.11520599545877264756135261804`30*ep
  +59.7799198446739180157589609808`29.8*ep^2
  -42.1984553244037974874378020975`29.8*ep^3
  +O[ep]^4
);

T[9, 6] = (
  +6.72847056008568105547188977521`30
  +26.0876465999666155389659770717`30*ep
  +192.511937028270755521957567384`29.9*ep^2
  +538.674291843537418086527085151`29.8*ep^3
  +O[ep]^4
);

T[9, 7] = (
  +3.4497511317390349922288`23
  -6.3127694885459812824115`23*ep
  +52.319545026695977054431`22.8*ep^2
  -170.04678406774887449161`22.8*ep^3
  +O[ep]^4
);

T[5,11] = (
  +(-1/4)*ep^-4
  +(-9/8)*ep^-3
  +(-5/2)*ep^-2
  +(15/8-3*Sqrt[3]*s2)*ep^-1
  +38.0704533519002103608024213485152329518935853379088`51
  +267.065844145517258452660246302380022925782281076783`51*ep
  +1686.90593105977560007153674285737381182591850405985`51*ep^2
  +8096.92980953150627103098678549127739733507884419312`51*ep^3
  +38379.7470702103703365291772316398184133613715150439`51*ep^4
  +167951.409808094415211034120864238042955059580672281`51*ep^5
  +O[ep]^6
);

T[7,17] = (
  +(-1/8)*ep^-4
  +(-3/4)*ep^-3
  +((-75+2*Pi^2)/24)*ep^-2
  +((5*Pi^2)/6-(3*(15+2*Zeta[3]))/4)*ep^-1
  +((212*Pi^2+Pi^4-2*(903+24*(4*Sqrt[3]-7*s2)*s2+464*Zeta[3]))/48)
  -1.4853010853959284159243952694`28.9*ep
  -121.68443497337153030056792706`29*ep^2
  +318.90697286781027948586655666`29*ep^3
  -1657.9303779095409850678955260`28.8*ep^4
  +O[ep]^5
);

(* ====== ep-finite ======*)

Tf[5, 2] = (
  +((3*Zeta[3])/2)*ep^-2
  +((-91*Pi^4)/360-(4*Pi^2*Log[2]^2)/3+(4*Log[2]^4)/3+32*PolyLog[4,1/2]+(15*Zeta[3])/2)*ep^-1
  +27.300681`8
  +O[ep]^1
);

Tf[5, 4] = (
  +((5*Zeta[5])/4)*ep^-1
  -6.7622402385`11
  +O[ep]^1
);

Tf[5, 8] = (
  +((3*Zeta[3])/2)*ep^-2
  +((-91*Pi^4)/360-(4*Pi^2*Log[2]^2)/3+(4*Log[2]^4)/3+32*PolyLog[4,1/2]+(15*Zeta[3])/2)*ep^-1
  +37.52435523`9
  +O[ep]^1
);

Tf[5, 9] = (
  +6.8481671`8
  +O[ep]^1
);

Tf[5,10] = (
  +((-3*Zeta[3])/2+(5*Zeta[5])/2)*ep^-1
  -3.272794`7
  +O[ep]^1
);

Tf[6, 1] = (
  +(5*Zeta[5])*ep^-1
  -29.7034624278`12
  +O[ep]^1
);

Tf[6, 2] = (
  +2.4403458233508`14
  +O[ep]^1
);

Tf[6, 4] = (
  +(5*Zeta[5])*ep^-1
  -18.0262459787292`15
  +O[ep]^1
);

Tf[6, 9] = (
  +3.2095107`8
  +O[ep]^1
);

Tf[6,10] = (
  +((-5*Zeta[5])/4)*ep^-1
  +5.68333946`9
  +O[ep]^1
);

Tf[6,12] = (
  +((-3*Zeta[3])/2+(5*Zeta[5])/2)*ep^-1
  -3.67204`6
  +O[ep]^1
);

Tf[6,13] = (
  +5.1732238`8
  +O[ep]^1
);

Tf[6,14] = (
  +(5*Zeta[5])*ep^-1
  -23.2862341`9
  +O[ep]^1
);

Tf[6,15] = (
  +8.816973`7
  +O[ep]^1
);

Tf[6,16] = (
  +9.5955369`8
  +O[ep]^1
);

Tf[6,17] = (
  +((-5*Zeta[5])/4)*ep^-1
  -4.872849`7
  +O[ep]^1
);

Tf[6,18] = (
  +((-3*Zeta[3])/2+(5*Zeta[5])/2)*ep^-1
  -1.589373`7
  +O[ep]^1
);

Tf[6,19] = (
  +((-3*Zeta[3])/2+(5*Zeta[5])/2)*ep^-1
  -0.80223`5
  +O[ep]^1
);

Tf[7, 1] = (
  +(5*Zeta[5])*ep^-1
  -5.198314`7
  +O[ep]^1
);

Tf[7, 2] = (
  +(5*Zeta[5])*ep^-1
  -26.49647940444747`16
  +O[ep]^1
);

Tf[7, 6] = (
  +(5*Zeta[5])*ep^-1
  -18.9284641`9
  +O[ep]^1
);

Tf[7,11] = (
  +5.126697`7
  +O[ep]^1
);

Tf[7,14] = (
  +((-3*Zeta[3])/2+(5*Zeta[5])/2)*ep^-1
  +0.0114944`6
  +O[ep]^1
);

Tf[7,15] = (
  +(5*Zeta[5])*ep^-1
  -24.8172810`9
  +O[ep]^1
);

Tf[9, 4] = (
  +((-5*Zeta[5])/4)*ep^-1
  +6.2801043`8
  +O[ep]^1
);

(* ====== full analytical results ======= *)

Ta[4, 1] = 256/((-4+d)^4*(-2+d)^4);
Ta[5, 1] = (12*Gamma[3-(3*d)/2]*Gamma[3-d]^2*Gamma[-2+d/2])/((-4+d)*(-2+d)*Gamma[6-2*d]*Gamma[3-d/2]^2);
Ta[5, 3] = (-2*Gamma[5-2*d]*Gamma[4-(3*d)/2]^2*Gamma[2-d]*Gamma[-1+d/2]^2)/(Gamma[8-3*d]*Gamma[3-d/2]^4);
Ta[5, 5] = (-3*Gamma[5-2*d]*Gamma[3-(3*d)/2]*Gamma[-1+d/2]^3)/Gamma[3-d/2]^4;
Ta[5, 6] = (-8*Gamma[4-(3*d)/2]*Gamma[2-d]*Gamma[-1+d/2]^2)/((-4+d)*(-2+d)*Gamma[3-d/2]^3);
Ta[5, 7] = (32*Gamma[2-d]*Gamma[-2+d/2])/((-4+d)^2*(-2+d)^2*Gamma[3-d/2]);
Ta[6, 3] = (8*(-2+d)*Gamma[6-2*d]*Gamma[5-(3*d)/2]^2*Gamma[-1+d/2]^3*Gamma[-4+(3*d)/2])/((-4+d)^2*Gamma[10-3*d]*Gamma[3-d/2]^2*Gamma[-1+d]^2);
Ta[6, 5] = (-8*Gamma[6-2*d]*Gamma[6-(3*d)/2]*Gamma[-1+d/2]^4*Gamma[-5+(3*d)/2])/((-4+d)^2*(-2+d)*Gamma[4-d]*Gamma[3-d/2]^2*Gamma[-2+d]^2);
Ta[6, 6] = (12*Gamma[8-2*d]*Gamma[4-d]*Gamma[-1+d/2]^4*Gamma[-8+2*d])/((-3+d)*Gamma[3-d/2]^3*Gamma[-2+d]*Gamma[-2+(3*d)/2]);
Ta[6, 7] = (-256*Gamma[6-(3*d)/2]*Gamma[d/2]^3*Gamma[-5+(3*d)/2])/((-4+d)^3*(-2+d)^3*Gamma[3-d/2]*Gamma[-1+d]^2);
Ta[6, 8] = (4*Gamma[2-d]^2*Gamma[-2+d/2]^2)/Gamma[3-d/2]^2;
Ta[7, 3] = (32*Gamma[8-2*d]*Gamma[-1+d/2]^5*Gamma[-8+2*d])/((-4+d)^2*(-2+d)*Gamma[3-d/2]*Gamma[-2+d]^3);

(* ====== Relations ====== *)

eq[80] = (
  T[6,13,0]-T[6,18,0] ==
  -11561/128-(209*Pi^2)/48-(11*Pi^4)/60+(9*Sqrt[3]*s2)/2+(205*Zeta[3])/12
);

eq[81] = (
  T[7,10,0]+T[7,11,0] ==
  -567/4-(19*Pi^2)/6-(11*Pi^4)/24-4*Sqrt[3]*s2-2*s2^2+42*Zeta[3]
);

eq[82] = (
  2*T[6,13,0]+3*T[6,14,0] ==
  89/4-(5*Pi^2)/3+(2*Pi^4)/9+12*Sqrt[3]*s2+6*s2^2-(35*Zeta[3])/6
);

eq[83] = (
  2*T[5,10,0]-3*T[6,10,0] ==
  -3651/16-(9*Pi^4)/40+9*Sqrt[3]*s2-12*s2^2+(115*Zeta[3])/2
);

eq[84] = (
  2*T[5,10,0]+T[6,9,0] ==
  -1731/16-(3*Pi^4)/10+19*Sqrt[3]*s2-16*s2^2+(185*Zeta[3])/6
);

eq[85] = (
  6*T[5,9,1]+T[5,10,0]-T[6,11,0] ==
  -2097/64+(59*Pi^4)/120-20*Sqrt[3]*s2+20*s2^2-(23*Zeta[3])/6
);

eq[86] = (
  6*T[5,9,1]-T[5,10,0]+3*T[6,12,0] ==
  5435/64+(2*Pi^4)/5-9*Sqrt[3]*s2+24*s2^2+Zeta[3]
);

eq[87] = (
  4*T[6,13,0]-2*T[6,13,1]-3*T[6,14,1]-T[8,10,0] ==
  -323/4+(41*Pi^2)/6-(3*Pi^4)/5+(32*Zeta[3])/3+(7*Pi^2*Zeta[3])/3+(197*Zeta[5])/3
);

eq[88] = (
  T[5,10,0]-T[6,17,0]+T[7,5,0] ==
  -50447/192-(37*Pi^2)/6-Pi^4/6+9*Sqrt[3]*s2+9*Zeta[3]+(2*Pi^2*Zeta[3])/9+(437*Zeta[5])/9
);

eq[89] = (
  5*T[5,8,1]-21*T[5,9,1]-2*T[7,14,1] ==
  2716223/2592+5*Pi^2+(19*Pi^4)/9+(467*s2)/Sqrt[3]+(4*s2^2)/3-(908*Zeta[3])/3+(Pi^2*Zeta[3])/3-(160*Zeta[5])/3
);

eq[90] = (
  T[5,8,1]-9*T[5,9,1]-2*T[7,10,1]-2*T[7,11,1] ==
  3471265/2592+(22*Pi^2)/3+(163*Pi^4)/36+(217*s2)/Sqrt[3]-(40*s2^2)/3-(890*Zeta[3])/3+(17*Pi^2*Zeta[3])/3-(1280*Zeta[5])/3
);

eq[91] = (
  T[5,8,1]-15*T[5,9,1]+2*T[6,16,1] ==
  -8415131/10368+(1255*Pi^2)/48-(2569*Pi^4)/360+(349*s2)/Sqrt[3]-(298*s2^2)/3-(64*Pi^2*Log[2]^2)/3+(64*Log[2]^4)/3+512*PolyLog[4,1/2]+(7777*Zeta[3])/12-Pi^2*Zeta[3]-9*Zeta[5]
);

eq[92] = (
  T[5,8,1]+15*T[5,9,1]+2*T[5,10,0]+4*T[5,10,1]-2*T[6,11,1]-6*T[6,12,1]+2*T[7,16,0] ==
  -957005/2592+(347*Pi^4)/180-(53*s2)/Sqrt[3]+(140*s2^2)/3-(157*Zeta[3])/3+(13*Pi^2*Zeta[3])/3-(203*Zeta[5])/3
);

eq[93] = (
  T[5,8,1]+15*T[5,9,1]-2*T[5,10,0]+4*T[5,10,1]-2*T[6,11,1]-6*T[6,12,1]-12*T[6,13,0]+4*T[6,13,1]+6*T[6,14,1]+6*T[7,8,0] ==
  -2290913/2592-(31*Pi^2)/3+(59*Pi^4)/20-(179*s2)/Sqrt[3]+(212*s2^2)/3-87*Zeta[3]+Pi^2*Zeta[3]-(389*Zeta[5])/3
);

eq[94] = (
  4*T[5,8,1]-12*T[5,9,1]-5*T[5,10,0]-2*T[5,10,1]+T[6,11,1]+3*T[6,12,1]-12*T[6,13,0]+4*T[6,13,1]+6*T[6,14,1]-3*T[7,9,0] ==
  964973/1296-(97*Pi^2)/12+Pi^4/4+(184*s2)/Sqrt[3]-(52*s2^2)/3-138*Zeta[3]+Pi^2*Zeta[3]-(47*Zeta[5])/3
);

eq[95] = (
  T[5,8,1]+15*T[5,9,1]+T[5,10,0]+4*T[5,10,1]-2*T[6,11,1]-6*T[6,12,1]-6*T[6,13,0]+4*T[6,13,1]+6*T[6,14,1]-3*T[6,17,0]-3*T[7,7,0] ==
  -1284721/5184-(101*Pi^2)/6+(16*Pi^4)/5-(8*s2)/Sqrt[3]+(176*s2^2)/3+62*Zeta[3]-(Pi^2*Zeta[3])/3-361*Zeta[5]
);

eq[96] = (
  7*T[5,8,1]-31*T[5,9,1]-3*T[5,10,0]-2*T[5,10,1]+T[6,11,1]+3*T[6,12,1]-20*T[6,13,0]+8*T[6,13,1]+12*T[6,14,1]+2*T[9,3,0]-6*T[9,4,0] ==
  590641/2592-23*Pi^2-(169*Pi^4)/120+(508*s2)/Sqrt[3]-(376*s2^2)/3-(53*Pi^4*Log[2])/15-(4*Pi^2*Log[2]^2)/3-(16*Pi^2*Log[2]^3)/3+(4*Log[2]^4)/3+(16*Log[2]^5)/5+32*PolyLog[4,1/2]-384*PolyLog[5,1/2]-(647*Zeta[3])/6+(7*Pi^2*Zeta[3])/3+425*Zeta[5]
);
