(* Old master integrals *)

Equation34 = ( -1039/6 - 1/(6*ep^4) - 3/(2*ep^3) - 26/(3*ep^2) - 41/ep
    - (7*Pi^4)/30 - Pi^4/(60*ep) + (5 - (2*(9 + ep^(-1)))/9 - ep^(-2) -
    28/(9*ep))*Zeta[3] + 53*Zeta[5] ); 

Equation36 = ( -635/24 - 2/ep^4 - 29/(3*ep^3) - 163/(6*ep^2) -
    601/(12*ep) + (4549*ep)/16 + (136*ep*Pi^4)/45 +
    (32*ep*Pi^2*Log[2]^2)/3 - (32*ep*Log[2]^4)/3 - 256*ep*PolyLog[4, 1/2]
    + (-220/9 + 8/(3*ep) - (2530*ep)/9 - (2*(12 + 58*ep +
    163*ep^2))/(9*ep))*Zeta[3] ); 

Equation38 = ( -89383/1296 - ep^(-4) - 9/(2*ep^3) - 415/(36*ep^2) -
    4991/(216*ep) - (3367679*ep)/7776 - (137735095*ep^2)/46656 -
    (4352*ep^2*Pi^4)/135 - (1024*ep^2*Pi^2*Log[2]^2)/9 +
    (1024*ep^2*Log[2]^4)/9 + (8192*ep^2*PolyLog[4, 1/2])/3 + ((1792*ep)/9
    + (68992*ep^2)/27)*Zeta[3] + ((4*ep)/5 + (18*ep^2)/5 - (2*ep*(2 +
    9*ep))/5)*Zeta[5] ); 

Equation40 = ( -237/2 + 3/(2*ep^4) + 19/(2*ep^3) + 67/(2*ep^2) +
     127/(2*ep) - (3969*ep)/2 - Pi^4/20 - (559*ep*Pi^4)/60 -
     32*ep*Pi^2*Log[2]^2 + 32*ep*Log[2]^4 + 768*ep*PolyLog[4, 1/2] + (253/3
     - 5/ep + (3649*ep)/3 + (2*(3 + 19*ep + 67*ep^2))/(3*ep))*Zeta[3] +
     189*ep*Zeta[5] ); 

Equation42 = ( -118 + 2/(3*ep^4) + 4/ep^3 + 38/(3*ep^2) + 44/(3*ep) -
    1156*ep - (4*Pi^4)/15 - (374*ep*Pi^4)/45 - (64*ep*Pi^2*Log[2]^2)/3 +
    (64*ep*Log[2]^4)/3 + 512*ep*PolyLog[4, 1/2] + (248/3 + 40/(9*ep) +
    (6304*ep)/9 + (8*(1 + 6*ep + 19*ep^2))/(9*ep))*Zeta[3] + 96*ep*Zeta[5]
    ); 

Equation44 = ( -147/2 + 1/(4*ep^4) + 3/(2*ep^3) + 9/(2*ep^2) + 2/ep -
    624*ep - (5*Pi^4)/8 - (27*ep*Pi^4)/10 + (28 + 37/(6*ep) + 39*ep + (1 +
    6*ep + 18*ep^2)/(3*ep))*Zeta[3] + (693*ep*Zeta[5])/2 ); 

Equation46 = ( -1039/6 - 1/(6*ep^4) - 3/(2*ep^3) - 26/(3*ep^2) - 41/ep
    - (29*Pi^4)/120 - Pi^4/(120*ep) + (21/2 - (2*(9 + ep^(-1)))/9 -
    1/(2*ep^2) + 7/(18*ep))*Zeta[3] + (43*Zeta[5])/2 ); 


(* equations (34) to (46) in the SeriesData[..] format *)

SeriesDataEquation34 = ( SeriesData[ep, 0, {-1/6, -3/2, -26/3 -
    Zeta[3], -41 -  Pi^4/60 - (10*Zeta[3])/3, -299/3 + (49*Pi^2)/18 -
    (7*Pi^4)/30 -  49*(3/2 + Pi^2/18) + 3*Zeta[3] + 53*Zeta[5]},
    -4, 1, 1] );  

