' ****************************************
' EoT Macro
Function EoT(The_Year, The_Month, The_Day, The_Hour, The_Minute, The_Second)
bbb = 367 * The_Year - 730531.5
ccc = Int((7# * Int(The_Year + (The_Month + 9) / 12)) / 4)
ddd = Int(275 * The_Month / 9) + The_Day
D2000 = bbb - ccc + ddd + (The_Hour + The_Minute / 60 + The_Second / 3600) / 24
Cycle = Int(D2000 / 365.25)
Theta = 0.0172024 * (D2000 - 365.25 * Cycle)
Average = 0.00526
Amp1 = 7.36303 - Cycle * 9e-05
Amp2 = 9.92465 - Cycle * 0.00014
Phi1 = 3.07892 + Cycle * -0.00019
Phi2 = -1.38995 + Cycle * 0.00013

EoT1 = Amp1 * Sin(1 * (Theta + Phi1))
EoT2 = Amp2 * Sin(2 * (Theta + Phi2))
EoT3 = 0.3173 * Sin(3 * (Theta - 0.94686))
EoT4 = 0.21922 * Sin(4 * (Theta - 0.60716))

EoT = Average + EoT1 + EoT2 + EoT3 + EoT4
End Function

' ****************************************
' Declination Macro
Function Decl(The_Year, The_Month, The_Day, The_Hour, The_Minute, The_Second)
bbb = 367 * The_Year - 730531.5
ccc = Int((7# * Int(The_Year + (The_Month + 9) / 12)) / 4)
ddd = Int(275 * The_Month / 9) + The_Day
D2000 = bbb - ccc + ddd + (The_Hour + The_Minute / 60 + The_Second / 3600) / 24
Cycle = Int(D2000 / 365.25)
Theta = 0.0172024 * (D2000 - 365.25 * Cycle)
Amp1 = 23.2639 - Cycle * 0.000131 + 0.0024 * Sin(Cycle * 0.335103 - 0.4)
Amp1 = 23.2639 - Cycle * 0.000131 + 0.0024 * Sin(Cycle * 0.335 - 0.4)
Phi1 = -1.38819 + Cycle * 0.000135
Decl1 = Amp1 * Sin(1 * (Theta + Phi1))
Decl2 = 0.380897 * Sin(2 * (Theta - 0.720483))
Decl3 = 0.171178 * Sin(3 * (Theta - 0.347175))
Decl4 = 0.008067 * Sin(4 * (Theta - 0.272216))
Decl = 0.37657 + Decl1 + Decl2 + Decl3 + Decl4
End Function

' ****************************************
' Right Ascension Function
Function RA(The_Year, The_Month, The_Day, The_Hour, The_Minute, The_Second)
bbb = 367 * The_Year - 730531.5
ccc = Int((7# * Int(The_Year + (The_Month + 9) / 12)) / 4)
ddd = Int(275 * The_Month / 9) + The_Day
D2000 = bbb - ccc + ddd + (The_Hour + The_Minute / 60 + The_Second / 3600) / 24
Cycle = Int(D2000 / 365.25)
Theta = 0.0172024 * (D2000 - 365.25 * Cycle)
Average = 0.00526
Amp1 = 7.36303 - Cycle * 9e-05
Amp2 = 9.92465 - Cycle * 0.00014
Phi1 = 3.07892 + Cycle * -0.00019
Phi2 = -1.38995 + Cycle * 0.00013

EoT1 = Amp1 * Sin(1 * (Theta + Phi1))
EoT2 = Amp2 * Sin(2 * (Theta + Phi2))
EoT3 = 0.3173 * Sin(3 * (Theta - 0.94686))
EoT4 = 0.21922 * Sin(4 * (Theta - 0.60716))

EOT_hrs = (Average + EoT1 + EoT2 + EoT3 + EoT4) / 60
RA = (18.6974 + 3.8198 * Theta + 24.00051 * Cycle - EOT_hrs)
RA = RA - (24 * (RA \ 24))
If RA < 0 Then RA = RA + 24
End Function