不知原始数据有问题，还是我的程序没编好，呵呵

依原始资料如下：

由儒略日推算历日

将JD加上0.5，令 Z 为其整数部分，F 为尾数（小数）部分。

若 Z < 2299161，取A = Z

若 Z 大于等于2299 161，计算 α=INT((Z-1867216.25)/36524.25) ，A=Z+1+α-INT(α/4)

然后计算

B = A+1524

C = INT((B-122.1)/365.25)

D = INT(365.25C)

E = INT((B-D)/30.6001)

该月日期（带小数部分）则为：d = B - D - INT(30.6001E) + F

月份m为：
IF E < 14 THEN m = E – 1
IF E=14 or E=15 THEN m = E – 13

年份为y：
IF m>2 THEN y = C – 4716 IF m =1 or m=2 THEN y = C – 4715

这个公式里求E时用的数30.6001不能代之以30.6，哪怕计算机没有先前所说的问题。否则，你得到的结果会是2月0日而不是1月31日，或者4月0日而不是3月31日。