SeriesDataEquation36 = ( SeriesData[ep, 0, {-2, -29/3, -163/6,
    -601/12, 733/24 -     (125*Pi^2)/18 - (2*Pi^4)/9 - (Pi^2*(-163/6 -
    (2*Pi^2)/3))/3 -     38*(3/2 + Pi^2/18) - (112*Zeta[3])/3, 4549/16
    + (136*Pi^4)/45 +     (32*Pi^2*Log[2]^2)/3 - (32*Log[2]^4)/3 -
    256*PolyLog[4, 1/2] -     (952*Zeta[3])/3}, -4, 2, 1] ); 

SeriesDataEquation38 = ( SeriesData[ep, 0, {-1, -9/2, -415/36, -4991/216,
    -6439/1296 - (53*Pi^2)/36 - Pi^4/9 - (Pi^2*(-415/36 - Pi^2/3))/3 -
    (128*(3/2 + Pi^2/18))/3, -3367679/7776 + (1792*Zeta[3])/9,
    -137735095/46656 - (4352*Pi^4)/135 - (1024*Pi^2*Log[2]^2)/9 +
    (1024*Log[2]^4)/9 + (8192*PolyLog[4, 1/2])/3 +
    (68992*Zeta[3])/27}, -4, 3, 1] ); 

SeriesDataEquation40 = ( SeriesData[ep, 0, {3/2, 19/2, 67/2, 127/2 -
    3*Zeta[3], -687/2 + (17*Pi^2)/6 + (7*Pi^4)/60 + 150*(3/2 +
    Pi^2/18) - (Pi^2*(67/2 + Pi^2/2))/3 + 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]}, -4, 2, 1] ); 

SeriesDataEquation42 = ( SeriesData[ep, 0, {2/3, 4, 38/3, 44/3 +
    (16*Zeta[3])/3,  -5419/12 - (59*Pi^2)/9 - (26*Pi^4)/135 -
    73*(-43/12 - Pi^2/9) +     48*(3/2 + Pi^2/18) - (Pi^2*(38/3 +
    (2*Pi^2)/9))/3 + 88*Zeta[3],     -1156 - (374*Pi^4)/45 -
    (64*Pi^2*Log[2]^2)/3 + (64*Log[2]^4)/3 +     512*PolyLog[4, 1/2] +
    (2152*Zeta[3])/3 + 96*Zeta[5]}, -4, 2, 1] );  

SeriesDataEquation44 = ( SeriesData[ep, 0, {1/4, 3/2, 9/2, 2 +
    (13*Zeta[3])/2,  -78 + (4*Pi^2)/3 - (43*Pi^4)/72 + 3*(3/2 +
    Pi^2/18) - (Pi^2*(9/2 +  Pi^2/12))/3 + 30*Zeta[3], -624 -
    (27*Pi^4)/10 + 45*Zeta[3] +  (693*Zeta[5])/2}, -4, 2, 1] ); 

SeriesDataEquation46 = ( SeriesData[ep, 0, {-1/6, -3/2, -26/3 -
    Zeta[3]/2, -41 - Pi^4/120 + Zeta[3]/6, -479/3 + Pi^2/2 -
    (29*Pi^4)/120 - 9*(3/2 + Pi^2/18) + (17*Zeta[3])/2 +
    (43*Zeta[5])/2}, -4, 1, 1] );  

(* Relation between different constants *)

Equation50 = ( 5511907345/31850496 + (89*Pi^4)/1440 + Pi^6/56 +
    (3*T62ep2)/64 + (11*Log[2]^4)/6 + 44*PolyLog[4, 1/2] ==
    (9*T54ep3)/1024 + T91ep0/4 + (11*Pi^2*Log[2]^2)/6 + (4103*Zeta[3])/144
    + Zeta[3]^2 + (273*Zeta[5])/8 );



