(*  This file contains the DeWitt-Schwinger approximation for the stress-energy
    tensor for free quantized massive scalar fields in a general static spherically
    symmetric spacetime.  The quantity phidsmm2 is equal to the approximation for
    <\phi^2> to order 1/m^2 with m the mass of the field.  The quantity phidsmm4
    contains the order 1/m^4 terms for this approximation.  Each component of the 
    stress-energy tensor has one index up and the other down.  The variables are 
    as follows:  si is the coupling to the scalar curvature, g_{tt} = f[r], 
    g_{rr} = h[r], the temperature is equal to k/(2*Pi).  The units are 
    G = c = hbar = k_B = 1 and conventions are those of Misner, Thorne, and Wheeler.*)

(*  Please feel free to use these results in whatever way you like.  We only
    ask that you reference us.  The appropriate reference for this approximation
    is Paul R. Anderson, William A. Hiscock, and David A. Samuel, Physical
    Review D 51, 4337 (1995).  *)

phidsmm2 = (1/m^2)*((7*fp^4)/(1280*f^4*h^2*Pi^2) - (13*fp^2*fpp)/(960*f^3*h^2*Pi^2) + 
   fpp^2/(192*f^2*h^2*Pi^2) + (fp*f3p)/(192*f^2*h^2*Pi^2) - 
   f4p/(480*f*h^2*Pi^2) + (13*fp^3*hp)/(1920*f^3*h^3*Pi^2) - 
   (5*fp*fpp*hp)/(384*f^2*h^3*Pi^2) + (f3p*hp)/(160*f*h^3*Pi^2) + 
   (5*fp^2*hp^2)/(768*f^2*h^4*Pi^2) - (19*fpp*hp^2)/(1920*f*h^4*Pi^2) + 
   (7*fp*hp^3)/(960*f*h^5*Pi^2) - (fp^2*hpp)/(384*f^2*h^3*Pi^2) + 
   (fpp*hpp)/(240*f*h^3*Pi^2) - (13*fp*hp*hpp)/(1920*f*h^4*Pi^2) + 
   (fp*h3p)/(960*f*h^3*Pi^2) + 1/(240*Pi^2*r^4) - 1/(240*h^2*Pi^2*r^4) + 
   fp/(96*f*h^2*Pi^2*r^3) - fp/(96*f*h*Pi^2*r^3) - 
   (3*hp)/(160*h^3*Pi^2*r^3) + hp/(96*h^2*Pi^2*r^3) + 
   (23*fp^2)/(5760*f^2*h^2*Pi^2*r^2) + fpp/(288*f*h^2*Pi^2*r^2) + 
   fp^2/(576*f^2*h*Pi^2*r^2) - fpp/(288*f*h*Pi^2*r^2) - 
   (fp*hp)/(96*f*h^3*Pi^2*r^2) + (fp*hp)/(576*f*h^2*Pi^2*r^2) - 
   (3*hp^2)/(640*h^4*Pi^2*r^2) + hpp/(240*h^3*Pi^2*r^2) - 
   (23*fp^3)/(1920*f^3*h^2*Pi^2*r) + (7*fp*fpp)/(320*f^2*h^2*Pi^2*r) - 
   f3p/(120*f*h^2*Pi^2*r) - (fp^2*hp)/(80*f^2*h^3*Pi^2*r) + 
   (13*fpp*hp)/(960*f*h^3*Pi^2*r) - (11*fp*hp^2)/(640*f*h^4*Pi^2*r) + 
   (7*hp^3)/(240*h^5*Pi^2*r) + (fp*hpp)/(120*f*h^3*Pi^2*r) - 
   (13*hp*hpp)/(480*h^4*Pi^2*r) + h3p/(240*h^3*Pi^2*r) - 
   (11*fp^4*si)/(384*f^4*h^2*Pi^2) + (7*fp^2*fpp*si)/(96*f^3*h^2*Pi^2) - 
   (fpp^2*si)/(32*f^2*h^2*Pi^2) - (5*fp*f3p*si)/(192*f^2*h^2*Pi^2) + 
   (f4p*si)/(96*f*h^2*Pi^2) - (7*fp^3*hp*si)/(192*f^3*h^3*Pi^2) + 
   (9*fp*fpp*hp*si)/(128*f^2*h^3*Pi^2) - (f3p*hp*si)/(32*f*h^3*Pi^2) - 
   (13*fp^2*hp^2*si)/(384*f^2*h^4*Pi^2) + (19*fpp*hp^2*si)/(384*f*h^4*Pi^2) - 
   (7*fp*hp^3*si)/(192*f*h^5*Pi^2) + (5*fp^2*hpp*si)/(384*f^2*h^3*Pi^2) - 
   (fpp*hpp*si)/(48*f*h^3*Pi^2) + (13*fp*hp*hpp*si)/(384*f*h^4*Pi^2) - 
   (fp*h3p*si)/(192*f*h^3*Pi^2) - si/(24*Pi^2*r^4) + si/(24*h*Pi^2*r^4) - 
   (5*fp*si)/(48*f*h^2*Pi^2*r^3) + (5*fp*si)/(48*f*h*Pi^2*r^3) + 
   (7*hp*si)/(48*h^3*Pi^2*r^3) - (5*hp*si)/(48*h^2*Pi^2*r^3) - 
   (fp^2*si)/(32*f^2*h^2*Pi^2*r^2) - (fpp*si)/(24*f*h^2*Pi^2*r^2) - 
   (fp^2*si)/(48*f^2*h*Pi^2*r^2) + (fpp*si)/(24*f*h*Pi^2*r^2) + 
   (11*fp*hp*si)/(96*f*h^3*Pi^2*r^2) - (fp*hp*si)/(48*f*h^2*Pi^2*r^2) - 
   (hpp*si)/(48*h^3*Pi^2*r^2) + (7*fp^3*si)/(96*f^3*h^2*Pi^2*r) - 
   (13*fp*fpp*si)/(96*f^2*h^2*Pi^2*r) + (f3p*si)/(24*f*h^2*Pi^2*r) + 
   (fp^2*hp*si)/(16*f^2*h^3*Pi^2*r) - (fpp*hp*si)/(24*f*h^3*Pi^2*r) + 
   (7*fp*hp^2*si)/(96*f*h^4*Pi^2*r) - (7*hp^3*si)/(48*h^5*Pi^2*r) - 
   (fp*hpp*si)/(24*f*h^3*Pi^2*r) + (13*hp*hpp*si)/(96*h^4*Pi^2*r) - 
   (h3p*si)/(48*h^3*Pi^2*r) + (fp^4*si^2)/(128*f^4*h^2*Pi^2) - 
   (fp^2*fpp*si^2)/(32*f^3*h^2*Pi^2) + (fpp^2*si^2)/(32*f^2*h^2*Pi^2) + 
   (fp^3*hp*si^2)/(64*f^3*h^3*Pi^2) - (fp*fpp*hp*si^2)/(32*f^2*h^3*Pi^2) + 
   (fp^2*hp^2*si^2)/(128*f^2*h^4*Pi^2) + si^2/(8*Pi^2*r^4) + 
   si^2/(8*h^2*Pi^2*r^4) - si^2/(4*h*Pi^2*r^4) + 
   (fp*si^2)/(4*f*h^2*Pi^2*r^3) - (fp*si^2)/(4*f*h*Pi^2*r^3) - 
   (hp*si^2)/(4*h^3*Pi^2*r^3) + (hp*si^2)/(4*h^2*Pi^2*r^3) + 
   (fp^2*si^2)/(16*f^2*h^2*Pi^2*r^2) + (fpp*si^2)/(8*f*h^2*Pi^2*r^2) + 
   (fp^2*si^2)/(16*f^2*h*Pi^2*r^2) - (fpp*si^2)/(8*f*h*Pi^2*r^2) - 
   (5*fp*hp*si^2)/(16*f*h^3*Pi^2*r^2) + (fp*hp*si^2)/(16*f*h^2*Pi^2*r^2) + 
   (hp^2*si^2)/(8*h^4*Pi^2*r^2) - (fp^3*si^2)/(16*f^3*h^2*Pi^2*r) + 
   (fp*fpp*si^2)/(8*f^2*h^2*Pi^2*r) - (fpp*hp*si^2)/(8*f*h^3*Pi^2*r) + 
   (fp*hp^2*si^2)/(16*f*h^4*Pi^2*r))
 
