DETERMINING ALTITUDE
AND AZIMUTH
Disp "M,D,Y"
Prompt A,B,C
If A\<=\2
Then
A+12\->\A
C-1\->\C
End
19\->\D
\(-)\13\->\E
iPart((365.25*C)\->\F
iPart(((A+1)*30.6001)\->\G
E+F+G+B+1720994.5\->\X
X-2451545\->\S
S/36525\->\T
6.697374558+(2400.051336*T)+(0.000025862*T\^2\)\->\V
If V>24
Goto M
If V<0
Goto H
Lbl M
While V>24
V-24\->\V
End
Lbl H
While V<0
V+24\->\V
End
Disp "MST"
Prompt O,P,Q
O+(P/60)+(Q/3600)\->\W
W+7\->\W
(W*1.0027379094)+V\->\Z
Z-7.461944\->\Y
If Y<0
24+Y\->\Y
If Y>24
Then
Y-24\->\Y
End
Disp "RTASC"
Prompt E,F,I
E+(F/60)+(I/3600)\->\Q
If Y>Q
Y-Q\->\N
If Y<Q
Q-Y\->\N
If Y=Q
N=0
N*15\->\Z
Disp "DECLIN"
Prompt A,B,C
If A<0
B*\(-)\1\->\B
C*\(-)\1\->\C
A+(B/60)+(C/3600)\->\D
(sin(D)*.55614251)+((cos(D)*.83108694*cos(Z))\->\E
sin\^-1\(E)\->\E
iPart(E)\->\\theta\
Fix 0
Disp "ALTITUDE:DEG",\theta\
Fix 2
60*fPart((E)\->\K
iPart(K)\->\W
Fix 0
Disp "MIN",W
Fix 2
60*fPart(K)\->\G
Fix 0
Disp "SEC",G
Float
Pause
(sin(D)-(.55614251*sin(E)))/(.83108694*cos(E))\->\F
cos\^-1\(F)\->\F
iPart(F)\->\W
Fix 0
Disp "AZIMUTH:DEGREE",W
Fix 2
60*fPart((F)\->\M
iPart(M)\->\W
Fix 0
Disp "MIN",W
Fix 2
60*fPart(M)\->\G
Fix 0
Disp "SEC",G
Float
All
comments and suggestions are always welcomed: gar8@earthlink.net