IMPLICIT REAL*8 (A-K, O-Z)
      Open(9,file='fit.csv')
      DATA 
     :  pi, kB, eVmc2,
     &  e, c, h
     =    / 3.141592653589793238D0, 1.3806503D-23, 0.510998902D6,
     &      1.602176462D-19, 2.99792458D8, 6.62606876D-34 /


c      Jmc2 = eVmc2 * e
	TH = 0.03

	DO X=0.D0, 15.1D0, 0.2D0

	WRITE(9,*) X,',', DSZ(TH,X)

	ENDDO

      END

C===== Double SZ Subroutine=======================================
      REAL FUNCTION DSZ*8(TH,X)
      IMPLICIT REAL*8 (A-K, O-Z)

	DSZ = 1.D0/2.D0*TH**2 * FFF(X) * SUMZ(TH,X)
      RETURN
	END

C===== X^3/(e^X-1) * Xe^X/(e^X-1) =====================
C
      REAL FUNCTION FFF*8(X)
	IMPLICIT REAL*8 (A-K, O-Z)

      IF(X.EQ.0.D0)THEN
	FFF2 = EXP(+X)
	FFF1 = X**2
	ELSEIF(X.GT.0.D0)THEN
      FFF2 = X*EXP(+X)/(EXP(+X) - 1.D0)
      FFF1 = X**3/(EXP(+X) - 1.D0)
	ENDIF
      FFF = FFF1*FFF2

	RETURN
	END

C===== Xe^X/(e^X-1) =====================
C
      REAL FUNCTION FF*8(X)
	IMPLICIT REAL*8 (A-K, O-Z)

      IF(X.EQ.0.D0)THEN
	F = EXP(+X)
	ELSEIF(X.GT.0.D0)THEN
      F = X*EXP(+X)/(EXP(+X) - 1.D0)
	ENDIF
      FF = F

	RETURN
	END


C===== ƒ°(\theta_e^k Z_k) ==================================	

      REAL FUNCTION SUMZ*8(TH,X)
      IMPLICIT REAL*8 (A-K, O-Z)
	DIMENSION Z(0:10)
      A=2689.D0/4320.D0
	B=529.D0/2160.D0
	C=89.D0/1440.D0
	D=1.D0/120.D0
	IF(X.EQ.0.D0)THEN
	XX = 2.D0
	SS = 2.D0
	ELSEIF(X.GT.0.D0)THEN
      XX = X*(EXP(+X) + 1.D0)/(EXP(+X) - 1.D0)
      SS = 2.D0*X/(EXP(+X/2.D0)-EXP(-X/2.D0))
	ENDIF
