というわけで、このプログラムの動作原理の解説をしておきます。 このプログラムでは標準時との差を示す経度（U=9ならば東経135度）の指定した日付における南中時刻と南中高度を求め実際の計測値との差分によって現在地を求めています。 南中高度＝90度－観測場所の緯度＋太陽の視赤緯 太陽の視赤緯は下記の式で求めることができます。 0.006918-0.399912*cos(khi)+0.070257*sin(khi)-0.006758*cos(2*khi)+0.000907*sin(2*khi)-0.002697*cos(3*khi)+0.001480*sin(3*khi) 　※khi=元旦からの経過日数*2π/365