如何使用Delphi计算当前日期的日出和日落时间

Bri*_*ost 9 delphi delphi-xe2

给定特定纬度和经度并使用Delphi(我使用的是XE2),我需要确定当天的日出和日落时间(当地时间).我发现了这篇文章并使用了SysTools库,但结果与英国公布的表格不一致.

有没有人看过这个pascal算法?