Run Code  | Code Wall  | Users  | Misc  | Feedback  | About  | Login  | Theme  | Privacy 

1568-1644 Solar eclipses, make 1617(8) and 1620(12) Lunar eclipses the...

Language: Editor: Layout:
+ ] Compiler args + ] Show input
Compilation time: 1.24 sec, absolute running time: 0.15 sec, cpu time: 0.01 sec, memory peak: 3 Mb, absolute service time: 1.42 sec 
edit mode |  history  | discussion
	Solar eclipse time is in 24 hours
	Beijing local noon time is 12 as reference
	month means total months from winter solitice
	ReqMDay=ReqMoonDay
	t	t	t	t	t	t	t	t	t	t	t	t	t	t	t	t
 N year(month)
 1 1568(10)				[2]MinSun=6.730724
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242213)104824.52(59.58)11.11	  164.78(48.46) 8.82(17.85)	   94.84(43.77) 1.03(13.82)	   (184.77)  94.84  87.78(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	11.62/23997.74(2.40)-1(182.37)	  1.03/ 0.98(12.54)0.04(2)	(12.54)10.21  1.34(1.34)1.1984	   -2.401/-1.343/-0.256(-6.147h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  18.3859h(12.2387h) 0.24h 	   7.31m(12.3605h)  2456/7544	  2456.06	150.22		   138.35/94.59(1)88.03

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	   5.06/24009.36(2.40)		  92.19(0)  92.19/-0.08/(-0.08)	 92.19/4.46/0.27(0.27)			 1/-0.08/ 0.27

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  188.24/182.37/-5.87		    5.87(6)0.0221		  0.77(2) 9.41(1.20)(1.12)	   107.29(12.1030h)(12.6180h)

----------------------------------------------
 1 1572(2)				[2]MinSun=6.824532
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242209)106285.48(20.54)25.08	  150.71(55.46)22.99(22.60)	   209.77(54.52)26.94( 0.03)	   ( 0.34)  27.15  27.15(2)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	370.99/12005.28(1.21)1(1.55)	 13.17/13.12(160.63)0.05(1)	( 6.37)10.68  0.80(0.80)1.2031	   1.206/0.801/0.137(3.284h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  12.4894h(15.7731h) 3.77h 	  80.84m(17.1204h)  2119/7881	  2118.62	2133.49		   264.49/209.97(2)27.35

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  370.99/12005.28(1.21)		  211.18(2)  28.56/4.02/(-0.03)	 28.56/2.35/2.01(2.01)			 1/ 0.03/ 2.01

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  359.68/365.34/ 5.66		    5.66(6)0.0562		 27.08(1) 4.70(1.20)(1.12)	   169.11(16.7145h)(17.5262h)

----------------------------------------------
 2 1572(8)				[2]MinSun=6.824532
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242209)106285.48(20.54)25.08	  150.71(55.46)22.99(22.60)	   21.72(51.70)11.25(13.94)	   (186.30)  21.72  21.72(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	388.29/9228.65(0.95)-1(185.35)	 11.25/11.24(137.21)0.01(3)	(29.79) 8.51  3.06(3.06)1.0112	   -0.951/-3.059/-0.325(-7.804h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  16.8948h(9.0908h)-2.91h 	  -68.87m(7.9429h)  2099/7901	  2900.54	1690.45		   73.05/21.34(0)21.34

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  388.29/9228.65(0.94)		  20.41(0)  20.41/4.24/(1.77)	 20.41/1.77/1.20(1.20)			 -1/-1.77/-1.20

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  185.09/185.35/ 0.27		   -0.27(8)0.9668		 10.92(3)33.75(1.02)(0.93)	   614.28(6.4686h)(9.4172h)

----------------------------------------------
 2 1575(0)				[2]MinSun=6.894888
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242206)107381.21(36.27)28.18	  147.55( 8.09)13.45(26.74)	   147.55( 8.09)13.45(26.74)	   (357.53)  147.55  35.07(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	326.94/14819.79(1.48)1(356.05)	 13.45/13.37(164.07)0.08(3)	( 3.00)10.90  0.44(0.44)0.9873	   -1.484/-0.439/-0.160(-3.833h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  2.1876h(22.3545h)10.35h 	  44.37m(23.0940h)  2150/7850	  2150.05	4622.51		   155.51/147.42(1)35.20

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  326.94/14819.79(1.49)		  145.93(1)  36.69/3.74/(-4.30)	 145.93/2.86/5.29(0.43)			 1/ 4.30/-0.43

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  361.51/356.05/-5.46		   -5.46(8)0.3174		 13.29(3) 4.87(0.99)(0.91)	   462.76(21.9834h)(24.2046h)

----------------------------------------------
 3 1575(6)				[2]MinSun=6.894888
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242206)107381.21(36.27)28.18	  147.55( 8.09)13.45(26.74)	   324.73( 5.27)25.31(13.44)	   (179.70)  142.11  40.51(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	294.76/15773.60(1.59)-1(181.29)	 11.53/11.48(140.65)0.05(1)	(26.35) 8.90  2.76(2.76)1.1853	   1.592/2.761/0.301(7.229h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  6.5929h(13.8215h) 1.82h 	  48.28m(14.6262h)  2173/7827	  2827.25	1094.25		   330.34/325.07(3)40.18

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  294.76/15773.60(1.58)		  326.65(3)  38.59/3.66/(2.25)	 144.03/2.97/1.30(1.30)			 1/-2.25/-1.30

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  184.50/181.29/-3.21		    3.21(6)0.4650		 25.61(1)22.67(1.19)(1.11)	   438.02(13.5750h)(15.6774h)

----------------------------------------------
 4 1582(7)				[2]MinSun=7.059052
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242199)109937.90(12.96)15.71	  159.85(57.25)20.05(10.74)	   1.32(23.97) 6.33(26.97)	   (360.59)   1.32   1.32(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	480.41/484.85(0.06)1(360.53)	  6.33/ 6.31(77.27)0.02(2)	(77.27) 0.84  5.41(5.41)1.1047	   -0.064/-5.408/-0.406(-9.747h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  23.2054h(13.4579h) 1.46h 	  40.02m(14.1250h)  2068/7932	  2931.64	885.41		   24.91/ 0.95(0) 0.95

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  484.85/ 0.00(0.05)		   0.90(0)  0.90/4.46/(3.11)	 0.90/0.09/0.03(0.03)			 1/ 3.11/-0.03

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  360.72/360.53/-0.19		   -0.19(8)0.9759		  5.93(2)72.31(1.11)(1.03)	   555.62(12.7915h)(15.4585h)

----------------------------------------------
 5 1583(13)				[2]MinSun=7.082504
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242198)110303.14(18.20)26.58	  148.95(51.62)16.21(11.35)	   167.61(15.52)14.35(14.28)	   (190.89)  167.61  15.01(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	434.82/7135.84(0.71)-1(190.18)	  0.57/ 0.49( 7.00)0.08(0)	( 7.00)10.63  0.76(0.76)0.9900	   -0.714/0.764/0.004(0.099h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  12.4215h(12.5206h) 0.52h 	  15.56m(12.7799h)  2085/7915	  2084.53	324.97		   183.14/167.63(1)15.00

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  440.02/6695.81(0.71)		  166.91(1)  15.71/4.33/(3.65)	 166.91/1.40/0.18(0.18)			 1/ 3.65/ 0.18

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  191.89/190.18/-1.71		    1.71(6)0.7154		 14.36(0) 7.05(0.99)(0.91)	   606.00(11.3255h)(14.2343h)

----------------------------------------------
 6 1587(10)				[2]MinSun=7.176312
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242194)111764.11(39.17)11.02	  164.42(28.15) 4.81(18.42)	   94.49(23.46)24.57(14.39)	   (192.42)  94.49  88.13(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	 5.06/24009.36(2.40)-1(190.02)	 10.79/10.74(131.68)0.05(1)	(35.32) 7.85  3.51(3.51)1.1747	   -2.401/3.508/0.077(1.855h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  10.9304h(12.7857h) 0.79h 	  22.95m(13.1681h)  2463/7537	  2462.69	486.72		   118.04/94.58(1)88.04

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	   5.06/24009.36(2.40)		  92.18(0)  92.18/-0.08/(-0.07)	 92.18/4.46/0.87(0.87)			 1/-0.07/ 0.87

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  188.85/190.02/ 1.17		   -1.17(8)0.8537		 24.65(1)34.38(1.18)(1.09)	   519.07(11.9223h)(14.4139h)

----------------------------------------------
 7 1589(2)				[2]MinSun=7.223216
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242192)112494.60(49.66) 3.24	  172.16(46.41)26.66(21.96)	   231.22(45.48) 3.06(26.59)	   (355.52)  48.60  48.60(2)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	252.82/18628.68(1.88)1(357.40)	  3.06/ 3.03(37.30)0.02(2)	(37.30) 7.59  3.58(3.58)1.1722	   1.878/-3.585/-0.119(-2.866h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  11.4124h(8.5468h)-3.45h 	  -76.86m(7.2658h)  2213/7787	  2213.27	1972.56		   276.52/231.05(2)48.42

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  252.82/18628.68(1.87)		  232.92(2)  50.30/3.11/(0.34)	 50.30/3.56/2.81(2.81)			 -1/-0.34/-2.81

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  354.49/357.40/ 2.91		    2.91(6)0.5157		  2.94(2)35.85(1.17)(1.09)	   459.85(6.1622h)(8.3695h)

----------------------------------------------
 7 1590(2)				[2]MinSun=7.246668
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242191)112859.84(54.90)14.12	  161.26(40.78)22.82(22.56)	   220.32(39.84)26.77(27.20)	   (363.65)  37.70  37.70(2)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	315.74/15468.08(1.57)1(365.22)	 13.00/12.96(158.54)0.04(1)	( 8.46)10.53  1.02(1.02)1.2016	   1.569/1.023/0.177(4.245h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  20.2232h(0.4677h)-11.53h 	  -14.05m(0.2336h)  2159/7841	  2158.87	4902.66		   260.33/220.49(2)37.87

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  315.74/15468.08(1.57)		  222.06(2)  39.44/3.63/(-4.61)	 39.44/3.02/5.92(0.12)			 -1/ 4.61/-0.12

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  362.13/365.22/ 3.09		    3.09(6)0.4853		 26.95(1) 6.30(1.20)(1.12)	   438.98(-0.8199h)(1.2872h)

----------------------------------------------
 8 1590(8)				[2]MinSun=7.246668
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242191)112859.84(54.90)14.12	  161.26(40.78)22.82(22.56)	   32.26(37.03)11.08(13.90)	   (185.82)  32.26  32.26(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	334.85/13234.41(1.33)-1(184.49)	 11.08/11.07(135.12)0.00(3)	(31.88) 8.27  3.23(3.23)1.0136	   -1.332/-3.234/-0.369(-8.865h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  0.6286h(15.7636h) 3.76h 	  80.73m(17.1090h)  2138/7862	  2862.47	2128.76		   68.97/31.95(0)31.95

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  339.79/12894.61(1.32)		  30.63(0)  30.63/3.96/(1.01)	 30.63/2.49/2.12(2.12)			 1/-1.01/ 2.12

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  189.16/184.49/-4.67		    4.67(6)0.2223		 10.71(3)36.39(1.02)(0.94)	   385.00(16.1850h)(18.0330h)

----------------------------------------------
 8 1594(0)				[2]MinSun=7.340476
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242187)114320.80(15.86)28.08	  147.20(47.78) 9.44( 0.10)	   147.20(47.78) 9.44( 0.10)	   ( 1.39)  147.20  35.42(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	326.94/14819.79(1.50)1(-0.11)	  9.44/ 9.43(115.14)0.01(3)	(51.86) 5.49  4.61(4.61)1.0414	   -1.496/-4.608/-0.481(-11.535h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  18.7321h(7.1975h)-4.80h 	  -90.02m(5.6972h)  2150/7850	  2150.05	2626.16		   194.43/146.65(1)35.97

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  326.94/14819.79(1.51)		  145.14(1)  37.48/3.71/(-0.82)	 145.14/2.91/3.06(2.76)			 -1/ 0.82/ 2.76

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  361.22/363.69/ 2.46		    2.46(6)0.5896		  8.96(3)57.73(1.05)(0.97)	   540.18(4.4008h)(6.9937h)

----------------------------------------------
 9 1594(6)				[2]MinSun=7.340476
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242187)114320.80(15.86)28.08	  147.20(47.78) 9.44( 0.10)	   324.38(44.96)21.29(14.01)	   (187.35)  141.76  40.86(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	294.76/15773.60(1.60)-1(188.95)	  7.52/ 7.46(91.71)0.05(1)	(75.29) 1.25  5.40(5.40)1.1088	   1.603/5.397/0.518(12.424h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  23.1375h(11.5617h)-0.44h 	  -13.20m(11.3417h)  2173/7827	  2827.25	274.28		   369.85/324.89(3)40.35

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  294.76/15773.60(1.59)		  326.48(3)  38.77/3.66/(3.30)	 143.86/2.98/0.33(0.33)			 -1/-3.30/ 0.33

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  185.08/188.95/ 3.88		   -3.88(8)0.5152		 21.81(1)69.00(1.12)(1.04)	   483.09(10.1823h)(12.5011h)

----------------------------------------------
10 1596(10)				[2]MinSun=7.387380
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242185)115051.29(26.35)20.30	  154.93( 6.05) 3.73( 3.64)	   84.99( 1.35)23.50(26.82)	   (358.58)  84.99  84.99(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	55.73/23718.98(2.38)1(356.20)	  9.72/ 9.68(118.58)0.04(1)	(48.42) 6.03  4.41(4.41)1.1566	   -2.377/4.412/0.144(3.462h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  8.4302h(11.8921h)-0.11h 	  -3.34m(11.8365h)  2436/7564	  2563.84	68.14		   86.49/85.14(0)85.14

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  49.94/23774.71(2.38)		  82.76(0)  82.76/0.80/(0.78)	 82.76/4.42/0.12(0.12)			 -1/ 0.78/ 0.12

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  358.54/356.20/-2.33		   -2.33(8)0.7084		 23.64(1)46.66(1.16)(1.08)	   509.72(10.6131h)(13.0598h)

----------------------------------------------
11 1603(5)				[2]MinSun=7.551544
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242178)117607.98( 3.04) 7.83	  167.24(55.21)10.34(14.85)	   314.89(22.86)20.22(26.44)	   (353.51)  132.27  50.35(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	243.19/18490.50(1.86)1(355.37)	  6.44/ 6.40(78.61)0.05(0)	(78.61) 0.56  5.42(5.42)1.0907	   1.858/5.418/0.547(13.129h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  20.6372h(9.7658h)-2.23h 	  -56.82m(8.8188h)  2224/7776	  2775.88	1325.52		   338.26/315.40(3)49.85

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  248.42/18242.08(1.85)		  317.24(3)  48.00/3.23/(1.69)	 134.62/3.46/1.83(1.83)			 -1/ 1.69/-1.83

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  357.49/355.37/-2.13		   -2.13(8)0.7343		 20.77(1)81.72(1.10)(1.01)	   546.11(7.5081h)(10.1294h)

----------------------------------------------
12 1604(5)				[2]MinSun=7.574996
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242177)117973.22( 8.28)18.71	  156.34(49.57) 6.50(15.46)	   303.99(17.23)16.38(27.05)	   (361.64)  121.37  61.25(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	184.59/20874.40(2.09)1(363.73)	  2.60/ 2.54(31.77)0.06(0)	(31.77) 8.28  3.14(3.14)1.0135	   2.092/3.142/0.423(10.164h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  5.4480h(15.6119h) 3.61h 	  78.90m(16.9268h)  2288/7712	  2711.98	2052.85		   321.70/304.47(3)60.77

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  190.00/20684.40(2.08)		  306.55(3)  58.69/2.62/(0.64)	 123.93/3.89/3.19(3.19)			 1/ 0.64/ 3.19

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  361.47/363.73/ 2.26		    2.26(6)0.6227		 16.80(0)36.93(1.02)(0.94)	   566.77(15.5666h)(18.2871h)

----------------------------------------------
12 1605( 3) ReqEcl(交泛)=25.341045	LunEcl(食月望交泛)=12.8941175

	SumDay=118338.465	[LunStep 3]: LunDay(食月經望)=56.8288195		LunSun=278.32644	LunMoon=25.33046
	Winter_solitice=13.525	[LunStep 4]: LunEclA(交常度)=172.37866026
	LeapDay= 0.0532800	[LunStep 5]: LunChSun(食月盈縮曆)=95.70534750
	SumDay + Q=118393.525	[LunStep 6]: LUseSun=86.91574050(3)
	RefMoonDay=13.471744	[LunStep 7]: LSunX(盈縮差)=2.386527878
	RefSun=174.969360	[LunStep 8]: LSunXS=2.386527878 	LunEclA_S(交定度)=174.765188134
	RefMoon= 4.636492	[LunStep 9]: LUseMoon(食月望遲疾曆)=11.5532795(1) 	Lmoon_x12_2(遲疾限)=140.9500099	LMoonCh=26.04999
	RefEclipse=18.385938	[LunStep 10]:(+0.08306=>)LMoonCh_=27.1330501	LMoonX(遲疾差)=2.74268889	LMoonSp=1.1844087
						古今律厯考卷五十(787-563)原計算推遲末限=總限(168.08306限)減遲限(138.3628779限)=29.7201821限
						有別大統以168限或167限減遲限計算,因此 LMoonCh_ 加回0.08306
						該遲行度以29限整數計算,得出1.1823。現改回29.7201821限計算,於是 LMoonSp=1.1814954

				[LunStep 11]:LMoonXS=-2.74268889	LSunMoonX=5.12921676 	LSunMoonT(加減差)=0.35511034
	ReqMoonDay=263.561139	[LunStep 12]:LConjuncT(食月定望)=57.18392984	LConjuncTDay=57		LConjuncTime=0.18392984
	ReqSun=80.940051	[LunStep 13]:LunSeeT(時差)=0.0070774	LunMaxT(食甚定分)=0.19100728  =>	LunMaxT24=4.5841748h
	ReqMoon=10.565167

				[LunStep 16]:LunEclSun=278.68862 	LEcl_ChSun(食甚入盈縮曆)=96.06753528 	LEclUseSun=86.55355272(3)
				[LunStep 17]:LEcl_SunX(食甚盈縮差)=2.38497926
				[LunStep 18]:LEcl_SunX_=2.38497926 	LEclNetSun(食甚入盈縮曆定度)=98.45251454(3)
					     原計算=103.57495347(食甚入盈縮曆)+2.33579615(食甚盈縮差)=105.91074962
						     但古今律厯考卷五十(787-565)給出錯誤食甚入盈曆定度=105.93074962

						LEclNetSun_= LEclNetSun
						if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875}
							As LEclSunStatus=3	LEclNetSun_=288.67971254

						total_year=calculated_year-Base_year=324 => Added_drift=0.015*total_year=4.86000
						Ref_drift(天正加時赤道度)=10-Added_drift=5.14000(邢雲路計算方法)
						  (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=5.03502)
						Next_Ref_drift=Ref_drift_check-4.3445=0.69052

						Sun_drift=Next_Ref_drift/1.0849=0.63649
						Net_drift=4+Sun_drift=4.63649
						LEclNetSun_All=LEclNetSun_ + Net_drift=293.31620
						LEclNetSun_From(食甚月躔黃道軫宿)=LEclNetSun_All-277.9775(翼宿)=15.33870
				[LunStep 22]:
					求交常交定度:(refer to LunStep 4 and LunStep 8)
						置經朔望入交泛日及分秒(LunEcl),以月平行度(13.36878311)乘之,為交常度(LunEclA);
						以盈縮差(LSunXS)盈加縮減之,為交定度(LunEclA_S)。

					求月食入陰陽曆去交前後度:
						視交定度,在交中度(181.8967)已下,為陽曆;已上,減去交中,為陰曆。
						視入陰陽曆,在後準十五度半已下,為交後度;
						前準一百六十六度三十九分六十八秒已上,覆減交中,餘為交前度及分。

						交中度=181.8967
						SexLunEclA=LunEclA_S=174.765188134
						if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S}
						Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=174.765188134

						LunEclFront=-1(負數[-1]表示交前度);
						if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);}

						SexLunEclA_=SexLunEclA=174.765188134
						if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1}
						if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1}
						Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= 7.131511866(-1)
				[LunStep 23]:
					求月食分秒:
						視去交前後度,
						不用南北東西差者。用減食限,不及減者不食。
						餘如定法而一,為月食之分秒。
						LunFraction=(13.05-abs(SexLunEclA_))/8.7=0.680285992
				[LunStep 24]:
					LEclM=25.68557384(1)	LEclM_Used(定入遲疾曆)=11.90838984
					LEclM_x12_2(定入遲疾限)=145.28235601	LEclM_Ch=21.71764399	LEclM_Sp=1.18899494
					LEclM_SunSp(定限行度)=1.10699	    LHalfEclT(定用分)=0.05583191(邢雲路用 4920,不是 5740)
				[LunStep 25]:
					LEclBeginT(初虧)=0.135175	   LEclEndT(復圓)=0.246839
					LEclBeginT24=3.24420901h	LEclEndT24=5.92414049h
----------------------------------------------------------
13 1607(3)				[2]MinSun=7.645352
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242174)119068.95(24.01)21.80	  153.17( 2.21)24.51(19.60)	   241.76(30.80) 2.89(26.56)	   (355.04)  59.14  59.14(2)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	187.65/21085.95(2.11)1(357.15)	  2.89/ 2.87(35.21)0.02(2)	(35.21) 7.86  3.42(3.42)1.1749	   2.111/-3.422/-0.091(-2.195h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  19.1461h(16.9514h) 4.95h 	  90.89m(18.4662h)  2276/7724	  2275.93	2694.25		   272.53/241.74(2)59.12

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  187.65/21085.95(2.11)		  243.85(2)  61.23/2.46/(-0.45)	 61.23/3.97/4.28(3.67)			 1/ 0.45/ 3.67

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  361.76/357.15/-4.60		   -4.60(8)0.4244		  2.79(2)34.09(1.18)(1.09)	   428.93(17.4368h)(19.4956h)

----------------------------------------------
13 1609(8)				[2]MinSun=7.692256
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242172)119799.43(34.49)14.02	  160.91(20.47)18.81(23.14)	   31.91(16.72) 7.06(14.47)	   (193.47)  31.91  31.91(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	339.79/12894.61(1.32)-1(192.15)	  7.06/ 7.05(86.18)0.01(3)	(80.82) 0.09  5.43(5.43)1.0954	   -1.320/-5.429/-0.505(-12.126h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  17.1731h(5.0473h)-6.95h 	  -91.39m(3.5242h)  2134/7866	  2866.46	3531.58		   48.06/31.34(0)31.34

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  339.79/12894.61(1.30)		  30.04(0)  30.04/3.98/(-0.92)	 30.04/2.45/3.46(1.44)			 -1/ 0.92/-1.44

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  187.53/192.15/ 4.62		   -4.62(8)0.4228		  6.56(2)80.00(1.10)(1.02)	   460.92(2.4180h)(4.6304h)

----------------------------------------------
15 1610(13)				[2]MinSun=7.715708
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242171)120164.67(39.73)24.90	  150.01(14.84)14.97(23.75)	   168.66(38.74)13.11(26.67)	   (356.57)  168.66  13.96(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	445.20/6250.61(0.67)1(355.90)	 13.11/13.04(159.88)0.07(3)	( 7.12)10.63  0.88(0.88)0.9901	   -0.668/-0.881/-0.128(-3.077h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  17.6550h(14.5777h) 2.58h 	  63.25m(15.6318h)  2081/7919	  2080.57	1513.26		   207.32/168.58(1)14.04

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  440.02/6695.81(0.67)		  167.91(1)  14.71/4.34/(1.18)	 167.91/1.32/0.80(0.80)			 1/-1.18/-0.80

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  355.66/355.90/ 0.25		    0.25(6)0.9588		 12.98(3) 8.68(0.99)(0.91)	   630.79(14.1179h)(17.1457h)

----------------------------------------------
15 1612(6)				[2]MinSun=7.762612
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242169)120895.16(50.22)17.12	  157.74(33.10) 9.27( 0.07)	   334.93(30.29)21.12(13.98)	   (186.87)  152.31  30.32(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	344.71/12549.90(1.27)-1(188.14)	  7.35/ 7.30(89.62)0.05(1)	(77.38) 0.82  5.42(5.42)1.1045	   1.266/5.417/0.496(11.907h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  6.8712h(18.7785h) 6.78h 	  92.17m(20.3147h)  2130/7870	  2870.34	3464.46		   365.77/335.49(3)29.76

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  349.62/12200.28(1.25)		  336.73(3)  28.51/4.03/(-0.83)	 154.11/2.35/3.26(1.44)			 1/ 0.83/-1.44

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  187.44/188.14/ 0.70		   -0.70(8)0.9130		 21.62(1)71.33(1.12)(1.03)	   552.64(18.9884h)(21.6410h)

----------------------------------------------
16 1613(6)				[2]MinSun=7.786064
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242168)121260.40(55.46)27.99	  146.85(27.47) 5.43( 0.68)	   324.03(24.65)17.28(14.59)	   (195.00)  141.41  41.21(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	289.68/16068.36(1.61)-1(196.61)	  3.51/ 3.44(42.77)0.06(0)	(42.77) 6.85  3.98(3.98)1.0278	   1.613/3.984/0.447(10.718h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  15.6820h(2.3996h)-9.60h 	  -59.99m(1.3997h)  2178/7822	  2822.46	4416.78		   349.09/324.43(3)40.81

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  294.76/15773.60(1.60)		  326.03(3)  39.21/3.64/(-2.05)	 143.41/3.01/5.31(0.70)			 -1/ 2.05/ 0.70

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  190.81/196.61/ 5.81		   -5.81(8)0.2740		 17.73(0)48.22(1.04)(0.95)	   413.91(0.4063h)(2.3931h)

----------------------------------------------
17 1615(4)				[2]MinSun=7.832968
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242166)121990.88( 5.94)20.21	  154.58(45.73)27.28( 4.21)	   272.70(43.86) 7.63(13.48)	   (180.27)  90.08  92.54(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	 8.90/24001.62(2.40)-1(182.67)	  7.63/ 7.63(93.06)0.00(3)	(73.94) 1.52  5.38(5.38)1.0811	   2.401/-5.380/-0.226(-5.424h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  20.5693h(15.1458h) 3.15h 	  72.54m(16.3547h)  2489/7511	  2510.86	1814.47		   316.38/272.53(3)92.72

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	   8.90/24001.62(2.40)		  274.93(3)  90.32/0.10/(0.03)	 92.31/4.46/3.24(3.24)			 1/-0.03/-3.24

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  184.79/182.67/-2.12		    2.12(6)0.6471		  7.40(3)76.70(1.09)(1.00)	   534.55(15.0718h)(17.6376h)

----------------------------------------------
17 1617( 8) ReqEcl(交泛)=26.292590	LunEcl(食月望交泛)=13.8456625
		《大明神宗顯皇帝實錄卷之五百四十七》
				邢雲路獻七政真數..(望交泛分:13.845663)
		以此法布算今歲七月十六日戊寅夜望月食,推得是月:
			望交泛分一十三日八十四刻五十六分六十三秒(13.845663)
			陰歷交前度一度二十二分零三秒九十五微(1.220395)
			月食分一十三分五十九秒六十九微(1.35969)
			定用分六刻六十八分九十五秒二十一微(0.06689521)
			既內分二刻八十六分六十五秒
			既外分三刻八十一分七十九秒
			初虧分六刻四十三分二十秒八十九微(0.064[3]2089,應修正為:復圓分-2*定用分=0.19801131-0.06689521*2=0.064[2]2089)
			食既分一十刻二十五分零四十微
			食甚分一十三刻一十一分六十六秒一十微(0.1311[6]610,應修正為:復圓分-定用分=0.19801131-0.06689521=0.1311[1]610)
			生光分一十五刻九十八分三十一秒八十微
			復圓分一十九刻八十分一十一秒三十一微(0.19801131)

				   初虧丑初二刻(1.48h<=>0.064[2]2089*24=1.54130136,仍然算是丑初二刻,所以該記錄是可靠的)
			食既丑正一刻,食甚寅初二刻(3.48h? <=>0.1311[1]610*24=3.1467864h=>應修正為寅初一刻)
			生光寅初三刻,復明寅正三刻(4.72h<=>0.19801131*24=4.75227144h=寅正三刻,所以該記錄亦是可靠的)

	SumDay=122721.367	[LunStep 3]: LunDay(食月經望)=15.0095485		LunSun=48.08149	LunMoon=24.59484
	Winter_solitice=16.427	[LunStep 4]: LunEclA(交常度)=185.09965897
	LeapDay=12.4275960	[LunStep 5]: LunChSun(食月盈縮曆)=48.08149050
	SumDay + Q=122776.427	[LunStep 6]: LUseSun=48.08149050(0)
	RefMoonDay= 3.999508	[LunStep 7]: LSunX(盈縮差)=1.800930527
	RefSun=162.313614	[LunStep 8]: LSunXS=-1.800930527 	LunEclA_S(交定度)=183.298728443
	RefMoon=21.574112	[LunStep 9]: LUseMoon(食月望遲疾曆)=10.8176565(1) 	Lmoon_x12_2(遲疾限)=131.9754093	LMoonCh=35.02459
	RefEclipse= 7.745638	[LunStep 10]:(+0.08306=>)LMoonCh_=36.1076507	LMoonX(遲疾差)=3.49220640	LMoonSp=1.1737438
				[LunStep 11]:LMoonXS=-3.49220640	LSunMoonX=1.69127587 	LSunMoonT(加減差)=0.11815579
	ReqMoonDay=33.316194	[LunStep 12]:LConjuncT(食月定望)=15.12770429	LConjuncTDay=15		LConjuncTime=0.12770429
	ReqSun=33.316194	[LunStep 13]:LunSeeT(時差)=0.0034118	LunMaxT(食甚定分)=0.13111609  =>	LunMaxT24=3.1467861h
	ReqMoon=9.829544
											邢雲路食甚分 0.13111610*24=3.1467864h
													   理論值:3.219

				[LunStep 16]:LunEclSun=48.20306 	LEcl_ChSun(食甚入盈縮曆)=48.20305809 	LEclUseSun=48.20305809(0)
				[LunStep 17]:LEcl_SunX(食甚盈縮差)=1.80403655
				[LunStep 18]:LEcl_SunX_=-1.80403655 	LEclNetSun(食甚入盈縮曆定度)=46.39902153(0)
						LEclNetSun_= LEclNetSun
						if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875}
							As LEclSunStatus=0	LEclNetSun_=54.27889353

						total_year=calculated_year-Base_year=336 => Added_drift=0.015*total_year=5.04000
						Ref_drift(天正加時赤道度)=10-Added_drift=4.96000(邢雲路計算方法)
						  (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.84710)
						Next_Ref_drift=Ref_drift_check-4.3445=0.50260

						Sun_drift=Next_Ref_drift/1.0849=0.46327
						Net_drift=4+Sun_drift=4.46327
						LEclNetSun_All=LEclNetSun_ + Net_drift=58.74217
						LEclNetSun_From(食甚月躔黃道虚宿)=LEclNetSun_All-51.08(女宿)=7.66217
				[LunStep 22]:
					求交常交定度:(refer to LunStep 4 and LunStep 8)
						置經朔望入交泛日及分秒(LunEcl),以月平行度(13.36878311)乘之,為交常度(LunEclA);
						以盈縮差(LSunXS)盈加縮減之,為交定度(LunEclA_S)。

					求月食入陰陽曆去交前後度:
						視交定度,在交中度(181.8967)已下,為陽曆;已上,減去交中,為陰曆。
						視入陰陽曆,在後準十五度半已下,為交後度;
						前準一百六十六度三十九分六十八秒已上,覆減交中,餘為交前度及分。

						交中度=181.8967
						SexLunEclA=LunEclA_S=183.298728443
						if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S}
						Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-1.402028443

						LunEclFront=-1(負數[-1]表示交前度);
						if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);}

						SexLunEclA_=SexLunEclA=-1.402028443
						if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1}
						if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1}
						Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -1.402028443(-1)
						邢雲路獻七政真數..以此法布算..陰歷交前度一度二十二分零三秒九十五微=>-1.220395(-1)
				[LunStep 23]:
					求月食分秒:
						視去交前後度,
						不用南北東西差者。用減食限,不及減者不食。
						餘如定法而一,為月食之分秒。
						LunFraction=(13.05-abs(SexLunEclA_))/8.7=1.338847305
					邢雲路獻七政真數..以此法布算..月食分一十三分五十九秒六十九微(1.35969)
											 理論值:1.403
				[LunStep 24]:
					LEclM=24.71299629(1)	LEclM_Used(定入遲疾曆)=10.93581229
					LEclM_x12_2(定入遲疾限)=133.41690993	LEclM_Ch=33.58309007	LEclM_Sp=1.17556264
					LEclM_SunSp(定限行度)=1.09356	    LHalfEclT(定用分)=0.06709524(邢雲路用 4920,不是 5740)
					邢雲路獻七政真數..以此法布算..定用分六刻六十八分九十五秒二十一微(0.06689521)
				[LunStep 25]:
					LEclBeginT(初虧)=0.064021	   LEclEndT(復圓)=0.198211
					   (邢雲路的初虧: 0.064[2]2089		    復圓: 0.19801131)

					LEclBeginT24=1.53650019h	LEclEndT24=4.75707192h
					(邢雲路的初虧*24=1.54130136h=>算是丑初二刻	   4.75227144h=寅正三刻,兩項記錄都可靠)
						理論值: 1.49h      		   4.95h

				邢雲路原有的陰陽曆 SexLunEclA_1617(不分陰陽)=13.05-1.35969*8.7=1.22070,與記錄所指的 1.220395 基本一致

				邢雲路變法後的食分 NewFraction=(13.125-1.220395)/7.5=1.58728,大於理論值
									     理論值=1.403
								   邢雲路原有的食分: 1.35969

		  邢雲路原有的定限行度 LEclM_SunSp_1617=(sqrt((30-13.5969)*13.5969)*49.20)/(0.06689521*10000)=1.09838
		  邢雲路變法後改大的定用分 NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1617*10000)=0.07805
		  邢雲路變法後的 NewBeginT=0.13111610-NewHalfT=0.05307	 NewEndT=0.13111610+NewHalfT=0.20916

			變法後的初虧=NewBeginT*24=1.27362		食甚不變=0.13111610*24=3.14679		復圓=NewEndT*24=5.01995
				     理論值: 初虧=1.49h				  食甚=3.219h			   復圓=4.95h
				 邢雲路原有的初虧=1.54130136h							   復圓=4.75227144h
----------------------------------------------------------
18 1617(8)				[2]MinSun=7.879872
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242164)122721.37(16.43)12.43	  162.31( 4.00)21.57( 7.75)	   33.32( 0.24) 9.83(26.29)	   (351.50)  33.32  33.32(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	329.90/13569.26(1.37)1(350.13)	  9.83/ 9.76(119.92)0.07(3)	(47.08) 6.23  4.33(4.33)1.0340	   -1.367/-4.330/-0.452(-10.844h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  5.8620h(19.0180h) 7.02h 	  91.05m(20.5355h)  2142/7858	  2858.39	3556.48		   34.17/33.93(0)33.93

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  329.90/13569.26(1.39)		  32.54(0)  32.54/3.89/(-0.95)	 32.54/2.61/3.72(1.51)			 1/-0.95/-1.51

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  355.18/350.13/-5.05		   -5.05(8)0.3689		  9.38(3)52.59(1.04)(0.96)	   463.56(19.4230h)(21.6481h)

----------------------------------------------
18 1620(12) ReqEcl(交泛)=12.496823	LunEcl(食月望交泛)= 0.0498955
			《大明熹宗哲皇帝實錄卷之七》
				按新法推泰昌元年庚申歲十一月十六日己丑夜望月食:
					初虧漏下二百七十三籌五十三分,計九十七刻二十四分(0.9724)
					食既漏下二百九十四籌六十分,計八十八分
					生光漏下三百三十二籌四十分,計七刻四十分
					復圓漏下三百五十三籌四十八分,計一十一刻(0.11)
				各以發斂求之,得:
					初虧夜子初一刻(23.24h)
					食既子正初刻
					食甚子正四刻(0.96h)
					生光丑初三刻
					復圓丑正二刻(2.48h,但見有改自三刻的痕跡)
					月食一十四分九十九秒(1.499)
					食甚月離黃道畢宿一十四度三十分一十秒(14.3010)
				而授時則推:
					初虧子正一刻(0.24h)
					食既丑初一刻
					食甚丑正初刻(2.00h?)
					生光丑正三刻
					復圓寅初三刻(3.72h?)
					月食一十三分三十一秒(1.331)
					月離黃道畢宿一十四度三十四分四十四秒(14.3444),與天不合

	SumDay=123817.093	[LunStep 3]: LunDay(食月經望)=25.7638615		LunSun=163.03997	LunMoon=22.95696
	Winter_solitice=32.153	[LunStep 4]: LunEclA(交常度)=0.66704211
	LeapDay=15.5211300	[LunStep 5]: LunChSun(食月盈縮曆)=163.03997400
	SumDay + Q=123872.153	[LunStep 6]: LUseSun=19.58110650(1)
	RefMoonDay=16.631449	[LunStep 7]: LSunX(盈縮差)=0.908488692
	RefSun=159.149722	[LunStep 8]: LSunXS=-0.908488692 	LunEclA_S(交定度)=-0.241446579
	RefMoon=12.031333	[LunStep 9]: LUseMoon(食月望遲疾曆)= 9.1797775(1) 	Lmoon_x12_2(遲疾限)=111.9932855	LMoonCh=55.00671
	RefEclipse=11.888619	[LunStep 10]:(+0.08306=>)LMoonCh_=56.0897745	LMoonX(遲疾差)=4.77402895	LMoonSp=1.1443566
				[LunStep 11]:LMoonXS=-4.77402895	LSunMoonX=3.86554026 	LSunMoonT(加減差)=0.27698911
	ReqMoonDay=148.274677	[LunStep 12]:LConjuncT(食月定望)=26.04085061	LConjuncTDay=26		LConjuncTime=0.04085061
	ReqSun=34.346403	[LunStep 13]:LunSeeT(時差)=0.0003491	LunMaxT(食甚定分)=0.04119972  =>	LunMaxT24=0.9887933h[但丑正初刻略遲]
	ReqMoon=8.191665
   邢雲路..按新法推..初虧..計九十七刻二十四分..復圓..計一十一刻=>食甚分=0.11-(0.11+(1-0.9724))/2=0.0412=>食甚=0.0412*24=0.9888h=子正四刻,與記錄一致
													    理論值:1.17h 

				[LunStep 16]:LunEclSun=163.31731 	LEcl_ChSun(食甚入盈縮曆)=163.31731222 	LEclUseSun=19.30376828(1)
				[LunStep 17]:LEcl_SunX(食甚盈縮差)=0.89700279
				[LunStep 18]:LEcl_SunX_=-0.89700279 	LEclNetSun(食甚入盈縮曆定度)=162.42030943(1)
						LEclNetSun_= LEclNetSun
						if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875}
							As LEclSunStatus=1	LEclNetSun_=170.37053743

						total_year=calculated_year-Base_year=339 => Added_drift=0.015*total_year=5.08500
						Ref_drift(天正加時赤道度)=10-Added_drift=4.91500(邢雲路計算方法)
						  (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.80008)
						Next_Ref_drift=Ref_drift_check-4.3445=0.45558

						Sun_drift=Next_Ref_drift/1.0849=0.41993
						Net_drift=4+Sun_drift=4.41993
						LEclNetSun_All=LEclNetSun_ + Net_drift=174.79046
					LEclNetSun_From(食甚月躔黃道畢宿)=LEclNetSun_All-160.8175(昴宿)=13.97296(幾乎=14.3444)
							而授時則推....月離黃道畢宿一十四度三十四分四十四秒(14.3444 與計算值符合=>記錄可靠)
						  邢雲路...按新法推...食甚月離黃道畢宿一十四度三十分一十秒(14.3010)
						按日行每天一度,估計授時比新法遲(14.3444-14.3010)*24=1.0416h
						與食甚遲的時刻比較,差距亦大致相同 LunMaxT24-0.9888h=-0.0000067h,表示該些記錄亦可靠
				[LunStep 22]:
					求交常交定度:(refer to LunStep 4 and LunStep 8)
						置經朔望入交泛日及分秒(LunEcl),以月平行度(13.36878311)乘之,為交常度(LunEclA);
						以盈縮差(LSunXS)盈加縮減之,為交定度(LunEclA_S)。

					求月食入陰陽曆去交前後度:
						視交定度,在交中度(181.8967)已下,為陽曆;已上,減去交中,為陰曆。
						視入陰陽曆,在後準十五度半已下,為交後度;
						前準一百六十六度三十九分六十八秒已上,覆減交中,餘為交前度及分。

						交中度=181.8967
						SexLunEclA=LunEclA_S=-0.241446579
						if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S}
						Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-0.241446579

						LunEclFront=-1(負數[-1]表示交前度);
						if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);}

						SexLunEclA_=SexLunEclA=-0.241446579
						if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1}
						if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1}
						Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -0.241446579(1)
				[LunStep 23]:
					求月食分秒:
						視去交前後度,
						不用南北東西差者。用減食限,不及減者不食。
						餘如定法而一,為月食之分秒。
						LunFraction=(13.05-abs(SexLunEclA_))/8.7=1.472247520
							    授時則推..月食一十三分三十一秒(1.331)
						邢雲路.....按新法推....月食一十四分九十九秒(1.499)
										   理論值:1.6
				[LunStep 24]:
					LEclM=23.23395061(1)	LEclM_Used(定入遲疾曆)=9.45676661
					LEclM_x12_2(定入遲疾限)=115.37255260	LEclM_Ch=51.62744740	LEclM_Sp=1.14987339
					LEclM_SunSp(定限行度)=1.06787	    LHalfEclT(定用分)=0.06909749(邢雲路用 4920,不是 5740)
				[LunStep 25]:
					LEclBeginT(初虧)=0.972102	   LEclEndT(復圓)=0.110297
		邢雲路...按新法推..初虧..計九十七刻二十四分(0.9724)		複圓....計一十一刻(0.11 )

					LEclBeginT24=23.33045365h[初虧子正一刻]是可靠的	LEclEndT24=2.64713302h[復圓寅初三刻]略為稍遲
		 (邢雲路...按新法推..初虧=0.9724*24 = 23.3376h				復圓=0.11*24=2.64h)
					理論值:初虧= 23.28h					復圓=3.05h

				邢雲路原有的陰陽曆 SexLunEclA_1620(不分陰陽)=13.05-1.499*8.7=0.00870
				邢雲路變法後的食分 NewFraction=(13.125-SexLunEclA_1620)/7.5=1.74884

				邢雲路原有的定用分 HalfT_1620=(0.11+(1-0.9724))/2=0.06880
				邢雲路原有的定限行度 LEclM_SunSp_1620=(sqrt((30-14.99)*14.99)*49.20)/(HalfT_1620*10000)=1.07267
		  邢雲路變法後的定用分 NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1620*10000)=0.08027
				邢雲路變法後的 NewBeginT=1+0.0412-NewHalfT=0.96093	NewEndT=0.0412+NewHalfT=0.12147

		邢雲路變法後的初虧=NewBeginT*24=23.06240		食甚不變=0.0412*24=0.98880	     復圓=NewEndT*24=2.91520
				  理論值:初虧= 23.28h			     食甚=1.17h				復圓=3.05h
		    邢雲路原有的初虧=0.9724*24=23.33760							復圓=0.11*24=2.64000h
