10 REM ================================
20 REM CASIO FX-850P ASTRONOMY UTILITY
30 REM RAD/DEC TO ALT/AZ CONVERTER
40 REM ================================
50 CLEAR
60 PRINT "ASTRONOMY CALC - FX850P"
70 INPUT "LAT (DEG, S=-): ";LA
80 INPUT "LON (DEG, W=-): ";LO
90 INPUT "YEAR (YYYY): ";Y
100 INPUT "MONTH (1-12): ";M
110 INPUT "DAY (1-31): ";D
120 INPUT "LOCAL TIME (HR.DEC): ";UT
130 INPUT "OBJECT RA (HRS.DEC): ";RA
140 INPUT "OBJECT DEC (DEG): ";DC
150 REM --- JULIAN DAY CALCULATION ---
160 IF M>2 THEN GOTO 190
170 Y=Y-1 180 M=M+12 190 A=INT(Y/100)
200 B=2-A+INT(A/4)
210 JD=INT(365.25*(Y+4716)) +INT(30.6001*(M+1))+D+B-1524.5+(UT/24)
220 REM --- GREENWICH SIDEREAL TIME ---
230 T=(JD-2451545)/36525
240 GST=280.46061837+360.98564736629*(JD-245 1545)+T T(0.000387933-T/38710000)
250 REM --- NORMALIZE GST TO 0-360 DEGREES ---
260 GST=GST-360 INT(GST/360)
270 IF GST<0 THEN GST=GST+360
280 REM --- LOCAL SIDEREAL TIME & HOUR ANGLE ---
290 LST=GST+LO
300 HA=LST-(RA15)
310 REM --- CONVERT DEGREES TO RADIANS FOR TRIG FUNCTIONS ---
320 R=3.1415926535/180
330 H=HA R
340 L=LAR
350 D1=DC*R
360 REM --- ALTITUDE CALCULATION ---
370 SA=SIN(L)*SIN(D1)+COS(L)*COS(D1)*COS(H)
380 ALT=ASN(SA)/R
390 REM --- AZIMUTH CALCULATION ---
400 CA=(SIN(D1)-SIN(L)*SA)/ (COS(L) COS(ALTR))
410 AZ=ACS(CA)/R
420 IF SIN(H)>0 THEN AZ=360-AZ
430 REM --- DISPLAY RESULTS ---
440 PRINT "--------------------------------"
450 PRINT "ALTITUDE: ";USING "###.##";ALT;" DEG"
460 PRINT "AZIMUTH : ";USING "###.##";AZ;" DEG"
470 END