6/1

lat<-34.147643*pi/180
long<- -131.468115*pi/180
year<-55347 #2010年6月に対応
day<-26 #observe day(UT)
MJD<-year+day+hour/24
D<-MJD-40000
a1<- 0.67239+1.00273781*D
a2<-a1-floor(a1)
thitaG<-2*pi*(a2) #2000年分点でのグリニッジ恒星時
thita<-thitaG-long
thitah<-(thitaG-long)*12/pi
(thitah-4.7)/1.002738
http://milkyway.sci.kagoshima-u.ac.jp/~imai/J-Net-guide.html#%8A%CF%91%AA%83X%83P%83W%83%85%81[%83%8B%82%CC%93]%91%97
http://milkyway.sci.kagoshima-u.ac.jp/users/kameno/weblog/fa7dd/Web.html