----------------------------------------------------------
18 1620(13) LunEcl=2.3683	LunMoon=24.9332	LunFraction=-2.1950	LunMaxT24=13.1410h	LEclNetSun_=201.2793
18 1621( 0) LunEcl=0.0499	LunMoon=22.9570	LunFraction=1.4721	LunMaxT24=0.9871h	LEclNetSun_=170.3698
18 1621( 1) LunEcl=2.3683	LunMoon=24.9332	LunFraction=-2.1949	LunMaxT24=13.1392h	LEclNetSun_=201.2785
18 1621( 2) LunEcl=4.6866	LunMoon=26.9094	LunFraction=-5.8883	LunMaxT24=23.8593h	LEclNetSun_=231.8570
18 1621( 3) LunEcl=7.0050	LunMoon=1.3313	LunFraction=-9.5249	LunMaxT24=9.4888h	LEclNetSun_=261.8989
18 1621( 4) LunEcl=9.3234	LunMoon=3.3075	LunFraction=-13.0994	LunMaxT24=18.6015h	LEclNetSun_=291.3826
18 1621( 5) LunEcl=11.6417	LunMoon=5.2837	LunFraction=-16.6140	LunMaxT24=4.6552h	LEclNetSun_=320.3892
18 1621( 6) LunEcl=13.9601	LunMoon=7.2599	LunFraction=0.8322	LunMaxT24=14.8351h	LEclNetSun_=348.9415
18 1621( 6) ReqEcl(交泛)=26.407037	LunEcl(食月望交泛)=13.9601095

	SumDay=124182.334	[LunStep 3]: LunDay(食月經望)=22.9474195		LunSun=340.20042	LunMoon=7.25994
	Winter_solitice=37.394	[LunStep 4]: LunEclA(交常度)=186.62967609
	LeapDay=26.3958350	[LunStep 5]: LunChSun(食月盈縮曆)=157.57933950
	SumDay + Q=124237.394	[LunStep 6]: LUseSun=25.04174050(3)
	RefMoonDay=10.998565	[LunStep 7]: LSunX(盈縮差)=1.076856479
	RefSun=148.251565	[LunStep 8]: LSunXS=1.076856479 	LunEclA_S(交定度)=187.706532570
	RefMoon= 8.191665	[LunStep 9]: LUseMoon(食月望遲疾曆)= 7.2599435(3) 	Lmoon_x12_2(遲疾限)=88.5713107	LMoonCh=78.42869
	RefEclipse=12.496823	[LunStep 10]:(+0.08306=>)LMoonCh_=79.5117493	LMoonX(遲疾差)=5.42352139	LMoonSp=1.0926015
				[LunStep 11]:LMoonXS=5.42352139	LSunMoonX=-4.34666491 	LSunMoonT(加減差)=-0.32621823
	ReqMoonDay=325.435123	[LunStep 12]:LConjuncT(食月定望)=22.62120127	LConjuncTDay=22		LConjuncTime=0.62120127
	ReqSun=39.807037	[LunStep 13]:LunSeeT(時差)=-0.0030732	LunMaxT(食甚定分)=0.61812810  =>	LunMaxT24=14.8350745h
	ReqMoon=20.049015

				[LunStep 16]:LunEclSun=339.87113 	LEcl_ChSun(食甚入盈縮曆)=157.25004810 	LEclUseSun=25.37103190(3)
				[LunStep 17]:LEcl_SunX(食甚盈縮差)=1.08905677
				[LunStep 18]:LEcl_SunX_=1.08905677 	LEclNetSun(食甚入盈縮曆定度)=158.33910487(3)
						LEclNetSun_= LEclNetSun
						if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875}
							As LEclSunStatus=3	LEclNetSun_=348.94153487

						total_year=calculated_year-Base_year=340 => Added_drift=0.015*total_year=5.10000
						Ref_drift(天正加時赤道度)=10-Added_drift=4.90000(邢雲路計算方法)
						  (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.78440)
						Next_Ref_drift=Ref_drift_check-4.3445=0.43990

						Sun_drift=Next_Ref_drift/1.0849=0.40548
						Net_drift=4+Sun_drift=4.40548
						LEclNetSun_All=LEclNetSun_ + Net_drift=353.34701
				[LunStep 22]:
					求交常交定度:(refer to LunStep 4 and LunStep 8)
						置經朔望入交泛日及分秒(LunEcl),以月平行度(13.36878311)乘之,為交常度(LunEclA);
						以盈縮差(LSunXS)盈加縮減之,為交定度(LunEclA_S)。

					求月食入陰陽曆去交前後度:
						視交定度,在交中度(181.8967)已下,為陽曆;已上,減去交中,為陰曆。
						視入陰陽曆,在後準十五度半已下,為交後度;
						前準一百六十六度三十九分六十八秒已上,覆減交中,餘為交前度及分。

						交中度=181.8967
						SexLunEclA=LunEclA_S=187.706532570
						if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S}
						Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-5.809832570

						LunEclFront=-1(負數[-1]表示交前度);
						if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);}

						SexLunEclA_=SexLunEclA=-5.809832570
						if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1}
						if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1}
						Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -5.809832570(-1)
				[LunStep 23]:
					求月食分秒:
						視去交前後度,
						不用南北東西差者。用減食限,不及減者不食。
						餘如定法而一,為月食之分秒。
						LunFraction=(13.05-abs(SexLunEclA_))/8.7=0.832203153
				[LunStep 24]:
					LEclM=6.93372527(3)	LEclM_Used(定入遲疾曆)=6.93372527
					LEclM_x12_2(定入遲疾限)=84.59144833	LEclM_Ch=82.40855167	LEclM_Sp=1.10120211
					LEclM_SunSp(定限行度)=1.01920	    LHalfEclT(定用分)=0.06483787(邢雲路用 4920,不是 5740)
				[LunStep 25]:
					LEclBeginT(初虧)=0.553290	   LEclEndT(復圓)=0.682966
					LEclBeginT24=13.27896554h	LEclEndT24=16.39118343h
