一部訂正です。 2017年10月28日の東京における南中時刻は11時25分、南中高度は41.2度なのでT=685にプログラムを書き直して (誤)南中高度は40.9を入力してみてください。    ↓ (正)南中高度は41.2を入力してみてください。