phidsmm4 = (1/m^4) ((355*fp^6)/(32256*f^6*h^3*Pi^2) - 
   (333*fp^4*fpp)/(8960*f^5*h^3*Pi^2) + 
   (439*fp^2*fpp^2)/(13440*f^4*h^3*Pi^2) - (19*fpp^3)/(4032*f^3*h^3*Pi^2) + 
   (733*fp^3*f3p)/(53760*f^4*h^3*Pi^2) - 
   (899*fp*fpp*f3p)/(53760*f^3*h^3*Pi^2) + (47*f3p^2)/(26880*f^2*h^3*Pi^2) - 
   (209*fp^2*f4p)/(53760*f^3*h^3*Pi^2) + (37*fpp*f4p)/(13440*f^2*h^3*Pi^2) + 
   (fp*f5p)/(1120*f^2*h^3*Pi^2) - f6p/(4480*f*h^3*Pi^2) + 
   (333*fp^5*hp)/(17920*f^5*h^4*Pi^2) - 
   (5711*fp^3*fpp*hp)/(107520*f^4*h^4*Pi^2) + 
   (3457*fp*fpp^2*hp)/(107520*f^3*h^4*Pi^2) + 
   (2153*fp^2*f3p*hp)/(107520*f^3*h^4*Pi^2) - 
   (121*fpp*f3p*hp)/(8960*f^2*h^4*Pi^2) - 
   (157*fp*f4p*hp)/(26880*f^2*h^4*Pi^2) + (3*f5p*hp)/(1792*f*h^4*Pi^2) + 
   (293*fp^4*hp^2)/(13440*f^4*h^5*Pi^2) - 
   (689*fp^2*fpp*hp^2)/(13440*f^3*h^5*Pi^2) + 
   (1829*fpp^2*hp^2)/(107520*f^2*h^5*Pi^2) + 
   (107*fp*f3p*hp^2)/(5376*f^2*h^5*Pi^2) - (25*f4p*hp^2)/(3584*f*h^5*Pi^2) + 
   (1819*fp^3*hp^3)/(80640*f^3*h^6*Pi^2) - 
   (4463*fp*fpp*hp^3)/(107520*f^2*h^6*Pi^2) + 
   (135*f3p*hp^3)/(7168*f*h^6*Pi^2) + (1193*fp^2*hp^4)/(53760*f^2*h^7*Pi^2) - 
   (113*fpp*hp^4)/(3584*f*h^7*Pi^2) + (13*fp*hp^5)/(512*f*h^8*Pi^2) - 
   (733*fp^4*hpp)/(107520*f^4*h^4*Pi^2) + 
   (347*fp^2*fpp*hpp)/(21504*f^3*h^4*Pi^2) - 
   (37*fpp^2*hpp)/(6720*f^2*h^4*Pi^2) - 
   (167*fp*f3p*hpp)/(26880*f^2*h^4*Pi^2) + (f4p*hpp)/(448*f*h^4*Pi^2) - 
   (113*fp^3*hp*hpp)/(6720*f^3*h^5*Pi^2) + 
   (167*fp*fpp*hp*hpp)/(5376*f^2*h^5*Pi^2) - 
   (51*f3p*hp*hpp)/(3584*f*h^5*Pi^2) - 
   (195*fp^2*hp^2*hpp)/(7168*f^2*h^6*Pi^2) + 
   (1403*fpp*hp^2*hpp)/(35840*f*h^6*Pi^2) - 
   (797*fp*hp^3*hpp)/(17920*f*h^7*Pi^2) + 
   (359*fp^2*hpp^2)/(107520*f^2*h^5*Pi^2) - 
   (11*fpp*hpp^2)/(2240*f*h^5*Pi^2) + (493*fp*hp*hpp^2)/(35840*f*h^6*Pi^2) + 
   (209*fp^3*h3p)/(107520*f^3*h^4*Pi^2) - 
   (97*fp*fpp*h3p)/(26880*f^2*h^4*Pi^2) + (3*f3p*h3p)/(1792*f*h^4*Pi^2) + 
   (53*fp^2*hp*h3p)/(10752*f^2*h^5*Pi^2) - 
   (129*fpp*hp*h3p)/(17920*f*h^5*Pi^2) + 
   (359*fp*hp^2*h3p)/(35840*f*h^6*Pi^2) - (9*fp*hpp*h3p)/(3584*f*h^5*Pi^2) - 
   (fp^2*h4p)/(2240*f^2*h^4*Pi^2) + (3*fpp*h4p)/(4480*f*h^4*Pi^2) - 
   (13*fp*hp*h4p)/(8960*f*h^5*Pi^2) + (fp*h5p)/(8960*f*h^4*Pi^2) + 
   1/(1260*Pi^2*r^6) - 5/(1008*h^3*Pi^2*r^6) + 1/(240*h*Pi^2*r^6) + 
   fp/(160*f*h^3*Pi^2*r^5) - fp/(288*f*h^2*Pi^2*r^5) - 
   fp/(360*f*h*Pi^2*r^5) - (19*hp)/(1120*h^4*Pi^2*r^5) + 
   hp/(288*h^3*Pi^2*r^5) + hp/(360*h^2*Pi^2*r^5) + 
   (11*fp^2)/(16128*f^2*h^3*Pi^2*r^4) - fpp/(360*f*h^3*Pi^2*r^4) + 
   fpp/(288*f*h^2*Pi^2*r^4) + fp^2/(2880*f^2*h*Pi^2*r^4) - 
   fpp/(1440*f*h*Pi^2*r^4) + (89*fp*hp)/(5760*f*h^4*Pi^2*r^4) - 
   (fp*hp)/(96*f*h^3*Pi^2*r^4) + (fp*hp)/(2880*f*h^2*Pi^2*r^4) - 
   (529*hp^2)/(16128*h^5*Pi^2*r^4) + hp^2/(96*h^4*Pi^2*r^4) + 
   (89*hpp)/(10080*h^4*Pi^2*r^4) - hpp/(288*h^3*Pi^2*r^4) + 
   fp^3/(240*f^3*h^3*Pi^2*r^3) - (41*fp*fpp)/(5040*f^2*h^3*Pi^2*r^3) + 
   f3p/(576*f*h^3*Pi^2*r^3) - fp^3/(288*f^3*h^2*Pi^2*r^3) + 
   (7*fp*fpp)/(1152*f^2*h^2*Pi^2*r^3) - f3p/(576*f*h^2*Pi^2*r^3) + 
   (13*fp^2*hp)/(1680*f^2*h^4*Pi^2*r^3) - (fpp*hp)/(960*f*h^4*Pi^2*r^3) - 
   (fp^2*hp)/(288*f^2*h^3*Pi^2*r^3) + (fpp*hp)/(384*f*h^3*Pi^2*r^3) + 
   (7*fp*hp^2)/(1440*f*h^5*Pi^2*r^3) - (fp*hp^2)/(192*f*h^4*Pi^2*r^3) - 
   (367*hp^3)/(10080*h^6*Pi^2*r^3) + (7*hp^3)/(576*h^5*Pi^2*r^3) - 
   (fp*hpp)/(288*f*h^4*Pi^2*r^3) + (fp*hpp)/(384*f*h^3*Pi^2*r^3) + 
   (193*hp*hpp)/(6720*h^5*Pi^2*r^3) - (13*hp*hpp)/(1152*h^4*Pi^2*r^3) - 
   (71*h3p)/(20160*h^4*Pi^2*r^3) + h3p/(576*h^3*Pi^2*r^3) + 
   (617*fp^4)/(161280*f^4*h^3*Pi^2*r^2) - 
   (229*fp^2*fpp)/(40320*f^3*h^3*Pi^2*r^2) - 
   (19*fpp^2)/(26880*f^2*h^3*Pi^2*r^2) + 
   (53*fp*f3p)/(40320*f^2*h^3*Pi^2*r^2) + f4p/(1440*f*h^3*Pi^2*r^2) + 
   (7*fp^4)/(3840*f^4*h^2*Pi^2*r^2) - (13*fp^2*fpp)/(2880*f^3*h^2*Pi^2*r^2) + 
   fpp^2/(576*f^2*h^2*Pi^2*r^2) + (fp*f3p)/(576*f^2*h^2*Pi^2*r^2) - 
   f4p/(1440*f*h^2*Pi^2*r^2) + (fp^3*hp)/(32256*f^3*h^4*Pi^2*r^2) + 
   (439*fp*fpp*hp)/(80640*f^2*h^4*Pi^2*r^2) - 
   (17*f3p*hp)/(2880*f*h^4*Pi^2*r^2) + (13*fp^3*hp)/(5760*f^3*h^3*Pi^2*r^2) - 
   (5*fp*fpp*hp)/(1152*f^2*h^3*Pi^2*r^2) + (f3p*hp)/(480*f*h^3*Pi^2*r^2) - 
   (601*fp^2*hp^2)/(161280*f^2*h^5*Pi^2*r^2) + 
   (61*fpp*hp^2)/(3840*f*h^5*Pi^2*r^2) + 
   (5*fp^2*hp^2)/(2304*f^2*h^4*Pi^2*r^2) - 
   (19*fpp*hp^2)/(5760*f*h^4*Pi^2*r^2) - 
   (607*fp*hp^3)/(23040*f*h^6*Pi^2*r^2) + (7*fp*hp^3)/(2880*f*h^5*Pi^2*r^2) - 
   (53*hp^4)/(40320*h^7*Pi^2*r^2) - (fp^2*hpp)/(2240*f^2*h^4*Pi^2*r^2) - 
   (23*fpp*hpp)/(5760*f*h^4*Pi^2*r^2) - (fp^2*hpp)/(1152*f^2*h^3*Pi^2*r^2) + 
   (fpp*hpp)/(720*f*h^3*Pi^2*r^2) + (71*fp*hp*hpp)/(3840*f*h^5*Pi^2*r^2) - 
   (13*fp*hp*hpp)/(5760*f*h^4*Pi^2*r^2) + (67*hp^2*hpp)/(8960*h^6*Pi^2*r^2) - 
   (173*hpp^2)/(80640*h^5*Pi^2*r^2) - (11*fp*h3p)/(5760*f*h^4*Pi^2*r^2) + 
   (fp*h3p)/(2880*f*h^3*Pi^2*r^2) - (53*hp*h3p)/(20160*h^5*Pi^2*r^2) + 
   h4p/(2240*h^4*Pi^2*r^2) - (113*fp^5)/(5760*f^5*h^3*Pi^2*r) + 
   (1469*fp^3*fpp)/(26880*f^4*h^3*Pi^2*r) - 
   (2567*fp*fpp^2)/(80640*f^3*h^3*Pi^2*r) - 
   (1597*fp^2*f3p)/(80640*f^3*h^3*Pi^2*r) + 
   (71*fpp*f3p)/(5760*f^2*h^3*Pi^2*r) + (fp*f4p)/(180*f^2*h^3*Pi^2*r) - 
   (3*f5p)/(2240*f*h^3*Pi^2*r) - (247*fp^4*hp)/(8064*f^4*h^4*Pi^2*r) + 
   (5539*fp^2*fpp*hp)/(80640*f^3*h^4*Pi^2*r) - 
   (79*fpp^2*hp)/(3840*f^2*h^4*Pi^2*r) - 
   (29*fp*f3p*hp)/(1152*f^2*h^4*Pi^2*r) + (73*f4p*hp)/(10080*f*h^4*Pi^2*r) - 
   (1427*fp^3*hp^2)/(40320*f^3*h^5*Pi^2*r) + 
   (139*fp*fpp*hp^2)/(2304*f^2*h^5*Pi^2*r) - 
   (1781*f3p*hp^2)/(80640*f*h^5*Pi^2*r) - 
   (13*fp^2*hp^3)/(360*f^2*h^6*Pi^2*r) + 
   (3193*fpp*hp^3)/(80640*f*h^6*Pi^2*r) - 
   (1063*fp*hp^4)/(20160*f*h^7*Pi^2*r) + (13*hp^5)/(128*h^8*Pi^2*r) + 
   (481*fp^3*hpp)/(40320*f^3*h^4*Pi^2*r) - 
   (119*fp*fpp*hpp)/(5760*f^2*h^4*Pi^2*r) + 
   (311*f3p*hpp)/(40320*f*h^4*Pi^2*r) + 
   (161*fp^2*hp*hpp)/(5760*f^2*h^5*Pi^2*r) - 
   (1261*fpp*hp*hpp)/(40320*f*h^5*Pi^2*r) + 
   (307*fp*hp^2*hpp)/(4480*f*h^6*Pi^2*r) - (797*hp^3*hpp)/(4480*h^7*Pi^2*r) - 
   (743*fp*hpp^2)/(80640*f*h^5*Pi^2*r) + (493*hp*hpp^2)/(8960*h^6*Pi^2*r) - 
   (13*fp^2*h3p)/(3840*f^2*h^4*Pi^2*r) + (157*fpp*h3p)/(40320*f*h^4*Pi^2*r) - 
   (107*fp*hp*h3p)/(8064*f*h^5*Pi^2*r) + (359*hp^2*h3p)/(8960*h^6*Pi^2*r) - 
   (9*hpp*h3p)/(896*h^5*Pi^2*r) + (3*fp*h4p)/(2240*f*h^4*Pi^2*r) - 
   (13*hp*h4p)/(2240*h^5*Pi^2*r) + h5p/(2240*h^4*Pi^2*r) - 
   (139*fp^6*si)/(2560*f^6*h^3*Pi^2) + (143*fp^4*fpp*si)/(768*f^5*h^3*Pi^2) - 
   (323*fp^2*fpp^2*si)/(1920*f^4*h^3*Pi^2) + 
   (5*fpp^3*si)/(192*f^3*h^3*Pi^2) - (43*fp^3*f3p*si)/(640*f^4*h^3*Pi^2) + 
   (109*fp*fpp*f3p*si)/(1280*f^3*h^3*Pi^2) - 
   (17*f3p^2*si)/(1920*f^2*h^3*Pi^2) + (73*fp^2*f4p*si)/(3840*f^3*h^3*Pi^2) - 
   (7*fpp*f4p*si)/(480*f^2*h^3*Pi^2) - (fp*f5p*si)/(240*f^2*h^3*Pi^2) + 
   (f6p*si)/(960*f*h^3*Pi^2) - (143*fp^5*hp*si)/(1536*f^5*h^4*Pi^2) + 
   (1033*fp^3*fpp*hp*si)/(3840*f^4*h^4*Pi^2) - 
   (427*fp*fpp^2*hp*si)/(2560*f^3*h^4*Pi^2) - 
   (51*fp^2*f3p*hp*si)/(512*f^3*h^4*Pi^2) + 
   (9*fpp*f3p*hp*si)/(128*f^2*h^4*Pi^2) + 
   (9*fp*f4p*hp*si)/(320*f^2*h^4*Pi^2) - (f5p*hp*si)/(128*f*h^4*Pi^2) - 
   (839*fp^4*hp^2*si)/(7680*f^4*h^5*Pi^2) + 
   (497*fp^2*fpp*hp^2*si)/(1920*f^3*h^5*Pi^2) - 
   (137*fpp^2*hp^2*si)/(1536*f^2*h^5*Pi^2) - 
   (31*fp*f3p*hp^2*si)/(320*f^2*h^5*Pi^2) + 
   (25*f4p*hp^2*si)/(768*f*h^5*Pi^2) - 
   (863*fp^3*hp^3*si)/(7680*f^3*h^6*Pi^2) + 
   (791*fp*fpp*hp^3*si)/(3840*f^2*h^6*Pi^2) - 
   (45*f3p*hp^3*si)/(512*f*h^6*Pi^2) - 
   (103*fp^2*hp^4*si)/(960*f^2*h^7*Pi^2) + 
   (113*fpp*hp^4*si)/(768*f*h^7*Pi^2) - (91*fp*hp^5*si)/(768*f*h^8*Pi^2) + 
   (43*fp^4*hpp*si)/(1280*f^4*h^4*Pi^2) - 
   (619*fp^2*fpp*hpp*si)/(7680*f^3*h^4*Pi^2) + 
   (7*fpp^2*hpp*si)/(240*f^2*h^4*Pi^2) + 
   (19*fp*f3p*hpp*si)/(640*f^2*h^4*Pi^2) - (f4p*hpp*si)/(96*f*h^4*Pi^2) + 
   (319*fp^3*hp*hpp*si)/(3840*f^3*h^5*Pi^2) - 
   (589*fp*fpp*hp*hpp*si)/(3840*f^2*h^5*Pi^2) + 
   (17*f3p*hp*hpp*si)/(256*f*h^5*Pi^2) + 
   (167*fp^2*hp^2*hpp*si)/(1280*f^2*h^6*Pi^2) - 
   (1403*fpp*hp^2*hpp*si)/(7680*f*h^6*Pi^2) + 
   (797*fp*hp^3*hpp*si)/(3840*f*h^7*Pi^2) - 
   (121*fp^2*hpp^2*si)/(7680*f^2*h^5*Pi^2) + 
   (11*fpp*hpp^2*si)/(480*f*h^5*Pi^2) - 
   (493*fp*hp*hpp^2*si)/(7680*f*h^6*Pi^2) - 
   (73*fp^3*h3p*si)/(7680*f^3*h^4*Pi^2) + 
   (17*fp*fpp*h3p*si)/(960*f^2*h^4*Pi^2) - (f3p*h3p*si)/(128*f*h^4*Pi^2) - 
   (3*fp^2*hp*h3p*si)/(128*f^2*h^5*Pi^2) + 
   (43*fpp*hp*h3p*si)/(1280*f*h^5*Pi^2) - 
   (359*fp*hp^2*h3p*si)/(7680*f*h^6*Pi^2) + 
   (3*fp*hpp*h3p*si)/(256*f*h^5*Pi^2) + (fp^2*h4p*si)/(480*f^2*h^4*Pi^2) - 
   (fpp*h4p*si)/(320*f*h^4*Pi^2) + (13*fp*hp*h4p*si)/(1920*f*h^5*Pi^2) - 
   (fp*h5p*si)/(1920*f*h^4*Pi^2) - si/(120*Pi^2*r^6) + si/(24*h^2*Pi^2*r^6) - 
   si/(30*h*Pi^2*r^6) - (31*fp*si)/(720*f*h^3*Pi^2*r^5) + 
   (fp*si)/(144*f*h^2*Pi^2*r^5) + (13*fp*si)/(360*f*h*Pi^2*r^5) + 
   (47*hp*si)/(720*h^4*Pi^2*r^5) + (hp*si)/(48*h^3*Pi^2*r^5) - 
   (13*hp*si)/(360*h^2*Pi^2*r^5) + (fp^2*si)/(288*f^2*h^3*Pi^2*r^4) + 
   (11*fpp*si)/(360*f*h^3*Pi^2*r^4) - (fp^2*si)/(144*f^2*h^2*Pi^2*r^4) - 
   (fpp*si)/(24*f*h^2*Pi^2*r^4) - (fp^2*si)/(180*f^2*h*Pi^2*r^4) + 
   (fpp*si)/(90*f*h*Pi^2*r^4) - (47*fp*hp*si)/(288*f*h^4*Pi^2*r^4) + 
   (fp*hp*si)/(8*f*h^3*Pi^2*r^4) - (fp*hp*si)/(180*f*h^2*Pi^2*r^4) + 
   (2*hp^2*si)/(9*h^5*Pi^2*r^4) - (7*hp^2*si)/(72*h^4*Pi^2*r^4) - 
   (19*hpp*si)/(360*h^4*Pi^2*r^4) + (hpp*si)/(36*h^3*Pi^2*r^4) - 
   (121*fp^3*si)/(2880*f^3*h^3*Pi^2*r^3) + 
   (257*fp*fpp*si)/(2880*f^2*h^3*Pi^2*r^3) - 
   (5*f3p*si)/(288*f*h^3*Pi^2*r^3) + (23*fp^3*si)/(576*f^3*h^2*Pi^2*r^3) - 
   (41*fp*fpp*si)/(576*f^2*h^2*Pi^2*r^3) + (5*f3p*si)/(288*f*h^2*Pi^2*r^3) - 
   (119*fp^2*hp*si)/(1440*f^2*h^4*Pi^2*r^3) - 
   (fpp*hp*si)/(2880*f*h^4*Pi^2*r^3) + 
   (5*fp^2*hp*si)/(144*f^2*h^3*Pi^2*r^3) - (fpp*hp*si)/(64*f*h^3*Pi^2*r^3) - 
   (7*fp*hp^2*si)/(320*f*h^5*Pi^2*r^3) + (3*fp*hp^2*si)/(64*f*h^4*Pi^2*r^3) + 
   (121*hp^3*si)/(480*h^6*Pi^2*r^3) - (35*hp^3*si)/(288*h^5*Pi^2*r^3) + 
   (29*fp*hpp*si)/(960*f*h^4*Pi^2*r^3) - (5*fp*hpp*si)/(192*f*h^3*Pi^2*r^3) - 
   (197*hp*hpp*si)/(960*h^5*Pi^2*r^3) + (65*hp*hpp*si)/(576*h^4*Pi^2*r^3) + 
   (37*h3p*si)/(1440*h^4*Pi^2*r^3) - (5*h3p*si)/(288*h^3*Pi^2*r^3) - 
   (367*fp^4*si)/(11520*f^4*h^3*Pi^2*r^2) + 
   (127*fp^2*fpp*si)/(2880*f^3*h^3*Pi^2*r^2) + 
   (17*fpp^2*si)/(1440*f^2*h^3*Pi^2*r^2) - 
   (23*fp*f3p*si)/(2880*f^2*h^3*Pi^2*r^2) - 
   (11*f4p*si)/(1440*f*h^3*Pi^2*r^2) - (59*fp^4*si)/(2880*f^4*h^2*Pi^2*r^2) + 
   (37*fp^2*fpp*si)/(720*f^3*h^2*Pi^2*r^2) - 
   (fpp^2*si)/(48*f^2*h^2*Pi^2*r^2) - (11*fp*f3p*si)/(576*f^2*h^2*Pi^2*r^2) + 
   (11*f4p*si)/(1440*f*h^2*Pi^2*r^2) + 
   (137*fp^3*hp*si)/(11520*f^3*h^4*Pi^2*r^2) - 
   (49*fp*fpp*hp*si)/(640*f^2*h^4*Pi^2*r^2) + 
   (169*f3p*hp*si)/(2880*f*h^4*Pi^2*r^2) - 
   (37*fp^3*hp*si)/(1440*f^3*h^3*Pi^2*r^2) + 
   (19*fp*fpp*hp*si)/(384*f^2*h^3*Pi^2*r^2) - 
   (11*f3p*hp*si)/(480*f*h^3*Pi^2*r^2) + 
   (547*fp^2*hp^2*si)/(11520*f^2*h^5*Pi^2*r^2) - 
   (173*fpp*hp^2*si)/(1152*f*h^5*Pi^2*r^2) - 
   (7*fp^2*hp^2*si)/(288*f^2*h^4*Pi^2*r^2) + 
   (209*fpp*hp^2*si)/(5760*f*h^4*Pi^2*r^2) + 
   (997*fp*hp^3*si)/(3840*f*h^6*Pi^2*r^2) - 
   (77*fp*hp^3*si)/(2880*f*h^5*Pi^2*r^2) - (21*hp^4*si)/(320*h^7*Pi^2*r^2) + 
   (fp^2*hpp*si)/(2880*f^2*h^4*Pi^2*r^2) + 
   (59*fpp*hpp*si)/(1440*f*h^4*Pi^2*r^2) + 
   (11*fp^2*hpp*si)/(1152*f^2*h^3*Pi^2*r^2) - 
   (11*fpp*hpp*si)/(720*f*h^3*Pi^2*r^2) - 
   (361*fp*hp*hpp*si)/(1920*f*h^5*Pi^2*r^2) + 
   (143*fp*hp*hpp*si)/(5760*f*h^4*Pi^2*r^2) + 
   (21*hp^2*hpp*si)/(640*h^6*Pi^2*r^2) + (hpp^2*si)/(160*h^5*Pi^2*r^2) + 
   (29*fp*h3p*si)/(1440*f*h^4*Pi^2*r^2) - 
   (11*fp*h3p*si)/(2880*f*h^3*Pi^2*r^2) + (hp*h3p*si)/(240*h^5*Pi^2*r^2) - 
   (h4p*si)/(480*h^4*Pi^2*r^2) + (1313*fp^5*si)/(11520*f^5*h^3*Pi^2*r) - 
   (41*fp^3*fpp*si)/(128*f^4*h^3*Pi^2*r) + 
   (1091*fp*fpp^2*si)/(5760*f^3*h^3*Pi^2*r) + 
   (649*fp^2*f3p*si)/(5760*f^3*h^3*Pi^2*r) - 
   (101*fpp*f3p*si)/(1440*f^2*h^3*Pi^2*r) - 
   (11*fp*f4p*si)/(360*f^2*h^3*Pi^2*r) + (f5p*si)/(160*f*h^3*Pi^2*r) + 
   (385*fp^4*hp*si)/(2304*f^4*h^4*Pi^2*r) - 
   (1067*fp^2*fpp*hp*si)/(2880*f^3*h^4*Pi^2*r) + 
   (17*fpp^2*hp*si)/(160*f^2*h^4*Pi^2*r) + 
   (37*fp*f3p*hp*si)/(288*f^2*h^4*Pi^2*r) - 
   (7*f4p*hp*si)/(240*f*h^4*Pi^2*r) + 
   (415*fp^3*hp^2*si)/(2304*f^3*h^5*Pi^2*r) - 
   (169*fp*fpp*hp^2*si)/(576*f^2*h^5*Pi^2*r) + 
   (51*f3p*hp^2*si)/(640*f*h^5*Pi^2*r) + 
   (203*fp^2*hp^3*si)/(1280*f^2*h^6*Pi^2*r) - 
   (223*fpp*hp^3*si)/(1920*f*h^6*Pi^2*r) + 
   (49*fp*hp^4*si)/(240*f*h^7*Pi^2*r) - (91*hp^5*si)/(192*h^8*Pi^2*r) - 
   (73*fp^3*hpp*si)/(1152*f^3*h^4*Pi^2*r) + 
   (617*fp*fpp*hpp*si)/(5760*f^2*h^4*Pi^2*r) - 
   (f3p*hpp*si)/(32*f*h^4*Pi^2*r) - 
   (81*fp^2*hp*hpp*si)/(640*f^2*h^5*Pi^2*r) + 
   (191*fpp*hp*hpp*si)/(1920*f*h^5*Pi^2*r) - 
   (269*fp*hp^2*hpp*si)/(960*f*h^6*Pi^2*r) + 
   (797*hp^3*hpp*si)/(960*h^7*Pi^2*r) + (13*fp*hpp^2*si)/(320*f*h^5*Pi^2*r) - 
   (493*hp*hpp^2*si)/(1920*h^6*Pi^2*r) + 
   (91*fp^2*h3p*si)/(5760*f^2*h^4*Pi^2*r) - 
   (13*fpp*h3p*si)/(960*f*h^4*Pi^2*r) + 
   (11*fp*hp*h3p*si)/(192*f*h^5*Pi^2*r) - 
   (359*hp^2*h3p*si)/(1920*h^6*Pi^2*r) + (3*hpp*h3p*si)/(64*h^5*Pi^2*r) - 
   (fp*h4p*si)/(160*f*h^4*Pi^2*r) + (13*hp*h4p*si)/(480*h^5*Pi^2*r) - 
   (h5p*si)/(480*h^4*Pi^2*r) + (29*fp^6*si^2)/(1536*f^6*h^3*Pi^2) - 
   (21*fp^4*fpp*si^2)/(256*f^5*h^3*Pi^2) + 
   (13*fp^2*fpp^2*si^2)/(128*f^4*h^3*Pi^2) - 
   (5*fpp^3*si^2)/(192*f^3*h^3*Pi^2) + (3*fp^3*f3p*si^2)/(128*f^4*h^3*Pi^2) - 
   (3*fp*fpp*f3p*si^2)/(64*f^3*h^3*Pi^2) + (f3p^2*si^2)/(192*f^2*h^3*Pi^2) - 
   (fp^2*f4p*si^2)/(192*f^3*h^3*Pi^2) + (fpp*f4p*si^2)/(96*f^2*h^3*Pi^2) + 
   (21*fp^5*hp*si^2)/(512*f^5*h^4*Pi^2) - 
   (35*fp^3*fpp*hp*si^2)/(256*f^4*h^4*Pi^2) + 
   (7*fp*fpp^2*hp*si^2)/(64*f^3*h^4*Pi^2) + 
   (5*fp^2*f3p*hp*si^2)/(128*f^3*h^4*Pi^2) - 
   (3*fpp*f3p*hp*si^2)/(64*f^2*h^4*Pi^2) - 
   (fp*f4p*hp*si^2)/(192*f^2*h^4*Pi^2) + 
   (25*fp^4*hp^2*si^2)/(512*f^4*h^5*Pi^2) - 
   (97*fp^2*fpp*hp^2*si^2)/(768*f^3*h^5*Pi^2) + 
   (47*fpp^2*hp^2*si^2)/(768*f^2*h^5*Pi^2) + 
   (5*fp*f3p*hp^2*si^2)/(192*f^2*h^5*Pi^2) + 
   (23*fp^3*hp^3*si^2)/(512*f^3*h^6*Pi^2) - 
   (59*fp*fpp*hp^3*si^2)/(768*f^2*h^6*Pi^2) + 
   (3*fp^2*hp^4*si^2)/(128*f^2*h^7*Pi^2) - 
   (3*fp^4*hpp*si^2)/(256*f^4*h^4*Pi^2) + 
   (13*fp^2*fpp*hpp*si^2)/(384*f^3*h^4*Pi^2) - 
   (fpp^2*hpp*si^2)/(48*f^2*h^4*Pi^2) - 
   (fp*f3p*hpp*si^2)/(192*f^2*h^4*Pi^2) - 
   (11*fp^3*hp*hpp*si^2)/(384*f^3*h^5*Pi^2) + 
   (5*fp*fpp*hp*hpp*si^2)/(96*f^2*h^5*Pi^2) - 
   (17*fp^2*hp^2*hpp*si^2)/(768*f^2*h^6*Pi^2) + 
   (fp^2*hpp^2*si^2)/(768*f^2*h^5*Pi^2) + 
   (fp^3*h3p*si^2)/(384*f^3*h^4*Pi^2) - 
   (fp*fpp*h3p*si^2)/(192*f^2*h^4*Pi^2) + 
   (fp^2*hp*h3p*si^2)/(384*f^2*h^5*Pi^2) + si^2/(24*Pi^2*r^6) + 
   si^2/(8*h^3*Pi^2*r^6) - (5*si^2)/(24*h^2*Pi^2*r^6) + 
   si^2/(24*h*Pi^2*r^6) + (fp*si^2)/(6*f*h^2*Pi^2*r^5) - 
   (fp*si^2)/(6*f*h*Pi^2*r^5) + (hp*si^2)/(6*h^4*Pi^2*r^5) - 
   (hp*si^2)/(3*h^3*Pi^2*r^5) + (hp*si^2)/(6*h^2*Pi^2*r^5) - 
   (7*fp^2*si^2)/(96*f^2*h^3*Pi^2*r^4) - (5*fpp*si^2)/(48*f*h^3*Pi^2*r^4) + 
   (fp^2*si^2)/(16*f^2*h^2*Pi^2*r^4) + (fpp*si^2)/(6*f*h^2*Pi^2*r^4) + 
   (fp^2*si^2)/(32*f^2*h*Pi^2*r^4) - (fpp*si^2)/(16*f*h*Pi^2*r^4) + 
   (17*fp*hp*si^2)/(32*f*h^4*Pi^2*r^4) - 
   (23*fp*hp*si^2)/(48*f*h^3*Pi^2*r^4) + (fp*hp*si^2)/(32*f*h^2*Pi^2*r^4) - 
   (hp^2*si^2)/(3*h^5*Pi^2*r^4) + (hp^2*si^2)/(4*h^4*Pi^2*r^4) + 
   (hpp*si^2)/(24*h^4*Pi^2*r^4) - (hpp*si^2)/(24*h^3*Pi^2*r^4) + 
   (11*fp^3*si^2)/(96*f^3*h^3*Pi^2*r^3) - 
   (7*fp*fpp*si^2)/(24*f^2*h^3*Pi^2*r^3) + (f3p*si^2)/(24*f*h^3*Pi^2*r^3) - 
   (13*fp^3*si^2)/(96*f^3*h^2*Pi^2*r^3) + 
   (fp*fpp*si^2)/(4*f^2*h^2*Pi^2*r^3) - (f3p*si^2)/(24*f*h^2*Pi^2*r^3) + 
   (13*fp^2*hp*si^2)/(48*f^2*h^4*Pi^2*r^3) + 
   (fpp*hp*si^2)/(12*f*h^4*Pi^2*r^3) - (fp^2*hp*si^2)/(12*f^2*h^3*Pi^2*r^3) - 
   (fpp*hp*si^2)/(24*f*h^3*Pi^2*r^3) - (3*fp*hp^2*si^2)/(32*f*h^5*Pi^2*r^3) - 
   (7*fp*hp^2*si^2)/(96*f*h^4*Pi^2*r^3) - (hp^3*si^2)/(3*h^6*Pi^2*r^3) + 
   (7*hp^3*si^2)/(24*h^5*Pi^2*r^3) - (fp*hpp*si^2)/(16*f*h^4*Pi^2*r^3) + 
   (fp*hpp*si^2)/(16*f*h^3*Pi^2*r^3) + (5*hp*hpp*si^2)/(16*h^5*Pi^2*r^3) - 
   (13*hp*hpp*si^2)/(48*h^4*Pi^2*r^3) - (h3p*si^2)/(24*h^4*Pi^2*r^3) + 
   (h3p*si^2)/(24*h^3*Pi^2*r^3) + (31*fp^4*si^2)/(384*f^4*h^3*Pi^2*r^2) - 
   (5*fp^2*fpp*si^2)/(48*f^3*h^3*Pi^2*r^2) - 
   (5*fpp^2*si^2)/(96*f^2*h^3*Pi^2*r^2) + 
   (fp*f3p*si^2)/(96*f^2*h^3*Pi^2*r^2) + (f4p*si^2)/(48*f*h^3*Pi^2*r^2) + 
   (23*fp^4*si^2)/(384*f^4*h^2*Pi^2*r^2) - 
   (5*fp^2*fpp*si^2)/(32*f^3*h^2*Pi^2*r^2) + 
   (7*fpp^2*si^2)/(96*f^2*h^2*Pi^2*r^2) + 
   (5*fp*f3p*si^2)/(96*f^2*h^2*Pi^2*r^2) - (f4p*si^2)/(48*f*h^2*Pi^2*r^2) - 
   (13*fp^3*hp*si^2)/(192*f^3*h^4*Pi^2*r^2) + 
   (19*fp*fpp*hp*si^2)/(64*f^2*h^4*Pi^2*r^2) - 
   (7*f3p*hp*si^2)/(48*f*h^4*Pi^2*r^2) + 
   (5*fp^3*hp*si^2)/(64*f^3*h^3*Pi^2*r^2) - 
   (29*fp*fpp*hp*si^2)/(192*f^2*h^3*Pi^2*r^2) + 
   (f3p*hp*si^2)/(16*f*h^3*Pi^2*r^2) - 
   (19*fp^2*hp^2*si^2)/(128*f^2*h^5*Pi^2*r^2) + 
   (21*fpp*hp^2*si^2)/(64*f*h^5*Pi^2*r^2) + 
   (9*fp^2*hp^2*si^2)/(128*f^2*h^4*Pi^2*r^2) - 
   (19*fpp*hp^2*si^2)/(192*f*h^4*Pi^2*r^2) - 
   (5*fp*hp^3*si^2)/(8*f*h^6*Pi^2*r^2) + 
   (7*fp*hp^3*si^2)/(96*f*h^5*Pi^2*r^2) + (3*hp^4*si^2)/(8*h^7*Pi^2*r^2) + 
   (fp^2*hpp*si^2)/(192*f^2*h^4*Pi^2*r^2) - 
   (5*fpp*hpp*si^2)/(48*f*h^4*Pi^2*r^2) - 
   (5*fp^2*hpp*si^2)/(192*f^2*h^3*Pi^2*r^2) + 
   (fpp*hpp*si^2)/(24*f*h^3*Pi^2*r^2) + 
   (91*fp*hp*hpp*si^2)/(192*f*h^5*Pi^2*r^2) - 
   (13*fp*hp*hpp*si^2)/(192*f*h^4*Pi^2*r^2) - 
   (17*hp^2*hpp*si^2)/(48*h^6*Pi^2*r^2) + (hpp^2*si^2)/(48*h^5*Pi^2*r^2) - 
   (5*fp*h3p*si^2)/(96*f*h^4*Pi^2*r^2) + (fp*h3p*si^2)/(96*f*h^3*Pi^2*r^2) + 
   (hp*h3p*si^2)/(24*h^5*Pi^2*r^2) - (41*fp^5*si^2)/(384*f^5*h^3*Pi^2*r) + 
   (61*fp^3*fpp*si^2)/(192*f^4*h^3*Pi^2*r) - 
   (5*fp*fpp^2*si^2)/(24*f^3*h^3*Pi^2*r) - 
   (3*fp^2*f3p*si^2)/(32*f^3*h^3*Pi^2*r) + 
   (fpp*f3p*si^2)/(16*f^2*h^3*Pi^2*r) + (fp*f4p*si^2)/(48*f^2*h^3*Pi^2*r) - 
   (15*fp^4*hp*si^2)/(128*f^4*h^4*Pi^2*r) + 
   (23*fp^2*fpp*hp*si^2)/(96*f^3*h^4*Pi^2*r) - 
   (fpp^2*hp*si^2)/(24*f^2*h^4*Pi^2*r) - 
   (5*fp*f3p*hp*si^2)/(96*f^2*h^4*Pi^2*r) - (f4p*hp*si^2)/(48*f*h^4*Pi^2*r) - 
   (9*fp^3*hp^2*si^2)/(128*f^3*h^5*Pi^2*r) + 
   (5*fp*fpp*hp^2*si^2)/(96*f^2*h^5*Pi^2*r) + 
   (5*f3p*hp^2*si^2)/(48*f*h^5*Pi^2*r) + 
   (17*fp^2*hp^3*si^2)/(384*f^2*h^6*Pi^2*r) - 
   (59*fpp*hp^3*si^2)/(192*f*h^6*Pi^2*r) + 
   (3*fp*hp^4*si^2)/(16*f*h^7*Pi^2*r) + 
   (7*fp^3*hpp*si^2)/(192*f^3*h^4*Pi^2*r) - 
   (5*fp*fpp*hpp*si^2)/(96*f^2*h^4*Pi^2*r) - 
   (f3p*hpp*si^2)/(48*f*h^4*Pi^2*r) - 
   (fp^2*hp*hpp*si^2)/(64*f^2*h^5*Pi^2*r) + 
   (5*fpp*hp*hpp*si^2)/(24*f*h^5*Pi^2*r) - 
   (17*fp*hp^2*hpp*si^2)/(96*f*h^6*Pi^2*r) + 
   (fp*hpp^2*si^2)/(96*f*h^5*Pi^2*r) - (fpp*h3p*si^2)/(48*f*h^4*Pi^2*r) + 
   (fp*hp*h3p*si^2)/(48*f*h^5*Pi^2*r) - (fp^6*si^3)/(768*f^6*h^3*Pi^2) + 
   (fp^4*fpp*si^3)/(128*f^5*h^3*Pi^2) - (fp^2*fpp^2*si^3)/(64*f^4*h^3*Pi^2) + 
   (fpp^3*si^3)/(96*f^3*h^3*Pi^2) - (fp^5*hp*si^3)/(256*f^5*h^4*Pi^2) + 
   (fp^3*fpp*hp*si^3)/(64*f^4*h^4*Pi^2) - 
   (fp*fpp^2*hp*si^3)/(64*f^3*h^4*Pi^2) - 
   (fp^4*hp^2*si^3)/(256*f^4*h^5*Pi^2) + 
   (fp^2*fpp*hp^2*si^3)/(128*f^3*h^5*Pi^2) - 
   (fp^3*hp^3*si^3)/(768*f^3*h^6*Pi^2) - si^3/(12*Pi^2*r^6) + 
   si^3/(12*h^3*Pi^2*r^6) - si^3/(4*h^2*Pi^2*r^6) + si^3/(4*h*Pi^2*r^6) + 
   (fp*si^3)/(4*f*h^3*Pi^2*r^5) - (fp*si^3)/(2*f*h^2*Pi^2*r^5) + 
   (fp*si^3)/(4*f*h*Pi^2*r^5) - (hp*si^3)/(4*h^4*Pi^2*r^5) + 
   (hp*si^3)/(2*h^3*Pi^2*r^5) - (hp*si^3)/(4*h^2*Pi^2*r^5) + 
   (3*fp^2*si^3)/(16*f^2*h^3*Pi^2*r^4) + (fpp*si^3)/(8*f*h^3*Pi^2*r^4) - 
   (fp^2*si^3)/(8*f^2*h^2*Pi^2*r^4) - (fpp*si^3)/(4*f*h^2*Pi^2*r^4) - 
   (fp^2*si^3)/(16*f^2*h*Pi^2*r^4) + (fpp*si^3)/(8*f*h*Pi^2*r^4) - 
   (9*fp*hp*si^3)/(16*f*h^4*Pi^2*r^4) + (5*fp*hp*si^3)/(8*f*h^3*Pi^2*r^4) - 
   (fp*hp*si^3)/(16*f*h^2*Pi^2*r^4) + (hp^2*si^3)/(4*h^5*Pi^2*r^4) - 
   (hp^2*si^3)/(4*h^4*Pi^2*r^4) - (fp^3*si^3)/(24*f^3*h^3*Pi^2*r^3) + 
   (fp*fpp*si^3)/(4*f^2*h^3*Pi^2*r^3) + (fp^3*si^3)/(8*f^3*h^2*Pi^2*r^3) - 
   (fp*fpp*si^3)/(4*f^2*h^2*Pi^2*r^3) - (fp^2*hp*si^3)/(4*f^2*h^4*Pi^2*r^3) - 
   (fpp*hp*si^3)/(4*f*h^4*Pi^2*r^3) + (fpp*hp*si^3)/(4*f*h^3*Pi^2*r^3) + 
   (3*fp*hp^2*si^3)/(8*f*h^5*Pi^2*r^3) - (fp*hp^2*si^3)/(8*f*h^4*Pi^2*r^3) - 
   (hp^3*si^3)/(12*h^6*Pi^2*r^3) - (3*fp^4*si^3)/(64*f^4*h^3*Pi^2*r^2) + 
   (fp^2*fpp*si^3)/(16*f^3*h^3*Pi^2*r^2) + 
   (fpp^2*si^3)/(16*f^2*h^3*Pi^2*r^2) - (fp^4*si^3)/(64*f^4*h^2*Pi^2*r^2) + 
   (fp^2*fpp*si^3)/(16*f^3*h^2*Pi^2*r^2) - 
   (fpp^2*si^3)/(16*f^2*h^2*Pi^2*r^2) + 
   (3*fp^3*hp*si^3)/(32*f^3*h^4*Pi^2*r^2) - 
   (5*fp*fpp*hp*si^3)/(16*f^2*h^4*Pi^2*r^2) - 
   (fp^3*hp*si^3)/(32*f^3*h^3*Pi^2*r^2) + 
   (fp*fpp*hp*si^3)/(16*f^2*h^3*Pi^2*r^2) + 
   (5*fp^2*hp^2*si^3)/(64*f^2*h^5*Pi^2*r^2) + 
   (fpp*hp^2*si^3)/(8*f*h^5*Pi^2*r^2) - 
   (fp^2*hp^2*si^3)/(64*f^2*h^4*Pi^2*r^2) - 
   (fp*hp^3*si^3)/(16*f*h^6*Pi^2*r^2) + (fp^5*si^3)/(64*f^5*h^3*Pi^2*r) - 
   (fp^3*fpp*si^3)/(16*f^4*h^3*Pi^2*r) + 
   (fp*fpp^2*si^3)/(16*f^3*h^3*Pi^2*r) + (fp^4*hp*si^3)/(64*f^4*h^4*Pi^2*r) - 
   (fpp^2*hp*si^3)/(16*f^2*h^4*Pi^2*r) - 
   (fp^3*hp^2*si^3)/(64*f^3*h^5*Pi^2*r) + 
   (fp*fpp*hp^2*si^3)/(16*f^2*h^5*Pi^2*r) - 
   (fp^2*hp^3*si^3)/(64*f^2*h^6*Pi^2*r))
 