----------------------------------------------------------
		《大明熹宗哲皇帝實錄卷之七》
		原任陝西按察使邢雲路獻七政真數.....復以此法推天啟元年辛酉歲四月壬申朔日食:
				初虧申正一刻(16.24h)
				食甚酉初刻(17.00h)
				復圓酉正初刻(18.00h)
				日食一分八十六秒(0.186)
				復圓日在天未入地,食不及三分不救(以上記錄缺乏其它數據可供證驗,是否可靠,甚有疑問?)
		而授時則推是日食:
				初虧申正三刻(16.72h)
				食甚酉正初刻(18.00h
				復圓戍初刻(19.00h)
				日食三分九十一秒(0.391)
				日未入已復光三分一十一秒,日已入未復光八十秒,與天不合(以上授時數據可按計算方法證驗)

		《明熹宗悊皇帝實錄卷之十六》
		庚申原任陝西按察使閒住邢雲路奏:
				臣本年春言四月朔日食歷數,臨期策驗分刻乃有不合
				臣始悟歷元玄法原自無差,而其差者,乃今大統歷沿郭守敬授時歷之誤也
				蓋日食者,月食之也,故正歷一準於月食,而月之分數與限度定法,皆其綱領最大者
				授時定月食分三十分及十五分,月食限十三度五分,定限八十七分
				此守敬率南北日官以儀器測驗于天,自丙子至庚辰又經甲午,再驗幾二十年,自謂推步之精,可永久無弊
				以故大統歷沿用之,臣今所步歷亦用之,及四月朔食測驗,方見其差,乃悟以異乘同
				除另變其法,變月食分三十分為三十五分(30->35)
				         變十五分為十七分五十秒
				         變定限十三度五分為十三度十二分五十秒(13.05->13.1250)
				         變定法八十七分為七十五分(8.7->7.5)
				諸數既變,乃依舊歷元玄法布歷,一一方合,然後知歷元法原自無差,而其差者,大統相沿授時法誤之也,故必如臣所推:
				如五千七百四十為日食成限,以四千九百二十為月食成限,以三十五分為月食分,而黜郭守敬之三十分,方為正法,下禮部

				  LunFraction=(13.05 - abs(SexLunEclA_))/8.7
			月食食分=定限(13.05->13.1250) - 視入陰陽曆/定法(8.7->7.5)

				LHalfEclT=(sqrt((30-LunFraction*10)*LunFraction*10)*49.20)/(LEclM_SunSp*10000)
				定用分=(sqrt(三十分(30->35)-月食分)*月食分)*49.20)/(定限行度*10000)