C-----------------------------------------------------

      Z(0)=-16.D0 + 34.D0*XX - 12.D0*XX ** 2 + XX ** 3
     * + SS ** 2*(-6.D0 + 2*XX)

	Z(1)=-80.D0 + 590.D0*XX - (3492.D0*XX ** 2)/5.D0
     * + (1271.D0*XX ** 3)/5.D0 - (168.D0*XX ** 4)/5.D0
     * + (7.D0*XX ** 5)/5.D0 + SS ** 4*(-168.D0/5.D0
     * + (119.D0*XX)/10.D0) + 
     * SS ** 2*(-1746.D0/5.D0 + (2542.D0*XX)/5.D0
     * - (924.D0*XX ** 2)/5.D0 + (91.D0*XX ** 3)/5.D0)

	Z(2)=-160.D0 + 4792.D0*XX - (357144.D0*XX ** 2)/25.D0
     * + (312912.D0*XX ** 3)/25.D0 - (110196.D0*XX ** 4)/25.D0
     * + (34873.D0*XX ** 5)/50.D0 - (734.D0*XX ** 6)/15.D0
     * + (367.D0*XX ** 7)/300.D0 + 
     *SS ** 6*(-6239.D0/30.D0 + (11377.D0*XX)/150.D0)  + 
     *SS ** 4*(-110196.D0/25.D0 + (592841.D0*XX)/100.D0
     * - 2202.D0*XX ** 2 + (5872.D0*XX ** 3)/25.D0)  + 
     *SS ** 2*(-178572.D0/25.D0 + (625824.D0*XX)/25.D0
     * - (606078.D0*XX ** 2)/25.D0 + (453349.D0*XX ** 3)/50.D0
     * - (6973.D0*XX ** 4)/5.D0 + (367.D0*XX ** 5)/5.D0) 

	Z(3)=-90.D0 + (96651.D0*XX)/4.D0 - (8659449.D0*XX ** 2)/50.D0
     * + (62384943.D0*XX ** 3)/200.D0 - (38586081.D0*XX ** 4)/175.D0
     * + (103117227.D0*XX ** 5)/1400.D0 - (1325008.D0*XX ** 6)/105.D0
     * + (590831.D0*XX ** 7)/525.D0 - (1718.D0*XX ** 8)/35.D0
     * + (859.D0*XX ** 9)/1050.D0 + 
     *SS ** 8*(-53258.D0/35.D0 + (593569.D0*XX)/1050.D0)  + 
     *SS ** 6*(-5631284.D0/105.D0 + (36631522.D0*XX)/525.D0
     * - (920848.D0*XX ** 2)/35.D0 + (1544482.D0*XX ** 3)/525.D0)  + 
     *SS ** 4*(-38586081.D0/175.D0 + (1752992859.D0*XX)/2800.D0
     * - (3975024.D0*XX ** 2)/7.D0 + (37813184.D0*XX ** 3)/175.D0
     * - (1247268.D0*XX ** 4)/35.D0 + (1459441.D0*XX ** 5)/700.D0)  + 
     *SS ** 2*(-8659449.D0/100.D0 + (62384943.D0*XX)/100.D0
     * - (424446891.D0*XX ** 2)/350.D0 + (1340523951.D0*XX ** 3)/1400.D0
     * - (12587576.D0*XX ** 4)/35.D0 + (2363324.D0*XX ** 5)/35.D0
     * - (212173.D0*XX ** 6)/35.D0 + (215609.D0*XX ** 7)/1050.D0) 

	Z(4)=60.D0 + 82497.D0*XX - (36883086.D0*XX ** 2)/25.D0
     * + (129233103.D0*XX ** 3)/25.D0 - (1154992263.D0*XX ** 4)/175.D0
     * + (5504779501.D0*XX ** 5)/1400.D0 - (129898756.D0*XX ** 6)/105.D0
     * + (114929504.D0*XX ** 7)/525.D0 - (2337809.D0*XX ** 8)/105.D0
     * + (2681837.D0*XX ** 9)/2100.D0 - (2851.D0*XX ** 10)/75.D0
     * + (2851.D0*XX ** 11)/6300.D0 + 
     *SS ** 10*(-1970041.D0/150.D0 + (15569311.D0*XX)/3150.D0)  + 
     *SS ** 8*(-72472079.D0/105.D0 + (1853149367.D0*XX)/2100.D0
     * - (25228499.D0*XX ** 2)/75.D0 + (488977861.D0*XX ** 3)/12600.D0)
     *  + SS ** 6*(-552069713.D0/105.D0 + (7125629248.D0*XX)/525.D0
     * - (1253065624.D0*XX ** 2)/105.D0 + (2410971463.D0*XX ** 3)/525.D0
     * - (236692871.D0*XX ** 4)/300.D0 + (309598643.D0*XX ** 5)/6300.D0)
     *  + SS ** 4*(-1154992263.D0/175.D0 + (93581251517.D0*XX)/2800.D0
     * - (389696268.D0*XX ** 2)/7.D0 + (7355488256.D0*XX ** 3)/175.D0
     * - (565749778.D0*XX ** 4)/35.D0 + (4556441063.D0*XX ** 5)/1400.D0
     * - (24301924.D0*XX ** 6)/75.D0 + (39340949.D0*XX ** 7)/3150.D0)  + 
     *SS ** 2*(-18441543.D0/25.D0 + (258466206.D0*XX)/25.D0
     * - (12704914893.D0*XX ** 2)/350.D0
     * + (71562133513.D0*XX ** 3)/1400.D0
     * - (1234038182.D0*XX ** 4)/35.D0 + (459718016.D0*XX ** 5)/35.D0
     * - (577438823.D0*XX ** 6)/210.D0 + (673141087.D0*XX ** 7)/2100.D0
     * - (2888063.D0*XX ** 8)/150.D0 + (1451159.D0*XX ** 9)/3150.D0) 

	Z(5)=-135.D0/8.D0 + (12368565.D0*XX)/64.D0
     * - (1523246139.D0*XX ** 2)/160.D0
     * + (41024053941.D0*XX ** 3)/640.D0
     * - (78913341669.D0*XX ** 4)/560.D0
     * + (124274226315.D0*XX ** 5)/896.D0
     * - (2505515368.D0*XX ** 6)/35.D0
     * + (29643451897.D0*XX ** 7)/1400.D0 - (105635617.D0*XX ** 8)/28.D0
     * + (10436409287.D0*XX ** 9)/25200.D0
     * - (6284921.D0*XX ** 10)/225.D0 + (2347649.D0*XX ** 11)/2100.D0
     * - (16312.D0*XX ** 12)/675.D0 + (2039.D0*XX ** 13)/9450.D0
     * + SS ** 12*(-89079832.D0/675.D0 + (1895391191.D0*XX)/37800.D0)  + 
     *SS ** 10*(-4342880411.D0/450.D0 + (12820511189.D0*XX)/1050.D0
     * - (1055741186.D0*XX ** 2)/225.D0
     * + (3487385299.D0*XX ** 3)/6300.D0)  + 
     *SS ** 8*(-3274704127.D0/28.D0 + (7211558817317.D0*XX)/25200.D0
     * - (55615265929.D0*XX ** 2)/225.D0
     * + (402647627639.D0*XX ** 3)/4200.D0
     * - (761509408.D0*XX ** 4)/45.D0 + (692470907.D0*XX ** 5)/630.D0)
     *  + SS ** 6*(-10648440314.D0/35.D0 + (918947008807.D0*XX)/700.D0
     * - (14155172678.D0*XX ** 2)/7.D0
     * + (9382331949013.D0*XX ** 3)/6300.D0
     * - (521780426341.D0*XX ** 4)/900.D0
     * + (254938247857.D0*XX ** 5)/2100.D0
     * - (1732220216.D0*XX ** 6)/135.D0
     * + (504318104.D0*XX ** 7)/945.D0)  + 
     *SS ** 4*(-78913341669.D0/560.D0 + (2112661847355.D0*XX)/1792.D0
     * - (22549638312.D0*XX ** 2)/7.D0
     * + (711442845528.D0*XX ** 3)/175.D0
     * - (38345728971.D0*XX ** 4)/14.D0
     * + (17731459378613.D0*XX ** 5)/16800.D0
     * - (53572666604.D0*XX ** 6)/225.D0
     * + (32395208551.D0*XX ** 7)/1050.D0 - (19003480.D0*XX ** 8)/9.D0
     * + (74038129.D0*XX ** 9)/1260.D0)  + 
     *SS ** 2*(-1523246139.D0/320.D0 + (41024053941.D0*XX)/320.D0
     * - (868046758359.D0*XX ** 2)/1120.D0
     * + (1615564942095.D0*XX ** 3)/896.D0
     * - (71407187988.D0*XX ** 4)/35.D0 + (88930355691.D0*XX ** 5)/70.D0
     * - (26091997399.D0*XX ** 6)/56.D0
     * + (2619538731037.D0*XX ** 7)/25200.D0
     * - (6366624973.D0*XX ** 8)/450.D0
     * + (1194953341.D0*XX ** 9)/1050.D0 - (11100316.D0*XX ** 10)/225.D0
     * + (2779157.D0*XX ** 11)/3150.D0) 

	Z(6)=(2310525.D0*XX)/8.D0 - (4808540583.D0*XX ** 2)/100.D0
     * + (252517854951.D0*XX ** 3)/400.D0
     * - (11473454766573.D0*XX ** 4)/4900.D0
     * + (143434835467311.D0*XX ** 5)/39200.D0
     * - (429688246765.D0*XX ** 6)/147.D0
     * + (9794517932561.D0*XX ** 7)/7350.D0
     * - (72661274793.D0*XX ** 8)/196.D0
     * + (2312186142587.D0*XX ** 9)/35280.D0
     * - (11845630792.D0*XX ** 10)/1575.D0
     * + (2067628712.D0*XX ** 11)/3675.D0
     * - (127687796.D0*XX ** 12)/4725.D0
     * + (105557789.D0*XX ** 13)/132300.D0 - (48128.D0*XX ** 14)/3675.D0
     * + (3008.D0*XX ** 15)/33075.D0 + 
     *SS ** 14*(-5592287104.D0/3675.D0 + (19264982656.D0*XX)/33075.D0)
     *  + SS ** 12*(-697303053956.D0/4725.D0
     * + (98123248362941.D0*XX)/529200.D0
     * - (263501425664.D0*XX ** 2)/3675.D0
     * + (284503269632.D0*XX ** 3)/33075.D0)  + 
     *SS ** 10*(-4092665438636.D0/1575.D0
     * + (22582640792464.D0*XX)/3675.D0
     * - (8264177610763.D0*XX ** 2)/1575.D0
     * + (180539814396049.D0*XX ** 3)/88200.D0
     * - (451257055744.D0*XX ** 4)/1225.D0
     * + (271877622784.D0*XX ** 5)/11025.D0)  + 
     *SS ** 8*(-2252499518583.D0/196.D0
     * + (1597720624527617.D0*XX)/35280.D0
     * - (104821986878408.D0*XX ** 2)/1575.D0
     * + (177310534011916.D0*XX ** 3)/3675.D0
     * - (5960977068464.D0*XX ** 4)/315.D0
     * + (35848797395657.D0*XX ** 5)/8820.D0
     * - (328499926016.D0*XX ** 6)/735.D0
     * + (129412835072.D0*XX ** 7)/6615.D0)  + 
     *SS ** 6*(-7304700195005.D0/588.D0+ (303630055909391.D0*XX)/3675.D0
     * - (9736610822262.D0*XX ** 2)/49.D0
     * + (2078655342185713.D0*XX ** 3)/8820.D0
     * - (245859028495658.D0*XX ** 4)/1575.D0
     * + (224530004722216.D0*XX ** 5)/3675.D0
     * - (13559550120628.D0*XX ** 6)/945.D0
     * + (13054120650052.D0*XX ** 7)/6615.D0
     * - (107026432768.D0*XX ** 8)/735.D0
     * + (29281785088.D0*XX ** 9)/6615.D0)  + 
     *SS ** 4*(-11473454766573.D0/4900.D0
     * + (2438392202944287.D0*XX)/78400.D0
     * - (6445323701475.D0*XX ** 2)/49.D0
     * + (313424573841952.D0*XX ** 3)/1225.D0
     * - (26376042749859.D0*XX ** 4)/98.D0
     * + (3928404256255313.D0*XX ** 5)/23520.D0
     * - (100972156871008.D0*XX ** 6)/1575.D0
     * + (57062417193776.D0*XX ** 7)/3675.D0
     * - (148756282340.D0*XX ** 8)/63.D0
     * + (3832908876379.D0*XX ** 9)/17640.D0
     * - (13467610112.D0*XX ** 10)/1225.D0
     * + (2574414848.D0*XX ** 11)/11025.D0)  + 
     *SS ** 2*(-4808540583.D0/200.D0 + (252517854951.D0*XX)/200.D0
     * - (126208002432303.D0*XX ** 2)/9800.D0
     * + (1864652861075043.D0*XX ** 3)/39200.D0
     * - (8164076688535.D0*XX ** 4)/98.D0
     * + (19589035865122.D0*XX ** 5)/245.D0
     * - (17947334873871.D0*XX ** 6)/392.D0
     * + (580358721789337.D0*XX ** 7)/35280.D0
     * - (5999811996148.D0*XX ** 8)/1575.D0
     * + (2104846028816.D0*XX ** 9)/3675.D0
     * - (86891545178.D0*XX ** 10)/1575.D0
     * + (143875266407.D0*XX ** 11)/44100.D0
     * - (393903616.D0*XX ** 12)/3675.D0
     * + (49259008.D0*XX ** 13)/33075.D0)

	Z(7) = -7425.D0/64.D0 + 97088895.D0 *XX/512.D0 -
     *   1222337267253.D0 *XX**2/6400.D0 + 130359634505091.D0
     *   *XX**3/25600.D0 - 711605106982737.D0 *XX**4/22400.D0 +
     *   13836548832678459.D0 *XX**5/179200.D0 - 15517594943657.D0
     *   *XX**6/168.D0 + 166674106414085.D0 *XX**7/2688.D0 -
     *   3121309975423931.D0 *XX**8/123200.D0 + 29449145431967693.D0
     *   *XX**9/4435200.D0 - 31998166697201.D0 *XX**10/27720.D0 +
     *   175594056521981.D0 *XX**11/1293600.D0 - 2826080185913.D0
     *   *XX**12/259875.D0 + 17178303350393.D0 *XX**13/29106000.D0 -
     *   1429785592.D0 *XX**14/67375.D0 +
     *   250186931.D0 *XX**15/519750.D0 -
     *   752408.D0 *XX**16/121275.D0 + 188102.D0 *XX**17/5457375.D0 +
     *   (-1222337267253.D0/12800.D0 + 130359634505091.D0 *XX/12800.D0 -
     *   7827656176810107.D0 *XX**2/44800.D0 + 179875134824819967.D0
     *   *XX**3/179200.D0 - 294834303929483.D0 *XX**4/112.D0 +
     *   833370532070425.D0 *XX**5/224.D0 - 770963563929710957.D0
     *   *XX**6/246400.D0 + 7391735503423890943.D0 *XX**7/4435200.D0 -
     *   32414142864264613.D0 *XX**8/55440.D0 + 89377374769688329.D0
     *   *XX**9/646800.D0 - 3846295133027593.D0 *XX**10/173250.D0 +
     *   23414027466585659.D0 *XX**11/9702000.D0 -
     *   11702080177724.D0 *XX**12
     *   /67375.D0 + 2048530591028.D0 *XX**13/259875.D0 -
     *   24648509876.D0 *XX
     *   **14/121275.D0 + 160074802.D0 *XX**15/70875.D0 ) *SS**2 +
     *   (-711605106982737.D0/22400.D0 +
     *   235221330155533803.D0 *XX/358400.D0 -
     *   232763924154855.D0 *XX**2/56.D0 +
     *   166674106414085.D0 *XX**3/14.D0 -
     *   103003229188989723.D0 *XX**4/5600.D0 + 50034098088913110407.D0
     *   *XX**5/2956800.D0 - 68188093231735331.D0 *XX**6/6930.D0 +
     *   2423022385946815819.D0 *XX**7/646800.D0 -
     *   658476683317729.D0 *XX**8
     *   /693.D0 + 56705579359647293.D0 *XX**9/352800.D0 -
     *   109116947024664.D0*XX
     *   **10/6125.D0 + 107061993008968.D0 *XX**11/86625.D0 -
     *   5875665428384.D0 *XX**12/121275.D0 + 4449883963571.D0
     *   *XX**13/5457375.D0 ) *SS**4 +
     *   ( -263799114042169.D0/672.D0 +
     *   5166897298836635.D0 *XX/1344.D0 - 209127768353403377.D0
     *   *XX**2/15400.D0 + 26474781743338956007.D0 *XX**3/1108800.D0 -
     *   2656519797368324221.D0 *XX**4/110880.D0 +
     *   19068285379891482733.D0*XX
     *   **5/1293600.D0 - 300109933182659209.D0 *XX**6/51975.D0 +
     *   531101604684100381.D0 *XX**7/363825.D0 -
     *   3179538969723302.D0 *XX**8
     *   /13475.D0 + 1217739352621058.D0 *XX**9/51975.D0 -
     *   157113582940352.D0*XX
     *   **10/121275.D0 + 15084280098448.D0 *XX**11/496125.D0 ) *SS**6 + 
     *   ( -96760609238141861.D0/123200.D0 + 20349359493489675863.D0
     *   *XX/4435200.D0 - 283151777103531649.D0 *XX**2/27720.D0 +
     *   30116312228141483291.D0 *XX**3/2587200.D0 -
     *   11993884309014772.D0
     *   *XX**4/1575.D0 + 5833975135737017909.D0 *XX**5/1940400.D0 -
     *   9759068758118824.D0 *XX**6/13475.D0 + 5381881655364502.D0
     *   *XX**7/51975.D0 - 194813691539056.D0 *XX**8/24255.D0 +
     *   282795510158908.D0 *XX**9/1091475.D0 ) *SS**8 +
     *   (-22110733187765891.D0/55440.D0 + 958919142666538241.D0
     *   *XX/646800.D0 - 731635421090458831.D0 *XX**2/346500.D0 +
     *   29380756530614514013.D0 *XX**3/19404000.D0 -
     *   3656164073405268.D0
     *   *XX**4/6125.D0 + 11306553865010744.D0 *XX**5/86625.D0 -
     *   1790601425683472.D0 *XX**6/121275.D0 + 332724599656408.D0
     *   *XX**7/496125.D0 ) *SS**10 +
     *   (-15433223895270893.D0/259875.D0 + 15968418267121470617.D0
     *   *XX/116424000.D0 - 7828094703412696.D0 *XX**2/67375.D0 +
     *   11831615672987962.D0 *XX**3/259875.D0 - 91654360816864.D0
     *   *XX**4/11025.D0 + 3107297514987236.D0 *XX**5/5457375.D0 ) 
     *   *SS**12 + 
     *   ( -166135545371231.D0/67375.D0 + 801171357458921.D0*XX
     *   /259875.D0 - 145297480288496.D0 *XX**2/121275.D0 +
     *   794693945069744.D0
     *   *XX**3/5457375.D0 ) *SS**14 + 
     *   ( -2409429366728.D0/121275.D0 +
     *   41745586147262.D0 *XX/5457375.D0 ) *SS**16

      Z(8) = 675.D0 - 387045.D0 *XX/4.D0 - 29211251373.D0 *XX**2/50.D0 +
     *   852417952857.D0 *XX**3/25.D0 -
     *   1772124191636301.D0 *XX**4/4900.D0 +
     *   53149849031963007.D0 *XX**5/39200.D0 -
     *   1401315611345291.D0 *XX**6/588.D0 +
     *   26948037029579767.D0 *XX**7/11760.D0 - 285771169184667599.D0
     *   *XX**8/215600.D0 + 3810572956513200337.D0 *XX**9/7761600.D0 -
     *   240515303839073.D0 *XX**10/1980.D0 +
     *   13358379764254099.D0 *XX**11/646800.D0 -
     *   638244896394439.D0 *XX**12/259875.D0 + 5995664444857159.D0
     *   *XX**13/29106000.D0 - 820428095376.D0 *XX**14/67375.D0 +
     *   304077691466.D0
     *   *XX**15/606375.D0 - 1698615488.D0 *XX**16/121275.D0 +
     *   915950053.D0 *XX**17
     *   /3638250.D0 - 259541.D0 *XX**18/99225.D0 +
     *   259541.D0 *XX**19/21829500.D0 +
     *   ( -29211251373.D0/100.D0 + 1704835905714.D0 *XX/25.D0 -
     *   19493366107999311.D0 *XX**2/9800.D0 +
     *   690948037415519091.D0 *XX**3/39200.D0 -
     *   26624996615560529.D0 *XX**4/392.D0 +
     *   26948037029579767.D0 *XX**5/196.D0 -
     *   70585478788612896953.D0 *XX**6/431200.D0 +
     *   956453812084813284587.D0
     *   *XX**7/7761600.D0 - 243642002788980949.D0 *XX**8/3960.D0 +
     *   6799415300005336391.D0 *XX**9/323400.D0 - 868651303992831479.D0
     *   *XX**10/173250.D0 + 8172090638340307717.D0 *XX**11/9702000.D0 -
     *   6714793746604872.D0 *XX**12/67375.D0 + 4979576275447216.D0
     *   *XX**13/606375.D0 - 55645794079136.D0 *XX**14/121275.D0 +
     *   779473495103.D0 *XX
     *   **15/47250.D0 - 302365265.D0 *XX**16/882.D0 + 3779695583.D0
     *   *XX**17/1212750.D0 ) *SS**2 + 
     *   ( -1772124191636301.D0/4900.D0 +
     *   903547433543371119.D0 *XX/78400.D0 -
     *   21019734170179365.D0 *XX**2/196.D0 +
     *   107792148118319068.D0 *XX**3/245.D0 -
     *   9430448583094030767.D0 *XX**4/9800.D0 +
     *   6474163453115927372563.D0 *XX**5/5174400.D0 -
     *   512538112481064563.D0
     *   *XX**6/495.D0 + 184332282366942312101.D0 *XX**7/323400.D0 -
     *   148711060859904287.D0 *XX**8/693.D0 + 19791688332473481859.D0
     *   *XX**9/352800.D0 - 62612610954810192.D0 *XX**10/6125.D0 +
     *   260246716707324896.D0 *XX**11/202125.D0 - 13264739740884224.D0
     *   *XX**12/121275.D0 + 43336822067566613.D0 *XX**13/7276500.D0 -
     *   6182416634698.D0 *XX**14/33075.D0 +
     *   94131108503.D0 *XX**15/36750.D0 )
     *   *SS**4 + 
     *   ( -23822365392869947.D0/2352.D0 + 835389147916972777.D0
     *   *XX/5880.D0 - 19146668335372729133.D0 *XX**2/26950.D0 +
     *   3425705087905367102963.D0 *XX**3/1940400.D0 -
     *   19967821040023679533.D0
     *   *XX**4/7920.D0 + 1450626533739645372707.D0 *XX**5/646800.D0 -
     *   67777140282814660727.D0 *XX**6/51975.D0 +
     *   185367957641648784803.D0 *XX**7
     *   /363825.D0 - 1824457538038932756.D0 *XX**8/13475.D0 +
     *   2960085642141816376.D0 *XX
     *   **9/121275.D0 - 354695278834959872.D0 *XX**10/121275.D0 +
     *   36725944316488636.D0 *XX**11/165375.D0 -
     *   60678510586879.D0 *XX**12/6300.D0 +
     *   37643589679067.D0 *XX**13/207900.D0 ) *SS**6 +
     *   (-8858906244724695569.D0/215600.D0 + 2633105912950621432867.D0
     *   *XX/7761600.D0 - 2128319923671956977.D0 *XX**2/1980.D0 +
     *   2291109071746984773589.D0 *XX**3/1293600.D0 -
     *   2708711340297999116.D0 *XX**4
     *   /1575.D0 + 2036205589111274339467.D0 *XX**5/1940400.D0 -
     *   5599870525109370672.D0 *XX**6/13475.D0 +
     *   13082299246929507944.D0
     *   *XX**7/121275.D0 - 439806001162527616.D0 *XX**8/24255.D0 +
     *   688526869778667481.D0 *XX**9/363825.D0 -
     *   174827688100514.D0 *XX**10/1575.D0 +
     *   1930312412765777.D0 *XX**11/693000.D0 ) *SS**8 +
     *   (-166196074952799443.D0/3960.D0 +
     *   72950111892591634639.D0 *XX/323400.D0 -
     *   165233306492867129393.D0 *XX**2/346500.D0 +
     *   10254630722281438181219.D0
     *   *XX**3/19404000.D0 - 2097950730451927704.D0 *XX**4/6125.D0 +
     *   27484015923345307168.D0 *XX**5/202125.D0 -
     *   4042412247744344192.D0
     *   *XX**6/121275.D0 + 810090043406478106.D0 *XX**7/165375.D0 -
     *   265846057283887.D0 *XX**8/675.D0 +
     *   382998603539749.D0 *XX**9/28875.D0 )
     *    *SS**10 + 
     *    ( -3485455379210031379.D0/259875.D0 +
     *   5573383802341424434471.D0
     *   *XX/116424000.D0 - 4491854487748839888.D0 *XX**2/67375.D0 +
     *   28760338246086271064.D0 *XX**3/606375.D0 -
     *   206916349674997504.D0
     *   *XX**4/11025.D0 + 7565388256210266527.D0 *XX**5/1819125.D0 -
     *   454160444366756.D0 *XX**6/945.D0 +
     *   833546612742332.D0 *XX**7/37125.D0 )
     *    *SS**12 + 
     *    ( -95330565523821618.D0/67375.D0 + 1947490509364697212.D0
     *   *XX/606375.D0 - 328019572340291456.D0 *XX**2/121275.D0 +
     *   1934854390451486708.D0 *XX**3/1819125.D0 - 721874509066993.D0
     *   *XX**4/3675.D0 +
     *   8305150942611073.D0 *XX**5/606375.D0 ) *SS**14 + 
     *    ( -5439461089683008.D0/121275.D0 +
     *   203277327429270793.D0 *XX/3638250.D0 -
     *   721055087779039.D0 *XX**2/33075.D0 +
     *   28864054745507.D0 *XX**3/10780.D0 )
     *   *SS**16 + 
     *    ( -57600084923321.D0/198450.D0 + 1225582843976861.D0
     *   *XX/10914750.D0 ) *SS**18

      Z(9) = -5716575.D0/2048.D0 - 712307115.D0 *XX/16384.D0 -
     *   53478412581339.D0 *XX**2/40960.D0 +
     *   31100390895173733.D0 *XX**3/163840.D0 -
     *   502781754480920721.D0 *XX**4/143360.D0 + 3325681861227511821.D0
     *   *XX**5/163840.D0 - 11627885912123209.D0 *XX**6/224.D0 +
     *   12594490460272652453.D0 *XX**7/179200.D0 -
     *   7915277597603095461.D0*XX**8/140800.D0 +
     *   337682546305005247727.D0 *XX**9/11827200.D0 -
     *    92978013607990233689.D0 *XX**10/9609600.D0 +
     *   40634464745833234669.D0
     *   *XX**11/17937920.D0 -
     *   6763724664910402927.D0 *XX**12/18018000.D0 +
     *   90278540532324113647.D0 *XX**13/2018016000.D0 -
     *   3387831607554819.D0 *XX**14
     *   /875875.D0 + 625697507310091.D0 *XX**15/2574000.D0 -
     *   3481283602417.D0
     *   *XX**16/315315.D0 + 11251974691379.D0 *XX**17/31531500.D0 -
     *   1468025159.D0 *XX
     *   **18/184275.D0 + 4713758111.D0 *XX**19/40540500.D0 - 43636.D0
     *   *XX**20/43875.D0 + 10909.D0 *XX**21/2895750.D0 +
     *   (-53478412581339.D0/81920.D0 +
     *   31100390895173733.D0 *XX/81920.D0 -
     *   5530599299290127931.D0 *XX**2/286720.D0 +
     *   43233864195957653673.D0
     *   *XX**3/163840.D0 - 662789496991022913.D0 *XX**4/448.D0 +
     *   37783471380817957359.D0 *XX**5/8960.D0 -
     *   1955073566607964578867.D0 *XX**6
     *   /281600.D0 + 84758319122556317179477.D0 *XX**7/11827200.D0 -
     *   94186727784894106726957.D0 *XX**8/19219200.D0 +
     *   20682942555629116446521.D0 *XX
     *   **9/8968960.D0 -
     *   9205429268943058383647.D0 *XX**10/12012000.D0 +
     *   123049650745557766900861.D0 *XX**11/672672000.D0 -
     *   55455415584064832211.D0 *XX
     *   **12/1751750.D0 + 1280802797463756277.D0 *XX**13/321750.D0 -
     *   228090220346759423.D0 *XX**14/630630.D0 + 9575430462363529.D0
     *   *XX**15/409500.D0 - 1710249310235.D0 *XX**16/1638.D0 +
     *   68646459370493.D0
     *   *XX**17/2252250.D0 - 4575474598.D0 *XX**18/8775.D0 +
     *   2287759117.D0 *XX**19
     *   /579150.D0 ) *SS**2 + 
     *    ( -502781754480920721.D0/143360.D0 +
     *   56536591640867700957.D0 *XX/327680.D0 -
     *   523254866045544405.D0 *XX**2/224.D0 +
     *   37783471380817957359.D0 *XX**3/2800.D0 -
     *   261204160720902150213.D0
     *   *XX**4/6400.D0 +
     *   573722646172203915888173.D0 *XX**5/7884800.D0 -
     *   198136146998627187991259.D0 *XX**6/2402400.D0 +
     *   560714979027752805197531.D0
     *   *XX**7/8968960.D0 - 1575947846924123881991.D0 *XX**8/48048.D0 +
     *   298009462297201899148747.D0 *XX**9/24460800.D0 -
     *   258549144793761121623.D0
     *   *XX**10/79625.D0 + 33469185363524077681.D0 *XX**11/53625.D0 -
     *   27185858881247510716.D0 *XX**12/315315.D0 +
     *   532370540852028835459.D0
     *   *XX**13/63063000.D0 - 34969207805921902.D0 *XX**14/61425.D0 +
     *   1709599932971813.D0 *XX**15/68250.D0 -
     *   48374738692.D0 *XX**16/75.D0 +
     *   945595141793.D0 *XX**17/128700.D0 ) *SS**4 +
     *   ( -197674060506094553.D0/896.D0 +
     *   390429204268452226043.D0 *XX/89600.D0 -
     *   530323599039407395887.D0 *XX**2/17600.D0 +
     *   303576609128199717706573.D0
     *   *XX**3/2956800.D0 -
     *   7719127667748957191094469.D0 *XX**4/38438400.D0 +
     *   4412618430144268452410717.D0 *XX**5/17937920.D0 -
     *   718260213340830418026911.D0
     *   *XX**6/3603600.D0 +
     *   2791141637637864621624299.D0 *XX**7/25225200.D0 -
     *   30135266936109640673031.D0 *XX**8/700700.D0 +
     *   1522734369942768752969.D0 *XX
     *   **9/128700.D0 - 726941951716662669448.D0 *XX**10/315315.D0 +
     *   225579656124614340274.D0 *XX**11/716625.D0 -
     *   343211978655342421.D0
     *   *XX**12/11700.D0 + 683679172758284657.D0 *XX**13/386100.D0 -
     *   181910249319532.D0 *XX**14/2925.D0 +
     *   92841739468754.D0 *XX**15/96525.D0 )
     *   *SS**6 + 
     *    ( -245373605525695959291.D0/140800.D0 +
     *   233338639496758626179357.D0 *XX/11827200.D0 -
     *   822762442417105577913961.D0
     *   *XX**2/9609600.D0 +
     *   6969257683022603911314859.D0 *XX**3/35875840.D0 -
     *   7176311869469937505547.D0 *XX**4/27300.D0 +
     *   30659765985804189207998611.D0
     *   *XX**5/134534400.D0 -
     *   23123803865450972160993.D0 *XX**6/175175.D0 +
     *   6729827819022799280611.D0 *XX**7/128700.D0 -
     *   901374932059785521594.D0
     *   *XX**8/63063.D0 + 8458198007314238292383.D0 *XX**9/3153150.D0 -
     *   988866670858783286.D0 *XX**10/2925.D0 + 35058144156178258067.D0
     *   *XX**11/1287000.D0 - 3695623019335772.D0 *XX**12/2925.D0 +
     *   2468935855239448.D0 *XX**13/96525.D0 ) *SS**8 +
     *   (-64247807403121251479099.D0/19219200.D0 +
     *   221904811976995294527409.D0 *XX
     *   /8968960.D0 - 1751040387324659482562249.D0 *XX**2/24024000.D0 +
     *   154407089292595756859123627.D0 *XX**3/1345344000.D0 -
     *   17326329597380282379777.D0 *XX**4/159250.D0 +
     *   3534598380762591493823.D0 *XX
     *   **5/53625.D0 - 8284855266951406033078.D0 *XX**6/315315.D0 +
     *   4975769495450756418379.D0 *XX**7/716625.D0 -
     *   10525816358863953091.D0 *XX**8
     *   /8775.D0 + 6955982961987374479.D0 *XX**9/53625.D0 -
     *   38701255504228037.D0
     *   *XX**10/4875.D0 + 134581881609362971.D0 *XX**11/643500.D0 )
     *   *SS**10 +
     *   ( -36936700395075710384347.D0/18018000.D0 +
     *   83920132644091993998728143.D0
     *   *XX/8072064000.D0 -
     *   18548422093173532237647.D0 *XX**2/875875.D0 +
     *   14794962319673629555141.D0 *XX**3/643500.D0 -
     *   424071545493614949236.D0
     *   *XX**4/28665.D0 +
     *   92936898590183087947561.D0 *XX**5/15765750.D0 -
     *   2568838675018658444.D0 *XX**6/1755.D0 +
     *   105971541093604493804.D0
     *   *XX**7/482625.D0 - 53059260588421304.D0 *XX**8/2925.D0 +
     *   487518590631633239.D0 *XX**9/772200.D0 ) *SS**12 +
     *   (-3149223239603125543011.D0/7007000.D0 +
     *   2003665496381538618481.D0
     *   *XX/1287000.D0 - 672270544173969978154.D0 *XX**2/315315.D0 +
     *   11884344873149797348022.D0 *XX**3/7882875.D0 -
     *   4083092617163066107.D0
     *   *XX**4/6825.D0 + 150837334443545493283.D0 *XX**5/1126125.D0 -
     *   46036944777647392.D0 *XX**6/2925.D0 +
     *   72566575071858044.D0 *XX**7/96525.D0 )
     *    *SS**14 + 
     *    ( -11148083148467537347.D0/315315.D0 +
     *   2497157280655037161199.D0 *XX/31531500.D0 -
     *   313727521303017097.D0
     *   *XX**2/4725.D0 + 524226122936189897.D0 *XX**3/20020.D0 -
     *   4758051520483876.D0
     *   *XX**4/975.D0 + 22214852115500569.D0 *XX**5/64350.D0 )
     *   *SS**16 +
     *   ( -325799676459487379.D0/368550.D0 + 22258915051950851831.D0
     *   *XX/20270250.D0 - 3772319621333914.D0 *XX**2/8775.D0 +
     *   30842030895512227.D0
     *   *XX**3/579150.D0 ) *SS**18 + 
     *    ( -206054276510356.D0/43875.D0 +
     *   10564097574141743.D0 *XX/5791500.D0 ) *SS**20

	Z(10) = 91125.D0/8.D0 + 21438675.D0 *XX/128.D0 - 607296231621.D0
     *   *XX**2/320.D0 + 1115737995860937.D0 *XX**3/1280.D0 -
     *   65356894821275343.D0 *XX
     *   **4/2240.D0 + 4712475624988559301.D0 *XX**5/17920.D0 -
     *   27311138259540955.D0 *XX**6/28.D0 +
     *   18410040997415653.D0 *XX**7/10.D0 -
     *   2470816192013890991.D0 *XX**8/1232.D0 +
     *   101184393090441496963.D0
     *   *XX**9/73920.D0 - 1144746381575325625.D0 *XX**10/1848.D0 +
     *   10013612053934469175.D0 *XX**11/51744.D0 -
     *   425435768722970419.D0
     *   *XX**12/9900.D0 + 10751805350563854881.D0 *XX**13/1552320.D0 -
     *   99688344381420658.D0 *XX**14/121275.D0 + 23433617409470623.D0
     *   *XX**15/323400.D0 - 577302804852667.D0 *XX**16/121275.D0 +
     *   5075875341420931.D0 *XX**17/21829500.D0 -
     *   830418475216.D0 *XX**18/99225.D0 +
     *   108140875343.D0 *XX**19/496125.D0 -
     *   5921877194.D0 *XX**20/1488375.D0 +
     *   3149587499.D0 *XX**21/65488500.D0 -
     *   5641586.D0 *XX**22/16372125.D0 +
     *   2820793.D0 *XX**23/2554051500.D0 + 
     *    ( -607296231621.D0/640.D0 +
     *   1115737995860937.D0 *XX/640.D0 -
     *   718925843034028773.D0 *XX**2/4480.D0 +
     *   61262183124851270913.D0 *XX**3/17920.D0 -
     *   1556734880793834435.D0*XX**4/56.D0 + 
     *   110460245984493918.D0 *XX**5 - 610291599427431074777.D0
     *   *XX**6/2464.D0 + 25397282665700815737713.D0 *XX**7/73920.D0 -
     *   1159628084535804858125.D0 *XX**8/3696.D0 +
     *   5096928535452644810075.D0 *XX**9
     *   /25872.D0 - 579018081231962740259.D0 *XX**10/6600.D0 +
     *   14654710692818534202803.D0 *XX**11/517440.D0 -
     *   815899254589737375401.D0 *XX
     *   **12/121275.D0 + 47968614837186365281.D0 *XX**13/40425.D0 -
     *   37824302471141889173.D0 *XX**14/242550.D0 +
     *   4319569915549212281.D0
     *   *XX**15/283500.D0 - 483718761813320.D0 *XX**16/441.D0 +
     *   3149711135240218.D0
     *   *XX**17/55125.D0 - 620941394115467.D0 *XX**18/297675.D0 +
     *   94358491882541.D0
     *   *XX**19/1871100.D0 - 1690171212119.D0 *XX**20/2338875.D0 +
     *   2957807378389.D0*XX**21/638512875.D0 ) *SS**2 + 
     *    ( -65356894821275343.D0/2240.D0 +
     *   80112085624805508117.D0 *XX/35840.D0 -
     *   1229001221679342975.D0 *XX**2/28.D0 +
     *   1767363935751902688.D0 *XX**3/5.D0 -
     *   81536934336458402703.D0 *XX**4/56.D0 +
     *   171912283860660103340137.D0 *XX**5/49280.D0 -
     *   2439454539137018906875.D0
     *   *XX**6/462.D0 + 138177832732241740145825.D0 *XX**7/25872.D0 -
     *   495632670562260538135.D0 *XX**8/132.D0 +
     *   35491709462211284962181.D0
     *   *XX**9/18816.D0 - 2535971792718960118862.D0 *XX**10/3675.D0 +
     *   2506975257699986189786.D0 *XX**11/13475.D0 -
     *   4508243043909594797716.D0*XX**12/121275.D0 +
     *   240157534559701388708051.D0 *XX**13/43659000.D0 -
     *   19781048061523794848.D0 *XX**14/33075.D0 +
     *   862858856024555918.D0
     *   *XX**15/18375.D0 - 12192096970182662.D0 *XX**16/4725.D0 +
     *   273007116849057223.D0 *XX**17/2910600.D0 - 6612931784476618.D0
     *   *XX**18/3274425.D0 + 496561443759988.D0 *XX**19/25540515.D0 )
     *   *SS**4 +
     *   ( -464289350412196235.D0/112.D0 +
     *   570711270919885243.D0 *XX/5.D0 -
     *   165544684864930696397.D0 *XX**2/154.D0 +
     *   90964769388306905769737.D0
     *   *XX**3/18480.D0 - 95037989344765108713125.D0 *XX**4/7392.D0 +
     *   1087408173772905811120775.D0 *XX**5/51744.D0 -
     *   45178300587998397704867.D0
     *   *XX**6/1980.D0 + 332413566023382701355877.D0 *XX**7/19404.D0 -
     *   443371338418024792294021.D0 *XX**8/48510.D0 +
     *   57029433911114866899557.D0
     *   *XX**9/16170.D0 -
     *   120549106484669676255448.D0 *XX**10/121275.D0 +
     *   101761179300058873213586.D0 *XX**11/496125.D0 -
     *   48536219942065003676.D0 *XX
     *   **12/1575.D0 + 15684653827129531841.D0 *XX**13/4725.D0 -
     *   24687188486575088678.D0 *XX**14/99225.D0 +
     *   13402382528930355647.D0
     *   *XX**15/1091475.D0 - 111525585134660881.D0 *XX**16/311850.D0 +
     *   395299246809670547.D0 *XX**17/85135050.D0 ) *SS**6 +
     *   ( -76595301952430620721.D0/1232.D0 + 69918415625495074401433.D0
     *   *XX/73920.D0 - 10129860730560056455625.D0 *XX**2/1848.D0 +
     *   1717444616982354742673425.D0 *XX**3/103488.D0 -
     *   451387350615071614559.D0
     *   *XX**4/15.D0 + 3651452870521042447701053.D0 *XX**5/103488.D0 -
     *   680427479927572627949926.D0 *XX**6/24255.D0 +
     *   252045450877008778684183.D0
     *   *XX**7/16170.D0 - 149475405031843508092094.D0 *XX**8/24255.D0 +
     *   3815575476816178205373487.D0 *XX**9/2182950.D0 -
     *   559372670129063474464.D0
     *   *XX**10/1575.D0 + 804287854334958455171.D0 *XX**11/15750.D0 -
     *   10235428936146118462.D0 *XX**12/2025.D0 +
     *   50915557963248163852.D0 *XX**13
     *   /155925.D0 - 2697658069891053532.D0 *XX**14/218295.D0 +
     *   8820521016380855371.D0
     *   *XX**15/42567525.D0 ) *SS**8 + 
     *    ( -791019749668550006875.D0/3696.D0 +
     *   54684335426536136164675.D0 *XX/25872.D0 -
     *   110139789857383642863653.D0
     *   *XX**2/13200.D0 +
     *   18389253515088734121024421.D0 *XX**3/1034880.D0 -
     *   84972400828007531176369.D0 *XX**4/3675.D0 +
     *   264755493455644502540038.D0 *XX
     *   **5/13475.D0 - 1373881225904366710466578.D0 *XX**6/121275.D0 +
     *   2244618067431649501723931.D0 *XX**7/496125.D0 -
     *   850591920089199136112.D0
     *   *XX**8/675.D0 + 638323875486879917308.D0 *XX**9/2625.D0 -
     *   214374652598823086429.D0 *XX**10/6750.D0 +
     *   792974555569633788269.D0 *XX**11
     *   /297000.D0 - 101340488710467827698.D0 *XX**12/779625.D0 +
     *   6495042153943108697.D0 *XX**13/2338875.D0 ) *SS**10 +
     *   ( -2323304732996141458159.D0/9900.D0 +
     *   9994544947918292017876289.D0
     *   *XX/6209280.D0 - 545794981436755061018554.D0 *XX**2/121275.D0 +
     *   554100795570125499148273.D0 *XX**3/80850.D0 -
     *   70323972600708613626236.D0
     *   *XX**4/11025.D0 +
     *   41924739861304629914696729.D0 *XX**5/10914750.D0 -
     *   1453116169370012885056.D0 *XX**6/945.D0 +
     *   1389228836155342241744.D0
     *   *XX**7/3375.D0 - 1028673654810167033188.D0 *XX**8/14175.D0 +
     *   2872525135739803987921.D0 *XX**9/356400.D0 -
     *   56739639943102900592.D0
     *   *XX**10/111375.D0 +
     *   847071322368712273799.D0 *XX**11/60810750.D0 )*SS**12 + 
     *    ( -46333597299146409818201.D0/485100.D0 +
     *   75041262127791090797293.D0 *XX/161700.D0 -
     *   111482922707591939968654.D0
     *   *XX**2/121275.D0 +
     *   5361143687674816604521558.D0 *XX**3/5457375.D0 -
     *   2309684901871807898768.D0 *XX**4/3675.D0 +
     *   13841763617920914109516.D0
     *   *XX**5/55125.D0 - 6247711370432383620368.D0 *XX**6/99225.D0 +
     *   10475514606818641584242.D0 *XX**7/1091475.D0 -
     *   3536318451369159685891.D0
     *   *XX**8/4365900.D0 +
     *   985287925995036150079.D0 *XX**9/34054020.D0 )
     *   *SS**14 + 
     *    ( -1848691576254451860097.D0/121275.D0 +
     *   1126491963605120582390911.D0 *XX/21829500.D0 -
     *   2307063109644366794864.D0
     *   *XX**2/33075.D0 + 12026555134361815961.D0 *XX**3/245.D0 -
     *   645719057361593389154.D0 *XX**4/33075.D0 +
     *   916250285021729060537.D0
     *   *XX**5/207900.D0 - 81840197248656023962.D0 *XX**6/155925.D0 +
     *   8692250878034489195497.D0 *XX**7/340540200.D0 ) *SS**16 +
     *   ( -92147627338910490248.D0/99225.D0 + 1021307628105163683406.D0
     *   *XX/496125.D0 - 73134939480036212483.D0 *XX**2/42525.D0 +
     *   8904544408495470286997.D0 *XX**3/13097700.D0 -
     *   83597711386408585033.D0
     *   *XX**4/654885.D0 +
     *   1166977249498219845157.D0 *XX**5/127702575.D0 )
     *   *SS**18 + 
     *    ( -27963794133120522074.D0/1488375.D0 +
     *   3050009135368325181673.D0 *XX/130977000.D0 -
     *   21381758576937593158.D0 *XX**2
     *   /2338875.D0 + 728285122468077990602.D0 *XX**3/638512875.D0 )
     *   *SS**20 +
     *   ( -2731609908191035811.D0/32744250.D0 + 41515492371927341743.D0
     *   *XX/1277025750.D0 ) *SS**22

      
	ZA=Z(0)+Z(1)*TH**1+Z(2)*TH**2+Z(3)*TH**3+Z(4)*TH**4+A*Z(5)*TH**5
     *+B*Z(6)*TH**6+C*Z(7)*TH**7+D*Z(8)*TH**8
	
	ZB=Z(0)+Z(1)*TH**1+Z(2)*TH**2+Z(3)*TH**3+Z(4)*TH**4+Z(5)*TH**5
     *+A*Z(6)*TH**6+B*Z(7)*TH**7+C*Z(8)*TH**8+D*Z(9)*TH**9

	ZC=Z(0)+Z(1)*TH**1+Z(2)*TH**2+Z(3)*TH**3+Z(4)*TH**4+Z(5)*TH**5
     *+Z(6)*TH**6+A*Z(7)*TH**7+B*Z(8)*TH**8+C*Z(9)*TH**9+D*Z(10)*TH**10
      SUMZ = (ZA+2.D0*ZB+ZC)/4.D0

      RETURN
      END