Tththdsmm2 = (1/m^2)((-6163*fp^6)/(322560*f^6*h^3*Pi^2) + 
   (1137*fp^4*fpp)/(17920*f^5*h^3*Pi^2) - 
   (1451*fp^2*fpp^2)/(26880*f^4*h^3*Pi^2) + 
   (283*fpp^3)/(40320*f^3*h^3*Pi^2) - (2557*fp^3*f3p)/(107520*f^4*h^3*Pi^2) + 
   (607*fp*fpp*f3p)/(21504*f^3*h^3*Pi^2) - (163*f3p^2)/(53760*f^2*h^3*Pi^2) + 
   (243*fp^2*f4p)/(35840*f^3*h^3*Pi^2) - (39*fpp*f4p)/(8960*f^2*h^3*Pi^2) - 
   (11*fp*f5p)/(6720*f^2*h^3*Pi^2) + (11*f6p)/(26880*f*h^3*Pi^2) - 
   (1137*fp^5*hp)/(35840*f^5*h^4*Pi^2) + 
   (19279*fp^3*fpp*hp)/(215040*f^4*h^4*Pi^2) - 
   (11369*fp*fpp^2*hp)/(215040*f^3*h^4*Pi^2) - 
   (7409*fp^2*f3p*hp)/(215040*f^3*h^4*Pi^2) + 
   (397*fpp*f3p*hp)/(17920*f^2*h^4*Pi^2) + 
   (557*fp*f4p*hp)/(53760*f^2*h^4*Pi^2) - (11*f5p*hp)/(3584*f*h^4*Pi^2) - 
   (167*fp^4*hp^2)/(4480*f^4*h^5*Pi^2) + 
   (467*fp^2*fpp*hp^2)/(5376*f^3*h^5*Pi^2) - 
   (1971*fpp^2*hp^2)/(71680*f^2*h^5*Pi^2) - 
   (629*fp*f3p*hp^2)/(17920*f^2*h^5*Pi^2) + 
   (275*f4p*hp^2)/(21504*f*h^5*Pi^2) - 
   (1249*fp^3*hp^3)/(32256*f^3*h^6*Pi^2) + 
   (1025*fp*fpp*hp^3)/(14336*f^2*h^6*Pi^2) - 
   (495*f3p*hp^3)/(14336*f*h^6*Pi^2) - (845*fp^2*hp^4)/(21504*f^2*h^7*Pi^2) + 
   (1243*fpp*hp^4)/(21504*f*h^7*Pi^2) - (143*fp*hp^5)/(3072*f*h^8*Pi^2) + 
   (2557*fp^4*hpp)/(215040*f^4*h^4*Pi^2) - 
   (5951*fp^2*fpp*hpp)/(215040*f^3*h^4*Pi^2) + 
   (39*fpp^2*hpp)/(4480*f^2*h^4*Pi^2) + 
   (201*fp*f3p*hpp)/(17920*f^2*h^4*Pi^2) - (11*f4p*hpp)/(2688*f*h^4*Pi^2) + 
   (391*fp^3*hp*hpp)/(13440*f^3*h^5*Pi^2) - 
   (963*fp*fpp*hp*hpp)/(17920*f^2*h^5*Pi^2) + 
   (187*f3p*hp*hpp)/(7168*f*h^5*Pi^2) + 
   (2089*fp^2*hp^2*hpp)/(43008*f^2*h^6*Pi^2) - 
   (15433*fpp*hp^2*hpp)/(215040*f*h^6*Pi^2) + 
   (8767*fp*hp^3*hpp)/(107520*f*h^7*Pi^2) - 
   (1307*fp^2*hpp^2)/(215040*f^2*h^5*Pi^2) + 
   (121*fpp*hpp^2)/(13440*f*h^5*Pi^2) - 
   (5423*fp*hp*hpp^2)/(215040*f*h^6*Pi^2) - 
   (243*fp^3*h3p)/(71680*f^3*h^4*Pi^2) + 
   (337*fp*fpp*h3p)/(53760*f^2*h^4*Pi^2) - (11*f3p*h3p)/(3584*f*h^4*Pi^2) - 
   (953*fp^2*hp*h3p)/(107520*f^2*h^5*Pi^2) + 
   (473*fpp*hp*h3p)/(35840*f*h^5*Pi^2) - 
   (3949*fp*hp^2*h3p)/(215040*f*h^6*Pi^2) + 
   (33*fp*hpp*h3p)/(7168*f*h^5*Pi^2) + (11*fp^2*h4p)/(13440*f^2*h^4*Pi^2) - 
   (11*fpp*h4p)/(8960*f*h^4*Pi^2) + (143*fp*hp*h4p)/(53760*f*h^5*Pi^2) - 
   (11*fp*h5p)/(53760*f*h^4*Pi^2) - 1/(1260*Pi^2*r^6) + 
   67/(5040*h^3*Pi^2*r^6) - 1/(80*h*Pi^2*r^6) - (3*fp)/(160*f*h^3*Pi^2*r^5) + 
   fp/(64*f*h^2*Pi^2*r^5) + fp/(320*f*h*Pi^2*r^5) + 
   (31*hp)/(672*h^4*Pi^2*r^5) - hp/(64*h^3*Pi^2*r^5) - 
   hp/(320*h^2*Pi^2*r^5) - (223*fp^2)/(23040*f^2*h^3*Pi^2*r^4) + 
   fpp/(120*f*h^3*Pi^2*r^4) + (13*fp^2)/(2304*f^2*h^2*Pi^2*r^4) - 
   (5*fpp)/(576*f*h^2*Pi^2*r^4) - fp^2/(5760*f^2*h*Pi^2*r^4) + 
   fpp/(2880*f*h*Pi^2*r^4) - (19*fp*hp)/(768*f*h^4*Pi^2*r^4) + 
   (5*fp*hp)/(288*f*h^3*Pi^2*r^4) - (fp*hp)/(5760*f*h^2*Pi^2*r^4) + 
   (3973*hp^2)/(53760*h^5*Pi^2*r^4) - (19*hp^2)/(768*h^4*Pi^2*r^4) - 
   (27*hpp)/(1120*h^4*Pi^2*r^4) + hpp/(96*h^3*Pi^2*r^4) - 
   (11*fp^3)/(1920*f^3*h^3*Pi^2*r^3) + (43*fp*fpp)/(6720*f^2*h^3*Pi^2*r^3) + 
   (31*f3p)/(40320*f*h^3*Pi^2*r^3) + fp^3/(384*f^3*h^2*Pi^2*r^3) - 
   (fp*fpp)/(256*f^2*h^2*Pi^2*r^3) + f3p/(1152*f*h^2*Pi^2*r^3) - 
   (659*fp^2*hp)/(80640*f^2*h^4*Pi^2*r^3) - (fpp*hp)/(448*f*h^4*Pi^2*r^3) + 
   (fp^2*hp)/(288*f^2*h^3*Pi^2*r^3) - (7*fpp*hp)/(2304*f*h^3*Pi^2*r^3) + 
   (17*fp*hp^2)/(40320*f*h^5*Pi^2*r^3) + (fp*hp^2)/(144*f*h^4*Pi^2*r^3) + 
   (4967*hp^3)/(80640*h^6*Pi^2*r^3) - (7*hp^3)/(384*h^5*Pi^2*r^3) + 
   (fp*hpp)/(5040*f*h^4*Pi^2*r^3) - (7*fp*hpp)/(2304*f*h^3*Pi^2*r^3) - 
   (65*hp*hpp)/(1344*h^5*Pi^2*r^3) + (13*hp*hpp)/(768*h^4*Pi^2*r^3) + 
   (27*h3p)/(4480*h^4*Pi^2*r^3) - h3p/(384*h^3*Pi^2*r^3) + 
   (37*fp^4)/(21504*f^4*h^3*Pi^2*r^2) - 
   (779*fp^2*fpp)/(107520*f^3*h^3*Pi^2*r^2) + 
   (557*fpp^2)/(161280*f^2*h^3*Pi^2*r^2) + 
   (253*fp*f3p)/(53760*f^2*h^3*Pi^2*r^2) - (101*f4p)/(40320*f*h^3*Pi^2*r^2) - 
   (5*fp^4)/(2304*f^4*h^2*Pi^2*r^2) + (fp^2*fpp)/(192*f^3*h^2*Pi^2*r^2) - 
   fpp^2/(576*f^2*h^2*Pi^2*r^2) - (5*fp*f3p)/(2304*f^2*h^2*Pi^2*r^2) + 
   f4p/(1152*f*h^2*Pi^2*r^2) + (683*fp^3*hp)/(107520*f^3*h^4*Pi^2*r^2) - 
   (2911*fp*fpp*hp)/(161280*f^2*h^4*Pi^2*r^2) + 
   (137*f3p*hp)/(10752*f*h^4*Pi^2*r^2) - (fp^3*hp)/(384*f^3*h^3*Pi^2*r^2) + 
   (23*fp*fpp*hp)/(4608*f^2*h^3*Pi^2*r^2) - (f3p*hp)/(384*f*h^3*Pi^2*r^2) + 
   (1031*fp^2*hp^2)/(64512*f^2*h^5*Pi^2*r^2) - 
   (3551*fpp*hp^2)/(107520*f*h^5*Pi^2*r^2) - 
   (fp^2*hp^2)/(384*f^2*h^4*Pi^2*r^2) + (19*fpp*hp^2)/(4608*f*h^4*Pi^2*r^2) + 
   (5161*fp*hp^3)/(107520*f*h^6*Pi^2*r^2) - 
   (7*fp*hp^3)/(2304*f*h^5*Pi^2*r^2) - (121*hp^4)/(7680*h^7*Pi^2*r^2) - 
   (1637*fp^2*hpp)/(322560*f^2*h^4*Pi^2*r^2) + 
   (281*fpp*hpp)/(26880*f*h^4*Pi^2*r^2) + 
   (5*fp^2*hpp)/(4608*f^2*h^3*Pi^2*r^2) - (fpp*hpp)/(576*f*h^3*Pi^2*r^2) - 
   (633*fp*hp*hpp)/(17920*f*h^5*Pi^2*r^2) + 
   (13*fp*hp*hpp)/(4608*f*h^4*Pi^2*r^2) + (77*hp^2*hpp)/(5120*h^6*Pi^2*r^2) - 
   (11*hpp^2)/(7680*h^5*Pi^2*r^2) + (641*fp*h3p)/(161280*f*h^4*Pi^2*r^2) - 
   (fp*h3p)/(2304*f*h^3*Pi^2*r^2) - (11*hp*h3p)/(7680*h^5*Pi^2*r^2) + 
   (1343*fp^5)/(46080*f^5*h^3*Pi^2*r) - 
   (4279*fp^3*fpp)/(53760*f^4*h^3*Pi^2*r) + 
   (2873*fp*fpp^2)/(64512*f^3*h^3*Pi^2*r) + 
   (9467*fp^2*f3p)/(322560*f^3*h^3*Pi^2*r) - 
   (397*fpp*f3p)/(23040*f^2*h^3*Pi^2*r) - (19*fp*f4p)/(2304*f^2*h^3*Pi^2*r) + 
   (11*f5p)/(5376*f*h^3*Pi^2*r) + (3085*fp^4*hp)/(64512*f^4*h^4*Pi^2*r) - 
   (34229*fp^2*fpp*hp)/(322560*f^3*h^4*Pi^2*r) + 
   (473*fpp^2*hp)/(15360*f^2*h^4*Pi^2*r) + 
   (917*fp*f3p*hp)/(23040*f^2*h^4*Pi^2*r) - (97*f4p*hp)/(8064*f*h^4*Pi^2*r) + 
   (18089*fp^3*hp^2)/(322560*f^3*h^5*Pi^2*r) - 
   (2189*fp*fpp*hp^2)/(23040*f^2*h^5*Pi^2*r) + 
   (12091*f3p*hp^2)/(322560*f*h^5*Pi^2*r) + 
   (575*fp^2*hp^3)/(9216*f^2*h^6*Pi^2*r) - 
   (24167*fpp*hp^3)/(322560*f*h^6*Pi^2*r) + 
   (3167*fp*hp^4)/(32256*f*h^7*Pi^2*r) - (299*hp^5)/(1536*h^8*Pi^2*r) - 
   (5833*fp^3*hpp)/(322560*f^3*h^4*Pi^2*r) + 
   (353*fp*fpp*hpp)/(11520*f^2*h^4*Pi^2*r) - 
   (1969*f3p*hpp)/(161280*f*h^4*Pi^2*r) - 
   (221*fp^2*hp*hpp)/(4608*f^2*h^5*Pi^2*r) + 
   (2351*fpp*hp*hpp)/(40320*f*h^5*Pi^2*r) - 
   (13429*fp*hp^2*hpp)/(107520*f*h^6*Pi^2*r) + 
   (18331*hp^3*hpp)/(53760*h^7*Pi^2*r) + 
   (1031*fp*hpp^2)/(64512*f*h^5*Pi^2*r) - 
   (11339*hp*hpp^2)/(107520*h^6*Pi^2*r) + 
   (89*fp^2*h3p)/(15360*f^2*h^4*Pi^2*r) - 
   (1157*fpp*h3p)/(161280*f*h^4*Pi^2*r) + 
   (769*fp*hp*h3p)/(32256*f*h^5*Pi^2*r) - 
   (8257*hp^2*h3p)/(107520*h^6*Pi^2*r) + (69*hpp*h3p)/(3584*h^5*Pi^2*r) - 
   (31*fp*h4p)/(13440*f*h^4*Pi^2*r) + (299*hp*h4p)/(26880*h^5*Pi^2*r) - 
   (23*h5p)/(26880*h^4*Pi^2*r) + (205*fp^6*si)/(1024*f^6*h^3*Pi^2) - 
   (343*fp^4*fpp*si)/(512*f^5*h^3*Pi^2) + 
   (739*fp^2*fpp^2*si)/(1280*f^4*h^3*Pi^2) - 
   (49*fpp^3*si)/(640*f^3*h^3*Pi^2) + (961*fp^3*f3p*si)/(3840*f^4*h^3*Pi^2) - 
   (583*fp*fpp*f3p*si)/(1920*f^3*h^3*Pi^2) + (f3p^2*si)/(30*f^2*h^3*Pi^2) - 
   (9*fp^2*f4p*si)/(128*f^3*h^3*Pi^2) + (3*fpp*f4p*si)/(64*f^2*h^3*Pi^2) + 
   (fp*f5p*si)/(60*f^2*h^3*Pi^2) - (f6p*si)/(240*f*h^3*Pi^2) + 
   (343*fp^5*hp*si)/(1024*f^5*h^4*Pi^2) - 
   (2439*fp^3*fpp*hp*si)/(2560*f^4*h^4*Pi^2) + 
   (73*fp*fpp^2*hp*si)/(128*f^3*h^4*Pi^2) + 
   (1393*fp^2*f3p*hp*si)/(3840*f^3*h^4*Pi^2) - 
   (77*fpp*f3p*hp*si)/(320*f^2*h^4*Pi^2) - 
   (41*fp*f4p*hp*si)/(384*f^2*h^4*Pi^2) + (f5p*hp*si)/(32*f*h^4*Pi^2) + 
   (6061*fp^4*hp^2*si)/(15360*f^4*h^5*Pi^2) - 
   (7087*fp^2*fpp*hp^2*si)/(7680*f^3*h^5*Pi^2) + 
   (381*fpp^2*hp^2*si)/(1280*f^2*h^5*Pi^2) + 
   (703*fp*f3p*hp^2*si)/(1920*f^2*h^5*Pi^2) - 
   (25*f4p*hp^2*si)/(192*f*h^5*Pi^2) + 
   (6259*fp^3*hp^3*si)/(15360*f^3*h^6*Pi^2) - 
   (1921*fp*fpp*hp^3*si)/(2560*f^2*h^6*Pi^2) + 
   (45*f3p*hp^3*si)/(128*f*h^6*Pi^2) + 
   (521*fp^2*hp^4*si)/(1280*f^2*h^7*Pi^2) - 
   (113*fpp*hp^4*si)/(192*f*h^7*Pi^2) + (91*fp*hp^5*si)/(192*f*h^8*Pi^2) - 
   (961*fp^4*hpp*si)/(7680*f^4*h^4*Pi^2) + 
   (1123*fp^2*fpp*hpp*si)/(3840*f^3*h^4*Pi^2) - 
   (3*fpp^2*hpp*si)/(32*f^2*h^4*Pi^2) - (7*fp*f3p*hpp*si)/(60*f^2*h^4*Pi^2) + 
   (f4p*hpp*si)/(24*f*h^4*Pi^2) - (1169*fp^3*hp*hpp*si)/(3840*f^3*h^5*Pi^2) + 
   (719*fp*fpp*hp*hpp*si)/(1280*f^2*h^5*Pi^2) - 
   (17*f3p*hp*hpp*si)/(64*f*h^5*Pi^2) - 
   (1283*fp^2*hp^2*hpp*si)/(2560*f^2*h^6*Pi^2) + 
   (1403*fpp*hp^2*hpp*si)/(1920*f*h^6*Pi^2) - 
   (797*fp*hp^3*hpp*si)/(960*f*h^7*Pi^2) + 
   (fp^2*hpp^2*si)/(16*f^2*h^5*Pi^2) - (11*fpp*hpp^2*si)/(120*f*h^5*Pi^2) + 
   (493*fp*hp*hpp^2*si)/(1920*f*h^6*Pi^2) + 
   (9*fp^3*h3p*si)/(256*f^3*h^4*Pi^2) - 
   (25*fp*fpp*h3p*si)/(384*f^2*h^4*Pi^2) + (f3p*h3p*si)/(32*f*h^4*Pi^2) + 
   (349*fp^2*hp*h3p*si)/(3840*f^2*h^5*Pi^2) - 
   (43*fpp*hp*h3p*si)/(320*f*h^5*Pi^2) + 
   (359*fp*hp^2*h3p*si)/(1920*f*h^6*Pi^2) - 
   (3*fp*hpp*h3p*si)/(64*f*h^5*Pi^2) - (fp^2*h4p*si)/(120*f^2*h^4*Pi^2) + 
   (fpp*h4p*si)/(80*f*h^4*Pi^2) - (13*fp*hp*h4p*si)/(480*f*h^5*Pi^2) + 
   (fp*h5p*si)/(480*f*h^4*Pi^2) + si/(120*Pi^2*r^6) - si/(15*h^3*Pi^2*r^6) - 
   si/(8*h^2*Pi^2*r^6) + (11*si)/(60*h*Pi^2*r^6) + 
   (127*fp*si)/(480*f*h^3*Pi^2*r^5) - (7*fp*si)/(32*f*h^2*Pi^2*r^5) - 
   (11*fp*si)/(240*f*h*Pi^2*r^5) - (15*hp*si)/(32*h^4*Pi^2*r^5) + 
   (5*hp*si)/(32*h^3*Pi^2*r^5) + (11*hp*si)/(240*h^2*Pi^2*r^5) + 
   (89*fp^2*si)/(640*f^2*h^3*Pi^2*r^4) - (43*fpp*si)/(360*f*h^3*Pi^2*r^4) - 
   (31*fp^2*si)/(384*f^2*h^2*Pi^2*r^4) + (fpp*si)/(8*f*h^2*Pi^2*r^4) + 
   (fp^2*si)/(360*f^2*h*Pi^2*r^4) - (fpp*si)/(180*f*h*Pi^2*r^4) + 
   (971*fp*hp*si)/(2880*f*h^4*Pi^2*r^4) - 
   (49*fp*hp*si)/(192*f*h^3*Pi^2*r^4) + (fp*hp*si)/(360*f*h^2*Pi^2*r^4) - 
   (1633*hp^2*si)/(1920*h^5*Pi^2*r^4) + (133*hp^2*si)/(384*h^4*Pi^2*r^4) + 
   (67*hpp*si)/(240*h^4*Pi^2*r^4) - (7*hpp*si)/(48*h^3*Pi^2*r^4) + 
   (943*fp^3*si)/(11520*f^3*h^3*Pi^2*r^3) - 
   (35*fp*fpp*si)/(384*f^2*h^3*Pi^2*r^3) - (11*f3p*si)/(960*f*h^3*Pi^2*r^3) - 
   (fp^3*si)/(32*f^3*h^2*Pi^2*r^3) + (17*fp*fpp*si)/(384*f^2*h^2*Pi^2*r^3) - 
   (f3p*si)/(192*f*h^2*Pi^2*r^3) + 
   (1289*fp^2*hp*si)/(11520*f^2*h^4*Pi^2*r^3) + 
   (83*fpp*hp*si)/(1920*f*h^4*Pi^2*r^3) - 
   (fp^2*hp*si)/(24*f^2*h^3*Pi^2*r^3) + (11*fpp*hp*si)/(384*f*h^3*Pi^2*r^3) - 
   (551*fp*hp^2*si)/(11520*f*h^5*Pi^2*r^3) - 
   (17*fp*hp^2*si)/(192*f*h^4*Pi^2*r^3) - 
   (2651*hp^3*si)/(3840*h^6*Pi^2*r^3) + (49*hp^3*si)/(192*h^5*Pi^2*r^3) + 
   (59*fp*hpp*si)/(5760*f*h^4*Pi^2*r^3) + 
   (5*fp*hpp*si)/(128*f*h^3*Pi^2*r^3) + (211*hp*hpp*si)/(384*h^5*Pi^2*r^3) - 
   (91*hp*hpp*si)/(384*h^4*Pi^2*r^3) - (67*h3p*si)/(960*h^4*Pi^2*r^3) + 
   (7*h3p*si)/(192*h^3*Pi^2*r^3) - (fp^4*si)/(64*f^4*h^3*Pi^2*r^2) + 
   (157*fp^2*fpp*si)/(1920*f^3*h^3*Pi^2*r^2) - 
   (121*fpp^2*si)/(2880*f^2*h^3*Pi^2*r^2) - 
   (11*fp*f3p*si)/(192*f^2*h^3*Pi^2*r^2) + (11*f4p*si)/(360*f*h^3*Pi^2*r^2) + 
   (5*fp^4*si)/(144*f^4*h^2*Pi^2*r^2) - (fp^2*fpp*si)/(12*f^3*h^2*Pi^2*r^2) + 
   (fpp^2*si)/(36*f^2*h^2*Pi^2*r^2) + (5*fp*f3p*si)/(144*f^2*h^2*Pi^2*r^2) - 
   (f4p*si)/(72*f*h^2*Pi^2*r^2) - (3*fp^3*hp*si)/(40*f^3*h^4*Pi^2*r^2) + 
   (65*fp*fpp*hp*si)/(288*f^2*h^4*Pi^2*r^2) - 
   (51*f3p*hp*si)/(320*f*h^4*Pi^2*r^2) + (fp^3*hp*si)/(24*f^3*h^3*Pi^2*r^2) - 
   (23*fp*fpp*hp*si)/(288*f^2*h^3*Pi^2*r^2) + 
   (f3p*hp*si)/(24*f*h^3*Pi^2*r^2) - 
   (299*fp^2*hp^2*si)/(1440*f^2*h^5*Pi^2*r^2) + 
   (823*fpp*hp^2*si)/(1920*f*h^5*Pi^2*r^2) + 
   (fp^2*hp^2*si)/(24*f^2*h^4*Pi^2*r^2) - 
   (19*fpp*hp^2*si)/(288*f*h^4*Pi^2*r^2) - 
   (103*fp*hp^3*si)/(160*f*h^6*Pi^2*r^2) + 
   (7*fp*hp^3*si)/(144*f*h^5*Pi^2*r^2) + (77*hp^4*si)/(320*h^7*Pi^2*r^2) + 
   (379*fp^2*hpp*si)/(5760*f^2*h^4*Pi^2*r^2) - 
   (13*fpp*hpp*si)/(96*f*h^4*Pi^2*r^2) - 
   (5*fp^2*hpp*si)/(288*f^2*h^3*Pi^2*r^2) + 
   (fpp*hpp*si)/(36*f*h^3*Pi^2*r^2) + 
   (113*fp*hp*hpp*si)/(240*f*h^5*Pi^2*r^2) - 
   (13*fp*hp*hpp*si)/(288*f*h^4*Pi^2*r^2) - 
   (147*hp^2*hpp*si)/(640*h^6*Pi^2*r^2) + (7*hpp^2*si)/(320*h^5*Pi^2*r^2) - 
   (151*fp*h3p*si)/(2880*f*h^4*Pi^2*r^2) + (fp*h3p*si)/(144*f*h^3*Pi^2*r^2) + 
   (7*hp*h3p*si)/(320*h^5*Pi^2*r^2) - (7427*fp^5*si)/(23040*f^5*h^3*Pi^2*r) + 
   (2267*fp^3*fpp*si)/(2560*f^4*h^3*Pi^2*r) - 
   (1153*fp*fpp^2*si)/(2304*f^3*h^3*Pi^2*r) - 
   (467*fp^2*f3p*si)/(1440*f^3*h^3*Pi^2*r) + 
   (1117*fpp*f3p*si)/(5760*f^2*h^3*Pi^2*r) + 
   (511*fp*f4p*si)/(5760*f^2*h^3*Pi^2*r) - (f5p*si)/(48*f*h^3*Pi^2*r) - 
   (2383*fp^4*hp*si)/(4608*f^4*h^4*Pi^2*r) + 
   (26471*fp^2*fpp*hp*si)/(23040*f^3*h^4*Pi^2*r) - 
   (1291*fpp^2*hp*si)/(3840*f^2*h^4*Pi^2*r) - 
   (2419*fp*f3p*hp*si)/(5760*f^2*h^4*Pi^2*r) + 
   (15*f4p*hp*si)/(128*f*h^4*Pi^2*r) - 
   (13409*fp^3*hp^2*si)/(23040*f^3*h^5*Pi^2*r) + 
   (22337*fp*fpp*hp^2*si)/(23040*f^2*h^5*Pi^2*r) - 
   (661*f3p*hp^2*si)/(1920*f*h^5*Pi^2*r) - 
   (313*fp^2*hp^3*si)/(512*f^2*h^6*Pi^2*r) + 
   (1701*fpp*hp^3*si)/(2560*f*h^6*Pi^2*r) - 
   (353*fp*hp^4*si)/(384*f*h^7*Pi^2*r) + (91*hp^5*si)/(48*h^8*Pi^2*r) + 
   (433*fp^3*hpp*si)/(2304*f^3*h^4*Pi^2*r) - 
   (3611*fp*fpp*hpp*si)/(11520*f^2*h^4*Pi^2*r) + 
   (109*f3p*hpp*si)/(960*f*h^4*Pi^2*r) + 
   (911*fp^2*hp*hpp*si)/(1920*f^2*h^5*Pi^2*r) - 
   (677*fpp*hp*hpp*si)/(1280*f*h^5*Pi^2*r) + 
   (4559*fp*hp^2*hpp*si)/(3840*f*h^6*Pi^2*r) - 
   (797*hp^3*hpp*si)/(240*h^7*Pi^2*r) - (59*fp*hpp^2*si)/(384*f*h^5*Pi^2*r) + 
   (493*hp*hpp^2*si)/(480*h^6*Pi^2*r) - 
   (67*fp^2*h3p*si)/(1152*f^2*h^4*Pi^2*r) + 
   (43*fpp*h3p*si)/(640*f*h^4*Pi^2*r) - 
   (89*fp*hp*h3p*si)/(384*f*h^5*Pi^2*r) + 
   (359*hp^2*h3p*si)/(480*h^6*Pi^2*r) - (3*hpp*h3p*si)/(16*h^5*Pi^2*r) + 
   (11*fp*h4p*si)/(480*f*h^4*Pi^2*r) - (13*hp*h4p*si)/(120*h^5*Pi^2*r) + 
   (h5p*si)/(120*h^4*Pi^2*r) - (1663*fp^6*si^2)/(3072*f^6*h^3*Pi^2) + 
   (2837*fp^4*fpp*si^2)/(1536*f^5*h^3*Pi^2) - 
   (1261*fp^2*fpp^2*si^2)/(768*f^4*h^3*Pi^2) + 
   (29*fpp^3*si^2)/(128*f^3*h^3*Pi^2) - 
   (263*fp^3*f3p*si^2)/(384*f^4*h^3*Pi^2) + 
   (337*fp*fpp*f3p*si^2)/(384*f^3*h^3*Pi^2) - 
   (13*f3p^2*si^2)/(128*f^2*h^3*Pi^2) + 
   (71*fp^2*f4p*si^2)/(384*f^3*h^3*Pi^2) - 
   (13*fpp*f4p*si^2)/(96*f^2*h^3*Pi^2) - (fp*f5p*si^2)/(24*f^2*h^3*Pi^2) + 
   (f6p*si^2)/(96*f*h^3*Pi^2) - (2837*fp^5*hp*si^2)/(3072*f^5*h^4*Pi^2) + 
   (1025*fp^3*fpp*hp*si^2)/(384*f^4*h^4*Pi^2) - 
   (53*fp*fpp^2*hp*si^2)/(32*f^3*h^4*Pi^2) - 
   (763*fp^2*f3p*hp*si^2)/(768*f^3*h^4*Pi^2) + 
   (91*fpp*f3p*hp*si^2)/(128*f^2*h^4*Pi^2) + 
   (53*fp*f4p*hp*si^2)/(192*f^2*h^4*Pi^2) - (5*f5p*hp*si^2)/(64*f*h^4*Pi^2) - 
   (3365*fp^4*hp^2*si^2)/(3072*f^4*h^5*Pi^2) + 
   (1323*fp^2*fpp*hp^2*si^2)/(512*f^3*h^5*Pi^2) - 
   (1339*fpp^2*hp^2*si^2)/(1536*f^2*h^5*Pi^2) - 
   (47*fp*f3p*hp^2*si^2)/(48*f^2*h^5*Pi^2) + 
   (125*f4p*hp^2*si^2)/(384*f*h^5*Pi^2) - 
   (1141*fp^3*hp^3*si^2)/(1024*f^3*h^6*Pi^2) + 
   (1565*fp*fpp*hp^3*si^2)/(768*f^2*h^6*Pi^2) - 
   (225*f3p*hp^3*si^2)/(256*f*h^6*Pi^2) - 
   (205*fp^2*hp^4*si^2)/(192*f^2*h^7*Pi^2) + 
   (565*fpp*hp^4*si^2)/(384*f*h^7*Pi^2) - 
   (455*fp*hp^5*si^2)/(384*f*h^8*Pi^2) + 
   (263*fp^4*hpp*si^2)/(768*f^4*h^4*Pi^2) - 
   (207*fp^2*fpp*hpp*si^2)/(256*f^3*h^4*Pi^2) + 
   (13*fpp^2*hpp*si^2)/(48*f^2*h^4*Pi^2) + 
   (119*fp*f3p*hpp*si^2)/(384*f^2*h^4*Pi^2) - 
   (5*f4p*hpp*si^2)/(48*f*h^4*Pi^2) + 
   (105*fp^3*hp*hpp*si^2)/(128*f^3*h^5*Pi^2) - 
   (1159*fp*fpp*hp*hpp*si^2)/(768*f^2*h^5*Pi^2) + 
   (85*f3p*hp*hpp*si^2)/(128*f*h^5*Pi^2) + 
   (333*fp^2*hp^2*hpp*si^2)/(256*f^2*h^6*Pi^2) - 
   (1403*fpp*hp^2*hpp*si^2)/(768*f*h^6*Pi^2) + 
   (797*fp*hp^3*hpp*si^2)/(384*f*h^7*Pi^2) - 
   (247*fp^2*hpp^2*si^2)/(1536*f^2*h^5*Pi^2) + 
   (11*fpp*hpp^2*si^2)/(48*f*h^5*Pi^2) - 
   (493*fp*hp*hpp^2*si^2)/(768*f*h^6*Pi^2) - 
   (71*fp^3*h3p*si^2)/(768*f^3*h^4*Pi^2) + 
   (11*fp*fpp*h3p*si^2)/(64*f^2*h^4*Pi^2) - 
   (5*f3p*h3p*si^2)/(64*f*h^4*Pi^2) - 
   (89*fp^2*hp*h3p*si^2)/(384*f^2*h^5*Pi^2) + 
   (43*fpp*hp*h3p*si^2)/(128*f*h^5*Pi^2) - 
   (359*fp*hp^2*h3p*si^2)/(768*f*h^6*Pi^2) + 
   (15*fp*hpp*h3p*si^2)/(128*f*h^5*Pi^2) + 
   (fp^2*h4p*si^2)/(48*f^2*h^4*Pi^2) - (fpp*h4p*si^2)/(32*f*h^4*Pi^2) + 
   (13*fp*hp*h4p*si^2)/(192*f*h^5*Pi^2) - (fp*h5p*si^2)/(192*f*h^4*Pi^2) - 
   si^2/(24*Pi^2*r^6) - (3*si^2)/(8*h^3*Pi^2*r^6) + 
   (11*si^2)/(8*h^2*Pi^2*r^6) - (23*si^2)/(24*h*Pi^2*r^6) - 
   (39*fp*si^2)/(32*f*h^3*Pi^2*r^5) + (47*fp*si^2)/(48*f*h^2*Pi^2*r^5) + 
   (23*fp*si^2)/(96*f*h*Pi^2*r^5) + (39*hp*si^2)/(32*h^4*Pi^2*r^5) - 
   (5*hp*si^2)/(16*h^3*Pi^2*r^5) - (23*hp*si^2)/(96*h^2*Pi^2*r^5) - 
   (127*fp^2*si^2)/(192*f^2*h^3*Pi^2*r^4) + 
   (53*fpp*si^2)/(96*f*h^3*Pi^2*r^4) + (3*fp^2*si^2)/(8*f^2*h^2*Pi^2*r^4) - 
   (7*fpp*si^2)/(12*f*h^2*Pi^2*r^4) - (fp^2*si^2)/(64*f^2*h*Pi^2*r^4) + 
   (fpp*si^2)/(32*f*h*Pi^2*r^4) - (289*fp*hp*si^2)/(192*f*h^4*Pi^2*r^4) + 
   (59*fp*hp*si^2)/(48*f*h^3*Pi^2*r^4) - (fp*hp*si^2)/(64*f*h^2*Pi^2*r^4) + 
   (3*hp^2*si^2)/(h^5*Pi^2*r^4) - (19*hp^2*si^2)/(12*h^4*Pi^2*r^4) - 
   (hpp*si^2)/(h^4*Pi^2*r^4) + (2*hpp*si^2)/(3*h^3*Pi^2*r^4) - 
   (19*fp^3*si^2)/(48*f^3*h^3*Pi^2*r^3) + 
   (7*fp*fpp*si^2)/(16*f^2*h^3*Pi^2*r^3) + (f3p*si^2)/(16*f*h^3*Pi^2*r^3) + 
   (7*fp^3*si^2)/(64*f^3*h^2*Pi^2*r^3) - 
   (13*fp*fpp*si^2)/(96*f^2*h^2*Pi^2*r^3) - (f3p*si^2)/(48*f*h^2*Pi^2*r^3) - 
   (49*fp^2*hp*si^2)/(96*f^2*h^4*Pi^2*r^3) - 
   (9*fpp*hp*si^2)/(32*f*h^4*Pi^2*r^3) + 
   (7*fp^2*hp*si^2)/(48*f^2*h^3*Pi^2*r^3) - 
   (fpp*hp*si^2)/(24*f*h^3*Pi^2*r^3) + 
   (43*fp*hp^2*si^2)/(96*f*h^5*Pi^2*r^3) + 
   (67*fp*hp^2*si^2)/(192*f*h^4*Pi^2*r^3) + 
   (109*hp^3*si^2)/(48*h^6*Pi^2*r^3) - (7*hp^3*si^2)/(6*h^5*Pi^2*r^3) - 
   (11*fp*hpp*si^2)/(96*f*h^4*Pi^2*r^3) - 
   (5*fp*hpp*si^2)/(32*f*h^3*Pi^2*r^3) - (15*hp*hpp*si^2)/(8*h^5*Pi^2*r^3) + 
   (13*hp*hpp*si^2)/(12*h^4*Pi^2*r^3) + (h3p*si^2)/(4*h^4*Pi^2*r^3) - 
   (h3p*si^2)/(6*h^3*Pi^2*r^3) - (7*fp^2*fpp*si^2)/(32*f^3*h^3*Pi^2*r^2) + 
   (7*fpp^2*si^2)/(48*f^2*h^3*Pi^2*r^2) + 
   (13*fp*f3p*si^2)/(64*f^2*h^3*Pi^2*r^2) - 
   (11*f4p*si^2)/(96*f*h^3*Pi^2*r^2) - 
   (35*fp^4*si^2)/(192*f^4*h^2*Pi^2*r^2) + 
   (7*fp^2*fpp*si^2)/(16*f^3*h^2*Pi^2*r^2) - 
   (7*fpp^2*si^2)/(48*f^2*h^2*Pi^2*r^2) - 
   (35*fp*f3p*si^2)/(192*f^2*h^2*Pi^2*r^2) + 
   (7*f4p*si^2)/(96*f*h^2*Pi^2*r^2) + 
   (15*fp^3*hp*si^2)/(64*f^3*h^4*Pi^2*r^2) - 
   (325*fp*fpp*hp*si^2)/(384*f^2*h^4*Pi^2*r^2) + 
   (5*f3p*hp*si^2)/(8*f*h^4*Pi^2*r^2) - 
   (7*fp^3*hp*si^2)/(32*f^3*h^3*Pi^2*r^2) + 
   (161*fp*fpp*hp*si^2)/(384*f^2*h^3*Pi^2*r^2) - 
   (7*f3p*hp*si^2)/(32*f*h^3*Pi^2*r^2) + 
   (163*fp^2*hp^2*si^2)/(192*f^2*h^5*Pi^2*r^2) - 
   (229*fpp*hp^2*si^2)/(128*f*h^5*Pi^2*r^2) - 
   (7*fp^2*hp^2*si^2)/(32*f^2*h^4*Pi^2*r^2) + 
   (133*fpp*hp^2*si^2)/(384*f*h^4*Pi^2*r^2) + 
   (91*fp*hp^3*si^2)/(32*f*h^6*Pi^2*r^2) - 
   (49*fp*hp^3*si^2)/(192*f*h^5*Pi^2*r^2) - (11*hp^4*si^2)/(8*h^7*Pi^2*r^2) - 
   (103*fp^2*hpp*si^2)/(384*f^2*h^4*Pi^2*r^2) + 
   (9*fpp*hpp*si^2)/(16*f*h^4*Pi^2*r^2) + 
   (35*fp^2*hpp*si^2)/(384*f^2*h^3*Pi^2*r^2) - 
   (7*fpp*hpp*si^2)/(48*f*h^3*Pi^2*r^2) - 
   (263*fp*hp*hpp*si^2)/(128*f*h^5*Pi^2*r^2) + 
   (91*fp*hp*hpp*si^2)/(384*f*h^4*Pi^2*r^2) + 
   (21*hp^2*hpp*si^2)/(16*h^6*Pi^2*r^2) - (hpp^2*si^2)/(8*h^5*Pi^2*r^2) + 
   (43*fp*h3p*si^2)/(192*f*h^4*Pi^2*r^2) - 
   (7*fp*h3p*si^2)/(192*f*h^3*Pi^2*r^2) - (hp*h3p*si^2)/(8*h^5*Pi^2*r^2) + 
   (1561*fp^5*si^2)/(1536*f^5*h^3*Pi^2*r) - 
   (1081*fp^3*fpp*si^2)/(384*f^4*h^3*Pi^2*r) + 
   (619*fp*fpp^2*si^2)/(384*f^3*h^3*Pi^2*r) + 
   (389*fp^2*f3p*si^2)/(384*f^3*h^3*Pi^2*r) - 
   (119*fpp*f3p*si^2)/(192*f^2*h^3*Pi^2*r) - 
   (25*fp*f4p*si^2)/(96*f^2*h^3*Pi^2*r) + (5*f5p*si^2)/(96*f*h^3*Pi^2*r) + 
   (795*fp^4*hp*si^2)/(512*f^4*h^4*Pi^2*r) - 
   (2641*fp^2*fpp*hp*si^2)/(768*f^3*h^4*Pi^2*r) + 
   (97*fpp^2*hp*si^2)/(96*f^2*h^4*Pi^2*r) + 
   (113*fp*f3p*hp*si^2)/(96*f^2*h^4*Pi^2*r) - 
   (49*f4p*hp*si^2)/(192*f*h^4*Pi^2*r) + 
   (2435*fp^3*hp^2*si^2)/(1536*f^3*h^5*Pi^2*r) - 
   (5*fp*fpp*hp^2*si^2)/(2*f^2*h^5*Pi^2*r) + 
   (229*f3p*hp^2*si^2)/(384*f*h^5*Pi^2*r) + 
   (2183*fp^2*hp^3*si^2)/(1536*f^2*h^6*Pi^2*r) - 
   (761*fpp*hp^3*si^2)/(768*f*h^6*Pi^2*r) + 
   (241*fp*hp^4*si^2)/(128*f*h^7*Pi^2*r) - (455*hp^5*si^2)/(96*h^8*Pi^2*r) - 
   (395*fp^3*hpp*si^2)/(768*f^3*h^4*Pi^2*r) + 
   (313*fp*fpp*hpp*si^2)/(384*f^2*h^4*Pi^2*r) - 
   (5*f3p*hpp*si^2)/(24*f*h^4*Pi^2*r) - 
   (437*fp^2*hp*hpp*si^2)/(384*f^2*h^5*Pi^2*r) + 
   (341*fpp*hp*hpp*si^2)/(384*f*h^5*Pi^2*r) - 
   (1975*fp*hp^2*hpp*si^2)/(768*f*h^6*Pi^2*r) + 
   (797*hp^3*hpp*si^2)/(96*h^7*Pi^2*r) + 
   (133*fp*hpp^2*si^2)/(384*f*h^5*Pi^2*r) - 
   (493*hp*hpp^2*si^2)/(192*h^6*Pi^2*r) + 
   (7*fp^2*h3p*si^2)/(48*f^2*h^4*Pi^2*r) - 
   (25*fpp*h3p*si^2)/(192*f*h^4*Pi^2*r) + 
   (13*fp*hp*h3p*si^2)/(24*f*h^5*Pi^2*r) - 
   (359*hp^2*h3p*si^2)/(192*h^6*Pi^2*r) + (15*hpp*h3p*si^2)/(32*h^5*Pi^2*r) - 
   (11*fp*h4p*si^2)/(192*f*h^4*Pi^2*r) + (13*hp*h4p*si^2)/(48*h^5*Pi^2*r) - 
   (h5p*si^2)/(48*h^4*Pi^2*r) + (215*fp^6*si^3)/(1536*f^6*h^3*Pi^2) - 
   (151*fp^4*fpp*si^3)/(256*f^5*h^3*Pi^2) + 
   (87*fp^2*fpp^2*si^3)/(128*f^4*h^3*Pi^2) - 
   (23*fpp^3*si^3)/(192*f^3*h^3*Pi^2) + 
   (13*fp^3*f3p*si^3)/(64*f^4*h^3*Pi^2) - 
   (13*fp*fpp*f3p*si^3)/(32*f^3*h^3*Pi^2) + (f3p^2*si^3)/(16*f^2*h^3*Pi^2) - 
   (fp^2*f4p*si^3)/(32*f^3*h^3*Pi^2) + (fpp*f4p*si^3)/(16*f^2*h^3*Pi^2) + 
   (151*fp^5*hp*si^3)/(512*f^5*h^4*Pi^2) - 
   (63*fp^3*fpp*hp*si^3)/(64*f^4*h^4*Pi^2) + 
   (101*fp*fpp^2*hp*si^3)/(128*f^3*h^4*Pi^2) + 
   (19*fp^2*f3p*hp*si^3)/(64*f^3*h^4*Pi^2) - 
   (3*fpp*f3p*hp*si^3)/(8*f^2*h^4*Pi^2) - 
   (fp*f4p*hp*si^3)/(32*f^2*h^4*Pi^2) + 
   (191*fp^4*hp^2*si^3)/(512*f^4*h^5*Pi^2) - 
   (243*fp^2*fpp*hp^2*si^3)/(256*f^3*h^5*Pi^2) + 
   (7*fpp^2*hp^2*si^3)/(16*f^2*h^5*Pi^2) + 
   (7*fp*f3p*hp^2*si^3)/(32*f^2*h^5*Pi^2) + 
   (503*fp^3*hp^3*si^3)/(1536*f^3*h^6*Pi^2) - 
   (71*fp*fpp*hp^3*si^3)/(128*f^2*h^6*Pi^2) + 
   (11*fp^2*hp^4*si^3)/(64*f^2*h^7*Pi^2) - 
   (13*fp^4*hpp*si^3)/(128*f^4*h^4*Pi^2) + 
   (17*fp^2*fpp*hpp*si^3)/(64*f^3*h^4*Pi^2) - 
   (fpp^2*hpp*si^3)/(8*f^2*h^4*Pi^2) - (fp*f3p*hpp*si^3)/(16*f^2*h^4*Pi^2) - 
   (13*fp^3*hp*hpp*si^3)/(64*f^3*h^5*Pi^2) + 
   (23*fp*fpp*hp*hpp*si^3)/(64*f^2*h^5*Pi^2) - 
   (21*fp^2*hp^2*hpp*si^3)/(128*f^2*h^6*Pi^2) + 
   (fp^2*hpp^2*si^3)/(64*f^2*h^5*Pi^2) + (fp^3*h3p*si^3)/(64*f^3*h^4*Pi^2) - 
   (fp*fpp*h3p*si^3)/(32*f^2*h^4*Pi^2) + 
   (fp^2*hp*h3p*si^3)/(64*f^2*h^5*Pi^2) + si^3/(12*Pi^2*r^6) + 
   (23*si^3)/(12*h^3*Pi^2*r^6) - (15*si^3)/(4*h^2*Pi^2*r^6) + 
   (7*si^3)/(4*h*Pi^2*r^6) + (29*fp*si^3)/(16*f*h^3*Pi^2*r^5) - 
   (11*fp*si^3)/(8*f*h^2*Pi^2*r^5) - (7*fp*si^3)/(16*f*h*Pi^2*r^5) - 
   (hp*si^3)/(16*h^4*Pi^2*r^5) - (3*hp*si^3)/(8*h^3*Pi^2*r^5) + 
   (7*hp*si^3)/(16*h^2*Pi^2*r^5) + (33*fp^2*si^3)/(32*f^2*h^3*Pi^2*r^4) - 
   (13*fpp*si^3)/(16*f*h^3*Pi^2*r^4) - (9*fp^2*si^3)/(16*f^2*h^2*Pi^2*r^4) + 
   (7*fpp*si^3)/(8*f*h^2*Pi^2*r^4) + (fp^2*si^3)/(32*f^2*h*Pi^2*r^4) - 
   (fpp*si^3)/(16*f*h*Pi^2*r^4) + (69*fp*hp*si^3)/(32*f*h^4*Pi^2*r^4) - 
   (31*fp*hp*si^3)/(16*f*h^3*Pi^2*r^4) + (fp*hp*si^3)/(32*f*h^2*Pi^2*r^4) - 
   (23*hp^2*si^3)/(8*h^5*Pi^2*r^4) + (19*hp^2*si^3)/(8*h^4*Pi^2*r^4) + 
   (hpp*si^3)/(h^4*Pi^2*r^4) - (hpp*si^3)/(h^3*Pi^2*r^4) + 
   (61*fp^3*si^3)/(96*f^3*h^3*Pi^2*r^3) - 
   (11*fp*fpp*si^3)/(16*f^2*h^3*Pi^2*r^3) - (f3p*si^3)/(8*f*h^3*Pi^2*r^3) - 
   (3*fp^3*si^3)/(32*f^3*h^2*Pi^2*r^3) + 
   (fp*fpp*si^3)/(16*f^2*h^2*Pi^2*r^3) + (f3p*si^3)/(8*f*h^2*Pi^2*r^3) + 
   (3*fp^2*hp*si^3)/(4*f^2*h^4*Pi^2*r^3) + 
   (5*fpp*hp*si^3)/(8*f*h^4*Pi^2*r^3) - (fp^2*hp*si^3)/(8*f^2*h^3*Pi^2*r^3) - 
   (fpp*hp*si^3)/(8*f*h^3*Pi^2*r^3) - (37*fp*hp^2*si^3)/(32*f*h^5*Pi^2*r^3) - 
   (13*fp*hp^2*si^3)/(32*f*h^4*Pi^2*r^3) - (83*hp^3*si^3)/(48*h^6*Pi^2*r^3) + 
   (7*hp^3*si^3)/(4*h^5*Pi^2*r^3) + (5*fp*hpp*si^3)/(16*f*h^4*Pi^2*r^3) + 
   (3*fp*hpp*si^3)/(16*f*h^3*Pi^2*r^3) + (13*hp*hpp*si^3)/(8*h^5*Pi^2*r^3) - 
   (13*hp*hpp*si^3)/(8*h^4*Pi^2*r^3) - (h3p*si^3)/(4*h^4*Pi^2*r^3) + 
   (h3p*si^3)/(4*h^3*Pi^2*r^3) + (5*fp^4*si^3)/(32*f^4*h^3*Pi^2*r^2) - 
   (fpp^2*si^3)/(8*f^2*h^3*Pi^2*r^2) - 
   (3*fp*f3p*si^3)/(16*f^2*h^3*Pi^2*r^2) + (f4p*si^3)/(8*f*h^3*Pi^2*r^2) + 
   (5*fp^4*si^3)/(16*f^4*h^2*Pi^2*r^2) - 
   (3*fp^2*fpp*si^3)/(4*f^3*h^2*Pi^2*r^2) + 
   (fpp^2*si^3)/(4*f^2*h^2*Pi^2*r^2) + 
   (5*fp*f3p*si^3)/(16*f^2*h^2*Pi^2*r^2) - (f4p*si^3)/(8*f*h^2*Pi^2*r^2) - 
   (fp^3*hp*si^3)/(8*f^3*h^4*Pi^2*r^2) + 
   (29*fp*fpp*hp*si^3)/(32*f^2*h^4*Pi^2*r^2) - 
   (3*f3p*hp*si^3)/(4*f*h^4*Pi^2*r^2) + 
   (3*fp^3*hp*si^3)/(8*f^3*h^3*Pi^2*r^2) - 
   (23*fp*fpp*hp*si^3)/(32*f^2*h^3*Pi^2*r^2) + 
   (3*f3p*hp*si^3)/(8*f*h^3*Pi^2*r^2) - 
   (35*fp^2*hp^2*si^3)/(32*f^2*h^5*Pi^2*r^2) + 
   (77*fpp*hp^2*si^3)/(32*f*h^5*Pi^2*r^2) + 
   (3*fp^2*hp^2*si^3)/(8*f^2*h^4*Pi^2*r^2) - 
   (19*fpp*hp^2*si^3)/(32*f*h^4*Pi^2*r^2) - 
   (67*fp*hp^3*si^3)/(16*f*h^6*Pi^2*r^2) + 
   (7*fp*hp^3*si^3)/(16*f*h^5*Pi^2*r^2) + (11*hp^4*si^3)/(4*h^7*Pi^2*r^2) + 
   (11*fp^2*hpp*si^3)/(32*f^2*h^4*Pi^2*r^2) - 
   (3*fpp*hpp*si^3)/(4*f*h^4*Pi^2*r^2) - 
   (5*fp^2*hpp*si^3)/(32*f^2*h^3*Pi^2*r^2) + 
   (fpp*hpp*si^3)/(4*f*h^3*Pi^2*r^2) + 
   (95*fp*hp*hpp*si^3)/(32*f*h^5*Pi^2*r^2) - 
   (13*fp*hp*hpp*si^3)/(32*f*h^4*Pi^2*r^2) - 
   (21*hp^2*hpp*si^3)/(8*h^6*Pi^2*r^2) + (hpp^2*si^3)/(4*h^5*Pi^2*r^2) - 
   (5*fp*h3p*si^3)/(16*f*h^4*Pi^2*r^2) + (fp*h3p*si^3)/(16*f*h^3*Pi^2*r^2) + 
   (hp*h3p*si^3)/(4*h^5*Pi^2*r^2) - (175*fp^5*si^3)/(256*f^5*h^3*Pi^2*r) + 
   (125*fp^3*fpp*si^3)/(64*f^4*h^3*Pi^2*r) - 
   (75*fp*fpp^2*si^3)/(64*f^3*h^3*Pi^2*r) - 
   (21*fp^2*f3p*si^3)/(32*f^3*h^3*Pi^2*r) + 
   (7*fpp*f3p*si^3)/(16*f^2*h^3*Pi^2*r) + (fp*f4p*si^3)/(8*f^2*h^3*Pi^2*r) - 
   (215*fp^4*hp*si^3)/(256*f^4*h^4*Pi^2*r) + 
   (117*fp^2*fpp*hp*si^3)/(64*f^3*h^4*Pi^2*r) - 
   (35*fpp^2*hp*si^3)/(64*f^2*h^4*Pi^2*r) - 
   (13*fp*f3p*hp*si^3)/(32*f^2*h^4*Pi^2*r) - (f4p*hp*si^3)/(8*f*h^4*Pi^2*r) - 
   (105*fp^3*hp^2*si^3)/(256*f^3*h^5*Pi^2*r) + 
   (7*fp*fpp*hp^2*si^3)/(32*f^2*h^5*Pi^2*r) + 
   (7*f3p*hp^2*si^3)/(8*f*h^5*Pi^2*r) + 
   (95*fp^2*hp^3*si^3)/(256*f^2*h^6*Pi^2*r) - 
   (71*fpp*hp^3*si^3)/(32*f*h^6*Pi^2*r) + 
   (11*fp*hp^4*si^3)/(8*f*h^7*Pi^2*r) + 
   (9*fp^3*hpp*si^3)/(64*f^3*h^4*Pi^2*r) - 
   (3*fp*fpp*hpp*si^3)/(32*f^2*h^4*Pi^2*r) - 
   (f3p*hpp*si^3)/(4*f*h^4*Pi^2*r) - 
   (11*fp^2*hp*hpp*si^3)/(64*f^2*h^5*Pi^2*r) + 
   (23*fpp*hp*hpp*si^3)/(16*f*h^5*Pi^2*r) - 
   (21*fp*hp^2*hpp*si^3)/(16*f*h^6*Pi^2*r) + 
   (fp*hpp^2*si^3)/(8*f*h^5*Pi^2*r) - (fpp*h3p*si^3)/(8*f*h^4*Pi^2*r) + 
   (fp*hp*h3p*si^3)/(8*f*h^5*Pi^2*r))

