日心赤道直交座標

 日心赤道直交座標(x', y', z')を日心黄道直交座標(x , y , z)と黄道傾斜角εから求める式は、次のようになります。

   x'=x
   y'=y cos ε−z sin ε
   z'=y sin ε+z cos ε

 時間引数Tから軌道傾斜角εを求める式は、次のようになります。

   ε=23+26/60+21.448/3600
    −46.8150/3600×T−0.59e-3/3600×pow(T ,2)+1.8137e-3/3600×pow(T ,3)

   T=(JD−2451545)/36525 (J2000.0)