def compute_obliquity(jc)
      return (84381.406        + \
             (  -46.836769     + \
             (   -0.0001831    + \
             (    0.00200340   + \
             (   -5.76 * 10e-7 + \
             (   -4.34 * 10e-8 ) \
             * jc) * jc) * jc) * jc) * jc) * Const::AS2R
    rescue => e
      raise
    end
def compute_obliquity(jc)
      return (84381.406        + \
             (  -46.836769     + \
             (   -0.0001831    + \
             (    0.00200340   + \
             (   -5.76 * 10e-7 + \
             (   -4.34 * 10e-8 ) \
             * jc) * jc) * jc) * jc) * jc) * Const::AS2R
    rescue => e
      raise
    end
=========================================================================
 黄道傾斜角計算

 * 黄道傾斜角 ε （単位: rad）を計算する。
   以下の計算式により求める。
     ε = 84381.406 - 46.836769 * T - 0.0001831 T^2 + 0.00200340 T^3
       - 5.76 * 10^(-7) * T^4 - 4.34 * 10^(-8) * T^5
   ここで、 T は J2000.0 からの経過時間を 36525 日単位で表したユリウス
   世紀数で、 T = (JD - 2451545) / 36525 である。

 @param:  jc   (ユリウス世紀)
 @return: eps  (平均黄道傾斜角)
=========================================================================