19 1621(6)				[2]MinSun=7.973680
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242160)124182.33(37.39)26.40	  148.25(11.00) 8.19(12.50)	   325.44( 8.18)20.05(26.41)	   (353.03)  142.81  39.81(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	299.83/15473.77(1.57)1(354.60)	  6.27/ 6.23(76.52)0.04(0)	(76.52) 1.00  5.40(5.40)1.0863	   1.572/5.400/0.526(12.629h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  4.3710h(17.0004h) 5.00h 	  91.15m(18.5196h)  2168/7832	  2831.96	2716.49		   334.21/326.02(3)39.22

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  299.83/15473.77(1.55)		  327.58(3)  37.66/3.70/(0.15)	 144.96/2.92/3.17(2.67)			 1/ 0.15/ 2.67

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  360.46/354.60/-5.86		   -5.86(8)0.2680		 20.58(0)82.94(1.10)(1.02)	   384.15(17.5976h)(19.4415h)

----------------------------------------------
19 1621( 7) LunEcl=16.2785	LunMoon=9.2362	LunFraction=-2.5819	LunMaxT24=2.2018h	LEclNetSun_=11.8802
19 1621( 8) LunEcl=18.5968	LunMoon=11.2124	LunFraction=-6.0092	LunMaxT24=15.3730h	LEclNetSun_=40.2466
19 1621( 9) LunEcl=20.9152	LunMoon=13.1886	LunFraction=-9.4863	LunMaxT24=7.3509h	LEclNetSun_=69.1614
19 1621(10) LunEcl=23.2336	LunMoon=15.1648	LunFraction=-13.0180	LunMaxT24=0.4774h	LEclNetSun_=98.6033
19 1621(11) LunEcl=25.5520	LunMoon=17.1411	LunFraction=-16.6120	LunMaxT24=17.4603h	LEclNetSun_=128.5868
19 1621(12) LunEcl=0.6581	LunMoon=19.1173	LunFraction=0.6413	LunMaxT24=10.2391h	LEclNetSun_=159.0944
19 1621(13) LunEcl=2.9765	LunMoon=21.0935	LunFraction=-3.0682	LunMaxT24=1.4936h	LEclNetSun_=190.0187
19 1623(10) ReqEcl(交泛)=12.003066	LunEcl(食月望交泛)=26.7683625
			《大明熹宗哲皇帝實錄卷之七》
				以新法推天啟三年癸亥歲九月十五壬寅夜望月食:
					初虧申正三刻
					食甚酉初三刻(17.72?,該記錄不可靠)
					復圓戍初初刻
					月食五分三十八秒
					月未出巳食五分二十五秒,月已出見食一十三秒,不及一分不救
				而授時則推月食:
					初虧酉初一刻
					食甚酉正二刻(18.48h,該記錄可靠)
					復圓戌正一刻
					月食七分八秒
					月未出已食一分九十八秒,月已出見食五分一十秒,與天不合

	SumDay=124912.818	[LunStep 3]: LunDay(食月經望)=39.3346165		LunSun=100.81492	LunMoon=9.46173
	Winter_solitice=47.878	[LunStep 4]: LunEclA(交常度)=357.86043247
	LeapDay=18.6146460	[LunStep 5]: LunChSun(食月盈縮曆)=100.81491750
	SumDay + Q=124967.878	[LunStep 6]: LUseSun=81.80616150(1)
	RefMoonDay=29.263390	[LunStep 7]: LSunX(盈縮差)=2.383266051
	RefSun=155.985849	[LunStep 8]: LSunXS=-2.383266051 	LunEclA_S(交定度)=355.477166424
	RefMoon= 2.488554	[LunStep 9]: LUseMoon(食月望遲疾曆)= 9.4617325(3) 	Lmoon_x12_2(遲疾限)=115.4331365	LMoonCh=51.56686
	RefEclipse=16.031600	[LunStep 10]:(+0.08306=>)LMoonCh_=52.6499235	LMoonX(遲疾差)=4.59614481	LMoonSp=1.0426297
				[LunStep 11]:LMoonXS=4.59614481	LSunMoonX=-6.97941086 	LSunMoonT(加減差)=-0.54891173
	ReqMoonDay=86.049621	[LunStep 12]:LConjuncT(食月定望)=38.78570477	LConjuncTDay=38		LConjuncTime=0.78570477
	ReqSun=86.049621	[LunStep 13]:LunSeeT(時差)=-0.0096072	LunMaxT(食甚定分)=0.77609756  =>	LunMaxT24=18.6263416h
	ReqMoon=22.250804
											  而授時則推...食甚酉正二刻(18.48h),該記錄可靠
													   理論值:19.23h 

				[LunStep 16]:LunEclSun=100.25640 	LEcl_ChSun(食甚入盈縮曆)=100.25639857 	LEclUseSun=82.36468043(1)
				[LunStep 17]:LEcl_SunX(食甚盈縮差)=2.38587969
				[LunStep 18]:LEcl_SunX_=-2.38587969 	LEclNetSun(食甚入盈縮曆定度)=97.87051887(1)
						LEclNetSun_= LEclNetSun
						if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875}
							As LEclSunStatus=1	LEclNetSun_=105.89110287

						total_year=calculated_year-Base_year=342 => Added_drift=0.015*total_year=5.13000
						Ref_drift(天正加時赤道度)=10-Added_drift=4.87000(邢雲路計算方法)
						  (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.75304)
						Next_Ref_drift=Ref_drift_check-4.3445=0.40854

						Sun_drift=Next_Ref_drift/1.0849=0.37657
						Net_drift=4+Sun_drift=4.37657
						LEclNetSun_All=LEclNetSun_ + Net_drift=110.26767
						LEclNetSun_From(食甚月躔黃道奎宿)=LEclNetSun_All-103.6975(壁宿)=6.57017
				[LunStep 22]:
					求交常交定度:(refer to LunStep 4 and LunStep 8)
						置經朔望入交泛日及分秒(LunEcl),以月平行度(13.36878311)乘之,為交常度(LunEclA);
						以盈縮差(LSunXS)盈加縮減之,為交定度(LunEclA_S)。

					求月食入陰陽曆去交前後度:
						視交定度,在交中度(181.8967)已下,為陽曆;已上,減去交中,為陰曆。
						視入陰陽曆,在後準十五度半已下,為交後度;
						前準一百六十六度三十九分六十八秒已上,覆減交中,餘為交前度及分。

						交中度=181.8967
						SexLunEclA=LunEclA_S=355.477166424
						if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S}
						Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-173.580466424

						LunEclFront=-1(負數[-1]表示交前度);
						if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);}

						SexLunEclA_=SexLunEclA=-173.580466424
						if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1}
						if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1}
						Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -8.316233576(-1)
				[LunStep 23]:
					求月食分秒:
						視去交前後度,
						不用南北東西差者。用減食限,不及減者不食。
						餘如定法而一,為月食之分秒。
						LunFraction=(13.05-abs(SexLunEclA_))/8.7=0.544111083
				以新法推天啟三年癸亥歲九月十五壬寅夜望月食...月食五分三十八秒(0.538)
											理論值:0.6
				[LunStep 24]:
					LEclM=8.91282077(3)	LEclM_Used(定入遲疾曆)=8.91282077
					LEclM_x12_2(定入遲疾限)=108.73641341	LEclM_Ch=58.26358659	LEclM_Sp=1.05377118
					LEclM_SunSp(定限行度)=0.97177	    LHalfEclT(定用分)=0.05852605(邢雲路用 4920,不是 5740)
				[LunStep 25]:
					LEclBeginT(初虧)=0.717572	   LEclEndT(復圓)=0.834624
					LEclBeginT24=17.22171630h	LEclEndT24=20.03096682h
    以新法推天啟三年癸亥歲九月十五壬寅夜望月食[初虧申正三刻](16.72)不可靠, 因為這應該是授時的計算值
			     而授時則推月食[初虧酉初一刻](17.24)不可靠, 因為與計算值不符合,很可能是新法推的初虧
			 以新法推天啟三年..[初虧申正三刻](16.72h ?)	 ....  復圓[戍初初刻](19h ?)
						理論值:17.7h				     20.79h

