昼长

已知日期和一个地方的纬度，可以用下列公式求昼长：
D=2/15°arc cos(-tgδtgφ)
式中D是昼长，δ是太阳赤纬（可以根据日期查得），φ是地理纬度
根据日期查太阳赤纬的办法，可以查天文年历，也可以用内插法查得近似值。
下面是我用内插法求得的近似值（其中的太阳直射点位置实际上就是太阳赤纬，太阳直射点日移动纬度数也就是太阳赤纬日变化数），可供参考：