def compute_lunisolar
      dp, de = 0.0, 0.0

      begin
        l  = compute_l_iers2003(@jc)
        lp = compute_lp_mhb2000(@jc)
        f  = compute_f_iers2003(@jc)
        d  = compute_d_mhb2000(@jc)
        om = compute_om_iers2003(@jc)
        Const::NUT_LS.map do |x|
          x[0, 5] + x[5..-1].map { |y| y.to_s.sub(/\./, "").to_i }
        end.reverse.each do |x|
          arg = (x[0] * l + x[1] * lp + x[2] * f +
                 x[3] * d + x[4] * om) % Const::PI2
          sarg, carg = Math.sin(arg), Math.cos(arg)
          dp += (x[5] + x[6] * @jc) * sarg + x[ 7] * carg
          de += (x[8] + x[9] * @jc) * carg + x[10] * sarg
        end
        return [dp * Const::U2R, de * Const::U2R]
      rescue => e
        raise
      end
    end
def compute_lunisolar
      dp, de = 0.0, 0.0

      begin
        l  = compute_l_iers2003(@jc)
        lp = compute_lp_mhb2000(@jc)
        f  = compute_f_iers2003(@jc)
        d  = compute_d_mhb2000(@jc)
        om = compute_om_iers2003(@jc)
        Const::NUT_LS.map do |x|
          x[0, 5] + x[5..-1].map { |y| y.to_s.sub(/\./, "").to_i }
        end.reverse.each do |x|
          arg = (x[0] * l + x[1] * lp + x[2] * f +
                 x[3] * d + x[4] * om) % Const::PI2
          sarg, carg = Math.sin(arg), Math.cos(arg)
          dp += (x[5] + x[6] * @jc) * sarg + x[ 7] * carg
          de += (x[8] + x[9] * @jc) * carg + x[10] * sarg
        end
        return [dp * Const::U2R, de * Const::U2R]
      rescue => e
        raise
      end
    end
=========================================================================
 日月章動(luni-solar nutation)の計算

 @param:  <none>
 @return: [dpsi_ls, deps_ls]
=========================================================================