----------------------------------------------------------
19 1629(1)				[2]MinSun=8.161296
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242152)127104.27(19.33)24.80	  149.66(54.53)10.96(24.32)	   179.19(24.06)12.93(26.64)	   (356.09)  179.19   3.43(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	495.99/1517.74(0.17)1(355.92)	 12.93/12.88(157.79)0.06(3)	( 9.21)10.47  1.10(1.10)0.9916	   -0.173/-1.101/-0.105(-2.529h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  1.3888h(22.8601h)10.86h 	  32.24m(23.3974h)  2069/7931	  2068.94	4748.92		   203.16/179.11(1) 3.52

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  495.99/1517.74(0.18)		  178.93(1)  3.69/4.45/(-5.77)	 178.93/0.35/0.67(0.04)			 1/ 5.77/-0.04

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  363.37/355.92/-7.46		   -7.46(8)0.0681		 12.83(3)10.50(0.99)(0.91)	   228.58(22.8488h)(23.9460h)

----------------------------------------------
20 1629(7)				[2]MinSun=8.161296
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242152)127104.27(19.33)24.80	  149.66(54.53)10.96(24.32)	   356.37(21.24)24.79(13.33)	   (178.26)  173.75   8.87(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	448.90/3753.66(0.41)-1(178.67)	 11.01/10.99(134.37)0.02(1)	(32.63) 8.18  3.30(3.29)1.1781	   0.414/3.295/0.258(6.196h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  5.7942h(11.9906h)-0.01h 	  -0.29m(11.9857h)  2073/7927	  2927.01	 5.98		   377.87/356.63(3) 8.61

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  448.90/3753.66(0.40)		  357.03(3)  8.21/4.42/(4.41)	 174.41/0.77/0.00(0.00)			 -1/-4.41/ 0.00

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  183.64/178.67/-4.96		    4.96(6)0.1726		 25.05(1)29.48(1.18)(1.10)	   293.15(11.2821h)(12.6892h)

----------------------------------------------
20 1631(6)				[2]MinSun=8.208200
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242150)127834.75(29.81)17.02	  157.39(12.79) 5.25( 0.64)	   334.58( 9.98)17.11(14.55)	   (194.52)  151.95  30.67(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	344.71/12549.90(1.28)-1(195.80)	  3.33/ 3.28(40.68)0.05(0)	(40.68) 7.14  3.84(3.84)1.0249	   1.278/3.836/0.409(9.820h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  23.4157h(9.2359h)-2.76h 	  -66.48m(8.1279h)  2130/7870	  2870.34	1613.37		   344.91/334.94(3)30.30

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  344.71/12549.90(1.27)		  336.20(3)  29.04/4.01/(1.76)	 153.58/2.38/1.54(1.54)			 -1/-1.76/ 1.54

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  187.83/195.80/ 7.96		   -7.96(8)0.0044		 17.52(0)45.67(1.03)(0.95)	   56.68(7.9919h)(8.2639h)

----------------------------------------------
21 1632(11)				[2]MinSun=8.231652
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242149)128199.99(35.05)27.90	  146.49( 7.16) 1.41( 1.25)	   106.09(32.00)23.15(26.75)	   (357.62)  106.09  76.53(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	82.51/23442.53(2.35)1(355.27)	  9.38/ 9.35(114.39)0.03(1)	(52.61) 5.37  4.65(4.65)1.1500	   -2.349/4.648/0.164(3.934h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  23.8977h(3.8321h)-8.17h 	  -81.51m(2.4736h)  2383/7617	  2383.35	3969.35		   138.19/106.20(1)76.43

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  82.51/23442.53(2.35)		  103.85(1)  78.77/1.14/(-0.76)	 103.85/4.37/6.95(1.80)			 -1/ 0.76/ 1.80

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  360.20/355.27/-4.93		   -4.93(8)0.3835		 23.32(1)50.61(1.15)(1.07)	   421.92(1.4610h)(3.4862h)

----------------------------------------------
22 1634(4)				[2]MinSun=8.278556
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242147)128930.48(45.54)20.11	  154.23(25.42)23.27( 4.78)	   272.35(23.55) 3.62(14.06)	   (187.92)  89.73  92.89(3)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	 8.90/24001.62(2.40)-1(190.32)	  3.62/ 3.61(44.12)0.01(2)	(44.12) 6.66  4.08(4.08)1.1629	   2.401/-4.076/-0.118(-2.834h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  13.1138h(10.2797h)-1.72h 	  -46.05m(9.5121h)  2489/7511	  2510.86	1036.62		   295.75/272.20(3)93.04

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	   2.98/24010.53(2.40)		  274.60(3)  90.64/0.07/(0.04)	 91.98/4.46/1.85(1.85)			 -1/-0.04/ 1.85

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  189.86/190.32/ 0.46		   -0.46(8)0.9424		  3.50(2)42.68(1.16)(1.08)	   529.16(8.2421h)(10.7821h)

----------------------------------------------
23 1636(8)				[2]MinSun=8.325460
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242145)129660.96(56.02)12.33	  161.96(43.69)17.56( 8.32)	   32.97(39.93) 5.82(26.86)	   (359.15)  32.97  32.97(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	334.85/13234.41(1.36)1(357.79)	  5.82/ 5.74(70.98)0.08(2)	(70.98) 2.11  5.31(5.31)1.1174	   -1.356/-5.308/-0.489(-11.737h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  22.4066h(10.6698h)-1.33h 	  -36.96m(10.0537h)  2138/7862	  2862.47	810.95		   72.38/32.45(0)32.45

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  334.85/13234.41(1.34)		  31.11(0)  31.11/3.94/(2.83)	 31.11/2.52/0.82(0.82)			 -1/ 2.83/ 0.82

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  361.28/357.79/-3.49		   -3.49(8)0.5639		  5.33(2)65.00(1.13)(1.05)	   493.47(8.8694h)(11.2380h)

----------------------------------------------
25 1638(1)				[2]MinSun=8.372364
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242143)130391.45( 6.51) 4.55	  169.70( 1.95)11.86(11.85)	   199.23(31.48)13.84(14.17)	   (189.45)  16.61  16.61(2)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	429.61/7570.66(0.78)-1(190.23)	  0.06/ 0.00( 0.72)0.06(0)	( 0.72)11.04  0.08(0.08)0.9859	   0.783/0.080/0.072(1.723h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  11.6227h(13.3460h) 1.35h 	  37.34m(13.9684h)  2087/7913	  2086.71	820.16		   230.81/199.33(2)16.70

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  429.61/7570.66(0.79)		  200.11(2)  17.49/4.30/(2.61)	 17.49/1.54/0.51(0.51)			 1/ 2.61/-0.51

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  190.15/190.23/ 0.08		   -0.08(8)0.9898		 13.91(0) 1.59(0.99)(0.90)	   634.62(12.4453h)(15.4915h)

----------------------------------------------
26 1641(11)				[2]MinSun=8.442720
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242140)131487.17(22.23) 7.64	  166.53(14.59) 2.32(16.00)	   126.13(39.42)24.06(14.29)	   (190.98)  126.13  56.49(1)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	205.65/20486.95(2.06)-1(188.92)	 10.28/10.25(125.39)0.03(1)	(41.61) 7.01  3.97(3.97)1.1664	   -2.059/3.972/0.134(3.228h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  10.1316h(13.3594h) 1.36h 	  37.67m(13.9872h)  2258/7742	  2258.16	827.98		   165.71/126.29(1)56.33

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  205.65/20486.95(2.06)		  124.23(1)  58.39/2.64/(1.67)	 124.23/3.88/1.28(1.28)			 1/ 1.67/ 1.28

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  191.00/188.92/-2.08		    2.08(6)0.6526		 24.19(1)40.00(1.17)(1.09)	   495.34(12.7983h)(15.1760h)

----------------------------------------------
27 1643(4)				[2]MinSun=8.489624
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242138)132217.65(32.71)29.39	  144.74( 3.32)22.19(17.21)	   262.86( 1.44) 2.54(26.49)	   (354.08)  80.24  80.24(2)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	56.99/23734.40(2.37)1(356.45)	  2.54/ 2.54(31.02)0.00(2)	(31.02) 8.37  3.08(3.08)1.1800	   2.375/-3.079/-0.049(-1.174h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  10.6136h(9.4396h)-2.56h 	  -62.94m(8.3906h)  2410/7590	  2409.68	1503.91		   264.21/262.77(2)80.15

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  56.99/23734.40(2.37)		  265.14(2)  82.52/0.82/(0.31)	 82.52/4.42/2.66(2.66)			 -1/-0.31/-2.66

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  354.68/356.45/ 1.78		    1.78(6)0.7035		  2.49(2)30.42(1.18)(1.10)	   498.95(7.1931h)(9.5881h)

----------------------------------------------
27 1644(3)				[2]MinSun=8.513076
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242137)132582.90(37.96)10.74	  163.37(27.22)20.33(20.14)	   251.96(55.81)26.26(27.09)	   (362.21)  69.34  69.34(2)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	126.46/22688.64(2.27)1(364.48)	 12.48/12.47(152.26)0.01(1)	(14.74)10.03  1.67(1.67)1.1965	   2.273/1.667/0.270(6.480h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  19.4244h(1.9045h)-10.10h 	  -50.07m(1.0700h)  2338/7662	  2338.01	4554.18		   308.01/252.20(2)69.58

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  126.46/22688.64(2.28)		  254.47(2)  71.85/1.70/(-1.61)	 71.85/4.26/7.75(0.76)			 -1/ 1.61/-0.76

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  358.49/364.48/ 5.99		    5.99(6)0.0014		 26.53(1)11.45(1.20)(1.12)	   27.23(1.0046h)(1.1353h)

----------------------------------------------
28 1644(9)				[2]MinSun=8.513076
	[1](365.)SumDay(WinSol)LeapDay	[2]RefS(RefMDay)RefM(RefEcl)	[3]ReqS(ReqMDay)ReqM(ReqEcl)	[4](EclA)[5]ChSun[6]UseSun(S)
	(0.242137)132582.90(37.96)10.74	  163.37(27.22)20.33(20.14)	   63.90(52.99)10.56(13.79)	   (184.38)  63.90  63.90(0)

	[7]Sun_Add/Sum(X)[8]Ft(EclA_S)	[9]UseM/SoD(M12_2)MSo(M)	(MoonCh)MAdd[10]Sum(X)Speed	[11]SunXS/MoonXS/SunMoonT
	173.73/21238.16(2.14)-1(182.24)	 10.56/10.50(128.83)0.06(3)	(38.17) 7.48  3.73(3.72)1.0215	   -2.140/-3.725/-0.471(-11.297h)

	[12]Conjunc0(T)  [13]TNoon	SeeT(MaxPhase) [14]Srise/Sset	HalfDayFan[15]NetTNoon		[16]Ecl)Ref/Sun(S)UseSun
	  23.8298h(12.5325h) 0.53h 	  15.90m(12.7976h)  2300/7700	  2699.72	332.32		   116.44/63.44(0)63.44

	[17]Ecl_)SAdd/SSum(SunX)	[18]NSun(S)[19]NS_R/NS_X(NS_)	[20]EW_R/EW_X/EW_H(EW_)		[21]A12_Sgn/NS_Sh/EW_Sh
	  173.73/21238.16(2.13)		  61.31(0)  61.31/2.45/(2.15)	 61.31/3.98/0.53(0.53)			 1/-2.15/ 0.53

	[22]EclRefA/NetEclA_S/PosA	[23]EclPos(Ref)Frac		[24]EclM(M)_Ch(_Sp)(_SunSp)	[25]HalfEclT(EclBegin)(EclEnd)
	  186.43/182.24/-4.19		    4.19(6)0.3014		 10.09(3)43.91(1.03)(0.95)	   433.50(11.7572h)(13.8380h)

----------------------------------------------
Reference:
(1)Accuracy of Solar Eclipse Calculations, 1644-1785, with the Chinese Shoushi Calandar(Table 2)
http://articles.adsabs.harvard.edu//full/1999JHA....30..161L/0000161.000.html

(2)Accuracy of Solar Eclipse Observations made by Jesuit Astronomers in China
http://adsabs.harvard.edu/full/1995JHA....26..227S
Observation data from Table 2
Observation error data from Table 3
Theorectical=Observation-Observation error

	 Shoushi(Observation)		Remark : 999 or 9999 = no data
	-)Theorectical(Obsevation-Shoushi)
	--------------
Shoushi-Theorectical(Observation-Theorectical)
==================================================================================================================
N  year(month)	  EclipseFraction	EclipseBegin		MaxPhase		EclipseEnd
      month/day
==================================================================================================================
 0 1568(10) 	   0.0221(0.110)	12.1030(11.137)   	12.3605(11.737)   	12.6180(12.343)
       9/21	-) 0.110(0.0878)   	11.137(-0.9663)   	11.737(-0.6239)   	12.343(-0.2752)
		----------------------------------------------------------------------------------------
		   -0.0878(0.000)   	0.9663(0.000)   	0.6239(0.000)   	0.2752(0.000)

 1 1572(8) 	   0.9668(0.916)	6.4686(6.660)   	7.9429(7.930)   	9.4172(9.930)
       6/10	-) 0.916(-0.0508)   	6.660(0.1914)   	7.930(-0.0129)   	9.930(0.5128)
		----------------------------------------------------------------------------------------
		   0.0508(0.000)   	-0.1914(0.000)   	0.0129(0.000)   	-0.5128(0.000)

 2 1575(6) 	   0.4650(0.701)	13.5750(13.050)   	14.6262(14.290)   	15.6774(15.450)
       5/10	-) 0.701(0.2360)   	13.050(-0.5250)   	14.290(-0.3362)   	15.450(-0.2274)
		----------------------------------------------------------------------------------------
		   -0.2360(0.000)   	0.5250(0.000)   	0.3362(0.000)   	0.2274(0.000)

 3 1582(7) 	   0.9759(0.723)	12.7915(12.206)   	14.1250(13.618)   	15.4585(14.915)
       6/20	-) 0.723(-0.2528)   	12.206(-0.5860)   	13.618(-0.5071)   	14.915(-0.5431)
		----------------------------------------------------------------------------------------
		   0.2528(0.000)   	0.5860(0.000)   	0.5071(0.000)   	0.5431(0.000)

 4 1583(13) 	   0.7154(0.825)	11.3255(11.230)   	12.7799(13.080)   	14.2343(14.770)
      12/14	-) 0.825(0.1096)   	11.230(-0.0955)   	13.080(0.3001)   	14.770(0.5357)
		----------------------------------------------------------------------------------------
		   -0.1096(0.000)   	0.0955(0.000)   	-0.3001(0.000)   	-0.5357(0.000)

 5 1587(10) 	   0.8537(0.736)	11.9223(11.322)   	13.1681(12.644)   	14.4139(13.934)
      10/2	-) 0.736(-0.1179)   	11.322(-0.6000)   	12.644(-0.5245)   	13.934(-0.4796)
		----------------------------------------------------------------------------------------
		   0.1179(0.000)   	0.6000(0.000)   	0.5245(0.000)   	0.4796(0.000)

 6 1589(2) 	   0.5157(0.301)	6.1622(6.450)   	7.2658(7.196)   	8.3695(7.999)
       2/15	-) 0.301(-0.2144)   	6.450(0.2877)   	7.196(-0.0697)   	7.999(-0.3703)
		----------------------------------------------------------------------------------------
		   0.2144(0.000)   	-0.2877(0.000)   	0.0697(0.000)   	0.3703(0.000)

 7 1590(8) 	   0.2223(0.431)	16.1850(16.080)   	17.1090(17.106)   	18.0330(18.037)
       7/31	-) 0.431(0.2086)   	16.080(-0.1047)   	17.106(-0.0031)   	18.037(0.0042)
		----------------------------------------------------------------------------------------
		   -0.2086(0.000)   	0.1047(0.000)   	0.0031(0.000)   	-0.0042(0.000)

 8 1594(6) 	   0.5152(0.319)	10.1823(9.830)   	11.3417(10.840)   	12.5011(11.900)
       5/20	-) 0.319(-0.1962)   	9.830(-0.3523)   	10.840(-0.5017)   	11.900(-0.6011)
		----------------------------------------------------------------------------------------
		   0.1962(0.000)   	0.3523(0.000)   	0.5017(0.000)   	0.6011(0.000)

 9 1596(10) 	   0.7084(0.949)	10.6131(9.870)   	11.8365(11.170)   	13.0598(12.500)
       9/22	-) 0.949(0.2406)   	9.870(-0.7431)   	11.170(-0.6665)   	12.500(-0.5598)
		----------------------------------------------------------------------------------------
		   -0.2406(0.000)   	0.7431(0.000)   	0.6665(0.000)   	0.5598(0.000)