Trrdsmm2 = (1/m^2)((71*fp^6)/(32256*f^6*h^3*Pi^2) - 
   (23*fp^4*fpp)/(3840*f^5*h^3*Pi^2) + (13*fp^2*fpp^2)/(4480*f^4*h^3*Pi^2) + 
   (11*fpp^3)/(20160*f^3*h^3*Pi^2) + (137*fp^3*f3p)/(53760*f^4*h^3*Pi^2) - 
   (11*fp*fpp*f3p)/(6720*f^3*h^3*Pi^2) + f3p^2/(8960*f^2*h^3*Pi^2) - 
   (3*fp^2*f4p)/(4480*f^3*h^3*Pi^2) - (fpp*f4p)/(4480*f^2*h^3*Pi^2) + 
   (fp*f5p)/(4480*f^2*h^3*Pi^2) + (23*fp^5*hp)/(7680*f^5*h^4*Pi^2) - 
   (241*fp^3*fpp*hp)/(35840*f^4*h^4*Pi^2) + 
   (11*fp*fpp^2*hp)/(6720*f^3*h^4*Pi^2) + 
   (19*fp^2*f3p*hp)/(6720*f^3*h^4*Pi^2) + 
   (3*fpp*f3p*hp)/(8960*f^2*h^4*Pi^2) - (9*fp*f4p*hp)/(8960*f^2*h^4*Pi^2) + 
   (11*fp^4*hp^2)/(3360*f^4*h^5*Pi^2) - 
   (101*fp^2*fpp*hp^2)/(17920*f^3*h^5*Pi^2) - 
   (29*fpp^2*hp^2)/(35840*f^2*h^5*Pi^2) + 
   (53*fp*f3p*hp^2)/(17920*f^2*h^5*Pi^2) + 
   (499*fp^3*hp^3)/(161280*f^3*h^6*Pi^2) - 
   (29*fp*fpp*hp^3)/(7168*f^2*h^6*Pi^2) + 
   (13*fp^2*hp^4)/(3584*f^2*h^7*Pi^2) - 
   (137*fp^4*hpp)/(107520*f^4*h^4*Pi^2) + 
   (29*fp^2*fpp*hpp)/(13440*f^3*h^4*Pi^2) + (fpp^2*hpp)/(2240*f^2*h^4*Pi^2) - 
   (11*fp*f3p*hpp)/(8960*f^2*h^4*Pi^2) - 
   (139*fp^3*hp*hpp)/(53760*f^3*h^5*Pi^2) + 
   (61*fp*fpp*hp*hpp)/(17920*f^2*h^5*Pi^2) - 
   (179*fp^2*hp^2*hpp)/(35840*f^2*h^6*Pi^2) + 
   (27*fp^2*hpp^2)/(35840*f^2*h^5*Pi^2) + (3*fp^3*h3p)/(8960*f^3*h^4*Pi^2) - 
   (fp*fpp*h3p)/(2240*f^2*h^4*Pi^2) + (9*fp^2*hp*h3p)/(8960*f^2*h^5*Pi^2) - 
   (fp^2*h4p)/(8960*f^2*h^4*Pi^2) + 1/(2520*Pi^2*r^6) - 
   67/(10080*h^3*Pi^2*r^6) + 1/(160*h*Pi^2*r^6) + 
   (7*fp)/(720*f*h^3*Pi^2*r^5) - fp/(96*f*h^2*Pi^2*r^5) + 
   fp/(1440*f*h*Pi^2*r^5) - (27*hp)/(1120*h^4*Pi^2*r^5) + 
   hp/(96*h^3*Pi^2*r^5) + (247*fp^2)/(23040*f^2*h^3*Pi^2*r^4) - 
   fpp/(288*f*h^3*Pi^2*r^4) - (5*fp^2)/(768*f^2*h^2*Pi^2*r^4) + 
   fpp/(288*f*h^2*Pi^2*r^4) - (17*fp*hp)/(3840*f*h^4*Pi^2*r^4) - 
   (fp*hp)/(1152*f*h^3*Pi^2*r^4) - (1381*hp^2)/(53760*h^5*Pi^2*r^4) + 
   (7*hp^2)/(768*h^4*Pi^2*r^4) + (27*hpp)/(2240*h^4*Pi^2*r^4) - 
   hpp/(192*h^3*Pi^2*r^4) - fp^3/(10752*f^3*h^3*Pi^2*r^3) + 
   (557*fp*fpp)/(80640*f^2*h^3*Pi^2*r^3) - (101*f3p)/(20160*f*h^3*Pi^2*r^3) + 
   fp^3/(768*f^3*h^2*Pi^2*r^3) - (fp*fpp)/(288*f^2*h^2*Pi^2*r^3) + 
   f3p/(576*f*h^2*Pi^2*r^3) - (281*fp^2*hp)/(53760*f^2*h^4*Pi^2*r^3) + 
   (281*fpp*hp)/(26880*f*h^4*Pi^2*r^3) + (fp^2*hp)/(1152*f^2*h^3*Pi^2*r^3) - 
   (fpp*hp)/(576*f*h^3*Pi^2*r^3) - (1303*fp*hp^2)/(53760*f*h^5*Pi^2*r^3) + 
   (7*fp*hp^2)/(2304*f*h^4*Pi^2*r^3) + (121*hp^3)/(23040*h^6*Pi^2*r^3) + 
   (281*fp*hpp)/(26880*f*h^4*Pi^2*r^3) - (fp*hpp)/(576*f*h^3*Pi^2*r^3) - 
   (11*hp*hpp)/(3840*h^5*Pi^2*r^3) - (517*fp^4)/(35840*f^4*h^3*Pi^2*r^2) + 
   (1039*fp^2*fpp)/(32256*f^3*h^3*Pi^2*r^2) - 
   (241*fpp^2)/(32256*f^2*h^3*Pi^2*r^2) - 
   (1067*fp*f3p)/(80640*f^2*h^3*Pi^2*r^2) + (11*f4p)/(3360*f*h^3*Pi^2*r^2) + 
   (7*fp^4)/(23040*f^4*h^2*Pi^2*r^2) - (fp^2*fpp)/(1920*f^3*h^2*Pi^2*r^2) - 
   fpp^2/(5760*f^2*h^2*Pi^2*r^2) + (fp*f3p)/(2880*f^2*h^2*Pi^2*r^2) - 
   (639*fp^3*hp)/(35840*f^3*h^4*Pi^2*r^2) + 
   (109*fp*fpp*hp)/(3584*f^2*h^4*Pi^2*r^2) - 
   (851*f3p*hp)/(80640*f*h^4*Pi^2*r^2) + (fp^3*hp)/(3840*f^3*h^3*Pi^2*r^2) - 
   (fp*fpp*hp)/(2880*f^2*h^3*Pi^2*r^2) - 
   (2669*fp^2*hp^2)/(107520*f^2*h^5*Pi^2*r^2) + 
   (589*fpp*hp^2)/(26880*f*h^5*Pi^2*r^2) + 
   (7*fp^2*hp^2)/(23040*f^2*h^4*Pi^2*r^2) - 
   (8209*fp*hp^3)/(322560*f*h^6*Pi^2*r^2) + (299*hp^4)/(5376*h^7*Pi^2*r^2) + 
   (871*fp^2*hpp)/(80640*f^2*h^4*Pi^2*r^2) - 
   (151*fpp*hpp)/(16128*f*h^4*Pi^2*r^2) - 
   (fp^2*hpp)/(5760*f^2*h^3*Pi^2*r^2) + 
   (1229*fp*hp*hpp)/(53760*f*h^5*Pi^2*r^2) - 
   (4117*hp^2*hpp)/(53760*h^6*Pi^2*r^2) + (207*hpp^2)/(17920*h^5*Pi^2*r^2) - 
   (3*fp*h3p)/(896*f*h^4*Pi^2*r^2) + (69*hp*h3p)/(4480*h^5*Pi^2*r^2) - 
   (23*h4p)/(13440*h^4*Pi^2*r^2) + (263*fp^5)/(53760*f^5*h^3*Pi^2*r) - 
   (1781*fp^3*fpp)/(107520*f^4*h^3*Pi^2*r) + 
   (29*fp*fpp^2)/(2240*f^3*h^3*Pi^2*r) + 
   (331*fp^2*f3p)/(53760*f^3*h^3*Pi^2*r) - 
   (169*fpp*f3p)/(26880*f^2*h^3*Pi^2*r) - 
   (53*fp*f4p)/(26880*f^2*h^3*Pi^2*r) + (11*f5p)/(13440*f*h^3*Pi^2*r) + 
   (111*fp^4*hp)/(17920*f^4*h^4*Pi^2*r) - 
   (1933*fp^2*fpp*hp)/(107520*f^3*h^4*Pi^2*r) + 
   (157*fpp^2*hp)/(17920*f^2*h^4*Pi^2*r) + 
   (13*fp*f3p*hp)/(1792*f^2*h^4*Pi^2*r) - (33*f4p*hp)/(8960*f*h^4*Pi^2*r) + 
   (79*fp^3*hp^2)/(10752*f^3*h^5*Pi^2*r) - 
   (657*fp*fpp*hp^2)/(35840*f^2*h^5*Pi^2*r) + 
   (583*f3p*hp^2)/(53760*f*h^5*Pi^2*r) + 
   (239*fp^2*hp^3)/(53760*f^2*h^6*Pi^2*r) - 
   (319*fpp*hp^3)/(21504*f*h^6*Pi^2*r) + (143*fp*hp^4)/(5376*f*h^7*Pi^2*r) - 
   (167*fp^3*hpp)/(53760*f^3*h^4*Pi^2*r) + 
   (139*fp*fpp*hpp)/(17920*f^2*h^4*Pi^2*r) - 
   (121*f3p*hpp)/(26880*f*h^4*Pi^2*r) - 
   (23*fp^2*hp*hpp)/(6720*f^2*h^5*Pi^2*r) + 
   (671*fpp*hp*hpp)/(53760*f*h^5*Pi^2*r) - 
   (1969*fp*hp^2*hpp)/(53760*f*h^6*Pi^2*r) + 
   (99*fp*hpp^2)/(17920*f*h^5*Pi^2*r) + (fp^2*h3p)/(2688*f^2*h^4*Pi^2*r) - 
   (11*fpp*h3p)/(6720*f*h^4*Pi^2*r) + (33*fp*hp*h3p)/(4480*f*h^5*Pi^2*r) - 
   (11*fp*h4p)/(13440*f*h^4*Pi^2*r) - (53*fp^6*si)/(2560*f^6*h^3*Pi^2) + 
   (217*fp^4*fpp*si)/(3840*f^5*h^3*Pi^2) - 
   (53*fp^2*fpp^2*si)/(1920*f^4*h^3*Pi^2) - (fpp^3*si)/(192*f^3*h^3*Pi^2) - 
   (23*fp^3*f3p*si)/(960*f^4*h^3*Pi^2) + (fp*fpp*f3p*si)/(64*f^3*h^3*Pi^2) - 
   (f3p^2*si)/(960*f^2*h^3*Pi^2) + (fp^2*f4p*si)/(160*f^3*h^3*Pi^2) + 
   (fpp*f4p*si)/(480*f^2*h^3*Pi^2) - (fp*f5p*si)/(480*f^2*h^3*Pi^2) - 
   (217*fp^5*hp*si)/(7680*f^5*h^4*Pi^2) + 
   (61*fp^3*fpp*hp*si)/(960*f^4*h^4*Pi^2) - 
   (fp*fpp^2*hp*si)/(64*f^3*h^4*Pi^2) - 
   (17*fp^2*f3p*hp*si)/(640*f^3*h^4*Pi^2) - 
   (fpp*f3p*hp*si)/(320*f^2*h^4*Pi^2) + (3*fp*f4p*hp*si)/(320*f^2*h^4*Pi^2) - 
   (79*fp^4*hp^2*si)/(2560*f^4*h^5*Pi^2) + 
   (17*fp^2*fpp*hp^2*si)/(320*f^3*h^5*Pi^2) + 
   (29*fpp^2*hp^2*si)/(3840*f^2*h^5*Pi^2) - 
   (53*fp*f3p*hp^2*si)/(1920*f^2*h^5*Pi^2) - 
   (223*fp^3*hp^3*si)/(7680*f^3*h^6*Pi^2) + 
   (29*fp*fpp*hp^3*si)/(768*f^2*h^6*Pi^2) - 
   (13*fp^2*hp^4*si)/(384*f^2*h^7*Pi^2) + 
   (23*fp^4*hpp*si)/(1920*f^4*h^4*Pi^2) - 
   (13*fp^2*fpp*hpp*si)/(640*f^3*h^4*Pi^2) - 
   (fpp^2*hpp*si)/(240*f^2*h^4*Pi^2) + 
   (11*fp*f3p*hpp*si)/(960*f^2*h^4*Pi^2) + 
   (31*fp^3*hp*hpp*si)/(1280*f^3*h^5*Pi^2) - 
   (61*fp*fpp*hp*hpp*si)/(1920*f^2*h^5*Pi^2) + 
   (179*fp^2*hp^2*hpp*si)/(3840*f^2*h^6*Pi^2) - 
   (9*fp^2*hpp^2*si)/(1280*f^2*h^5*Pi^2) - (fp^3*h3p*si)/(320*f^3*h^4*Pi^2) + 
   (fp*fpp*h3p*si)/(240*f^2*h^4*Pi^2) - 
   (3*fp^2*hp*h3p*si)/(320*f^2*h^5*Pi^2) + (fp^2*h4p*si)/(960*f^2*h^4*Pi^2) - 
   si/(240*Pi^2*r^6) + si/(30*h^3*Pi^2*r^6) + si/(16*h^2*Pi^2*r^6) - 
   (11*si)/(120*h*Pi^2*r^6) - (13*fp*si)/(80*f*h^3*Pi^2*r^5) + 
   (25*fp*si)/(144*f*h^2*Pi^2*r^5) - (fp*si)/(90*f*h*Pi^2*r^5) + 
   (67*hp*si)/(240*h^4*Pi^2*r^5) - (7*hp*si)/(48*h^3*Pi^2*r^5) - 
   (857*fp^2*si)/(5760*f^2*h^3*Pi^2*r^4) + (11*fpp*si)/(288*f*h^3*Pi^2*r^4) + 
   (101*fp^2*si)/(1152*f^2*h^2*Pi^2*r^4) - (11*fpp*si)/(288*f*h^2*Pi^2*r^4) + 
   (21*fp*hp*si)/(320*f*h^4*Pi^2*r^4) + (5*fp*hp*si)/(576*f*h^3*Pi^2*r^4) + 
   (187*hp^2*si)/(640*h^5*Pi^2*r^4) - (49*hp^2*si)/(384*h^4*Pi^2*r^4) - 
   (67*hpp*si)/(480*h^4*Pi^2*r^4) + (7*hpp*si)/(96*h^3*Pi^2*r^4) - 
   (17*fp^3*si)/(2880*f^3*h^3*Pi^2*r^3) - 
   (121*fp*fpp*si)/(1440*f^2*h^3*Pi^2*r^3) + 
   (11*f3p*si)/(180*f*h^3*Pi^2*r^3) - (fp^3*si)/(48*f^3*h^2*Pi^2*r^3) + 
   (fp*fpp*si)/(18*f^2*h^2*Pi^2*r^3) - (f3p*si)/(36*f*h^2*Pi^2*r^3) + 
   (13*fp^2*hp*si)/(192*f^2*h^4*Pi^2*r^3) - 
   (13*fpp*hp*si)/(96*f*h^4*Pi^2*r^3) - (fp^2*hp*si)/(72*f^2*h^3*Pi^2*r^3) + 
   (fpp*hp*si)/(36*f*h^3*Pi^2*r^3) + (101*fp*hp^2*si)/(320*f*h^5*Pi^2*r^3) - 
   (7*fp*hp^2*si)/(144*f*h^4*Pi^2*r^3) - (77*hp^3*si)/(960*h^6*Pi^2*r^3) - 
   (13*fp*hpp*si)/(96*f*h^4*Pi^2*r^3) + (fp*hpp*si)/(36*f*h^3*Pi^2*r^3) + 
   (7*hp*hpp*si)/(160*h^5*Pi^2*r^3) + (19*fp^4*si)/(120*f^4*h^3*Pi^2*r^2) - 
   (2083*fp^2*fpp*si)/(5760*f^3*h^3*Pi^2*r^2) + 
   (25*fpp^2*si)/(288*f^2*h^3*Pi^2*r^2) + 
   (419*fp*f3p*si)/(2880*f^2*h^3*Pi^2*r^2) - (f4p*si)/(30*f*h^3*Pi^2*r^2) - 
   (7*fp^4*si)/(1440*f^4*h^2*Pi^2*r^2) + 
   (fp^2*fpp*si)/(120*f^3*h^2*Pi^2*r^2) + (fpp^2*si)/(360*f^2*h^2*Pi^2*r^2) - 
   (fp*f3p*si)/(180*f^2*h^2*Pi^2*r^2) + 
   (727*fp^3*hp*si)/(3840*f^3*h^4*Pi^2*r^2) - 
   (611*fp*fpp*hp*si)/(1920*f^2*h^4*Pi^2*r^2) + 
   (31*f3p*hp*si)/(320*f*h^4*Pi^2*r^2) - 
   (fp^3*hp*si)/(240*f^3*h^3*Pi^2*r^2) + 
   (fp*fpp*hp*si)/(180*f^2*h^3*Pi^2*r^2) + 
   (499*fp^2*hp^2*si)/(1920*f^2*h^5*Pi^2*r^2) - 
   (61*fpp*hp^2*si)/(320*f*h^5*Pi^2*r^2) - 
   (7*fp^2*hp^2*si)/(1440*f^2*h^4*Pi^2*r^2) + 
   (863*fp*hp^3*si)/(3840*f*h^6*Pi^2*r^2) - (13*hp^4*si)/(24*h^7*Pi^2*r^2) - 
   (73*fp^2*hpp*si)/(640*f^2*h^4*Pi^2*r^2) + 
   (27*fpp*hpp*si)/(320*f*h^4*Pi^2*r^2) + 
   (fp^2*hpp*si)/(360*f^2*h^3*Pi^2*r^2) - 
   (407*fp*hp*hpp*si)/(1920*f*h^5*Pi^2*r^2) + 
   (179*hp^2*hpp*si)/(240*h^6*Pi^2*r^2) - (9*hpp^2*si)/(80*h^5*Pi^2*r^2) + 
   (fp*h3p*si)/(30*f*h^4*Pi^2*r^2) - (3*hp*h3p*si)/(20*h^5*Pi^2*r^2) + 
   (h4p*si)/(60*h^4*Pi^2*r^2) - (49*fp^5*si)/(960*f^5*h^3*Pi^2*r) + 
   (113*fp^3*fpp*si)/(640*f^4*h^3*Pi^2*r) - 
   (137*fp*fpp^2*si)/(960*f^3*h^3*Pi^2*r) - 
   (31*fp^2*f3p*si)/(480*f^3*h^3*Pi^2*r) + 
   (11*fpp*f3p*si)/(160*f^2*h^3*Pi^2*r) + (fp*f4p*si)/(48*f^2*h^3*Pi^2*r) - 
   (f5p*si)/(120*f*h^3*Pi^2*r) - (263*fp^4*hp*si)/(3840*f^4*h^4*Pi^2*r) + 
   (127*fp^2*fpp*hp*si)/(640*f^3*h^4*Pi^2*r) - 
   (3*fpp^2*hp*si)/(32*f^2*h^4*Pi^2*r) - 
   (13*fp*f3p*hp*si)/(160*f^2*h^4*Pi^2*r) + (3*f4p*hp*si)/(80*f*h^4*Pi^2*r) - 
   (53*fp^3*hp^2*si)/(640*f^3*h^5*Pi^2*r) + 
   (197*fp*fpp*hp^2*si)/(960*f^2*h^5*Pi^2*r) - 
   (53*f3p*hp^2*si)/(480*f*h^5*Pi^2*r) - 
   (211*fp^2*hp^3*si)/(3840*f^2*h^6*Pi^2*r) + 
   (29*fpp*hp^3*si)/(192*f*h^6*Pi^2*r) - (13*fp*hp^4*si)/(48*f*h^7*Pi^2*r) + 
   (13*fp^3*hpp*si)/(384*f^3*h^4*Pi^2*r) - 
   (41*fp*fpp*hpp*si)/(480*f^2*h^4*Pi^2*r) + 
   (11*f3p*hpp*si)/(240*f*h^4*Pi^2*r) + 
   (79*fp^2*hp*hpp*si)/(1920*f^2*h^5*Pi^2*r) - 
   (61*fpp*hp*hpp*si)/(480*f*h^5*Pi^2*r) + 
   (179*fp*hp^2*hpp*si)/(480*f*h^6*Pi^2*r) - 
   (9*fp*hpp^2*si)/(160*f*h^5*Pi^2*r) - (fp^2*h3p*si)/(240*f^2*h^4*Pi^2*r) + 
   (fpp*h3p*si)/(60*f*h^4*Pi^2*r) - (3*fp*hp*h3p*si)/(40*f*h^5*Pi^2*r) + 
   (fp*h4p*si)/(120*f*h^4*Pi^2*r) + (85*fp^6*si^2)/(1536*f^6*h^3*Pi^2) - 
   (119*fp^4*fpp*si^2)/(768*f^5*h^3*Pi^2) + 
   (31*fp^2*fpp^2*si^2)/(384*f^4*h^3*Pi^2) + (fpp^3*si^2)/(64*f^3*h^3*Pi^2) + 
   (49*fp^3*f3p*si^2)/(768*f^4*h^3*Pi^2) - 
   (3*fp*fpp*f3p*si^2)/(64*f^3*h^3*Pi^2) + (f3p^2*si^2)/(384*f^2*h^3*Pi^2) - 
   (fp^2*f4p*si^2)/(64*f^3*h^3*Pi^2) - (fpp*f4p*si^2)/(192*f^2*h^3*Pi^2) + 
   (fp*f5p*si^2)/(192*f^2*h^3*Pi^2) + 
   (119*fp^5*hp*si^2)/(1536*f^5*h^4*Pi^2) - 
   (271*fp^3*fpp*hp*si^2)/(1536*f^4*h^4*Pi^2) + 
   (3*fp*fpp^2*hp*si^2)/(64*f^3*h^4*Pi^2) + 
   (9*fp^2*f3p*hp*si^2)/(128*f^3*h^4*Pi^2) + 
   (fpp*f3p*hp*si^2)/(128*f^2*h^4*Pi^2) - 
   (3*fp*f4p*hp*si^2)/(128*f^2*h^4*Pi^2) + 
   (43*fp^4*hp^2*si^2)/(512*f^4*h^5*Pi^2) - 
   (37*fp^2*fpp*hp^2*si^2)/(256*f^3*h^5*Pi^2) - 
   (29*fpp^2*hp^2*si^2)/(1536*f^2*h^5*Pi^2) + 
   (53*fp*f3p*hp^2*si^2)/(768*f^2*h^5*Pi^2) + 
   (39*fp^3*hp^3*si^2)/(512*f^3*h^6*Pi^2) - 
   (145*fp*fpp*hp^3*si^2)/(1536*f^2*h^6*Pi^2) + 
   (65*fp^2*hp^4*si^2)/(768*f^2*h^7*Pi^2) - 
   (49*fp^4*hpp*si^2)/(1536*f^4*h^4*Pi^2) + 
   (7*fp^2*fpp*hpp*si^2)/(128*f^3*h^4*Pi^2) + 
   (fpp^2*hpp*si^2)/(96*f^2*h^4*Pi^2) - 
   (11*fp*f3p*hpp*si^2)/(384*f^2*h^4*Pi^2) - 
   (fp^3*hp*hpp*si^2)/(16*f^3*h^5*Pi^2) + 
   (61*fp*fpp*hp*hpp*si^2)/(768*f^2*h^5*Pi^2) - 
   (179*fp^2*hp^2*hpp*si^2)/(1536*f^2*h^6*Pi^2) + 
   (9*fp^2*hpp^2*si^2)/(512*f^2*h^5*Pi^2) + 
   (fp^3*h3p*si^2)/(128*f^3*h^4*Pi^2) - (fp*fpp*h3p*si^2)/(96*f^2*h^4*Pi^2) + 
   (3*fp^2*hp*h3p*si^2)/(128*f^2*h^5*Pi^2) - 
   (fp^2*h4p*si^2)/(384*f^2*h^4*Pi^2) + si^2/(48*Pi^2*r^6) + 
   (3*si^2)/(16*h^3*Pi^2*r^6) - (11*si^2)/(16*h^2*Pi^2*r^6) + 
   (23*si^2)/(48*h*Pi^2*r^6) + (43*fp*si^2)/(48*f*h^3*Pi^2*r^5) - 
   (23*fp*si^2)/(24*f*h^2*Pi^2*r^5) + (fp*si^2)/(16*f*h*Pi^2*r^5) - 
   (hp*si^2)/(h^4*Pi^2*r^5) + (2*hp*si^2)/(3*h^3*Pi^2*r^5) + 
   (65*fp^2*si^2)/(96*f^2*h^3*Pi^2*r^4) - (5*fpp*si^2)/(48*f*h^3*Pi^2*r^4) - 
   (3*fp^2*si^2)/(8*f^2*h^2*Pi^2*r^4) + (5*fpp*si^2)/(48*f*h^2*Pi^2*r^4) - 
   (5*fp*hp*si^2)/(16*f*h^4*Pi^2*r^4) - (fp*hp*si^2)/(48*f*h^3*Pi^2*r^4) - 
   (hp^2*si^2)/(h^5*Pi^2*r^4) + (7*hp^2*si^2)/(12*h^4*Pi^2*r^4) + 
   (hpp*si^2)/(2*h^4*Pi^2*r^4) - (hpp*si^2)/(3*h^3*Pi^2*r^4) + 
   (17*fp^3*si^2)/(192*f^3*h^3*Pi^2*r^3) + 
   (7*fp*fpp*si^2)/(24*f^2*h^3*Pi^2*r^3) - 
   (11*f3p*si^2)/(48*f*h^3*Pi^2*r^3) + (7*fp^3*si^2)/(64*f^3*h^2*Pi^2*r^3) - 
   (7*fp*fpp*si^2)/(24*f^2*h^2*Pi^2*r^3) + (7*f3p*si^2)/(48*f*h^2*Pi^2*r^3) - 
   (9*fp^2*hp*si^2)/(32*f^2*h^4*Pi^2*r^3) + 
   (9*fpp*hp*si^2)/(16*f*h^4*Pi^2*r^3) + 
   (7*fp^2*hp*si^2)/(96*f^2*h^3*Pi^2*r^3) - 
   (7*fpp*hp*si^2)/(48*f*h^3*Pi^2*r^3) - 
   (85*fp*hp^2*si^2)/(64*f*h^5*Pi^2*r^3) + 
   (49*fp*hp^2*si^2)/(192*f*h^4*Pi^2*r^3) + 
   (11*hp^3*si^2)/(24*h^6*Pi^2*r^3) + (9*fp*hpp*si^2)/(16*f*h^4*Pi^2*r^3) - 
   (7*fp*hpp*si^2)/(48*f*h^3*Pi^2*r^3) - (hp*hpp*si^2)/(4*h^5*Pi^2*r^3) - 
   (383*fp^4*si^2)/(768*f^4*h^3*Pi^2*r^2) + 
   (227*fp^2*fpp*si^2)/(192*f^3*h^3*Pi^2*r^2) - 
   (55*fpp^2*si^2)/(192*f^2*h^3*Pi^2*r^2) - 
   (11*fp*f3p*si^2)/(24*f^2*h^3*Pi^2*r^2) + (f4p*si^2)/(12*f*h^3*Pi^2*r^2) + 
   (7*fp^4*si^2)/(256*f^4*h^2*Pi^2*r^2) - 
   (3*fp^2*fpp*si^2)/(64*f^3*h^2*Pi^2*r^2) - 
   (fpp^2*si^2)/(64*f^2*h^2*Pi^2*r^2) + (fp*f3p*si^2)/(32*f^2*h^2*Pi^2*r^2) - 
   (71*fp^3*hp*si^2)/(128*f^3*h^4*Pi^2*r^2) + 
   (57*fp*fpp*hp*si^2)/(64*f^2*h^4*Pi^2*r^2) - 
   (f3p*hp*si^2)/(6*f*h^4*Pi^2*r^2) + 
   (3*fp^3*hp*si^2)/(128*f^3*h^3*Pi^2*r^2) - 
   (fp*fpp*hp*si^2)/(32*f^2*h^3*Pi^2*r^2) - 
   (571*fp^2*hp^2*si^2)/(768*f^2*h^5*Pi^2*r^2) + 
   (fpp*hp^2*si^2)/(4*f*h^5*Pi^2*r^2) + 
   (7*fp^2*hp^2*si^2)/(256*f^2*h^4*Pi^2*r^2) - 
   (17*fp*hp^3*si^2)/(48*f*h^6*Pi^2*r^2) + (65*hp^4*si^2)/(48*h^7*Pi^2*r^2) + 
   (fp^2*hpp*si^2)/(3*f^2*h^4*Pi^2*r^2) - 
   (13*fpp*hpp*si^2)/(96*f*h^4*Pi^2*r^2) - 
   (fp^2*hpp*si^2)/(64*f^2*h^3*Pi^2*r^2) + 
   (5*fp*hp*hpp*si^2)/(12*f*h^5*Pi^2*r^2) - 
   (179*hp^2*hpp*si^2)/(96*h^6*Pi^2*r^2) + (9*hpp^2*si^2)/(32*h^5*Pi^2*r^2) - 
   (fp*h3p*si^2)/(12*f*h^4*Pi^2*r^2) + (3*hp*h3p*si^2)/(8*h^5*Pi^2*r^2) - 
   (h4p*si^2)/(24*h^4*Pi^2*r^2) + (31*fp^5*si^2)/(256*f^5*h^3*Pi^2*r) - 
   (59*fp^3*fpp*si^2)/(128*f^4*h^3*Pi^2*r) + 
   (27*fp*fpp^2*si^2)/(64*f^3*h^3*Pi^2*r) + 
   (31*fp^2*f3p*si^2)/(192*f^3*h^3*Pi^2*r) - 
   (5*fpp*f3p*si^2)/(24*f^2*h^3*Pi^2*r) - 
   (5*fp*f4p*si^2)/(96*f^2*h^3*Pi^2*r) + (f5p*si^2)/(48*f*h^3*Pi^2*r) + 
   (73*fp^4*hp*si^2)/(384*f^4*h^4*Pi^2*r) - 
   (73*fp^2*fpp*hp*si^2)/(128*f^3*h^4*Pi^2*r) + 
   (13*fpp^2*hp*si^2)/(48*f^2*h^4*Pi^2*r) + 
   (23*fp*f3p*hp*si^2)/(96*f^2*h^4*Pi^2*r) - 
   (3*f4p*hp*si^2)/(32*f*h^4*Pi^2*r) + 
   (63*fp^3*hp^2*si^2)/(256*f^3*h^5*Pi^2*r) - 
   (239*fp*fpp*hp^2*si^2)/(384*f^2*h^5*Pi^2*r) + 
   (53*f3p*hp^2*si^2)/(192*f*h^5*Pi^2*r) + 
   (3*fp^2*hp^3*si^2)/(16*f^2*h^6*Pi^2*r) - 
   (145*fpp*hp^3*si^2)/(384*f*h^6*Pi^2*r) + 
   (65*fp*hp^4*si^2)/(96*f*h^7*Pi^2*r) - 
   (3*fp^3*hpp*si^2)/(32*f^3*h^4*Pi^2*r) + 
   (fp*fpp*hpp*si^2)/(4*f^2*h^4*Pi^2*r) - 
   (11*f3p*hpp*si^2)/(96*f*h^4*Pi^2*r) - 
   (25*fp^2*hp*hpp*si^2)/(192*f^2*h^5*Pi^2*r) + 
   (61*fpp*hp*hpp*si^2)/(192*f*h^5*Pi^2*r) - 
   (179*fp*hp^2*hpp*si^2)/(192*f*h^6*Pi^2*r) + 
   (9*fp*hpp^2*si^2)/(64*f*h^5*Pi^2*r) + 
   (fp^2*h3p*si^2)/(96*f^2*h^4*Pi^2*r) - (fpp*h3p*si^2)/(24*f*h^4*Pi^2*r) + 
   (3*fp*hp*h3p*si^2)/(16*f*h^5*Pi^2*r) - (fp*h4p*si^2)/(48*f*h^4*Pi^2*r) - 
   (11*fp^6*si^3)/(768*f^6*h^3*Pi^2) + (7*fp^4*fpp*si^3)/(128*f^5*h^3*Pi^2) - 
   (3*fp^2*fpp^2*si^3)/(64*f^4*h^3*Pi^2) - (fpp^3*si^3)/(96*f^3*h^3*Pi^2) - 
   (fp^3*f3p*si^3)/(64*f^4*h^3*Pi^2) + (fp*fpp*f3p*si^3)/(32*f^3*h^3*Pi^2) - 
   (7*fp^5*hp*si^3)/(256*f^5*h^4*Pi^2) + 
   (9*fp^3*fpp*hp*si^3)/(128*f^4*h^4*Pi^2) - 
   (fp*fpp^2*hp*si^3)/(32*f^3*h^4*Pi^2) - 
   (fp^2*f3p*hp*si^3)/(64*f^3*h^4*Pi^2) - 
   (7*fp^4*hp^2*si^3)/(256*f^4*h^5*Pi^2) + 
   (3*fp^2*fpp*hp^2*si^3)/(64*f^3*h^5*Pi^2) - 
   (11*fp^3*hp^3*si^3)/(768*f^3*h^6*Pi^2) + 
   (fp^4*hpp*si^3)/(128*f^4*h^4*Pi^2) - 
   (fp^2*fpp*hpp*si^3)/(64*f^3*h^4*Pi^2) + 
   (fp^3*hp*hpp*si^3)/(128*f^3*h^5*Pi^2) - si^3/(24*Pi^2*r^6) - 
   (23*si^3)/(24*h^3*Pi^2*r^6) + (15*si^3)/(8*h^2*Pi^2*r^6) - 
   (7*si^3)/(8*h*Pi^2*r^6) - (13*fp*si^3)/(8*f*h^3*Pi^2*r^5) + 
   (7*fp*si^3)/(4*f*h^2*Pi^2*r^5) - (fp*si^3)/(8*f*h*Pi^2*r^5) + 
   (hp*si^3)/(h^4*Pi^2*r^5) - (hp*si^3)/(h^3*Pi^2*r^5) - 
   (fp^2*si^3)/(f^2*h^3*Pi^2*r^4) + (fp^2*si^3)/(2*f^2*h^2*Pi^2*r^4) + 
   (fp*hp*si^3)/(2*f*h^4*Pi^2*r^4) + (7*hp^2*si^3)/(8*h^5*Pi^2*r^4) - 
   (7*hp^2*si^3)/(8*h^4*Pi^2*r^4) - (hpp*si^3)/(2*h^4*Pi^2*r^4) + 
   (hpp*si^3)/(2*h^3*Pi^2*r^4) - (13*fp^3*si^3)/(48*f^3*h^3*Pi^2*r^3) - 
   (fp*fpp*si^3)/(4*f^2*h^3*Pi^2*r^3) + (f3p*si^3)/(4*f*h^3*Pi^2*r^3) - 
   (3*fp^3*si^3)/(16*f^3*h^2*Pi^2*r^3) + (fp*fpp*si^3)/(2*f^2*h^2*Pi^2*r^3) - 
   (f3p*si^3)/(4*f*h^2*Pi^2*r^3) + (3*fp^2*hp*si^3)/(8*f^2*h^4*Pi^2*r^3) - 
   (3*fpp*hp*si^3)/(4*f*h^4*Pi^2*r^3) - (fp^2*hp*si^3)/(8*f^2*h^3*Pi^2*r^3) + 
   (fpp*hp*si^3)/(4*f*h^3*Pi^2*r^3) + (29*fp*hp^2*si^3)/(16*f*h^5*Pi^2*r^3) - 
   (7*fp*hp^2*si^3)/(16*f*h^4*Pi^2*r^3) - (11*hp^3*si^3)/(12*h^6*Pi^2*r^3) - 
   (3*fp*hpp*si^3)/(4*f*h^4*Pi^2*r^3) + (fp*hpp*si^3)/(4*f*h^3*Pi^2*r^3) + 
   (hp*hpp*si^3)/(2*h^5*Pi^2*r^3) + (43*fp^4*si^3)/(128*f^4*h^3*Pi^2*r^2) - 
   (29*fp^2*fpp*si^3)/(32*f^3*h^3*Pi^2*r^2) + 
   (7*fpp^2*si^3)/(32*f^2*h^3*Pi^2*r^2) + 
   (5*fp*f3p*si^3)/(16*f^2*h^3*Pi^2*r^2) - 
   (7*fp^4*si^3)/(128*f^4*h^2*Pi^2*r^2) + 
   (3*fp^2*fpp*si^3)/(32*f^3*h^2*Pi^2*r^2) + 
   (fpp^2*si^3)/(32*f^2*h^2*Pi^2*r^2) - (fp*f3p*si^3)/(16*f^2*h^2*Pi^2*r^2) + 
   (17*fp^3*hp*si^3)/(64*f^3*h^4*Pi^2*r^2) - 
   (5*fp*fpp*hp*si^3)/(16*f^2*h^4*Pi^2*r^2) - 
   (f3p*hp*si^3)/(4*f*h^4*Pi^2*r^2) - 
   (3*fp^3*hp*si^3)/(64*f^3*h^3*Pi^2*r^2) + 
   (fp*fpp*hp*si^3)/(16*f^2*h^3*Pi^2*r^2) + 
   (39*fp^2*hp^2*si^3)/(128*f^2*h^5*Pi^2*r^2) + 
   (3*fpp*hp^2*si^3)/(4*f*h^5*Pi^2*r^2) - 
   (7*fp^2*hp^2*si^3)/(128*f^2*h^4*Pi^2*r^2) - 
   (11*fp*hp^3*si^3)/(16*f*h^6*Pi^2*r^2) - 
   (5*fp^2*hpp*si^3)/(32*f^2*h^4*Pi^2*r^2) - 
   (fpp*hpp*si^3)/(4*f*h^4*Pi^2*r^2) + 
   (fp^2*hpp*si^3)/(32*f^2*h^3*Pi^2*r^2) + 
   (3*fp*hp*hpp*si^3)/(8*f*h^5*Pi^2*r^2) + 
   (3*fp^5*si^3)/(128*f^5*h^3*Pi^2*r) + (fp^3*fpp*si^3)/(16*f^4*h^3*Pi^2*r) - 
   (7*fp*fpp^2*si^3)/(32*f^3*h^3*Pi^2*r) + 
   (fpp*f3p*si^3)/(8*f^2*h^3*Pi^2*r) - (fp^4*hp*si^3)/(16*f^4*h^4*Pi^2*r) + 
   (fp^2*fpp*hp*si^3)/(4*f^3*h^4*Pi^2*r) - 
   (fpp^2*hp*si^3)/(8*f^2*h^4*Pi^2*r) - (fp*f3p*hp*si^3)/(8*f^2*h^4*Pi^2*r) - 
   (17*fp^3*hp^2*si^3)/(128*f^3*h^5*Pi^2*r) + 
   (3*fp*fpp*hp^2*si^3)/(8*f^2*h^5*Pi^2*r) - 
   (11*fp^2*hp^3*si^3)/(64*f^2*h^6*Pi^2*r) + 
   (fp^3*hpp*si^3)/(32*f^3*h^4*Pi^2*r) - 
   (fp*fpp*hpp*si^3)/(8*f^2*h^4*Pi^2*r) + 
   (3*fp^2*hp*hpp*si^3)/(32*f^2*h^5*Pi^2*r))

