【高手进！】日出日落时间【追加20分】

http://www.time.ac.cn/serve/sunriseset/default.asp
国家授时中心网站上有……应该很权威

网上查到一个算法……很麻烦，写程序算，根据的公式：

已知：日出日落时太阳的位置h＝-0.833°,要计算地的地理位置,经度Long,纬度Glat,时区zone,UTo为上次计算的日出日落时间,第一次计算时UTo＝180°。

(1)先计算出从格林威治时间公元2000年1月1日到计算日天数days;

(2)计算从格林威治时间公元2000年1月1日到计算日的世纪数t,

则 t＝(days+UTo／360)／36525;

(3)计算太阳的平黄径 L=280.460+36000.770×t;

(4)计算太阳的平近点角

G＝357.528+35999.050×t

(5)计算太阳的黄道经度

λ＝L+1.915×sinG+0.020xsin(2G);

(6)计算地球的倾角 ε＝23.4393-0.0130×t;

(7)计算太阳的偏差 δ＝arcsin(sinε×sinλ);

(8)计算格林威治时间的太阳时间角GHA：

GHA=UTo-180-1.915×sinG-0.020×sin(2G) +2.466×sin(2λ)-0.053×sin(4λ)

(9)计算修正值e：

e=arcos{[ sinh-sin(Glat)sin(δ)]/cos(Glat)cos(δ)}

(10)计算新的日出日落时间

UT＝UTo-(GHA+Long±e);

其中“+”表示计算日出时间,“-”表示计算