(* 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))