Real*8 TH,X,CO,BETA,THERMAL,KAINEMATIC Open(9,file='fit.csv') TH=4.D0/100.D0 Do X=2.D0/10.D0,15.D0,2.D0/10.D0 THERMAL=SZT(TH,X) Write(9,*)X,',',THERMAL End do Stop End REAL FUNCTION SZT*8(TH, X) real*8 X,XX,S2,Y0,Y1,Y2,Y3,Y4,C1,C2,C3,D0,D1,D2,D3,FFF2 real*8 TH,CO,BETA,FFF1,A,B,C,D,SZ1,SZ2,SZ3 A=2689.D0/4320.D0 B=529.D0/2160.D0 C=89.D0/1440.D0 D=1.D0/120.D0 XX = X*(EXP(+X) + 1.D0)/(EXP(+X) - 1.D0) S2 = 4.D0*X**2/(EXP(+X)+EXP(-X)-2.D0) Y0 = -4.D0 + XX Y1 = -10.D0 + 47.D0/2.D0*XX - 42.D0/5.D0*XX**2 * + 7.D0/10.D0*XX**3 * + 7.D0/5.D0*S2*(-3.D0 + XX) Y2 = -15.D0/2.D0 + 1023.D0/8.D0*XX * - 868.D0/5.D0*XX**2 + 329.D0/5.D0*XX**3 * - 44.D0/5.D0*XX**4 + 11.D0/30.D0*XX**5 * + S2 /30.D0*(-2604.D0 + 3948.D0*XX - 1452.D0*XX**2 * + 143.D0*XX**3) * + S2**2/60.D0*( -528.D0 + 187.D0*XX) Y3 = +15.D0/2.D0 + 2505.D0/8.D0*XX * - 7098.D0/5.D0*XX**2 + 14253.D0/10.D0*XX**3 * - 18594.D0/35.D0*XX**4 + 12059.D0/140.D0*XX**5 * - 128.D0/21.D0*XX**6 + 16.D0/105.D0*XX**7 * + S2*(-7098.D0/10.D0 + 14253.D0/5.D0*XX * - 102267.D0/35.D0*XX**2 + 156767.D0/140.D0*XX**3 * - 1216.D0/7.D0*XX**4 + 64.D0/7.D0*XX**5) * + S2**2*(-18594.D0/35.D0 + 205003.D0/280.D0*XX * - 1920.D0/7.D0*XX**2 + 1024.D0/35.D0*XX**3) * + S2**3*(-544.D0/21.D0 + 992.D0/105.D0*XX) Y4 = -135.D0/32.D0 + 30375.D0/128.D0*XX * - 62391.D0/10.D0*XX**2 + 614727.D0/40.D0*XX**3 * - 124389.D0/10.D0*XX**4 + 355703.D0/80.D0*XX**5 * - 16568.D0/21.D0*XX**6 + 7516.D0/105.D0*XX**7 * - 22.D0/7.D0*XX**8 + 11.D0/210.D0*XX**9 * + S2**4*(-682.D0/7.D0 + 7601.D0/210.D0*XX) * + S2**3*(-70414.D0/21.D0 + 465992.D0/105.D0*XX * - 11792.D0/7.D0*XX**2 + 19778.D0/105.D0*XX**3) * + S2**2*(-124389.D0/10.D0 + 6046951.D0/160.D0*XX * - 248520.D0/7.D0*XX**2 + 481024.D0/35.D0*XX**3 * - 15972.D0/7.D0*XX**4 + 18689.D0/140.D0*XX**5) * + S2 *(-62391.D0/20.D0 + 614727.D0/20.D0*XX * - 1368279.D0/20.D0*XX**2 + 4624139.D0/80.D0*XX**3 * - 157396.D0/7.D0*XX**4 + 30064.D0/ 7.D0*XX**5 * - 2717.D0/7.D0*XX**6 + 2761.D0/210.D0*XX**7) Y5=(-45.D0)/8.D0 - (7515.D0*XX)/32.D0 - (28917.D0*XX**2)/2.D0 + - (795429.D0*XX**3)/8.D0 - (2319993.D0*XX**4)/14.D0 + - (12667283.D0*XX**5)/112.D0 - (806524.D0*XX**6)/21.D0 + - (21310.D0*XX**7)/3.D0 - (46679.D0*XX**8)/63.D0 + - (10853.D0*XX**9)/252.D0 - (58.D0*XX**10)/45.D0 + - (29.D0*XX**11)/1890.D0 + - S2**5*((-20039.D0)/45.D0 + (158369.D0*XX)/945.D0) + - S2**4*((-1447049.D0)/63.D0 + (7499423.D0*XX)/252.D0 - - (513242.D0*XX**2)/45.D0 + (4973819.D0*XX**3)/3780.D0) + - S2**3*((-3427727.D0)/21.D0 + (1321220.D0*XX)/3.D0 - - (25019944.D0*XX**2)/63.D0 + (9756847.D0*XX**3)/63.D0 - - (2407609.D0*XX**4)/90.D0 + (3149197.D0*XX**5)/1890.D0) + - S2**2*((-2319993.D0)/14.D0 + (215343811.D0*XX)/224.D0 - - (12097860.D0*XX**2)/7.D0 + 1363840.D0*XX**3 - - (11296318.D0*XX**4)/21.D0 + (18439247.D0*XX**5)/168.D0 - - (494392.D0*XX**6)/45.D0 + (400171.D0*XX**7)/945.D0) + - S2*((-28917.D0)/4.D0 + (795429.D0*XX)/4.D0 - - (25519923.D0*XX**2)/28.D0 + (164674679.D0*XX**3)/112.D0 - - (7661978.D0*XX**4)/7.D0 + 426200.D0*XX**5 - - (11529713.D0*XX**6)/126.D0 + (2724103.D0*XX**7)/252.D0 - - (29377.D0*XX**8)/45.D0 + (14761.D0*XX**9)/945.D0) Y6=(7425.D0)/256.D0 + (128655.D0*XX)/1024.D0 - - (360675.D0*XX**2)/32.D0 + (50853555.D0*XX**3)/128.D0 - - (45719721.D0*XX**4)/32.D0 + (458203107.D0*XX**5)/256.D0 - - (22251961.D0*XX**6)/21.D0 + (71548297.D0*XX**7)/210.D0 - - (26865067.D0*XX**8)/420.D0 + (7313155.D0*XX**9)/1008.D0 - - (4492.D0*XX**10)/9.D0 + (6361.D0*XX**11)/315.D0 - - (296.D0*XX**12)/675.D0 + (37.D0*XX**13)/9450.D0 + - S2**6*((-1616456.D0)/675.D0 + (34394053.D0*XX)/37800.D0) + - S2**5*((-1551986.D0)/9.D0 + (69474842.D0*XX)/315.D0 - - (19157638.D0*XX**2)/225.D0 + (63282617.D0*XX**3)/6300.D0) - + S2**4*((-832817077.D0)/420.D0 + - (5053390105.D0*XX)/1008.D0 - (39749708.D0*XX**2)/9.D0 + - (1090981471.D0*XX**3)/630.D0 - - (13818464.D0*XX**4)/45.D0 + (12565681.D0*XX**5)/630.D0) + - S2**3*((-378283337.D0)/84.D0 + (2217997207.D0*XX)/105.D0 - - (3599918978.D0*XX**2)/105.D0 + - (6574526345.D0*XX**3)/252.D0 - (93232583.D0*XX**4)/9.D0 + - (690760073.D0*XX**5)/315.D0 - - (31433128.D0*XX**6)/135.D0 + (9151432.D0*XX**7)/945.D0) + - S2**2*((-45719721.D0)/32.D0 + (7789452819.D0*XX)/512.D0 - - (333779415.D0*XX**2)/7.D0 + (2289545504.D0*XX**3)/35.D0 - - (3250673107.D0*XX**4)/70.D0 + - (12425050345.D0*XX**5)/672.D0 - - (38289808.D0*XX**6)/9.D0 + (175550878.D0*XX**7)/315.D0 - - (344840.D0*XX**8)/9.D0 + (1343507.D0*XX**9)/1260.D0) + - S2*((-360675.D0)/64.D0 + (50853555.D0*XX)/64.D0 - - (502916931.D0*XX**2)/64.D0 + - (5956640391.D0*XX**3)/256.D0 - - (422787259.D0*XX**4)/14.D0 + (143096594.D0*XX**5)/7.D0 - - (6635671549.D0*XX**6)/840.D0 + - (1835601905.D0*XX**7)/1008.D0 - (2275198.D0*XX**8)/9.D0 + - (6475498.D0*XX**9)/315.D0 - (201428.D0*XX**10)/225.D0 + - (50431.D0*XX**11)/3150.D0) Y7=(-675.D0)/8.D0 + (6345.D0*XX)/32.D0 + (86751.D0*XX**2)/8.D0 + - (28579473.D0*XX**3)/32.D0 - (463090581.D0*XX**4)/56.D0 + - (8680356807.D0*XX**5)/448.D0 - (407333911.D0*XX**6)/21.D0 + - (304758409.D0*XX**7)/30.D0 - - (14281971623.D0*XX**8)/4620.D0 + - (32154229291.D0*XX**9)/55440.D0 - - (34276642.D0*XX**10)/495.D0 + (36841447.D0*XX**11)/6930.D0 - - (1927084.D0*XX**12)/7425.D0 + (229693.D0*XX**13)/29700.D0 - - (736.D0*XX**14)/5775.D0 + (46.D0*XX**15)/51975.D0 + - S2**7*((-85520348.D0)/5775.D0 + - (294610772.D0*XX)/51975.D0) + - S2**6*((-10523805724.D0)/7425.D0 + - (213515492317.D0*XX)/118800.D0 - - (4029609568.D0*XX**2)/5775.D0 + - (4350781384.D0*XX**3)/51975.D0) + - S2**5*((-11842579811.D0)/495.D0 + - (201191142067.D0*XX)/3465.D0 - - (124724248877.D0*XX**2)/2475.D0 + - (392853355313.D0*XX**3)/19800.D0 - - (627352048.D0*XX**4)/175.D0 + - (4157703008.D0*XX**5)/17325.D0) + - S2**4*((-442741120313.D0)/4620.D0 + - (22218572440081.D0*XX)/55440.D0 - - (303314005058.D0*XX**2)/495.D0 + - (6318713416417.D0*XX**3)/13860.D0 - - (8178544496.D0*XX**4)/45.D0 + - (78006728809.D0*XX**5)/1980.D0 - - (5023602592.D0*XX**6)/1155.D0 + - (1979052664.D0*XX**7)/10395.D0) + - S2**3*((-6924676487.D0)/84.D0 + (9447510679.D0*XX)/15.D0 - - (1913784197482.D0*XX**2)/1155.D0 + - (28906652132609.D0*XX**3)/13860.D0 - - (1422840547741.D0*XX**4)/990.D0 + - (4000723254071.D0*XX**5)/6930.D0 - - (204642831212.D0*XX**6)/1485.D0 + - (28405673924.D0*XX**7)/1485.D0 - - (1636707416.D0*XX**8)/1155.D0 + - (447793256.D0*XX**9)/10395.D0) + - S2**2*((-463090581.D0)/56.D0 + - (147566065719.D0*XX)/896.D0 - - (6110008665.D0*XX**2)/7.D0 + (9752269088.D0*XX**3)/5.D0 - - (157101687853.D0*XX**4)/70.D0 + - (54630035565409.D0*XX**5)/36960.D0 - - (292174096408.D0*XX**6)/495.D0 + - (508375127153.D0*XX**7)/3465.D0 - - (2245052860.D0*XX**8)/99.D0 + - (758216593.D0*XX**9)/360.D0 - - (18723104.D0*XX**10)/175.D0 + - (39369376.D0*XX**11)/17325.D0) + - S2*((86751.D0)/16.D0 + (28579473.D0*XX)/16.D0 - - (5093996391.D0*XX**2)/112.D0 + - (112844638491.D0*XX**3)/448.D0 - - (7739344309.D0*XX**4)/14.D0 + 609516818.D0*XX**5 - - (3527646990881.D0*XX**6)/9240.D0 + - (8070711552041.D0*XX**7)/55440.D0 - - (17361119173.D0*XX**8)/495.D0 + - (18752296523.D0*XX**9)/3465.D0 - - (1311380662.D0*XX**10)/2475.D0 + - (313071559.D0*XX**11)/9900.D0 - - (6023792.D0*XX**12)/5775.D0 + (753296.D0*XX**13)/51975.D0) Y8= (1905525.D0)/8192.D0 - (31843125.D0*XX)/32768.D0 - - (1274427.D0*XX**2)/256.D0 + (723764619.D0*XX**3)/1024.D0 - - (56131109271.D0*XX**4)/1792.D0 + - (2131228533597.D0*XX**5)/14336.D0 - - (28385005515.D0*XX**6)/112.D0 + - (33759855933.D0*XX**7)/160.D0 - - (7332664403233.D0*XX**8)/73920.D0 + - (8454102129551.D0*XX**9)/295680.D0 - - (159273899.D0*XX**10)/30.D0 + - (3623853049.D0*XX**11)/5544.D0 - - (133798003.D0*XX**12)/2475.D0 + - (118712629.D0*XX**13)/39600.D0 - - (631168.D0*XX**14)/5775.D0 + (1312.D0*XX**15)/525.D0 - - (16.D0*XX**16)/495.D0 + (4.D0*XX**17)/22275.D0 + - S2**8*((-51236656.D0)/495.D0 + - (887722324.D0*XX)/22275.D0) + - S2**7*((-73339275824.D0)/5775.D0 + - (8402811584.D0*XX)/525.D0 - - (3089759392.D0*XX**2)/495.D0 + - (16899213088.D0*XX**3)/22275.D0) + - S2**6*((-730670894383.D0)/2475.D0 + - (110351579826901.D0*XX)/158400.D0 - - (3455653005184.D0*XX**2)/5775.D0 + - (124091851648.D0*XX**3)/525.D0 - - (1949035328.D0*XX**4)/45.D0 + - (66076862872.D0*XX**5)/22275.D0) + - S2**5*((-110058264209.D0)/60.D0 + - (19789861500589.D0*XX)/2772.D0 - - (34638563602661.D0*XX**2)/3300.D0 + - (203039076596489.D0*XX**3)/26400.D0 - - (537995295424.D0*XX**4)/175.D0 + - (118584920576.D0*XX**5)/175.D0 - - (38077243744.D0*XX**6)/495.D0 + - (7075408016.D0*XX**7)/2025.D0) + - S2**4*((-227312596500223.D0)/73920.D0 + - (5841784571519741.D0*XX)/295680.D0 - - (1409414732251.D0*XX**2)/30.D0 + - (621530660287039.D0*XX**3)/11088.D0 - - (567838724732.D0*XX**4)/15.D0 + - (40316352072577.D0*XX**5)/2640.D0 - - (4308066848896.D0*XX**6)/1155.D0 + - (56446023808.D0*XX**7)/105.D0 - - (4142724512.D0*XX**8)/99.D0 + - (6013663016.D0*XX**9)/4455.D0) + - S2**3*((-482545093755.D0)/448.D0 + - (1046555533923.D0*XX)/80.D0 - - (491288515016611.D0*XX**2)/9240.D0 + - (7600237814466349.D0*XX**3)/73920.D0 - - (13223078368879.D0*XX**4)/120.D0 + - (393525074150057.D0*XX**5)/5544.D0 - - (14208411332579.D0*XX**6)/495.D0 + - (3670238350793.D0*XX**7)/495.D0 - - (1403583351008.D0*XX**8)/1155.D0 + - (12771842432.D0*XX**9)/105.D0 - - (3341029504.D0*XX**10)/495.D0 + - (320768096.D0*XX**11)/2025.D0) + - S2**2*((-56131109271.D0)/1792.D0 + - (36230885071149.D0*XX)/28672.D0 - - (1277325248175.D0*XX**2)/112.D0 + - (202559135598.D0*XX**3)/5.D0 - - (80659308435563.D0*XX**4)/1120.D0 + - (14363519518107149.D0*XX**5)/197120.D0 - - (678825357538.D0*XX**6)/15.D0 + - (50005548223151.D0*XX**7)/2772.D0 - - (155874673495.D0*XX**8)/33.D0 + - (391870388329.D0*XX**9)/480.D0 - - (16056282752.D0*XX**10)/175.D0 + - (1122883072.D0*XX**11)/175.D0 - - (124946368.D0*XX**12)/495.D0 + - (94627042.D0*XX**13)/22275.D0) + - S2*((-1274427.D0)/512.D0 + (723764619.D0*XX)/512.D0 - - (617442201981.D0*XX**2)/3584.D0 + - (27705970936761.D0*XX**3)/14336.D0 - - (1617945314355.D0*XX**4)/224.D0 + - (101279567799.D0*XX**5)/8.D0 - - (1811168107598551.D0*XX**6)/147840.D0 + - (2121979634517301.D0*XX**7)/295680.D0 - - (161344459687.D0*XX**8)/60.D0 + - (1844541201941.D0*XX**9)/2772.D0 - - (182099082083.D0*XX**10)/1650.D0 + - (161805313327.D0*XX**11)/13200.D0 - - (5165794496.D0*XX**12)/5775.D0 + - (21485312.D0*XX**13)/525.D0 - (524152.D0*XX**14)/495.D0 + - (23828.D0*XX**15)/2025.D0) Y9=(-91125.D0)/128.D0 + (1451925.D0*XX)/512.D0 - - (95823.D0*XX**2)/8.D0 - (21387969.D0*XX**3)/32.D0 - - (3867907059.D0*XX**4)/56.D0 + - (361018793313.D0*XX**5)/448.D0 - - (68346357865.D0*XX**6)/28.D0 + - (129419653687.D0*XX**7)/40.D0 - - (14191595238489.D0*XX**8)/6160.D0 + - (72688716977749.D0*XX**9)/73920.D0 - - (1153370108027.D0*XX**10)/4290.D0 + - (393223901251.D0*XX**11)/8008.D0 - - (197855054569.D0*XX**12)/32175.D0 + - (275619041167.D0*XX**13)/514800.D0 - - (813344128.D0*XX**14)/25025.D0 + - (102350176.D0*XX**15)/75075.D0 - - (248008.D0*XX**16)/6435.D0 + (22448.D0*XX**17)/32175.D0 - - (268.D0*XX**18)/36855.D0 + (67.D0*XX**19)/2.027025D6 + - S2**9*((-29738697854.D0)/36855.D0 + - (632763613814.D0*XX)/2.027025D6) + - S2**8*((-794193786328.D0)/6435.D0 + - (4981897682288.D0*XX)/32175.D0 - - (57273525044.D0*XX**2)/945.D0 + - (7451199109.D0*XX**3)/1001.D0) + - S2**7*((-94507435965104.D0)/25025.D0 + - (655510094906432.D0*XX)/75075.D0 - - (47892815455696.D0*XX**2)/6435.D0 + - (94838383849856.D0*XX**3)/32175.D0 - - (745401953564.D0*XX**4)/1365.D0 + - (8575833693404.D0*XX**5)/225225.D0) + - S2**6*((-1080486453001309.D0)/32175.D0 + - (256206916478567023.D0*XX)/2.0592D6 - - (4453069674273664.D0*XX**2)/25025.D0 + - (9680505225867904.D0*XX**3)/75075.D0 - - (30211022101664.D0*XX**4)/585.D0 + - (370823354437664.D0*XX**5)/32175.D0 - - (468962511088.D0*XX**6)/351.D0 + - (6024995840752.D0*XX**7)/96525.D0) + - S2**5*((-796978744646657.D0)/8580.D0 + - (2147395724731711.D0*XX)/4004.D0 - - (51222101512204703.D0*XX**2)/42900.D0 + - (471402546488607947.D0*XX**3)/343200.D0 - - (2079836023490112.D0*XX**4)/2275.D0 + - (9250905100533248.D0*XX**5)/25025.D0 - - (590216316653872.D0*XX**6)/6435.D0 + - (39707189785792.D0*XX**7)/2925.D0 - - (1921573868732.D0*XX**8)/1755.D0 + - (395481352652.D0*XX**9)/10725.D0) + - S2**4*((-439939452393159.D0)/6160.D0 + - (50227903431624559.D0*XX)/73920.D0 - - (10206172085930923.D0*XX**2)/4290.D0 + - (67442224527460261.D0*XX**3)/16016.D0 - - (839696851590836.D0*XX**4)/195.D0 + - (93603809427848371.D0*XX**5)/34320.D0 - - (5551518572838016.D0*XX**6)/5005.D0 + - (4403399749427584.D0*XX**7)/15015.D0 - - (64214301298256.D0*XX**8)/1287.D0 + - (33748676845792.D0*XX**9)/6435.D0 - - (180525698872.D0*XX**10)/585.D0 + - (498306362599.D0*XX**11)/64350.D0) + - S2**3*((-1161888083705.D0)/112.D0 + - (4012009264297.D0*XX)/20.D0 - - (950836880978763.D0*XX**2)/770.D0 + - (65347156562996351.D0*XX**3)/18480.D0 - - (95753939738509567.D0*XX**4)/17160.D0 + - (42701363108549843.D0*XX**5)/8008.D0 - - (21010821809845817.D0*XX**6)/6435.D0 + - (8521313895760139.D0*XX**7)/6435.D0 - - (1808704301708768.D0*XX**8)/5005.D0 + - (996341707895936.D0*XX**9)/15015.D0 - - (51787627826752.D0*XX**10)/6435.D0 + - (1800150554752.D0*XX**11)/2925.D0 - - (15664038473.D0*XX**12)/585.D0 + - (9717618829.D0*XX**13)/19305.D0) + - S2**2*((-3867907059.D0)/56.D0 + - (6137319486321.D0*XX)/896.D0 - - (3075586103925.D0*XX**2)/28.D0 + - (3106071688488.D0*XX**3)/5.D0 - - (468322642870137.D0*XX**4)/280.D0 + - (123498130145195551.D0*XX**5)/49280.D0 - - (4915663400411074.D0*XX**6)/2145.D0 + - (5426096613362549.D0*XX**7)/4004.D0 - - (230501138572885.D0*XX**8)/429.D0 + - (909818454892267.D0*XX**9)/6240.D0 - - (62071983816576.D0*XX**10)/2275.D0 + - (87597012230656.D0*XX**11)/25025.D0 - - (1936731177184.D0*XX**12)/6435.D0 + - (531046959704.D0*XX**13)/32175.D0 - - (6383914904.D0*XX**14)/12285.D0 + - (48599522.D0*XX**15)/6825.D0) + - S2*((-95823.D0)/16.D0 - (21387969.D0*XX)/16.D0 - - (42546977649.D0*XX**2)/112.D0 + - (4693244313069.D0*XX**3)/448.D0 - - (3895742398305.D0*XX**4)/56.D0 + - (388258961061.D0*XX**5)/2.D0 - - (3505324023906783.D0*XX**6)/12320.D0 + - (18244867961414999.D0*XX**7)/73920.D0 - - (1168363919431351.D0*XX**8)/8580.D0 + - (200150965736759.D0*XX**9)/4004.D0 - - (269280729268409.D0*XX**10)/21450.D0 + - (375668753110621.D0*XX**11)/171600.D0 - - (6656815015616.D0*XX**12)/25025.D0 + - (1676086482176.D0*XX**13)/75075.D0 - - (8124618076.D0*XX**14)/6435.D0 + - (133722736.D0*XX**15)/2925.D0 - - (780550.D0*XX**16)/819.D0 + (1951442.D0*XX**17)/225225.D0) Y10=(169255575.D0)/65536.D0 - (2090024775.D0*XX)/262144.D0 + - (433967625.D0*XX**2)/8192.D0 + - (8992650375.D0*XX**3)/32768.D0 - - (3162934444995.D0*XX**4)/57344.D0 + - (1353148643034945.D0*XX**5)/458752.D0 - - (15614127041155.D0*XX**6)/896.D0 + - (9649428040913.D0*XX**7)/256.D0 - - (1600697595147911.D0*XX**8)/39424.D0 + - (11942478518370683.D0*XX**9)/473088.D0 - - (135515512037513.D0*XX**10)/13728.D0 + - (985744536107759.D0*XX**11)/384384.D0 - - (330716455601941.D0*XX**12)/720720.D0 + - (667649235923203.D0*XX**13)/1.153152D7 - - (3734552464.D0*XX**14)/715.D0 + - (1173672524.D0*XX**15)/3465.D0 - - (101283718.D0*XX**16)/6435.D0 + - (46573313.D0*XX**17)/90090.D0 - - (231496.D0*XX**18)/19845.D0 + - (486254.D0*XX**19)/2.837835D6 - - (632.D0*XX**20)/429975.D0 + (79.D0*XX**21)/1.4189175D7 + - S2**10*((-2984377641272.D0)/429975.D0 + - (76502310785333.D0*XX)/2.837835D7) + - S2**9*((-25688020889588.D0)/19845.D0 + - (4592296093604668.D0*XX)/2.837835D6 - - (7805173973924.D0*XX**2)/12285.D0 + - (223349568314737.D0*XX**3)/2.837835D6) + - S2**8*((-324339938597938.D0)/6435.D0 + - (10336042413184853.D0*XX)/90090.D0 - - (643140654464984.D0*XX**2)/6615.D0 + - (270386221757290.D0*XX**3)/7007.D0 - - (68913020463512.D0*XX**4)/9555.D0 + - (22981984955077.D0*XX**5)/45045.D0) + - S2**7*((-433940524926002.D0)/715.D0 + - (7516881921104968.D0*XX)/3465.D0 - - (19558894934198716.D0*XX**2)/6435.D0 + - (98381542575140068.D0*XX**3)/45045.D0 - - (643871532247208.D0*XX**4)/735.D0 + - (62239305026156248.D0*XX**5)/315315.D0 - - (666773973312704.D0*XX**6)/28665.D0 + - (1051014654079528.D0*XX**7)/945945.D0) + - S2**6*((-1806042564042199801.D0)/720720.D0 + - (620626032587895889507.D0*XX)/4.612608D7 - - (20446723289582032.D0*XX**2)/715.D0 + - (111008534094162896.D0*XX**3)/3465.D0 - - (12337846533324344.D0*XX**4)/585.D0 + - (384677302074466867.D0*XX**5)/45045.D0 - - (405085617413536.D0*XX**6)/189.D0 + - (6246648886032032.D0*XX**7)/19305.D0 - - (109783051647664.D0*XX**8)/4095.D0 + - (72050541791741.D0*XX**9)/77220.D0) + - S2**5*((-93641218817921483.D0)/27456.D0 + - (5383150911684471899.D0*XX)/192192.D0 - - (85618191041419699667.D0*XX**2)/960960.D0 + - (1141907861818126942223.D0*XX**3)/7.68768D6 - - (9549779089621656.D0*XX**4)/65.D0 + - (106082212683516352.D0*XX**5)/1155.D0 - - (241037801099035012.D0*XX**6)/6435.D0 + - (41190649016484626.D0*XX**7)/4095.D0 - - (237119757098072.D0*XX**8)/135.D0 + - (2870214770932024.D0*XX**9)/15015.D0 - - (11439343978606.D0*XX**10)/975.D0 + - (19889903014249.D0*XX**11)/64350.D0) + - S2**4*((-49621625449585241.D0)/39424.D0 + - (8252252656194141953.D0*XX)/473088.D0 - - (1199176766019952537.D0*XX**2)/13728.D0 + - (169066031132377853849.D0*XX**3)/768768.D0 - - (350890159393659401.D0*XX**4)/1092.D0 + - (226742359959586740439.D0*XX**5)/768768.D0 - - (25490363366997808.D0*XX**6)/143.D0 + - (50494776854039216.D0*XX**7)/693.D0 - - (26224408826568476.D0*XX**8)/1287.D0 + - (35009526240086501.D0*XX**9)/9009.D0 - - (155936482037584.D0*XX**10)/315.D0 + - (1808234791337419.D0*XX**11)/45045.D0 - - (1092354818536.D0*XX**12)/585.D0 + - (5108388422768.D0*XX**13)/135135.D0) + - S2**3*((-265440159699635.D0)/3584.D0 + - (299132269268303.D0*XX)/128.D0 - - (107246738874910037.D0*XX**2)/4928.D0 + - (10736288188015244017.D0*XX**3)/118272.D0 - - (11250633324866366773.D0*XX**4)/54912.D0 + - (107044956409549873087.D0*XX**5)/384384.D0 - - (35119772569736920613.D0*XX**6)/144144.D0 + - (20641711427037667151.D0*XX**7)/144144.D0 - - (8304850153899284.D0*XX**8)/143.D0 + - (11425274804341264.D0*XX**9)/693.D0 - - (21149493132050992.D0*XX**10)/6435.D0 + - (1867404116927756.D0*XX**11)/4095.D0 - - (13530456158006.D0*XX**12)/315.D0 + - (70525836210098.D0*XX**13)/27027.D0 - - (2634688733384.D0*XX**14)/28665.D0 + - (1344669065548.D0*XX**15)/945945.D0) + - S2**2*((-3162934444995.D0)/57344.D0 + - (23003526931594065.D0*XX)/917504.D0 - - (702635716851975.D0*XX**2)/896.D0 + - (28948284122739.D0*XX**3)/4.D0 - - (52823020639881063.D0*XX**4)/1792.D0 + - (20290271002711790417.D0*XX**5)/315392.D0 - - (288783556151940203.D0*XX**6)/3432.D0 + - (13602288853750966441.D0*XX**7)/192192.D0 - - (1926423353881306325.D0*XX**8)/48048.D0 + - (2203910127782493103.D0*XX**9)/139776.D0 - - (285009840395088.D0*XX**10)/65.D0 + - (1004494671700544.D0*XX**11)/1155.D0 - - (790939543852264.D0*XX**12)/6435.D0 + - (2203547422665073.D0*XX**13)/180180.D0 - - (5514368524688.D0*XX**14)/6615.D0 + - (352712118964.D0*XX**15)/9555.D0 - - (100090472.D0*XX**16)/105.D0 + - (6847741883.D0*XX**17)/630630.D0) + - S2*((433967625.D0)/16384.D0 + - (8992650375.D0*XX)/16384.D0 - - (34792278894945.D0*XX**2)/114688.D0 + - (17590932359454285.D0*XX**3)/458752.D0 - - (890005241345835.D0*XX**4)/1792.D0 + - (144741420613695.D0*XX**5)/64.D0 - - (395372306001534017.D0*XX**6)/78848.D0 + - (2997562108111041433.D0*XX**7)/473088.D0 - - (137277213694000669.D0*XX**8)/27456.D0 + - (501743968878849331.D0*XX**9)/192192.D0 - - (450105096074241701.D0*XX**10)/480480.D0 + - (910005908563325689.D0*XX**11)/3.84384D6 - - (30565444641608.D0*XX**12)/715.D0 + - (19220061253024.D0*XX**13)/3465.D0 - - (3318003959821.D0*XX**14)/6435.D0 + - (39633889363.D0*XX**15)/1170.D0 - - (674232100.D0*XX**16)/441.D0 + - (14162634004.D0*XX**17)/315315.D0 - - (66268676.D0*XX**18)/85995.D0 + - (2366761.D0*XX**19)/405405.D0) FFF1 = X**3/(EXP(+X)-1.D0) FFF2 = X*EXP(+X)/(EXP(+X) - 1.D0) SZ1=FFF1*(TH * FFF2 * (Y0 + TH*Y1 + TH**2*Y2 + TH**3*Y3 +TH**4*Y4 *+A*TH**5*Y5 +B*TH**6*Y6+ C*TH**7*Y7 + D*TH**8*Y8) *) SZ2=FFF1*(TH * FFF2 * (Y0 + TH*Y1 + TH**2*Y2 + TH**3*Y3 +TH**4*Y4 *+TH**5*Y5 +A*TH**6*Y6+ B*TH**7*Y7 + C*TH**8*Y8 +D*TH**9*Y9) *) SZ3=FFF1*(TH * FFF2 * (Y0 + TH*Y1 + TH**2*Y2 + TH**3*Y3 +TH**4*Y4 *+TH**5*Y5 +TH**6*Y6+ A*TH**7*Y7 + B*TH**8*Y8 +C*TH**9*Y9 *+D*TH**10*Y10)) SZT=(SZ1+2.D0*SZ2+SZ3)/4.D0 END