Tttdsmm2 = (1/m^2)((-781*fp^6)/(32256*f^6*h^3*Pi^2) + 
   (2159*fp^4*fpp)/(26880*f^5*h^3*Pi^2) - 
   (131*fp^2*fpp^2)/(1920*f^4*h^3*Pi^2) + (179*fpp^3)/(20160*f^3*h^3*Pi^2) - 
   (229*fp^3*f3p)/(7680*f^4*h^3*Pi^2) + 
   (943*fp*fpp*f3p)/(26880*f^3*h^3*Pi^2) - (97*f3p^2)/(26880*f^2*h^3*Pi^2) + 
   (227*fp^2*f4p)/(26880*f^3*h^3*Pi^2) - (71*fpp*f4p)/(13440*f^2*h^3*Pi^2) - 
   (9*fp*f5p)/(4480*f^2*h^3*Pi^2) + f6p/(2240*f*h^3*Pi^2) - 
   (2159*fp^5*hp)/(53760*f^5*h^4*Pi^2) + 
   (347*fp^3*fpp*hp)/(3072*f^4*h^4*Pi^2) - 
   (709*fp*fpp^2*hp)/(10752*f^3*h^4*Pi^2) - 
   (461*fp^2*f3p*hp)/(10752*f^3*h^4*Pi^2) + 
   (239*fpp*f3p*hp)/(8960*f^2*h^4*Pi^2) + 
   (341*fp*f4p*hp)/(26880*f^2*h^4*Pi^2) - (3*f5p*hp)/(896*f*h^4*Pi^2) - 
   (3*fp^4*hp^2)/(64*f^4*h^5*Pi^2) + 
   (1163*fp^2*fpp*hp^2)/(10752*f^3*h^5*Pi^2) - 
   (3571*fpp^2*hp^2)/(107520*f^2*h^5*Pi^2) - 
   (2299*fp*f3p*hp^2)/(53760*f^2*h^5*Pi^2) + 
   (25*f4p*hp^2)/(1792*f*h^5*Pi^2) - (1555*fp^3*hp^3)/(32256*f^3*h^6*Pi^2) + 
   (9361*fp*fpp*hp^3)/(107520*f^2*h^6*Pi^2) - 
   (135*f3p*hp^3)/(3584*f*h^6*Pi^2) - (2581*fp^2*hp^4)/(53760*f^2*h^7*Pi^2) + 
   (113*fpp*hp^4)/(1792*f*h^7*Pi^2) - (13*fp*hp^5)/(256*f*h^8*Pi^2) + 
   (229*fp^4*hpp)/(15360*f^4*h^4*Pi^2) - 
   (617*fp^2*fpp*hpp)/(17920*f^3*h^4*Pi^2) + 
   (71*fpp^2*hpp)/(6720*f^2*h^4*Pi^2) + 
   (367*fp*f3p*hpp)/(26880*f^2*h^4*Pi^2) - (f4p*hpp)/(224*f*h^4*Pi^2) + 
   (649*fp^3*hp*hpp)/(17920*f^3*h^5*Pi^2) - 
   (3523*fp*fpp*hp*hpp)/(53760*f^2*h^5*Pi^2) + 
   (51*f3p*hp*hpp)/(1792*f*h^5*Pi^2) + 
   (2129*fp^2*hp^2*hpp)/(35840*f^2*h^6*Pi^2) - 
   (1403*fpp*hp^2*hpp)/(17920*f*h^6*Pi^2) + 
   (797*fp*hp^3*hpp)/(8960*f*h^7*Pi^2) - 
   (799*fp^2*hpp^2)/(107520*f^2*h^5*Pi^2) + 
   (11*fpp*hpp^2)/(1120*f*h^5*Pi^2) - (493*fp*hp*hpp^2)/(17920*f*h^6*Pi^2) - 
   (227*fp^3*h3p)/(53760*f^3*h^4*Pi^2) + 
   (103*fp*fpp*h3p)/(13440*f^2*h^4*Pi^2) - (3*f3p*h3p)/(896*f*h^4*Pi^2) - 
   (73*fp^2*hp*h3p)/(6720*f^2*h^5*Pi^2) + 
   (129*fpp*hp*h3p)/(8960*f*h^5*Pi^2) - 
   (359*fp*hp^2*h3p)/(17920*f*h^6*Pi^2) + (9*fp*hpp*h3p)/(1792*f*h^5*Pi^2) + 
   (9*fp^2*h4p)/(8960*f^2*h^4*Pi^2) - (3*fpp*h4p)/(2240*f*h^4*Pi^2) + 
   (13*fp*hp*h4p)/(4480*f*h^5*Pi^2) - (fp*h5p)/(4480*f*h^4*Pi^2) + 
   1/(2520*Pi^2*r^6) + 101/(10080*h^3*Pi^2*r^6) - 1/(96*h*Pi^2*r^6) - 
   fp/(72*f*h^3*Pi^2*r^5) + fp/(72*f*h^2*Pi^2*r^5) + 
   (23*hp)/(672*h^4*Pi^2*r^5) - hp/(72*h^3*Pi^2*r^5) - 
   hp/(1440*h^2*Pi^2*r^5) - (487*fp^2)/(53760*f^2*h^3*Pi^2*r^4) + 
   fpp/(144*f*h^3*Pi^2*r^4) + (5*fp^2)/(768*f^2*h^2*Pi^2*r^4) - 
   fpp/(144*f*h^2*Pi^2*r^4) - (193*fp*hp)/(11520*f*h^4*Pi^2*r^4) + 
   (13*fp*hp)/(1152*f*h^3*Pi^2*r^4) + (1823*hp^2)/(32256*h^5*Pi^2*r^4) - 
   (5*hp^2)/(256*h^4*Pi^2*r^4) - (373*hpp)/(20160*h^4*Pi^2*r^4) + 
   (5*hpp)/(576*h^3*Pi^2*r^4) - (23*fp^3)/(3584*f^3*h^3*Pi^2*r^3) + 
   (69*fp*fpp)/(8960*f^2*h^3*Pi^2*r^3) + (5*fp^3)/(2304*f^3*h^2*Pi^2*r^3) - 
   (fp*fpp)/(384*f^2*h^2*Pi^2*r^3) - (493*fp^2*hp)/(53760*f^2*h^4*Pi^2*r^3) - 
   (7*fpp*hp)/(11520*f*h^4*Pi^2*r^3) + (fp^2*hp)/(384*f^2*h^3*Pi^2*r^3) - 
   (fpp*hp)/(1152*f*h^3*Pi^2*r^3) + (49*fp*hp^2)/(23040*f*h^5*Pi^2*r^3) + 
   (5*fp*hp^2)/(2304*f*h^4*Pi^2*r^3) + (8971*hp^3)/(161280*h^6*Pi^2*r^3) - 
   (7*hp^3)/(576*h^5*Pi^2*r^3) - (7*fp*hpp)/(11520*f*h^4*Pi^2*r^3) - 
   (fp*hpp)/(1152*f*h^3*Pi^2*r^3) - (1133*hp*hpp)/(26880*h^5*Pi^2*r^3) + 
   (13*hp*hpp)/(1152*h^4*Pi^2*r^3) + (101*h3p)/(20160*h^4*Pi^2*r^3) - 
   h3p/(576*h^3*Pi^2*r^3) - (1877*fp^4)/(322560*f^4*h^3*Pi^2*r^2) + 
   (709*fp^2*fpp)/(80640*f^3*h^3*Pi^2*r^2) - 
   (53*fpp^2)/(53760*f^2*h^3*Pi^2*r^2) - 
   (13*fp*f3p)/(13440*f^2*h^3*Pi^2*r^2) - f4p/(1440*f*h^3*Pi^2*r^2) - 
   (49*fp^4)/(23040*f^4*h^2*Pi^2*r^2) + 
   (29*fp^2*fpp)/(5760*f^3*h^2*Pi^2*r^2) - fpp^2/(640*f^2*h^2*Pi^2*r^2) - 
   (fp*f3p)/(480*f^2*h^2*Pi^2*r^2) + f4p/(1440*f*h^2*Pi^2*r^2) - 
   (1031*fp^3*hp)/(322560*f^3*h^4*Pi^2*r^2) - 
   (479*fp*fpp*hp)/(161280*f^2*h^4*Pi^2*r^2) + 
   (7*f3p*hp)/(1152*f*h^4*Pi^2*r^2) - (29*fp^3*hp)/(11520*f^3*h^3*Pi^2*r^2) + 
   (3*fp*fpp*hp)/(640*f^2*h^3*Pi^2*r^2) - (f3p*hp)/(480*f*h^3*Pi^2*r^2) + 
   (2413*fp^2*hp^2)/(322560*f^2*h^5*Pi^2*r^2) - 
   (487*fpp*hp^2)/(23040*f*h^5*Pi^2*r^2) - 
   (19*fp^2*hp^2)/(7680*f^2*h^4*Pi^2*r^2) + 
   (19*fpp*hp^2)/(5760*f*h^4*Pi^2*r^2) + 
   (169*fp*hp^3)/(5120*f*h^6*Pi^2*r^2) - (7*fp*hp^3)/(2880*f*h^5*Pi^2*r^2) + 
   (433*hp^4)/(20160*h^7*Pi^2*r^2) - 
   (313*fp^2*hpp)/(161280*f^2*h^4*Pi^2*r^2) + 
   (73*fpp*hpp)/(11520*f*h^4*Pi^2*r^2) + (fp^2*hpp)/(960*f^2*h^3*Pi^2*r^2) - 
   (fpp*hpp)/(720*f*h^3*Pi^2*r^2) - (181*fp*hp*hpp)/(7680*f*h^5*Pi^2*r^2) + 
   (13*fp*hp*hpp)/(5760*f*h^4*Pi^2*r^2) - (53*hp^2*hpp)/(1680*h^6*Pi^2*r^2) + 
   (709*hpp^2)/(161280*h^5*Pi^2*r^2) + (29*fp*h3p)/(11520*f*h^4*Pi^2*r^2) - 
   (fp*h3p)/(2880*f*h^3*Pi^2*r^2) + (587*hp*h3p)/(80640*h^5*Pi^2*r^2) - 
   (11*h4p)/(13440*h^4*Pi^2*r^2) + (443*fp^5)/(10752*f^5*h^3*Pi^2*r) - 
   (12137*fp^3*fpp)/(107520*f^4*h^3*Pi^2*r) + 
   (487*fp*fpp^2)/(7680*f^3*h^3*Pi^2*r) + 
   (53*fp^2*f3p)/(1280*f^3*h^3*Pi^2*r) - (65*fpp*f3p)/(2688*f^2*h^3*Pi^2*r) - 
   (311*fp*f4p)/(26880*f^2*h^3*Pi^2*r) + (3*f5p)/(1120*f*h^3*Pi^2*r) + 
   (3557*fp^4*hp)/(53760*f^4*h^4*Pi^2*r) - 
   (751*fp^2*fpp*hp)/(5120*f^3*h^4*Pi^2*r) + 
   (1133*fpp^2*hp)/(26880*f^2*h^4*Pi^2*r) + 
   (61*fp*f3p*hp)/(1120*f^2*h^4*Pi^2*r) - (83*f4p*hp)/(5376*f*h^4*Pi^2*r) + 
   (583*fp^3*hp^2)/(7680*f^3*h^5*Pi^2*r) - 
   (13667*fp*fpp*hp^2)/(107520*f^2*h^5*Pi^2*r) + 
   (251*f3p*hp^2)/(5376*f*h^5*Pi^2*r) + 
   (865*fp^2*hp^3)/(10752*f^2*h^6*Pi^2*r) - 
   (3091*fpp*hp^3)/(35840*f*h^6*Pi^2*r) + (411*fp*hp^4)/(4480*f*h^7*Pi^2*r) - 
   (143*hp^5)/(768*h^8*Pi^2*r) - (19*fp^3*hpp)/(768*f^3*h^4*Pi^2*r) + 
   (743*fp*fpp*hpp)/(17920*f^2*h^4*Pi^2*r) - 
   (41*f3p*hpp)/(2688*f*h^4*Pi^2*r) - 
   (209*fp^2*hp*hpp)/(3360*f^2*h^5*Pi^2*r) + 
   (1201*fpp*hp*hpp)/(17920*f*h^5*Pi^2*r) - 
   (523*fp*hp^2*hpp)/(4480*f*h^6*Pi^2*r) + 
   (8767*hp^3*hpp)/(26880*h^7*Pi^2*r) + (67*fp*hpp^2)/(4480*f*h^5*Pi^2*r) - 
   (5423*hp*hpp^2)/(53760*h^6*Pi^2*r) + 
   (407*fp^2*h3p)/(53760*f^2*h^4*Pi^2*r) - 
   (221*fpp*h3p)/(26880*f*h^4*Pi^2*r) + 
   (599*fp*hp*h3p)/(26880*f*h^5*Pi^2*r) - 
   (3949*hp^2*h3p)/(53760*h^6*Pi^2*r) + (33*hpp*h3p)/(1792*h^5*Pi^2*r) - 
   (29*fp*h4p)/(13440*f*h^4*Pi^2*r) + (143*hp*h4p)/(13440*h^5*Pi^2*r) - 
   (11*h5p)/(13440*h^4*Pi^2*r) + (583*fp^6*si)/(2560*f^6*h^3*Pi^2) - 
   (969*fp^4*fpp*si)/(1280*f^5*h^3*Pi^2) + 
   (413*fp^2*fpp^2*si)/(640*f^4*h^3*Pi^2) - 
   (27*fpp^3*si)/(320*f^3*h^3*Pi^2) + (539*fp^3*f3p*si)/(1920*f^4*h^3*Pi^2) - 
   (319*fp*fpp*f3p*si)/(960*f^3*h^3*Pi^2) + 
   (11*f3p^2*si)/(320*f^2*h^3*Pi^2) - (19*fp^2*f4p*si)/(240*f^3*h^3*Pi^2) + 
   (fpp*f4p*si)/(20*f^2*h^3*Pi^2) + (3*fp*f5p*si)/(160*f^2*h^3*Pi^2) - 
   (f6p*si)/(240*f*h^3*Pi^2) + (969*fp^5*hp*si)/(2560*f^5*h^4*Pi^2) - 
   (273*fp^3*fpp*hp*si)/(256*f^4*h^4*Pi^2) + 
   (5*fp*fpp^2*hp*si)/(8*f^3*h^4*Pi^2) + 
   (155*fp^2*f3p*hp*si)/(384*f^3*h^4*Pi^2) - 
   (81*fpp*f3p*hp*si)/(320*f^2*h^4*Pi^2) - 
   (19*fp*f4p*hp*si)/(160*f^2*h^4*Pi^2) + (f5p*hp*si)/(32*f*h^4*Pi^2) + 
   (679*fp^4*hp^2*si)/(1536*f^4*h^5*Pi^2) - 
   (49*fp^2*fpp*hp^2*si)/(48*f^3*h^5*Pi^2) + 
   (403*fpp^2*hp^2*si)/(1280*f^2*h^5*Pi^2) + 
   (257*fp*f3p*hp^2*si)/(640*f^2*h^5*Pi^2) - 
   (25*f4p*hp^2*si)/(192*f*h^5*Pi^2) + 
   (697*fp^3*hp^3*si)/(1536*f^3*h^6*Pi^2) - 
   (131*fp*fpp*hp^3*si)/(160*f^2*h^6*Pi^2) + 
   (45*f3p*hp^3*si)/(128*f*h^6*Pi^2) + (9*fp^2*hp^4*si)/(20*f^2*h^7*Pi^2) - 
   (113*fpp*hp^4*si)/(192*f*h^7*Pi^2) + (91*fp*hp^5*si)/(192*f*h^8*Pi^2) - 
   (539*fp^4*hpp*si)/(3840*f^4*h^4*Pi^2) + 
   (623*fp^2*fpp*hpp*si)/(1920*f^3*h^4*Pi^2) - 
   (fpp^2*hpp*si)/(10*f^2*h^4*Pi^2) - (41*fp*f3p*hpp*si)/(320*f^2*h^4*Pi^2) + 
   (f4p*hpp*si)/(24*f*h^4*Pi^2) - (1307*fp^3*hp*hpp*si)/(3840*f^3*h^5*Pi^2) + 
   (197*fp*fpp*hp*hpp*si)/(320*f^2*h^5*Pi^2) - 
   (17*f3p*hp*hpp*si)/(64*f*h^5*Pi^2) - 
   (89*fp^2*hp^2*hpp*si)/(160*f^2*h^6*Pi^2) + 
   (1403*fpp*hp^2*hpp*si)/(1920*f*h^6*Pi^2) - 
   (797*fp*hp^3*hpp*si)/(960*f*h^7*Pi^2) + 
   (89*fp^2*hpp^2*si)/(1280*f^2*h^5*Pi^2) - 
   (11*fpp*hpp^2*si)/(120*f*h^5*Pi^2) + 
   (493*fp*hp*hpp^2*si)/(1920*f*h^6*Pi^2) + 
   (19*fp^3*h3p*si)/(480*f^3*h^4*Pi^2) - 
   (23*fp*fpp*h3p*si)/(320*f^2*h^4*Pi^2) + (f3p*h3p*si)/(32*f*h^4*Pi^2) + 
   (13*fp^2*hp*h3p*si)/(128*f^2*h^5*Pi^2) - 
   (43*fpp*hp*h3p*si)/(320*f*h^5*Pi^2) + 
   (359*fp*hp^2*h3p*si)/(1920*f*h^6*Pi^2) - 
   (3*fp*hpp*h3p*si)/(64*f*h^5*Pi^2) - (3*fp^2*h4p*si)/(320*f^2*h^4*Pi^2) + 
   (fpp*h4p*si)/(80*f*h^4*Pi^2) - (13*fp*hp*h4p*si)/(480*f*h^5*Pi^2) + 
   (fp*h5p*si)/(480*f*h^4*Pi^2) - si/(240*Pi^2*r^6) - si/(20*h^3*Pi^2*r^6) - 
   (5*si)/(48*h^2*Pi^2*r^6) + (19*si)/(120*h*Pi^2*r^6) + 
   (29*fp*si)/(144*f*h^3*Pi^2*r^5) - (29*fp*si)/(144*f*h^2*Pi^2*r^5) - 
   (257*hp*si)/(720*h^4*Pi^2*r^5) + (7*hp*si)/(48*h^3*Pi^2*r^5) + 
   (hp*si)/(90*h^2*Pi^2*r^5) + (251*fp^2*si)/(1920*f^2*h^3*Pi^2*r^4) - 
   (29*fpp*si)/(288*f*h^3*Pi^2*r^4) - (37*fp^2*si)/(384*f^2*h^2*Pi^2*r^4) + 
   (29*fpp*si)/(288*f*h^2*Pi^2*r^4) + (703*fp*hp*si)/(2880*f*h^4*Pi^2*r^4) - 
   (95*fp*hp*si)/(576*f*h^3*Pi^2*r^4) - (775*hp^2*si)/(1152*h^5*Pi^2*r^4) + 
   (313*hp^2*si)/(1152*h^4*Pi^2*r^4) + (319*hpp*si)/(1440*h^4*Pi^2*r^4) - 
   (35*hpp*si)/(288*h^3*Pi^2*r^4) + (53*fp^3*si)/(576*f^3*h^3*Pi^2*r^3) - 
   (53*fp*fpp*si)/(480*f^2*h^3*Pi^2*r^3) - 
   (5*fp^3*si)/(144*f^3*h^2*Pi^2*r^3) + (fp*fpp*si)/(24*f^2*h^2*Pi^2*r^3) + 
   (127*fp^2*hp*si)/(960*f^2*h^4*Pi^2*r^3) + 
   (11*fpp*hp*si)/(1440*f*h^4*Pi^2*r^3) - 
   (fp^2*hp*si)/(24*f^2*h^3*Pi^2*r^3) + (fpp*hp*si)/(72*f*h^3*Pi^2*r^3) - 
   (77*fp*hp^2*si)/(2880*f*h^5*Pi^2*r^3) - 
   (5*fp*hp^2*si)/(144*f*h^4*Pi^2*r^3) - (127*hp^3*si)/(192*h^6*Pi^2*r^3) + 
   (7*hp^3*si)/(36*h^5*Pi^2*r^3) + (11*fp*hpp*si)/(1440*f*h^4*Pi^2*r^3) + 
   (fp*hpp*si)/(72*f*h^3*Pi^2*r^3) + (81*hp*hpp*si)/(160*h^5*Pi^2*r^3) - 
   (13*hp*hpp*si)/(72*h^4*Pi^2*r^3) - (11*h3p*si)/(180*h^4*Pi^2*r^3) + 
   (h3p*si)/(36*h^3*Pi^2*r^3) + (91*fp^4*si)/(1152*f^4*h^3*Pi^2*r^2) - 
   (329*fp^2*fpp*si)/(2880*f^3*h^3*Pi^2*r^2) + 
   (3*fpp^2*si)/(320*f^2*h^3*Pi^2*r^2) + (fp*f3p*si)/(120*f^2*h^3*Pi^2*r^2) + 
   (f4p*si)/(90*f*h^3*Pi^2*r^2) + (49*fp^4*si)/(1440*f^4*h^2*Pi^2*r^2) - 
   (29*fp^2*fpp*si)/(360*f^3*h^2*Pi^2*r^2) + 
   (fpp^2*si)/(40*f^2*h^2*Pi^2*r^2) + (fp*f3p*si)/(30*f^2*h^2*Pi^2*r^2) - 
   (f4p*si)/(90*f*h^2*Pi^2*r^2) + (403*fp^3*hp*si)/(11520*f^3*h^4*Pi^2*r^2) + 
   (71*fp*fpp*hp*si)/(1152*f^2*h^4*Pi^2*r^2) - 
   (4*f3p*hp*si)/(45*f*h^4*Pi^2*r^2) + 
   (29*fp^3*hp*si)/(720*f^3*h^3*Pi^2*r^2) - 
   (3*fp*fpp*hp*si)/(40*f^2*h^3*Pi^2*r^2) + (f3p*hp*si)/(30*f*h^3*Pi^2*r^2) - 
   (19*fp^2*hp^2*si)/(144*f^2*h^5*Pi^2*r^2) + 
   (1777*fpp*hp^2*si)/(5760*f*h^5*Pi^2*r^2) + 
   (19*fp^2*hp^2*si)/(480*f^2*h^4*Pi^2*r^2) - 
   (19*fpp*hp^2*si)/(360*f*h^4*Pi^2*r^2) - 
   (125*fp*hp^3*si)/(256*f*h^6*Pi^2*r^2) + 
   (7*fp*hp^3*si)/(180*f*h^5*Pi^2*r^2) - (33*hp^4*si)/(320*h^7*Pi^2*r^2) + 
   (53*fp^2*hpp*si)/(1440*f^2*h^4*Pi^2*r^2) - 
   (271*fpp*hpp*si)/(2880*f*h^4*Pi^2*r^2) - 
   (fp^2*hpp*si)/(60*f^2*h^3*Pi^2*r^2) + (fpp*hpp*si)/(45*f*h^3*Pi^2*r^2) + 
   (223*fp*hp*hpp*si)/(640*f*h^5*Pi^2*r^2) - 
   (13*fp*hp*hpp*si)/(360*f*h^4*Pi^2*r^2) + 
   (27*hp^2*hpp*si)/(128*h^6*Pi^2*r^2) - (11*hpp^2*si)/(320*h^5*Pi^2*r^2) - 
   (107*fp*h3p*si)/(2880*f*h^4*Pi^2*r^2) + (fp*h3p*si)/(180*f*h^3*Pi^2*r^2) - 
   (61*hp*h3p*si)/(960*h^5*Pi^2*r^2) + (h4p*si)/(120*h^4*Pi^2*r^2) - 
   (543*fp^5*si)/(1280*f^5*h^3*Pi^2*r) + 
   (371*fp^3*fpp*si)/(320*f^4*h^3*Pi^2*r) - 
   (207*fp*fpp^2*si)/(320*f^3*h^3*Pi^2*r) - 
   (101*fp^2*f3p*si)/(240*f^3*h^3*Pi^2*r) + 
   (29*fpp*f3p*si)/(120*f^2*h^3*Pi^2*r) + 
   (11*fp*f4p*si)/(96*f^2*h^3*Pi^2*r) - (f5p*si)/(40*f*h^3*Pi^2*r) - 
   (427*fp^4*hp*si)/(640*f^4*h^4*Pi^2*r) + 
   (469*fp^2*fpp*hp*si)/(320*f^3*h^4*Pi^2*r) - 
   (33*fpp^2*hp*si)/(80*f^2*h^4*Pi^2*r) - 
   (127*fp*f3p*hp*si)/(240*f^2*h^4*Pi^2*r) + 
   (67*f4p*hp*si)/(480*f*h^4*Pi^2*r) - 
   (947*fp^3*hp^2*si)/(1280*f^3*h^5*Pi^2*r) + 
   (2311*fp*fpp*hp^2*si)/(1920*f^2*h^5*Pi^2*r) - 
   (97*f3p*hp^2*si)/(240*f*h^5*Pi^2*r) - 
   (241*fp^2*hp^3*si)/(320*f^2*h^6*Pi^2*r) + 
   (467*fpp*hp^3*si)/(640*f*h^6*Pi^2*r) - (5*fp*hp^4*si)/(6*f*h^7*Pi^2*r) + 
   (91*hp^5*si)/(48*h^8*Pi^2*r) + (31*fp^3*hpp*si)/(128*f^3*h^4*Pi^2*r) - 
   (379*fp*fpp*hpp*si)/(960*f^2*h^4*Pi^2*r) + 
   (2*f3p*hpp*si)/(15*f*h^4*Pi^2*r) + 
   (377*fp^2*hp*hpp*si)/(640*f^2*h^5*Pi^2*r) - 
   (37*fpp*hp*hpp*si)/(64*f*h^5*Pi^2*r) + 
   (517*fp*hp^2*hpp*si)/(480*f*h^6*Pi^2*r) - 
   (797*hp^3*hpp*si)/(240*h^7*Pi^2*r) - (67*fp*hpp^2*si)/(480*f*h^5*Pi^2*r) + 
   (493*hp*hpp^2*si)/(480*h^6*Pi^2*r) - (7*fp^2*h3p*si)/(96*f^2*h^4*Pi^2*r) + 
   (7*fpp*h3p*si)/(96*f*h^4*Pi^2*r) - (101*fp*hp*h3p*si)/(480*f*h^5*Pi^2*r) + 
   (359*hp^2*h3p*si)/(480*h^6*Pi^2*r) - (3*hpp*h3p*si)/(16*h^5*Pi^2*r) + 
   (fp*h4p*si)/(48*f*h^4*Pi^2*r) - (13*hp*h4p*si)/(120*h^5*Pi^2*r) + 
   (h5p*si)/(120*h^4*Pi^2*r) - (935*fp^6*si^2)/(1536*f^6*h^3*Pi^2) + 
   (527*fp^4*fpp*si^2)/(256*f^5*h^3*Pi^2) - 
   (231*fp^2*fpp^2*si^2)/(128*f^4*h^3*Pi^2) + 
   (47*fpp^3*si^2)/(192*f^3*h^3*Pi^2) - 
   (581*fp^3*f3p*si^2)/(768*f^4*h^3*Pi^2) + 
   (361*fp*fpp*f3p*si^2)/(384*f^3*h^3*Pi^2) - 
   (13*f3p^2*si^2)/(128*f^2*h^3*Pi^2) + 
   (79*fp^2*f4p*si^2)/(384*f^3*h^3*Pi^2) - 
   (9*fpp*f4p*si^2)/(64*f^2*h^3*Pi^2) - (3*fp*f5p*si^2)/(64*f^2*h^3*Pi^2) + 
   (f6p*si^2)/(96*f*h^3*Pi^2) - (527*fp^5*hp*si^2)/(512*f^5*h^4*Pi^2) + 
   (1505*fp^3*fpp*hp*si^2)/(512*f^4*h^4*Pi^2) - 
   (455*fp*fpp^2*hp*si^2)/(256*f^3*h^4*Pi^2) - 
   (835*fp^2*f3p*hp*si^2)/(768*f^3*h^4*Pi^2) + 
   (93*fpp*f3p*hp*si^2)/(128*f^2*h^4*Pi^2) + 
   (39*fp*f4p*hp*si^2)/(128*f^2*h^4*Pi^2) - (5*f5p*hp*si^2)/(64*f*h^4*Pi^2) - 
   (1855*fp^4*hp^2*si^2)/(1536*f^4*h^5*Pi^2) + 
   (2155*fp^2*fpp*hp^2*si^2)/(768*f^3*h^5*Pi^2) - 
   (459*fpp^2*hp^2*si^2)/(512*f^2*h^5*Pi^2) - 
   (271*fp*f3p*hp^2*si^2)/(256*f^2*h^5*Pi^2) + 
   (125*f4p*hp^2*si^2)/(384*f*h^5*Pi^2) - 
   (625*fp^3*hp^3*si^2)/(512*f^3*h^6*Pi^2) + 
   (1119*fp*fpp*hp^3*si^2)/(512*f^2*h^6*Pi^2) - 
   (225*f3p*hp^3*si^2)/(256*f*h^6*Pi^2) - 
   (299*fp^2*hp^4*si^2)/(256*f^2*h^7*Pi^2) + 
   (565*fpp*hp^4*si^2)/(384*f*h^7*Pi^2) - 
   (455*fp*hp^5*si^2)/(384*f*h^8*Pi^2) + 
   (581*fp^4*hpp*si^2)/(1536*f^4*h^4*Pi^2) - 
   (677*fp^2*fpp*hpp*si^2)/(768*f^3*h^4*Pi^2) + 
   (9*fpp^2*hpp*si^2)/(32*f^2*h^4*Pi^2) + 
   (43*fp*f3p*hpp*si^2)/(128*f^2*h^4*Pi^2) - 
   (5*f4p*hpp*si^2)/(48*f*h^4*Pi^2) + 
   (347*fp^3*hp*hpp*si^2)/(384*f^3*h^5*Pi^2) - 
   (417*fp*fpp*hp*hpp*si^2)/(256*f^2*h^5*Pi^2) + 
   (85*f3p*hp*hpp*si^2)/(128*f*h^5*Pi^2) + 
   (733*fp^2*hp^2*hpp*si^2)/(512*f^2*h^6*Pi^2) - 
   (1403*fpp*hp^2*hpp*si^2)/(768*f*h^6*Pi^2) + 
   (797*fp*hp^3*hpp*si^2)/(384*f*h^7*Pi^2) - 
   (91*fp^2*hpp^2*si^2)/(512*f^2*h^5*Pi^2) + 
   (11*fpp*hpp^2*si^2)/(48*f*h^5*Pi^2) - 
   (493*fp*hp*hpp^2*si^2)/(768*f*h^6*Pi^2) - 
   (79*fp^3*h3p*si^2)/(768*f^3*h^4*Pi^2) + 
   (3*fp*fpp*h3p*si^2)/(16*f^2*h^4*Pi^2) - (5*f3p*h3p*si^2)/(64*f*h^4*Pi^2) - 
   (33*fp^2*hp*h3p*si^2)/(128*f^2*h^5*Pi^2) + 
   (43*fpp*hp*h3p*si^2)/(128*f*h^5*Pi^2) - 
   (359*fp*hp^2*h3p*si^2)/(768*f*h^6*Pi^2) + 
   (15*fp*hpp*h3p*si^2)/(128*f*h^5*Pi^2) + 
   (3*fp^2*h4p*si^2)/(128*f^2*h^4*Pi^2) - (fpp*h4p*si^2)/(32*f*h^4*Pi^2) + 
   (13*fp*hp*h4p*si^2)/(192*f*h^5*Pi^2) - (fp*h5p*si^2)/(192*f*h^4*Pi^2) + 
   si^2/(48*Pi^2*r^6) - (5*si^2)/(16*h^3*Pi^2*r^6) + 
   (55*si^2)/(48*h^2*Pi^2*r^6) - (41*si^2)/(48*h*Pi^2*r^6) - 
   (23*fp*si^2)/(24*f*h^3*Pi^2*r^5) + (23*fp*si^2)/(24*f*h^2*Pi^2*r^5) + 
   (43*hp*si^2)/(48*h^4*Pi^2*r^5) - (hp*si^2)/(3*h^3*Pi^2*r^5) - 
   (hp*si^2)/(16*h^2*Pi^2*r^5) - (5*fp^2*si^2)/(8*f^2*h^3*Pi^2*r^4) + 
   (23*fpp*si^2)/(48*f*h^3*Pi^2*r^4) + (15*fp^2*si^2)/(32*f^2*h^2*Pi^2*r^4) - 
   (23*fpp*si^2)/(48*f*h^2*Pi^2*r^4) - (55*fp*hp*si^2)/(48*f*h^4*Pi^2*r^4) + 
   (19*fp*hp*si^2)/(24*f*h^3*Pi^2*r^4) + (227*hp^2*si^2)/(96*h^5*Pi^2*r^4) - 
   (115*hp^2*si^2)/(96*h^4*Pi^2*r^4) - (19*hpp*si^2)/(24*h^4*Pi^2*r^4) + 
   (13*hpp*si^2)/(24*h^3*Pi^2*r^4) - (85*fp^3*si^2)/(192*f^3*h^3*Pi^2*r^3) + 
   (17*fp*fpp*si^2)/(32*f^2*h^3*Pi^2*r^3) + 
   (35*fp^3*si^2)/(192*f^3*h^2*Pi^2*r^3) - 
   (7*fp*fpp*si^2)/(32*f^2*h^2*Pi^2*r^3) - 
   (5*fp^2*hp*si^2)/(8*f^2*h^4*Pi^2*r^3) - 
   (5*fpp*hp*si^2)/(96*f*h^4*Pi^2*r^3) + 
   (7*fp^2*hp*si^2)/(32*f^2*h^3*Pi^2*r^3) - 
   (7*fpp*hp*si^2)/(96*f*h^3*Pi^2*r^3) + 
   (35*fp*hp^2*si^2)/(192*f*h^5*Pi^2*r^3) + 
   (35*fp*hp^2*si^2)/(192*f*h^4*Pi^2*r^3) + 
   (227*hp^3*si^2)/(96*h^6*Pi^2*r^3) - (49*hp^3*si^2)/(48*h^5*Pi^2*r^3) - 
   (5*fp*hpp*si^2)/(96*f*h^4*Pi^2*r^3) - 
   (7*fp*hpp*si^2)/(96*f*h^3*Pi^2*r^3) - (59*hp*hpp*si^2)/(32*h^5*Pi^2*r^3) + 
   (91*hp*hpp*si^2)/(96*h^4*Pi^2*r^3) + (11*h3p*si^2)/(48*h^4*Pi^2*r^3) - 
   (7*h3p*si^2)/(48*h^3*Pi^2*r^3) - (85*fp^4*si^2)/(256*f^4*h^3*Pi^2*r^2) + 
   (7*fp^2*fpp*si^2)/(16*f^3*h^3*Pi^2*r^2) - 
   (fpp^2*si^2)/(64*f^2*h^3*Pi^2*r^2) + (fp*f3p*si^2)/(96*f^2*h^3*Pi^2*r^2) - 
   (f4p*si^2)/(16*f*h^3*Pi^2*r^2) - (49*fp^4*si^2)/(256*f^4*h^2*Pi^2*r^2) + 
   (29*fp^2*fpp*si^2)/(64*f^3*h^2*Pi^2*r^2) - 
   (9*fpp^2*si^2)/(64*f^2*h^2*Pi^2*r^2) - 
   (3*fp*f3p*si^2)/(16*f^2*h^2*Pi^2*r^2) + (f4p*si^2)/(16*f*h^2*Pi^2*r^2) - 
   (25*fp^3*hp*si^2)/(384*f^3*h^4*Pi^2*r^2) - 
   (23*fp*fpp*hp*si^2)/(48*f^2*h^4*Pi^2*r^2) + 
   (15*f3p*hp*si^2)/(32*f*h^4*Pi^2*r^2) - 
   (29*fp^3*hp*si^2)/(128*f^3*h^3*Pi^2*r^2) + 
   (27*fp*fpp*hp*si^2)/(64*f^2*h^3*Pi^2*r^2) - 
   (3*f3p*hp*si^2)/(16*f*h^3*Pi^2*r^2) + 
   (629*fp^2*hp^2*si^2)/(768*f^2*h^5*Pi^2*r^2) - 
   (51*fpp*hp^2*si^2)/(32*f*h^5*Pi^2*r^2) - 
   (57*fp^2*hp^2*si^2)/(256*f^2*h^4*Pi^2*r^2) + 
   (19*fpp*hp^2*si^2)/(64*f*h^4*Pi^2*r^2) + 
   (119*fp*hp^3*si^2)/(48*f*h^6*Pi^2*r^2) - 
   (7*fp*hp^3*si^2)/(32*f*h^5*Pi^2*r^2) - (55*hp^4*si^2)/(96*h^7*Pi^2*r^2) - 
   (23*fp^2*hpp*si^2)/(96*f^2*h^4*Pi^2*r^2) + 
   (47*fpp*hpp*si^2)/(96*f*h^4*Pi^2*r^2) + 
   (3*fp^2*hpp*si^2)/(32*f^2*h^3*Pi^2*r^2) - 
   (fpp*hpp*si^2)/(8*f*h^3*Pi^2*r^2) - 
   (113*fp*hp*hpp*si^2)/(64*f*h^5*Pi^2*r^2) + 
   (13*fp*hp*hpp*si^2)/(64*f*h^4*Pi^2*r^2) + 
   (17*hp^2*hpp*si^2)/(64*h^6*Pi^2*r^2) + (hpp^2*si^2)/(96*h^5*Pi^2*r^2) + 
   (3*fp*h3p*si^2)/(16*f*h^4*Pi^2*r^2) - (fp*h3p*si^2)/(32*f*h^3*Pi^2*r^2) + 
   (hp*h3p*si^2)/(12*h^5*Pi^2*r^2) - (h4p*si^2)/(48*h^4*Pi^2*r^2) + 
   (83*fp^5*si^2)/(64*f^5*h^3*Pi^2*r) - 
   (457*fp^3*fpp*si^2)/(128*f^4*h^3*Pi^2*r) + 
   (129*fp*fpp^2*si^2)/(64*f^3*h^3*Pi^2*r) + 
   (81*fp^2*f3p*si^2)/(64*f^3*h^3*Pi^2*r) - 
   (71*fpp*f3p*si^2)/(96*f^2*h^3*Pi^2*r) - 
   (31*fp*f4p*si^2)/(96*f^2*h^3*Pi^2*r) + (f5p*si^2)/(16*f*h^3*Pi^2*r) + 
   (1495*fp^4*hp*si^2)/(768*f^4*h^4*Pi^2*r) - 
   (815*fp^2*fpp*hp*si^2)/(192*f^3*h^4*Pi^2*r) + 
   (19*fpp^2*hp*si^2)/(16*f^2*h^4*Pi^2*r) + 
   (137*fp*f3p*hp*si^2)/(96*f^2*h^4*Pi^2*r) - 
   (5*f4p*hp*si^2)/(16*f*h^4*Pi^2*r) + 
   (755*fp^3*hp^2*si^2)/(384*f^3*h^5*Pi^2*r) - 
   (385*fp*fpp*hp^2*si^2)/(128*f^2*h^5*Pi^2*r) + 
   (145*f3p*hp^2*si^2)/(192*f*h^5*Pi^2*r) + 
   (445*fp^2*hp^3*si^2)/(256*f^2*h^6*Pi^2*r) - 
   (113*fpp*hp^3*si^2)/(96*f*h^6*Pi^2*r) + 
   (323*fp*hp^4*si^2)/(192*f*h^7*Pi^2*r) - (455*hp^5*si^2)/(96*h^8*Pi^2*r) - 
   (83*fp^3*hpp*si^2)/(128*f^3*h^4*Pi^2*r) + 
   (fp*fpp*hpp*si^2)/(f^2*h^4*Pi^2*r) - (25*f3p*hpp*si^2)/(96*f*h^4*Pi^2*r) - 
   (45*fp^2*hp*hpp*si^2)/(32*f^2*h^5*Pi^2*r) + 
   (197*fpp*hp*hpp*si^2)/(192*f*h^5*Pi^2*r) - 
   (887*fp*hp^2*hpp*si^2)/(384*f*h^6*Pi^2*r) + 
   (797*hp^3*hpp*si^2)/(96*h^7*Pi^2*r) + 
   (5*fp*hpp^2*si^2)/(16*f*h^5*Pi^2*r) - 
   (493*hp*hpp^2*si^2)/(192*h^6*Pi^2*r) + 
   (35*fp^2*h3p*si^2)/(192*f^2*h^4*Pi^2*r) - 
   (7*fpp*h3p*si^2)/(48*f*h^4*Pi^2*r) + 
   (47*fp*hp*h3p*si^2)/(96*f*h^5*Pi^2*r) - 
   (359*hp^2*h3p*si^2)/(192*h^6*Pi^2*r) + (15*hpp*h3p*si^2)/(32*h^5*Pi^2*r) - 
   (5*fp*h4p*si^2)/(96*f*h^4*Pi^2*r) + (13*hp*h4p*si^2)/(48*h^5*Pi^2*r) - 
   (h5p*si^2)/(48*h^4*Pi^2*r) + (121*fp^6*si^3)/(768*f^6*h^3*Pi^2) - 
   (85*fp^4*fpp*si^3)/(128*f^5*h^3*Pi^2) + 
   (49*fp^2*fpp^2*si^3)/(64*f^4*h^3*Pi^2) - 
   (13*fpp^3*si^3)/(96*f^3*h^3*Pi^2) + (7*fp^3*f3p*si^3)/(32*f^4*h^3*Pi^2) - 
   (7*fp*fpp*f3p*si^3)/(16*f^3*h^3*Pi^2) + (f3p^2*si^3)/(16*f^2*h^3*Pi^2) - 
   (fp^2*f4p*si^3)/(32*f^3*h^3*Pi^2) + (fpp*f4p*si^3)/(16*f^2*h^3*Pi^2) + 
   (85*fp^5*hp*si^3)/(256*f^5*h^4*Pi^2) - 
   (35*fp^3*fpp*hp*si^3)/(32*f^4*h^4*Pi^2) + 
   (55*fp*fpp^2*hp*si^3)/(64*f^3*h^4*Pi^2) + 
   (5*fp^2*f3p*hp*si^3)/(16*f^3*h^4*Pi^2) - 
   (3*fpp*f3p*hp*si^3)/(8*f^2*h^4*Pi^2) - 
   (fp*f4p*hp*si^3)/(32*f^2*h^4*Pi^2) + 
   (105*fp^4*hp^2*si^3)/(256*f^4*h^5*Pi^2) - 
   (65*fp^2*fpp*hp^2*si^3)/(64*f^3*h^5*Pi^2) + 
   (7*fpp^2*hp^2*si^3)/(16*f^2*h^5*Pi^2) + 
   (7*fp*f3p*hp^2*si^3)/(32*f^2*h^5*Pi^2) + 
   (265*fp^3*hp^3*si^3)/(768*f^3*h^6*Pi^2) - 
   (71*fp*fpp*hp^3*si^3)/(128*f^2*h^6*Pi^2) + 
   (11*fp^2*hp^4*si^3)/(64*f^2*h^7*Pi^2) - 
   (7*fp^4*hpp*si^3)/(64*f^4*h^4*Pi^2) + 
   (9*fp^2*fpp*hpp*si^3)/(32*f^3*h^4*Pi^2) - 
   (fpp^2*hpp*si^3)/(8*f^2*h^4*Pi^2) - (fp*f3p*hpp*si^3)/(16*f^2*h^4*Pi^2) - 
   (27*fp^3*hp*hpp*si^3)/(128*f^3*h^5*Pi^2) + 
   (23*fp*fpp*hp*hpp*si^3)/(64*f^2*h^5*Pi^2) - 
   (21*fp^2*hp^2*hpp*si^3)/(128*f^2*h^6*Pi^2) + 
   (fp^2*hpp^2*si^3)/(64*f^2*h^5*Pi^2) + (fp^3*h3p*si^3)/(64*f^3*h^4*Pi^2) - 
   (fp*fpp*h3p*si^3)/(32*f^2*h^4*Pi^2) + 
   (fp^2*hp*h3p*si^3)/(64*f^2*h^5*Pi^2) - si^3/(24*Pi^2*r^6) + 
   (37*si^3)/(24*h^3*Pi^2*r^6) - (25*si^3)/(8*h^2*Pi^2*r^6) + 
   (13*si^3)/(8*h*Pi^2*r^6) + (3*fp*si^3)/(2*f*h^3*Pi^2*r^5) - 
   (3*fp*si^3)/(2*f*h^2*Pi^2*r^5) + (hp*si^3)/(8*h^4*Pi^2*r^5) - 
   (hp*si^3)/(4*h^3*Pi^2*r^5) + (hp*si^3)/(8*h^2*Pi^2*r^5) + 
   (fp^2*si^3)/(f^2*h^3*Pi^2*r^4) - (3*fpp*si^3)/(4*f*h^3*Pi^2*r^4) - 
   (3*fp^2*si^3)/(4*f^2*h^2*Pi^2*r^4) + (3*fpp*si^3)/(4*f*h^2*Pi^2*r^4) + 
   (7*fp*hp*si^3)/(4*f*h^4*Pi^2*r^4) - (5*fp*hp*si^3)/(4*f*h^3*Pi^2*r^4) - 
   (17*hp^2*si^3)/(8*h^5*Pi^2*r^4) + (13*hp^2*si^3)/(8*h^4*Pi^2*r^4) + 
   (3*hpp*si^3)/(4*h^4*Pi^2*r^4) - (3*hpp*si^3)/(4*h^3*Pi^2*r^4) + 
   (35*fp^3*si^3)/(48*f^3*h^3*Pi^2*r^3) - 
   (7*fp*fpp*si^3)/(8*f^2*h^3*Pi^2*r^3) - 
   (5*fp^3*si^3)/(16*f^3*h^2*Pi^2*r^3) + 
   (3*fp*fpp*si^3)/(8*f^2*h^2*Pi^2*r^3) + (fp^2*hp*si^3)/(f^2*h^4*Pi^2*r^3) + 
   (fpp*hp*si^3)/(8*f*h^4*Pi^2*r^3) - (3*fp^2*hp*si^3)/(8*f^2*h^3*Pi^2*r^3) + 
   (fpp*hp*si^3)/(8*f*h^3*Pi^2*r^3) - (7*fp*hp^2*si^3)/(16*f*h^5*Pi^2*r^3) - 
   (5*fp*hp^2*si^3)/(16*f*h^4*Pi^2*r^3) - (55*hp^3*si^3)/(24*h^6*Pi^2*r^3) + 
   (7*hp^3*si^3)/(4*h^5*Pi^2*r^3) + (fp*hpp*si^3)/(8*f*h^4*Pi^2*r^3) + 
   (fp*hpp*si^3)/(8*f*h^3*Pi^2*r^3) + (15*hp*hpp*si^3)/(8*h^5*Pi^2*r^3) - 
   (13*hp*hpp*si^3)/(8*h^4*Pi^2*r^3) - (h3p*si^3)/(4*h^4*Pi^2*r^3) + 
   (h3p*si^3)/(4*h^3*Pi^2*r^3) + (55*fp^4*si^3)/(128*f^4*h^3*Pi^2*r^2) - 
   (15*fp^2*fpp*si^3)/(32*f^3*h^3*Pi^2*r^2) - 
   (fpp^2*si^3)/(32*f^2*h^3*Pi^2*r^2) - (fp*f3p*si^3)/(8*f^2*h^3*Pi^2*r^2) + 
   (f4p*si^3)/(8*f*h^3*Pi^2*r^2) + (49*fp^4*si^3)/(128*f^4*h^2*Pi^2*r^2) - 
   (29*fp^2*fpp*si^3)/(32*f^3*h^2*Pi^2*r^2) + 
   (9*fpp^2*si^3)/(32*f^2*h^2*Pi^2*r^2) + 
   (3*fp*f3p*si^3)/(8*f^2*h^2*Pi^2*r^2) - (f4p*si^3)/(8*f*h^2*Pi^2*r^2) - 
   (7*fp^3*hp*si^3)/(64*f^3*h^4*Pi^2*r^2) + 
   (37*fp*fpp*hp*si^3)/(32*f^2*h^4*Pi^2*r^2) - 
   (7*f3p*hp*si^3)/(8*f*h^4*Pi^2*r^2) + 
   (29*fp^3*hp*si^3)/(64*f^3*h^3*Pi^2*r^2) - 
   (27*fp*fpp*hp*si^3)/(32*f^2*h^3*Pi^2*r^2) + 
   (3*f3p*hp*si^3)/(8*f*h^3*Pi^2*r^2) - 
   (213*fp^2*hp^2*si^3)/(128*f^2*h^5*Pi^2*r^2) + 
   (91*fpp*hp^2*si^3)/(32*f*h^5*Pi^2*r^2) + 
   (57*fp^2*hp^2*si^3)/(128*f^2*h^4*Pi^2*r^2) - 
   (19*fpp*hp^2*si^3)/(32*f*h^4*Pi^2*r^2) - 
   (67*fp*hp^3*si^3)/(16*f*h^6*Pi^2*r^2) + 
   (7*fp*hp^3*si^3)/(16*f*h^5*Pi^2*r^2) + (11*hp^4*si^3)/(4*h^7*Pi^2*r^2) + 
   (fp^2*hpp*si^3)/(2*f^2*h^4*Pi^2*r^2) - 
   (7*fpp*hpp*si^3)/(8*f*h^4*Pi^2*r^2) - 
   (3*fp^2*hpp*si^3)/(16*f^2*h^3*Pi^2*r^2) + 
   (fpp*hpp*si^3)/(4*f*h^3*Pi^2*r^2) + 
   (95*fp*hp*hpp*si^3)/(32*f*h^5*Pi^2*r^2) - 
   (13*fp*hp*hpp*si^3)/(32*f*h^4*Pi^2*r^2) - 
   (21*hp^2*hpp*si^3)/(8*h^6*Pi^2*r^2) + (hpp^2*si^3)/(4*h^5*Pi^2*r^2) - 
   (5*fp*h3p*si^3)/(16*f*h^4*Pi^2*r^2) + (fp*h3p*si^3)/(16*f*h^3*Pi^2*r^2) + 
   (hp*h3p*si^3)/(4*h^5*Pi^2*r^2) - (53*fp^5*si^3)/(64*f^5*h^3*Pi^2*r) + 
   (19*fp^3*fpp*si^3)/(8*f^4*h^3*Pi^2*r) - 
   (23*fp*fpp^2*si^3)/(16*f^3*h^3*Pi^2*r) - 
   (3*fp^2*f3p*si^3)/(4*f^3*h^3*Pi^2*r) + (fpp*f3p*si^3)/(2*f^2*h^3*Pi^2*r) + 
   (fp*f4p*si^3)/(8*f^2*h^3*Pi^2*r) - 
   (127*fp^4*hp*si^3)/(128*f^4*h^4*Pi^2*r) + 
   (67*fp^2*fpp*hp*si^3)/(32*f^3*h^4*Pi^2*r) - 
   (19*fpp^2*hp*si^3)/(32*f^2*h^4*Pi^2*r) - 
   (3*fp*f3p*hp*si^3)/(8*f^2*h^4*Pi^2*r) - (f4p*hp*si^3)/(8*f*h^4*Pi^2*r) - 
   (7*fp^3*hp^2*si^3)/(16*f^3*h^5*Pi^2*r) + 
   (fp*fpp*hp^2*si^3)/(16*f^2*h^5*Pi^2*r) + 
   (7*f3p*hp^2*si^3)/(8*f*h^5*Pi^2*r) + 
   (61*fp^2*hp^3*si^3)/(128*f^2*h^6*Pi^2*r) - 
   (71*fpp*hp^3*si^3)/(32*f*h^6*Pi^2*r) + 
   (11*fp*hp^4*si^3)/(8*f*h^7*Pi^2*r) + 
   (5*fp^3*hpp*si^3)/(32*f^3*h^4*Pi^2*r) - 
   (fp*fpp*hpp*si^3)/(16*f^2*h^4*Pi^2*r) - (f3p*hpp*si^3)/(4*f*h^4*Pi^2*r) - 
   (7*fp^2*hp*hpp*si^3)/(32*f^2*h^5*Pi^2*r) + 
   (23*fpp*hp*hpp*si^3)/(16*f*h^5*Pi^2*r) - 
   (21*fp*hp^2*hpp*si^3)/(16*f*h^6*Pi^2*r) + 
   (fp*hpp^2*si^3)/(8*f*h^5*Pi^2*r) - (fpp*h3p*si^3)/(8*f*h^4*Pi^2*r) + 
   (fp*hp*h3p*si^3)/(8*f*h^5*Pi^2*r))
