反正古人定义的那四个星曜与现代的轨道理论不完全相同，建议全部用平轨道要素来计算吧。升交点平黄经、近点平黄经在《天文算法》中有现面的公式！降交点平黄经与长交点平黄经相差180度。近点平黄经与远点平黄经也相差180度。

F=W₁-W₃

l=W₁-W₂

l'=T-ω'

月球及月地参数：

W₁：月球平黄经

W₂：白道近地平黄经

W₃：白道升交点平黄经

T ：地月质心的日心平黄经

ω'：日月质心的近日平黄经

W₁、W₂及W₃是J2000之后的平黄道角度，详见5.1节。T和ω’是J2000.0起算的角度。这些参数随时间的变化规律可用多项式表示：A=a⁽⁰⁾+a⁽¹⁾t+a⁽²⁾t²+a⁽³⁾t³+a⁽⁴⁾t⁴。这些参数在修正前可由表1得到，表示日月运动的参数(系数W₁和T)也表示为ν和n’：v=W⁽¹⁾ ，n’=T⁽¹⁾。

W1 = 218º18'59.95571" +1732559343.73604"*t -6.8084"  *t2 +0.006604"  *t3 -0.00003169" *t4

W2 =  83º21'11.67475" +14643420.3171"   *t -38.2631" *t2 -0.045047"  *t3 +0.00021301" *t4

W3 = 125º02'40.39816" -6967919.5383"    *t +6.3590"  *t2 +0.007625"  *t3 -0.00003586" *t4

T  = 100º27'59.13885" +129597742.2930"  *t -0.0202"  *t2 +0.000009"  *t3 +0.00000015" *t4
ω' = 102º56'14.45766" +1161.24342"      *t +0.529265"*t2 -0.00011814"*t3 +0.000011379"*t4

t是J2000起算的地心时间TDB（力学时），单是儒略世纪

以上计算是超高精度计算，但不含与DE405/406拟合的修正量

已发布V4.01版

更新内容：

1、根据“逍遥子”的建议，修改了建国以来及民国时期的年号

2、改进了角度差j1_j2()函数，使误差控制在1毫角秒以内，原来的误差约0.05角秒。通过这个改进，可使食界验算时误差不被觉察。

3、改进视差修正算法。原来使用《天文算法》一书中的公式，现在直接转换到直角坐标中求解，简单、严格。

4、采用了《2008年中国天文年历》中的一些常数：地球赤道半径、1天文单位的长度

cs_rEar = 6378.1366; //地球赤道半径(千米)

cs_AU   = 1.49597870691e8; //天文单位长度(千米)

辛苦，程序是越做越完善了

谢谢

关于闰正月初一是不是该标注为"春节"?

许兄的万年历计算出来的结果,2262年2月20日为闰正月初一,标为"春节".按说一年只有一个春节,把闰正月初一标为"春节"似乎不一定合式.不知道到时候的政府部门是如何规定的?如果闰正月初一也放假的话,那么标为"春节"也没什么的.

2262年2月20日
星期四 双鱼座
4960年 闰正月小 初一日
壬寅年 壬寅月 壬寅日
回历[1690年11月30日]
朔 02:01:01 春节 七九第7天

3359年1月20日对应农历十二月小廿九日,第二天是初一,由于还有一个闰十二月,所以这一天没标注为"除夕",看来许兄是考虑过一年只有一个除夕的.那么一年也应当只有一个"春节"的.

建议许兄把2262年2月20日的闰正月初一的"春节"标注删除掉.闰正月初二的"大年初二"也删除掉.

3359年1月20日
星期六 水瓶座
6056年 十二月小 廿九日
戊午年 乙丑月 壬戌日
回历[2821年6月28日]
定大寒 20:34:32
大寒 四九第3天

3359年2月18日
星期日 水瓶座
6056年 闰十二月小 廿九日
戊午年 丙寅月 辛卯日
回历[2821年7月28日]
除夕 七九第5天