日の出と日の入りの時刻を求めるコードをPythonで作成した。 計算方法は、参考文献（下の本）で紹介されている逐次近似法による反復計算。このコードを実行すると、2016年（変数yearで指定）の毎月3日（変数dayで指定）における東京での（緯度latitudeと経度…