10 1603(5) 	   0.7343(0.903)	7.5081(7.400)   	8.8188(8.590)   	10.1294(9.930)
       5/11	-) 0.903(0.1687)   	7.400(-0.1081)   	8.590(-0.2288)   	9.930(-0.1994)
		----------------------------------------------------------------------------------------
		   -0.1687(0.000)   	0.1081(0.000)   	0.2288(0.000)   	0.1994(0.000)

11 1604(5) 	   0.6227(0.078)	15.5666(15.883)   	16.9268(16.408)   	18.2871(16.907)
       4/29	-) 0.078(-0.5446)   	15.883(0.3161)   	16.408(-0.5189)   	16.907(-1.3801)
		----------------------------------------------------------------------------------------
		   0.5446(0.000)   	-0.3161(0.000)   	0.5189(0.000)   	1.3801(0.000)

12 1607(3) 	   0.4244(999.000)	17.4368(17.320)   	18.4662(18.260)   	19.4956(999.000)
       2/26	-) 999.000(998.5756)   	17.320(-0.1168)   	18.260(-0.2062)   	999.000(979.5044)
		----------------------------------------------------------------------------------------
		   -998.5756(0.000)   	0.1168(0.000)   	0.2062(0.000)   	-979.5044(0.000)

13 1609(0) 	   -99.0000(0.044)	0.0000(16.197)   	0.0000(16.594)   	0.0000(16.975)
      12/26	-) 0.044(99.0444)   	16.197(16.1970)   	16.594(16.5940)   	16.975(16.9753)
		----------------------------------------------------------------------------------------
		   -99.0444(0.000)   	-16.1970(0.000)   	-16.5940(0.000)   	-16.9753(0.000)

