      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


     


    
