時間引数Tの元期を２０００年1月1日１２時UT（J2000.0）として、 黄道傾角εを、 ε＝23＋26/60＋21.448/3600−46.8150/3600×T −0.59e-3/3600×pow(T ,2)＋1.8137e-3/3600×pow(T ,3) で求める。 黄道傾角εを用いて、太陽の地心黄道座標（λ，β）、摂動補正を行っ…