14 1610(13) 	   0.9588(0.607)	14.1179(14.610)   	15.6318(16.070)   	17.1457(17.320)
      12/15	-) 0.607(-0.3518)   	14.610(0.4921)   	16.070(0.4382)   	17.320(0.1743)
		----------------------------------------------------------------------------------------
		   0.3518(0.000)   	-0.4921(0.000)   	-0.4382(0.000)   	-0.1743(0.000)

15 1612(0) 	   -99.0000(999.000)	0.0000(999.000)   	0.0000(999.000)   	0.0000(999.000)
      12/99	-) 999.000(1098.0000)   	999.000(999.0000)   	999.000(999.0000)   	999.000(999.0000)
		----------------------------------------------------------------------------------------
		   -1098.0000(0.000)   	-999.0000(0.000)   	-999.0000(0.000)   	-999.0000(0.000)

16 1615(4) 	   0.6471(0.890)	15.0718(15.162)   	16.3547(16.412)   	17.6376(17.546)
       3/29	-) 0.890(0.2428)   	15.162(0.0906)   	16.412(0.0571)   	17.546(-0.0913)
		----------------------------------------------------------------------------------------
		   -0.2428(0.000)   	-0.0906(0.000)   	-0.0571(0.000)   	0.0913(0.000)

17 1617(8) 	   0.3689(0.470)	19.4230(19.070)   	20.5355(19.904)   	21.6481(20.691)
       8/1	-) 0.470(0.1011)   	19.070(-0.3530)   	19.904(-0.6314)   	20.691(-0.9572)
		----------------------------------------------------------------------------------------
		   -0.1011(0.000)   	0.3530(0.000)   	0.6314(0.000)   	0.9572(0.000)

