Julian Day

Rの練習。
JDへの計算。

-----------
Y <- 2010                     # year
M <- 5                          # month
D <- 4.8474                 # day

if (M <= 2)
     {
      Y <- Y - 1
      M <- M + 12
     }

a <- trunc(Y / 100)
b <- trunc(2 - a + (a / 4))

ja <- trunc(365.25 * (Y + 4716))
jb <- trunc(30.6001 * (M + 1))
jc <- (D + b - 1524.5)

jd <- ja + jb + jc
     options(digits=14)     # 14桁表示
jd                                      # 結果
-----------

2010年2月1日で計算すると
2455228.5
てな感じの結果が出ました。合っていそうです。
時間・分・秒を日に直して計算してみると(ex. 2010年5月4.8474日UT)
2455321.3474
と出でた。
概ね良さそうかな。