在一網站中見到“黃赤大距”算法，全文如下：

黃赤交角（黃赤大距）算法

黃赤交角精確的算法，與太陽到山，日影測向等有著密切的關係，下列程式算法給需要的人做參考。

黃赤交角的常數是23度26分21.448秒（黃赤大距）若欲求得精確之數，可用下列公式計算：

ε黃赤交角＝23°27′8〞.26－46〞.845Ｔ－0〞.0059Ｔ2＋0〞.181Ｔ3

Ｔ為1900年1月0.5日後的儒略周日紀，如1991年12月22日的儒略周日數是2448613。

2448613－2415020＝33593（日）
33593÷36525＝0.919835紀（Ｔ）

則⊿ε（黃赤交角）＝46.845－0.059Ｔ2＋0.0018Ｔ3＝43.093366弧秒
43.093366弧秒÷3600＝0.011970（度）
23°27′8〞.26－0.011970（度）＝23°26′25〞.1

黃赤交角的數值每年都在變化。例如2002 年是23.4390（23°26′20.4〞） 度，到了2010 年是23.4380（23°26′16.8〞） 度。總的來說，其變化是週期性的。約37158 年為一週期，變化幅度約為1.612 度。2008 年到達中間點。

仔細算算，發現稍有差異，疑惑不解，請各位指點：
1、33593÷36525＝0.919835紀（Ｔ）按计算应等于：0.91972621492有点偏差。
2、⊿ε（黃赤交角）＝46.845－0.059Ｔ2＋0.0018Ｔ3＝43.093366弧秒 所得数值有差。
按上面所提供公式应为：
ε黃赤交角＝23°27′8〞.26－46〞.845Ｔ－0〞.0059Ｔ2＋0〞.181Ｔ3
⊿ε（黃赤交角）＝46.845×0.919835-0.059×0.919835×2+0.0018×0.919835×3=42.986097154弧秒 所得数值也稍有不同。

請教各位大大，以上所提供的計算方法正確嗎？

>ε黃赤交角＝23°27′8〞.26－46〞.845Ｔ－0〞.0059Ｔ2＋0〞.181Ｔ3
上方的公式有問題,T2,T3非乘法是指次方, 0.181也不對.

公式:
黃赤交角=23.452294-(46.845*T+0.0059*T*T-0.00181*T*T*T)/3600

答案:
T=0.919726214921287
黃赤交角=23.440325067(度)

謝謝天步大大指點！