(* Old masters in terms of ep-finite basis *)

T54 = (-((-2 + d)^3*(21458 - 20723*d + 7488*d^2 - 1200*d^3 + 72*d^4)*T41)/
      (2*(-3 + d)*(-7 + 2*d)*(-5 + 2*d)*(-10 + 3*d)*(-8 + 3*d)*
       (73 - 34*d + 4*d^2)*m^2) + (8*(-4 + d)^2*(-2 + d)*(25 - 22*d + 4*d^2)*
       T51)/((-7 + 2*d)^2*(-5 + 2*d)*(73 - 34*d + 4*d^2)) + 
     (256*(-5 + d)*(-4 + d)^3*(-3 + d)*m^6*T52f)/((-7 + 2*d)*(-5 + 2*d)*
       (-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) - 
     (32*(-4 + d)*(-33 + 10*d)*T53)/(3*(-11 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (512*(-4 + d)^4*m^6*T54f)/(3*(-5 + 2*d)*(-10 + 3*d)*(-8 + 3*d)*
       (73 - 34*d + 4*d^2)) + (128*(-4 + d)^4*(-37 + 11*d)*m^6*T61f)/
      (3*(-7 + 2*d)*(-5 + 2*d)*(-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (256*(-5 + d)*(-4 + d)^3*m^8*T62f)/((-5 + 2*d)*(-10 + 3*d)*(-8 + 3*d)*
       (73 - 34*d + 4*d^2)) - (256*(-4 + d)^2*(-3 + d)^3*m^2*T63)/
      ((-5 + 2*d)*(-11 + 3*d)*(-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (128*(-4 + d)^4*m^6*T64f)/((-7 + 2*d)*(-5 + 2*d)*(-10 + 3*d)*(-8 + 3*d)*
       (73 - 34*d + 4*d^2)) + (384*(-4 + d)^4*(-3 + d)*m^6*T71f)/
      ((-7 + 2*d)*(-5 + 2*d)*(-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) - 
     (128*(-4 + d)^4*(-89 + 25*d)*m^6*T72f)/(3*(-7 + 2*d)*(-5 + 2*d)*
       (-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (256*(-4 + d)^3*(-13 + 3*d)*m^8*T91)/(3*(-5 + 2*d)*(-10 + 3*d)*
       (-8 + 3*d)*(73 - 34*d + 4*d^2)));
 
T64 = ((3*(-4 + d)*(-2 + d)*(-8 + 3*d)*(82 - 48*d + 7*d^2)*T51)/
      (64*(-3 + d)^3*(-7 + 2*d)^2*m^2) - ((-5 + 2*d)*(-10 + 3*d)^2*(-8 + 3*d)*
       T53)/(4*(-4 + d)*(-3 + d)^3*m^2) - (3*(-10 + 3*d)*T63)/(4*(-3 + d)) + 
     (9*(-4 + d)^2*(-11 + 3*d)*(-10 + 3*d)*m^4*T71f)/
      (8*(-3 + d)^3*(-7 + 2*d)));
 
T52 = ((-2*(-2 + d)^2*(-1298 + 953*d - 233*d^2 + 19*d^3)*T41)/
      ((-3 + d)*(-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)*m^2) - 
     (4*(-5 + d)*(-4 + d)*(-17 + 4*d)*T51)/((-7 + 2*d)*(73 - 34*d + 4*d^2)) + 
     (32*(-5 + d)*(-4 + d)^2*(-3 + d)*m^6*T52f)/((-2 + d)*(-10 + 3*d)*
       (-8 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (16*(-7 + 2*d)*(-5 + 2*d)*(209 - 96*d + 11*d^2)*T53)/
      (3*(-2 + d)*(-11 + 3*d)*(73 - 34*d + 4*d^2)) - 
     (256*(-5 + d)*(-4 + d)^4*m^6*T54f)/(3*(-2 + d)*(-10 + 3*d)*(-8 + 3*d)*
       (73 - 34*d + 4*d^2)) - (16*(-4 + d)^3*(767 - 351*d + 40*d^2)*m^6*
       T61f)/(3*(-2 + d)*(-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) - 
     (128*(-5 + d)^2*(-4 + d)^3*m^8*T62f)/((-2 + d)*(-10 + 3*d)*(-8 + 3*d)*
       (73 - 34*d + 4*d^2)) + (32*(-4 + d)*(-3 + d)^2*(-7 + 2*d)*
       (149 - 69*d + 8*d^2)*m^2*T63)/((-2 + d)*(-11 + 3*d)*(-10 + 3*d)*
       (-8 + 3*d)*(73 - 34*d + 4*d^2)) - (32*(-4 + d)^3*(109 - 51*d + 6*d^2)*
       m^6*T64f)/((-2 + d)*(-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) - 
     (48*(-4 + d)^3*(149 - 69*d + 8*d^2)*m^6*T71f)/((-2 + d)*(-10 + 3*d)*
       (-8 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (16*(-4 + d)^3*(673 - 297*d + 32*d^2)*m^6*T72f)/
      (3*(-2 + d)*(-10 + 3*d)*(-8 + 3*d)*(73 - 34*d + 4*d^2)) - 
     (128*(-5 + d)*(-4 + d)^3*(-13 + 3*d)*m^8*T91)/(3*(-2 + d)*(-10 + 3*d)*
       (-8 + 3*d)*(73 - 34*d + 4*d^2)));
 
T61 = ((3*(-2 + d)^3*(11278 - 11551*d + 4411*d^2 - 745*d^3 + 47*d^4)*T41)/
      (16*(-3 + d)^3*(-7 + 2*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)*m^4) + 
     (3*(-5 + d)*(-4 + d)*(-2 + d)*(-8 + 3*d)*(-17 + 4*d)*(-16 + 5*d)*T51)/
      (8*(-3 + d)^2*(-7 + 2*d)^2*(73 - 34*d + 4*d^2)*m^2) - 
     (3*(-5 + d)*(-4 + d)^2*(-16 + 5*d)*m^4*T52f)/((-3 + d)*(-7 + 2*d)*
       (-10 + 3*d)*(73 - 34*d + 4*d^2)) - 
     ((-5 + 2*d)*(-8 + 3*d)*(21406 - 22007*d + 8444*d^2 - 1434*d^3 + 91*d^4)*
       T53)/(2*(-4 + d)*(-3 + d)^2*(-11 + 3*d)*(73 - 34*d + 4*d^2)*m^2) + 
     (8*(-5 + d)*(-4 + d)^4*(-16 + 5*d)*m^4*T54f)/((-3 + d)^2*(-7 + 2*d)*
       (-10 + 3*d)*(73 - 34*d + 4*d^2)) + 
     ((-4 + d)^2*(73178 - 75093*d + 28748*d^2 - 4869*d^3 + 308*d^4)*m^4*
       T61f)/(2*(-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (12*(-5 + d)^2*(-4 + d)^3*(-16 + 5*d)*m^6*T62f)/
      ((-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)) - 
     (3*(17566 - 18119*d + 6980*d^2 - 1191*d^3 + 76*d^4)*T63)/
      ((-11 + 3*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (3*(-4 + d)^2*(15006 - 15527*d + 6004*d^2 - 1029*d^3 + 66*d^4)*m^4*
       T64f)/((-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (9*(-4 + d)^2*(17566 - 18119*d + 6980*d^2 - 1191*d^3 + 76*d^4)*m^4*
       T71f)/(2*(-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)) - 
     ((-4 + d)^2*(18982 - 18219*d + 6388*d^2 - 963*d^3 + 52*d^4)*m^4*T72f)/
      (2*(-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (4*(-5 + d)*(-4 + d)^3*(-13 + 3*d)*(-16 + 5*d)*m^6*T91)/
      ((-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)));
 
T62 = (-((-4 + d)^2*(-2 + d)^3*(-1 + d)*T41)/(4*(-3 + d)^3*(-7 + 2*d)*
       (73 - 34*d + 4*d^2)*m^4) - ((-2 + d)*(-8 + 3*d)*
       (-499 + 409*d - 111*d^2 + 10*d^3)*T51)/(4*(-3 + d)^2*(-7 + 2*d)*
       (73 - 34*d + 4*d^2)*m^2) - (4*(-5 + d)*(-4 + d)^2*m^4*T52f)/
      ((-3 + d)*(-7 + 2*d)*(73 - 34*d + 4*d^2)) - 
     (2*(-5 + 2*d)*(-10 + 3*d)*(-8 + 3*d)*(63 - 28*d + 3*d^2)*T53)/
      (3*(-3 + d)^2*(-11 + 3*d)*(73 - 34*d + 4*d^2)*m^2) - 
     (16*(-4 + d)^3*(33 - 16*d + 2*d^2)*m^4*T54f)/(3*(-3 + d)^2*(-7 + 2*d)*
       (73 - 34*d + 4*d^2)) - (2*(-4 + d)^3*(255 - 125*d + 16*d^2)*m^4*T61f)/
      (3*(-3 + d)^2*(-7 + 2*d)*(73 - 34*d + 4*d^2)) - 
     (4*(-5 + d)*(-4 + d)^2*m^6*T62f)/((-3 + d)^2*(73 - 34*d + 4*d^2)) - 
     (4*(-4 + d)^2*(-19 + 4*d)*T63)/((-11 + 3*d)*(73 - 34*d + 4*d^2)) + 
     (4*(-4 + d)^4*(-9 + 2*d)*m^4*T64f)/((-3 + d)^2*(-7 + 2*d)*
       (73 - 34*d + 4*d^2)) + (6*(-4 + d)^4*(-19 + 4*d)*m^4*T71f)/
      ((-3 + d)^2*(-7 + 2*d)*(73 - 34*d + 4*d^2)) + 
     (2*(-4 + d)^3*(276 - 145*d + 20*d^2)*m^4*T72f)/(3*(-3 + d)^2*(-7 + 2*d)*
       (73 - 34*d + 4*d^2)) - (8*(-4 + d)^2*(-13 + 3*d)*(33 - 16*d + 2*d^2)*
       m^6*T91)/(3*(-3 + d)^2*(-7 + 2*d)*(73 - 34*d + 4*d^2)));
 
T72 = (-((-2 + d)^3*(1789 - 2111*d + 912*d^2 - 172*d^3 + 12*d^4)*T41)/
      (8*(-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)*m^6) + 
     ((-4 + d)*(-2 + d)*(-8 + 3*d)*(-914 + 525*d - 90*d^2 + 4*d^3)*T51)/
      (16*(-3 + d)*(-7 + 2*d)^2*(-10 + 3*d)*(73 - 34*d + 4*d^2)*m^4) - 
     (8*(-5 + d)*(-4 + d)^3*m^2*T52f)/((-7 + 2*d)*(-10 + 3*d)^2*
       (73 - 34*d + 4*d^2)) + ((-5 + 2*d)*(-8 + 3*d)*
       (9772 - 7470*d + 1165*d^2 + 399*d^3 - 142*d^4 + 12*d^5)*T53)/
      (6*(-4 + d)^2*(-3 + d)*(-11 + 3*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)*
       m^4) - (16*(-4 + d)^4*m^2*T54f)/(3*(-3 + d)*(-10 + 3*d)^2*
       (73 - 34*d + 4*d^2)) - (4*(-4 + d)^4*(-37 + 11*d)*m^2*T61f)/
      (3*(-3 + d)*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)) - 
     (8*(-5 + d)*(-4 + d)^3*m^4*T62f)/((-3 + d)*(-10 + 3*d)^2*
       (73 - 34*d + 4*d^2)) + ((-3 + d)^2*(12628 - 13328*d + 5255*d^2 - 
        918*d^3 + 60*d^4)*T63)/((-4 + d)*(-11 + 3*d)*(-10 + 3*d)^2*
       (73 - 34*d + 4*d^2)*m^2) - (2*(-4 + d)*(-35898 + 47927*d - 25554*d^2 + 
        6802*d^3 - 904*d^4 + 48*d^5)*m^2*T64f)/((-3 + d)*(-7 + 2*d)*
       (-10 + 3*d)^2*(73 - 34*d + 4*d^2)) - 
     (3*(-4 + d)*(12628 - 13328*d + 5255*d^2 - 918*d^3 + 60*d^4)*m^2*T71f)/
      (2*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)) + 
     (4*(-4 + d)^4*(-89 + 25*d)*m^2*T72f)/(3*(-3 + d)*(-7 + 2*d)*
       (-10 + 3*d)^2*(73 - 34*d + 4*d^2)) - 
     (8*(-4 + d)^3*(-13 + 3*d)*m^4*T91)/(3*(-3 + d)*(-10 + 3*d)^2*
       (73 - 34*d + 4*d^2)));
 
T71 = (-((-2 + d)^3*(5925 - 6293*d + 2499*d^2 - 440*d^3 + 29*d^4)*T41)/
      (4*(-3 + d)^2*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)*m^6) - 
     ((-2 + d)*(-8 + 3*d)*(9053 - 9587*d + 3799*d^2 - 668*d^3 + 44*d^4)*T51)/
      (4*(-3 + d)*(-7 + 2*d)^2*(-10 + 3*d)*(73 - 34*d + 4*d^2)*m^4) - 
     (4*(-5 + d)*(-4 + d)^2*(-19 + 5*d)*m^2*T52f)/((-7 + 2*d)*(-10 + 3*d)^2*
       (73 - 34*d + 4*d^2)) + (2*(-5 + 2*d)*(-8 + 3*d)*
       (-26546 + 35605*d - 19114*d^2 + 5132*d^3 - 689*d^4 + 37*d^5)*T53)/
      (3*(-4 + d)^2*(-3 + d)*(-11 + 3*d)*(-10 + 3*d)*(73 - 34*d + 4*d^2)*
       m^4) - (16*(-4 + d)^3*(-554 + 435*d - 114*d^2 + 10*d^3)*m^2*T54f)/
      (3*(-3 + d)*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)) - 
     (2*(-4 + d)^3*(-7327 + 5682*d - 1473*d^2 + 128*d^3)*m^2*T61f)/
      (3*(-3 + d)*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)) - 
     (4*(-5 + d)*(-4 + d)^2*(54 - 29*d + 4*d^2)*m^4*T62f)/
      ((-3 + d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)) + 
     (4*(-3 + d)*(-19090 + 25909*d - 14064*d^2 + 3815*d^3 - 517*d^4 + 28*d^5)*
       T63)/((-4 + d)*(-11 + 3*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)*m^2) - 
     (4*(-4 + d)*(-21906 + 29301*d - 15680*d^2 + 4195*d^3 - 561*d^4 + 30*d^5)*
       m^2*T64f)/((-3 + d)*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)) - 
     (6*(-4 + d)*(-19090 + 25909*d - 14064*d^2 + 3815*d^3 - 517*d^4 + 28*d^5)*
       m^2*T71f)/((-3 + d)*(-7 + 2*d)*(-10 + 3*d)^2*(73 - 34*d + 4*d^2)) + 
     (2*(-4 + d)*(-50762 + 63841*d - 32144*d^2 + 8107*d^3 - 1025*d^4 + 
        52*d^5)*m^2*T72f)/(3*(-3 + d)*(-7 + 2*d)*(-10 + 3*d)^2*
       (73 - 34*d + 4*d^2)) - (8*(-4 + d)^2*(-13 + 3*d)*
       (-554 + 435*d - 114*d^2 + 10*d^3)*m^4*T91)/(3*(-3 + d)*(-7 + 2*d)*
       (-10 + 3*d)^2*(73 - 34*d + 4*d^2)));