18 1621(6) 	   0.2680(0.800)	17.5976(17.240)   	18.5196(9999.000)   	19.4415(9999.000)
       5/21	-) 0.711(0.5320)   	17.270(-0.3576)   	18.230(9980.4804)   	19.190(9979.5585)
		----------------------------------------------------------------------------------------
		   -0.4430(0.089)   	0.3276(-0.030)   	0.2896(9980.770)   	0.2515(9979.810)

19 1629(7) 	   0.1726(0.166)	11.2821(10.980)   	11.9857(11.710)   	12.6892(12.450)
       6/21	-) 0.166(-0.0066)   	10.980(-0.3021)   	11.710(-0.2757)   	12.450(-0.2392)
		----------------------------------------------------------------------------------------
		   0.0066(0.000)   	0.3021(0.000)   	0.2757(0.000)   	0.2392(0.000)

20 1631(0) 	   -99.0000(0.108)	0.0000(12.340)   	0.0000(12.910)   	0.0000(13.480)
      10/25	-) 0.108(99.1080)   	12.340(12.3400)   	12.910(12.9100)   	13.480(13.4800)
		----------------------------------------------------------------------------------------
		   -99.1080(0.000)   	-12.3400(0.000)   	-12.9100(0.000)   	-13.4800(0.000)

21 1634(4) 	   0.9424(0.820)	8.2421(8.770)   	9.5121(9.944)   	10.7821(11.190)
       3/29	-) 0.820(-0.1228)   	8.770(0.5281)   	9.944(0.4315)   	11.190(0.4078)
		----------------------------------------------------------------------------------------
		   0.1228(0.000)   	-0.5281(0.000)   	-0.4315(0.000)   	-0.4078(0.000)

22 1636(8) 	   0.5639(0.633)	8.8694(8.380)   	10.0537(9.627)   	11.2380(10.964)
       8/1	-) 0.633(0.0689)   	8.380(-0.4898)   	9.627(-0.4272)   	10.964(-0.2737)
		----------------------------------------------------------------------------------------
		   -0.0689(0.000)   	0.4898(0.000)   	0.4272(0.000)   	0.2737(0.000)

23 1637(0) 	   -99.0000(0.199)	0.0000(12.710)   	0.0000(13.680)   	0.0000(14.580)
       1/26	-) 0.199(99.1990)   	12.710(12.7100)   	13.680(13.6800)   	14.580(14.5800)
		----------------------------------------------------------------------------------------
		   -99.1990(0.000)   	-12.7100(0.000)   	-13.6800(0.000)   	-14.5800(0.000)

24 1638(1) 	   0.9898(0.656)	12.4453(11.650)   	13.9684(13.320)   	15.4915(14.840)
       1/15	-) 0.656(-0.3338)   	11.650(-0.7953)   	13.320(-0.6484)   	14.840(-0.6515)
		----------------------------------------------------------------------------------------
		   0.3338(0.000)   	0.7953(0.000)   	0.6484(0.000)   	0.6515(0.000)

25 1641(11) 	   0.6526(0.827)	12.7983(12.650)   	13.9872(14.040)   	15.1760(15.330)
      11/3	-) 0.827(0.1744)   	12.650(-0.1483)   	14.040(0.0528)   	15.330(0.1540)
		----------------------------------------------------------------------------------------
		   -0.1744(0.000)   	0.1483(0.000)   	-0.0528(0.000)   	-0.1540(0.000)

26 1643(4) 	   0.7035(0.530)	7.1931(7.830)   	8.3906(8.850)   	9.5881(9.940)
       3/20	-) 0.530(-0.1735)   	7.830(0.6369)   	8.850(0.4594)   	9.940(0.3519)
		----------------------------------------------------------------------------------------
		   0.1735(0.000)   	-0.6369(0.000)   	-0.4594(0.000)   	-0.3519(0.000)

27 1644(9) 	   0.3014(0.303)	11.7572(10.930)   	12.7976(12.140)   	13.8380(13.340)
       9/1	-) 0.303(0.0016)   	10.930(-0.8272)   	12.140(-0.6576)   	13.340(-0.4980)
		----------------------------------------------------------------------------------------
		   -0.0016(0.000)   	0.8272(0.000)   	0.6576(0.000)   	0.4980(0.000)