home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #18 / NN_1992_18.iso / spool / comp / sources / hp48 / 240 < prev    next >
Encoding:
Text File  |  1992-08-12  |  27.2 KB  |  1,608 lines
  1. Newsgroups: comp.sources.hp48
  2. Path: sparky!uunet!seq!spell
  3. From: Kati Sinenmaa <sinenmaa@cc.helsinki.fi>
  4. Subject:  v06i018:  tyko_ks - Handheld Astronomy v3.0, Part02/03
  5. Message-ID: <1992Aug13.014546.699@seq.uncwil.edu>
  6. Followup-To: comp.sys.hp48
  7. Sender: spell@seq.uncwil.edu (Chris Spell)
  8. Organization: Univ. of North Carolina @ Wilmington
  9. References: <cshp48-v06i017=tyko_ks.214349@seq.uncwil.edu>
  10. Date: Thu, 13 Aug 1992 01:45:46 GMT
  11. Approved: spell@seq.uncwil.edu
  12. Lines: 1598
  13.  
  14. Checksum: 2881625009 (verify with brik -cv)
  15. Submitted-by: Kati Sinenmaa <sinenmaa@cc.helsinki.fi>
  16. Posting-number: Volume 6, Issue 18
  17. Archive-name: tyko_ks/part02
  18.  
  19. BEGIN_RPL tyko.rpl
  20. %%HP: T(3)A(D)F(.);
  21. DIR
  22.   START.
  23. \<< NOVA { { "SETUP"
  24.   \<< CLEAR SETUP
  25.   \>> } { "JULIA"
  26.   \<< CLEAR 'DN'
  27. PURGE
  28. "Julian Day.fraction"
  29. ":JDAY:" INPUT OBJ\->
  30. DUP DUP 'JD' STO FP
  31. 24 * 12 + TZ + DL +
  32. F24 'LT' STO LT TZ
  33. - DL - 24 / - 'JD0'
  34. STO CHUT JD0 CADA
  35. 'DT' STO NOVA
  36.   \>> } { "LT.ST"
  37.   \<< CLEAR JD0 GWUT
  38. LT TZ - DL -
  39. 1.00273790935 * +
  40. F24 DUP \->HMS "GMST"
  41. \->TAG SWAP \GlM 15 / +
  42. F24 DUP \GD\165 \Ge COS *
  43. 15 / + \->HMS "AST"
  44. \->TAG SWAP \->HMS
  45. "MST" \->TAG
  46.   \>> } { "ST.LT"
  47.   \<< CLEAR JD0 GWUT
  48. LT \GlM 15 / - F24
  49. DUP \->HMS "GMST"
  50. \->TAG SWAP ROT - F24
  51. .99726956633 * \->HMS
  52. TZ + DL + F24 "LT"
  53. \->TAG
  54.   \>> } { "SPACE"
  55.   \<< CLEAR SUDEP
  56.   \>> } { "END"
  57.   \<< WIPE PURGE 2
  58. MENU
  59.   \>> } } TMENU
  60. \>>
  61. ASCAL
  62. \<< -16 CF -19 SF
  63. " Type MONTH and YEAR"
  64. ":MM YYYY:" INPUT
  65. OBJ\-> \-> m y
  66.   \<<
  67.     IF m 0 == m 13
  68. \>= OR
  69.     THEN
  70. "Are you an Alien ?"
  71. DOERR
  72.     END m 1 y JDAY
  73. OVER 3 + 7 / FP 7 *
  74. 0 RND 1 + + 7 / 0
  75. RND ERASE -1 39
  76. XRNG 0 24 YRNG
  77. CLLCD
  78. "A calendar v1.0 by KS"
  79. 1 DISP (0,0) (0,21)
  80. (3,0) (3,21) (8,0)
  81. (8,3) (0,0) (13,0)
  82. (8,0) (8,3) (13,0)
  83. (13,3) 1 6
  84.     FOR j LINE
  85.     NEXT 1 7
  86.     FOR j 38 j 3 *
  87. \->V2 0 j 3 * \->V2
  88. LINE 43 j 5 * - 3
  89. \->V2 43 j 5 * - 21
  90. \->V2 LINE
  91.     NEXT PICT RCL 1
  92. 7
  93.     FOR j -1.1 j 5
  94. * + 20.3 \->V2 DAMO j
  95. 1 + GET 1 \->GROB GOR
  96.     NEXT SWAP DUP
  97.     IF 0 ==
  98.     THEN 52
  99.     ELSE DUP
  100.     END 1 \->GROB ROT
  101. (.5,17) ROT GOR
  102. (7,23.5) JG DAMO m
  103. 8 + GET y \->STR + +
  104. 1 \->GROB GOR SWAP 1
  105. + DUP DUP DUP DUP 6
  106. ROLL 0 4
  107.     FOR j .5 14 j 3
  108. * - \->V2 ROT j + 1
  109. \->GROB GOR
  110.     NEXT -42 CF
  111. SWAP DUP 1 + 7 / FP
  112. 7 * 0 RND 'c' STO 1
  113. 'w' STO 1 m + 1 y
  114. JDAY DROP SWAP -
  115.     IF y 1991 \<= y
  116. 2099 \>= OR
  117.     THEN 1 SWAP
  118. "Alarm detect disabled"
  119. 1 DISP
  120.     ELSE DUP y
  121. 1000000 / DUP .01 +
  122. m + SWAP ROT 100 /
  123. + m + 23.5959 2
  124. \->LIST FINDALARM
  125. SWAP FINDALARM -
  126. DUP
  127.       IF 0 ==
  128.       THEN DROP
  129. SWAP (19,2) "No"
  130.       ELSE ROT
  131. (19,2) ROT ABS \->STR
  132.       END
  133. " alarms" + 1
  134. \->GROB GOR 1 ROT
  135.     END
  136.     FOR j
  137.       IF c 0 \=/ c 7
  138. \=/ AND
  139.       THEN c 5 *
  140. 1.3 - 20 w 3 * -
  141. \->V2 j 1 \->GROB GOR c
  142. 1 + 'c' STO
  143.       ELSE 33.7 20
  144. w 3 * - \->V2 j 1
  145. \->GROB GOR 1 'c' STO
  146. w 1 + 'w' STO
  147.       END
  148.     NEXT PICT STO {
  149. } PVIEW { w c PPAR
  150. } PURGE
  151.   \>>
  152. \>>
  153. GSTR
  154. \<< STD \-> y
  155.   \<< y 19 MOD y 100
  156. MOD y 100 / IP DUP
  157. DUP 4 MOD 2 * SWAP
  158. 4 / IP ROT DUP DUP
  159. 8 + 25 / IP - 1 + 3
  160. / IP SWAP ROT -
  161. SWAP - 4 PICK 19 *
  162. + 15 + 30 MOD ROT
  163. DUP 4 / IP 2 * SWAP
  164. 4 MOD SWAP 3 PICK -
  165. SWAP - ROT + 32 + 7
  166. MOD DUP 22 * 3 PICK
  167. 11 * + 4 ROLL + 451
  168. / IP 7 * ROT ROT +
  169. SWAP - 114 + DUP 31
  170. / IP SWAP 31 MOD 1
  171. + \-> M D
  172.     \<< { "March "
  173. "April " } M 2 -
  174. GET D y ABS 10 y
  175. ABS \->STR SIZE DUP
  176. FIX ^ / + y SIGN *
  177. \->STR + STD M
  178.       IF 4 <
  179.       THEN D 22 -
  180. DUP
  181.         IF 0 ==
  182.         THEN DROP
  183. "April " 30
  184.         ELSE "May "
  185. SWAP
  186.         END
  187.       ELSE D 9 +
  188. DUP
  189.         IF 31 >
  190.         THEN
  191. "June " SWAP 31 -
  192.         ELSE "May "
  193. SWAP
  194.         END
  195.       END \->STR +
  196. "Ascension Day"
  197. \->TAG
  198.     \>>
  199.   \>>
  200. \>>
  201. JSTR
  202. \<< STD \-> y
  203.   \<< y 4 MOD 2 * y 7
  204. MOD 4 * + y 19 MOD
  205. 19 * 15 + 30 MOD
  206. SWAP OVER - 34 + 7
  207. MOD 114 + + DUP 31
  208. / IP { "March "
  209. "April " } SWAP 2 -
  210. GET SWAP 31 MOD 1 +
  211. y ABS 10 y ABS \->STR
  212. SIZE DUP FIX ^ / +
  213. y SIGN * \->STR + STD
  214. "Julian Easter"
  215. \->TAG
  216.   \>>
  217. \>>
  218. JDAY
  219. \<< STD \-> m D V
  220.   \<< D m
  221.     IF 3 \>=
  222.     THEN m 3 - V
  223.     ELSE m 9 + V 1
  224. -
  225.     END DUP 4712 +
  226. 365.25 * IP SWAP
  227. 100 / 49 + IP .75 *
  228. IP 38 - V \-> D m y g
  229. W
  230.     \<< y m 30.6 * .5
  231. + IP 59 D + + + DUP
  232. g - "GreCal" \->TAG
  233. SWAP "JulCal" \->TAG
  234.       IF JG SIZE 10
  235. <
  236.       THEN SWAP
  237. DROP DUP 0
  238.       ELSE DROP DUP
  239. W 1 - 100 / 49 + IP
  240. .75 * IP 38 -
  241.       END W 1 -
  242. 4712 + 365.25 * IP
  243. 365 + SWAP - -
  244. "DayNum" \->TAG
  245.     \>>
  246.   \>>
  247. \>>
  248. CADA
  249. \<< STD 0 RND \-> J
  250.   \<< J 1 + 7 / FP 7
  251. * 0 RND 1 + DAMO
  252. SWAP GET
  253.     IF JG SIZE 10
  254. ==
  255.     THEN J 4479.5 -
  256. 36524.25 / IP .75 *
  257. .5 + IP 37 - J +
  258.     ELSE J
  259.     END DUP 59.25 -
  260. 365.25 MOD IP DUP
  261. .5 + 30.6 / IP 2 +
  262. 12 MOD 9 + DAMO
  263. SWAP GET SWAP .5 +
  264. 30.6 MOD IP 1 + ROT
  265. 365.25 / IP 4712 -
  266. DUP ABS DUP 10 SWAP
  267. \->STR SIZE DUP FIX ^
  268. / ROT + SWAP SIGN *
  269. \->STR + + STD
  270.   \>>
  271. \>>
  272. EOFT
  273. \<< 2451545 -
  274. .00273790931 *
  275. .779072 + FP 360 *
  276. F36 \-> L
  277.   \<< '-106.2*SIN(L)+
  278. 596.2*SIN(2*L)+4.4*
  279. SIN(3*L)-12.7*SIN(4
  280. *L)-428.9*COS(L)-
  281. 2.1*COS(2*L)+19.3*
  282. COS(3*L)' EVAL 3600
  283. / 12 SWAP - TZ \GlM
  284. 15 / - + DL +
  285.   \>>
  286. \>>
  287. GWUT
  288. \<< 2451545 - 36525 /
  289. DUP DUP 3 ^
  290. -.0000062 * .093104
  291. ROT 2 ^ * +
  292. 8640184.81287 ROT *
  293. + 24110.54841 +
  294. 3600 / F24 DUP \->HMS
  295. "GMST0" \->TAG SWAP
  296. \>>
  297. EQSO
  298. \<<
  299. "Equinox and Solstice"
  300. ":MM YYYY:" INPUT
  301. OBJ\-> \-> M Y
  302.   \<< M 3 / FP
  303.     IF 0 ==
  304.     THEN Y 2000 -
  305. 1000 / 'P' STO DAES
  306. M 3 / GET EVAL 'P'
  307. PURGE DUP FP 24 *
  308. 12 + TZ + DL + F24
  309. DUO \->HMS "LT" \->TAG
  310. ROT ROT TZ - DL -
  311. 24 / - CADA
  312.     END
  313.   \>>
  314. \>>
  315. CONV
  316. \<< -16 CF \-> l b
  317.   \<< l COS l SIN \Ge
  318. COS * b TAN \Ge SIN *
  319. - \->V2 -16 SF V\-> F36
  320. SWAP DROP "RA" \->TAG
  321. b SIN \Ge COS * b COS
  322. \Ge SIN l SIN * * +
  323. ASIN "DECL" \->TAG
  324. -16 CF
  325.   \>>
  326. \>>
  327. EQEC
  328. \<< { { "MEAN"
  329.   \<< \Ge \GD\Ge - 'E' STO
  330. CONT
  331.   \>> } { "APPAR"
  332.   \<< \Ge 'E' STO CONT
  333.   \>> } } TMENU
  334. "Select conversion"
  335. PROMPT -16 CF STD
  336. HMS\-> '\Gd' STO HMS\->
  337. 15 * '\Ga' STO \Gd SIN
  338. E COS * \Gd COS E SIN
  339. \Ga SIN * * - ASIN
  340. "\Gb" \->TAG \Gd COS \Ga
  341. COS * \Gd SIN E SIN *
  342. \Gd COS E COS \Ga SIN *
  343. * + \->V2 -16 SF V\->
  344. SWAP DROP DUP
  345.   IF 0 <
  346.   THEN 360 +
  347.   END "\Gl" \->TAG { \Ga
  348. \Gd E } PURGE -16 CF
  349. SWAP 0 MENU
  350. \>>
  351. SETUP
  352. \<< { { "DATE"
  353.   \<<
  354. "DATE-form MM DD YYYY"
  355. ":DATE:" INPUT OBJ\->
  356. JDAY 'DN' STO .5 -
  357. 'JD0' STO JD0 LT TZ
  358. - DL - 24 / + 'JD'
  359. STO JD0 CADA 'DT'
  360. STO NOVA
  361.   \>> } { "TIME"
  362.   \<<
  363. "TIME-format HH.MMSS"
  364. ":L/ST:" INPUT OBJ\->
  365. HMS\-> DUP TZ - DL -
  366. 24 / JD0 + 'JD' STO
  367. 'LT' STO CHUT NOVA
  368.   \>> } { "LONG"
  369.   \<<
  370. "Longitude in decimals"
  371. "
  372. West negative
  373. " +
  374. ":\Gl:" INPUT OBJ\->
  375. '\GlM' STO NOVA
  376.   \>> } { "LATIT"
  377.   \<<
  378. "Latitude in decimals"
  379. "
  380. -90 < \O/ < 90
  381. " +
  382. ":\O/:" INPUT OBJ\->
  383. '\O/' STO NOVA
  384.   \>> } { "ZONE"
  385.   \<<
  386. "Integer  -12  to  12"
  387. "
  388. West=neg
  389. " +
  390. ":Time Zone:" INPUT
  391. OBJ\-> 'TZ' STO JD0
  392. LT TZ - DL - 24 / +
  393. 'JD' STO CHUT NOVA
  394.   \>> } { "BACK"
  395.   \<< 0 MENU NOVA
  396.   \>> } { "DST"
  397.   \<< 1 DL 1 * - 'DL'
  398. STO CHUT NOVA
  399.   \>> } { "CALE"
  400.   \<< { "Gregorian "
  401. "  Julian " } JG
  402. SIZE 8 - GET 'JG'
  403. STO JD0 CADA 'DT'
  404. STO NOVA
  405.   \>> } } TMENU NOVA
  406. \>>
  407. SOLAR
  408. \<< 1400 .12 BEEP
  409. 1000 .2 BEEP { {
  410. "Merc"
  411.   \<< CAL1 OBJ\-> DROP
  412. EVAL 'R' STO EVAL
  413. '\Gb' STO EVAL 3600 /
  414. L1 + '\Gl' STO 6.72 .5
  415. 87 G1 1 CONT
  416.   \>> } { "Venu"
  417.   \<< CAL2 OBJ\-> DROP
  418. EVAL 'R' STO EVAL
  419. '\Gb' STO EVAL 3600 /
  420. L2 + '\Gl' STO 16.68
  421. -4. 224 G2 2 CONT
  422.   \>> } { "Mars"
  423.   \<< CAL4 OBJ\-> DROP
  424. EVAL 'R' STO EVAL
  425. '\Gb' STO EVAL 3600 /
  426. L3 + '\Gl' STO 9.36
  427. -1.12 686 G3 3 CONT
  428.   \>> } { "Jupit"
  429.   \<< CAL5 OBJ\-> DROP
  430. EVAL 'R' STO EVAL
  431. '\Gb' STO EVAL 3600 /
  432. L4 + '\Gl' STO 196.88
  433. -9.33 4360 G4 4
  434. CONT
  435.   \>> } { "Satur"
  436.   \<< CAL6 OBJ\-> DROP
  437. EVAL 'R' STO EVAL
  438. '\Gb' STO EVAL 3600 /
  439. L5 + '\Gl' STO 165.46
  440. -8.5 10740 G5 5
  441. CONT
  442.   \>> } { "Uran"
  443.   \<< CAL7 OBJ\-> DROP
  444. EVAL 'R' STO EVAL
  445. '\Gb' STO EVAL 3600 /
  446. L6 + '\Gl' STO 70.04
  447. -7.2 30690 G6 6
  448. CONT
  449.   \>> } { "Neptu"
  450.   \<< CAL8 OBJ\-> DROP
  451. EVAL 'R' STO EVAL
  452. '\Gb' STO EVAL 3600 /
  453. L7 + '\Gl' STO 67
  454. -6.86 60181 G7 7
  455. CONT
  456.   \>> } { "Pluto"
  457.   \<< CAL9 OBJ\-> DROP
  458. EVAL 'R' STO EVAL
  459. '\Gb' STO EVAL 3600 /
  460. L8 + '\Gl' STO 4.14
  461. -1 90000 G8 8 CONT
  462.   \>> } { "Moon"
  463.   \<< CALM OBJ\-> DROP
  464. EVAL '\GD' STO r
  465. 23454.78 * 'R' STO
  466. EVAL 3600 / DUP '\Gb'
  467. STO 'b' STO EVAL
  468. 3600 / L9 \GD\165 + +
  469. F36 DUP '\Gl' STO 'l'
  470. STO \Gb COS \Gl \GlS -
  471. COS * ACOS '\165' STO
  472. \GD R \165 COS * - R \165
  473. SIN * -16 CF \->V2
  474. -16 SF V\-> F36 '\Gm'
  475. STO DROP 936.75 \GD /
  476. \->HMS 29 G9 9 CONT
  477.   \>> } { "CALE"
  478.   \<< { "Gregorian "
  479. "  Julian " } JG
  480. SIZE 8 - GET 'JG'
  481. STO JD0 CADA 'DT'
  482. STO NOVA
  483.   \>> } } TMENU
  484. "   Select  Object"
  485. PROMPT 'n' STO *
  486. 360 / NEG JD + CADA
  487. '\Gt' STO
  488.   IF n 9 <
  489.   THEN \Gl F36 '\Gl'
  490. STO \Gb 3600 / '\Gb'
  491. STO GEOC 5 R \GD * 1
  492. l \Gl - COS + .5 * \v/
  493. / LOG * + 'V' STO r
  494. r * \GD \GD * + R R * -
  495. 2 r \GD * * / ACOS
  496. '\165' STO R R * \GD \GD *
  497. + r r * - 2 R \GD * *
  498. / ACOS '\Gm' STO \GD /
  499.   END 'SD' STO l b
  500. CONV '\Gd' STO '\Ga'
  501. STO AARS 400 .2
  502. BEEP 800 .2 BEEP
  503. 1600 .2 BEEP { {
  504. "PICT"
  505.   \<< DRAV
  506.   \>> } { "FILE"
  507.   \<< CROB
  508.   \>> } { "DATA"
  509.   \<< PDRV
  510.   \>> } { "NEW"
  511.   \<< SOLAR
  512.   \>> } { "END"
  513.   \<< { D L1 G1 F1 L2
  514. G2 F2 L3 F3 L4 F4
  515. L5 G5 F5 L6 G6 F6
  516. L7 G7 F7 L8 G8 F8
  517. G9 F9 \Gt l b \165 \Gm \GD }
  518. PURGE WIPE PURGE 2
  519. MENU
  520.   \>> } } TMENU
  521. \>>
  522. SUDEP
  523. \<< DEG MISC 900 .13
  524. BEEP 500 .2 BEEP {
  525. { "SUN"
  526.   \<< 10 'n' STO r
  527. 'R' STO 15.99 R /
  528. \->HMS 'SD' STO .62
  529. 'BV' STO -26.78 'V'
  530. STO JD0 .5 + EOFT
  531. 'ET' STO \GlS '\Gl' STO
  532. \Gl \GbS '\Gb' STO \Gb CONV
  533. '\Gd' STO '\Ga' STO
  534. AARS SDRV
  535.   \>> } { "SOLAR"
  536.   \<<
  537. [ .374897 .259091 .827362 .700695 .485541 .566441 .505498 .292498 .987353
  538. .849694
  539. .089608 .814794 .133295 .882987 .821218 .870169 .400589 .664614 .846912
  540. .725368
  541. .480856 .663854 .04102 .357355 ]
  542. [ .03629164709 .0367481952 .03386319198 .011367714 .01136759566 .01136762384
  543. .00445046867 .00445040017 .00145575328 .00145569465 .00023080893
  544. .00023080893
  545. .00009294371 .00009294371 .00009294371 .00003269438 .00003269438
  546. .00003265562
  547. .00001672092 .00001672092 .00001663715 .00001115482 .00001104864
  548. .00001104864 ]
  549. C * + OBJ\-> DROP {
  550. F8 G8 L8 F7 G7 L7
  551. F6 G6 L6 F5 G5 L5
  552. F4 L4 F3 L3 F2 L2
  553. F1 G1 L1 D F9 G9 }
  554. \-> i
  555.     \<< 1 24
  556.       FOR j FP 360
  557. * i j GET STO
  558.       NEXT
  559.     \>> SOLAR
  560.   \>> } { "STARS"
  561.   \<< DEEP
  562.   \>> } { "CALE"
  563.   \<< { "Gregorian "
  564. "  Julian " } JG
  565. SIZE 8 - GET 'JG'
  566. STO JD0 CADA 'DT'
  567. STO NOVA
  568.   \>> } { "BACK"
  569.   \<< START.
  570.   \>> } { "END"
  571.   \<< WIPE PURGE 2
  572. MENU
  573.   \>> } } TMENU
  574. \>>
  575. MISC
  576. \<< STD CLLCD
  577. "Wait ..." 1 DISP
  578. JD0 GWUT LT TZ - DL
  579. - 1.00273790935 * +
  580. \GlM 15 / + F24 'ST'
  581. STO HMS\-> 'GST' STO
  582. JD 2451545 - DUP
  583. 'C' STO 36525 / 't'
  584. STO t 1 + 'T' STO
  585. .606434
  586. .03660110129 C * +
  587. .779072
  588. .00273790931 C * +
  589. .993126 .0027377785
  590. C * + .140023
  591. .00445036173 C * +
  592. .053856
  593. .00145561327 C * +
  594. .056531
  595. .00023080893 C * +
  596. { G4 G3 G2 GS LS L9
  597. } DEG \-> i
  598.   \<< 1 6
  599.     FOR j FP 360 *
  600. i j GET STO
  601.     NEXT 125.044555
  602. 1934.1361849 t * -
  603. .0020762 t 2 ^ * +
  604. t 3 ^ 467410 / + t
  605. 4 ^ 60616000 / +
  606. '\GW' STO -.00172 \GW
  607. SIN * .000132 2 LS
  608. * SIN * - .000023 2
  609. L9 * SIN * -
  610. .000021 2 \GW * SIN *
  611. + HMS\-> '\GD\165' STO
  612. 23.2621448 .0046815
  613. t * - .000000059 t
  614. 2 ^ * - .0000001813
  615. t 3 ^ * + HMS\->
  616. .00092 \GW COS *
  617. .000057 2 LS * COS
  618. * + .00001 2 L9 *
  619. COS * + .000009 2 \GW
  620. * COS * - HMS\-> DUP
  621. '\GD\Ge' STO + '\Ge' STO
  622. t 10 / 'p' STO CALS
  623. OBJ\-> DROP EVAL 'r'
  624. STO EVAL 3600 / LS
  625. + F36 \GD\165 .0052 r /
  626. - + '\GlS' STO EVAL
  627. SWAP EVAL p * +
  628. 100000000 / R\->D NEG
  629. '\GbS' STO
  630.   \>>
  631. \>>
  632. CROB
  633. \<< ERASE NM n GET 2
  634. \->GROB PICT STO PICT
  635. RCL { # 0d # 15d }
  636. DT 1 \->GROB GOR 6
  637. FIX { # 60d # 7d }
  638. "GMST0= " GST \->HMS
  639. ZZ { # 71d # 14d }
  640. "LST= " ST \->HMS ZZ
  641. { # 0d # 23d }
  642. "\Gl= " \Gl \->HMS ZZ {
  643. # 0d # 29d } "\Gb= "
  644. \Gb \->HMS ZZ { # 0d
  645. # 50d } "r= " R ZZ
  646. { # 0d # 57d }
  647. "\GD= " \GD ZZ { # 0d
  648. # 37d } "\Ga= " \Ga 15
  649. / \->HMS ZZ { # 0d
  650. # 43d } "\Gd= " \Gd
  651. \->HMS ZZ { # 60d
  652. # 21d } "AZ= " A
  653. \->HMS ZZ { # 60d
  654. # 27d } "a = " a
  655. \->HMS ZZ { # 53d
  656. # 58d } "\Gt= " \Gt + 1
  657. \->GROB GOR 3 FIX {
  658. # 66d # 40d }
  659. "RISE = " RI ZZ {
  660. # 65d # 46d }
  661. "  ET = " ET \->HMS ZZ
  662. { # 65d # 52d }
  663. " SET = " SE ZZ {
  664. # 57d # 33d } "D= "
  665. SD ZZ { # 94d # 33d
  666. } "M= " V ZZ {
  667. # 60d # 1d }
  668. "JDAY= " JD ZZ {
  669. # 0d # 9d } "CAL="
  670. JG + 1 \->GROB GOR
  671. STD PICT STO PICT
  672. RCL \->LCD { } PVIEW
  673. \>>
  674. DRAV
  675. \<< -16 CF -19 SF
  676. ERASE -35 35 XRNG
  677. -15 15 YRNG CLLCD
  678. "Drawing system ..."
  679. 1 DISP DEG
  680.   IF 'IMAG' VTYPE 0
  681. <
  682.   THEN 3 FIX (0,0)
  683. 1.5 0 360 ARC
  684. (-35,0) (0,0) LINE
  685. 1 RES CONIC 'X'
  686. INDEP 'Y' DEPND 'X^
  687. 2/196+Y^2/16=1'
  688. STEQ DRAW 'X^2/625+
  689. Y^2/81=1' STEQ DRAW
  690. PICT RCL { # 0d
  691. # 8d } "Earth" 1
  692. \->GROB GOR 'IMAG'
  693. STO
  694.   END IMAG { # 0d
  695. # 0d } \GlS 0 CONV ST
  696. 15 * ROT - \-> d h
  697.   \<< \GlS 180 + h COS
  698. d COS \O/ SIN * * d
  699. SIN \O/ COS * - h SIN
  700. d COS * \->V2 -16 SF
  701. V\-> SWAP DROP 180 -
  702. + F36 \-> x
  703.     \<<
  704.       IF n 2 \<=
  705.       THEN \Gl COS 14
  706. * \Gl SIN 4 * R\->C 1.2
  707. 0 360 ARC \GlS 180 -
  708. DUP COS 25 * SWAP
  709. SIN 9 * R\->C DUP 1.8
  710. 0 360 ARC x COS 2.3
  711. * x SIN 2.3 * R\->C +
  712. .6 0 360 ARC
  713.       ELSE
  714.         IF n 9 ==
  715.         THEN \GlS 180
  716. - DUP COS 25 * SWAP
  717. SIN 9 * R\->C DUP 1.5
  718. 0 360 ARC x COS 2.3
  719. * x SIN 2.3 * R\->C +
  720. .6 0 360 ARC PICT
  721. RCL GOR \GlS COS 25 *
  722. \GlS SIN 9 * R\->C
  723. CENTR { # 0d # 0d }
  724. 'X' INDEP 'Y' DEPND
  725. CONIC 'X^2/81+Y^2/9
  726. =1' STEQ DRAW \Gl COS
  727. 9 * \Gl SIN 3 * R\->C
  728. 1.1 0 360 ARC
  729.         ELSE \GlS 180
  730. - DUP COS 14 * SWAP
  731. SIN 4 * R\->C DUP 1.5
  732. 0 360 ARC x COS 2.3
  733. * x SIN 2.3 * R\->C +
  734. .6 0 360 ARC \Gl COS
  735. 25 * \Gl SIN 9 * R\->C
  736. 1.5 0 360 ARC
  737.         END
  738.       END 'EQ'
  739. PURGE PICT RCL GOR
  740. { # 0d # 0d } NM n
  741. GET 2 \->GROB GOR {
  742. # 70d # 0d } DT 1
  743. \->GROB GOR { # 1d
  744. # 58d } "\Gb= "
  745.       IF n 9 ==
  746.       THEN \GD 6378 '
  747. 1_km' * *
  748.       ELSE R '1_AU'
  749. *
  750.       END \Gb SIN * 5
  751. FIX ZZ { # 88d
  752. # 58d } "LT= " LT
  753. \->HMS ZZ STD -16 CF
  754. PICT STO PICT RCL
  755. \->LCD { } PVIEW
  756.     \>>
  757.   \>>
  758. \>>
  759. PDRV
  760. \<< STD -19 SF "G" n
  761. + STR\-> F36 "L" n +
  762. STR\-> F36 "F" n +
  763. STR\-> F36 \-> M L F
  764.   \<< 6 FIX "L = " L
  765. \->HMS + "  M = " M
  766. \->HMS + + "\pi = " L M
  767. - F36 \->HMS +
  768. "  \GW = " L F - F36
  769. \->HMS + + "T = " t +
  770. "  \Ge = " \Ge \->HMS + +
  771. "\GD\165 = " \GD\165 \->HMS
  772. 10000 * + "  \GD\Ge = "
  773. \GD\Ge \->HMS 10000 * + +
  774. "l = " l \->HMS +
  775. "  b = " b \->HMS + +
  776. "\165 = " \165 \->HMS +
  777. "  \GDt = " \GD 8.31675
  778. * \->HMS + + "i = " \Gm
  779. \->HMS + "  k = " 1 \Gm
  780. COS + 2 / + + STD
  781. -6.7 6.5 XRNG 0 8.2
  782. YRNG DT " " LT + +
  783.   \>> 2 \->GROB PICT
  784. STO PICT RCL 1 7
  785.   FOR j -6.5 j \->V2
  786. ROT 1 \->GROB GOR
  787.   NEXT PICT STO { }
  788. PVIEW
  789. \>>
  790. SDRV
  791. \<< -16 CF -19 SF 6
  792. FIX "          " DT
  793. + "R= " R + " UT= "
  794. UT + + "\Gl = " \Gl
  795. \->HMS + "  \Gb = " \Gb
  796. \->HMS + + "\Ga =  " \Ga
  797. 15 / \->HMS +
  798. "  \Gd = " \Gd \->HMS + +
  799. "AZI = " A +
  800. "  ALT = " a + + 3
  801. FIX "RISE " RI +
  802. " SET " SE + +
  803. " ET= " ET \->HMS + +
  804. "V= " V + " B-V= "
  805. BV + + " SD= " SD +
  806. + "\GD\165 " \GD\165 \->HMS
  807. 10000 * + " \GD\Ge " \GD\Ge
  808. \->HMS 10000 * + +
  809. STD " \Ge " \Ge \->HMS +
  810. + ERASE -6.7 6.5
  811. XRNG .3 8.3 YRNG
  812. STD NM n GET 2
  813. \->GROB PICT STO PICT
  814. RCL 1 8
  815.   FOR j -6.5 j \->V2
  816. ROT 1 \->GROB GOR
  817.   NEXT PICT STO { }
  818. PVIEW
  819. \>>
  820. CHUT
  821. \<< LT TZ - DL - \-> U
  822.   \<< U F24 \->HMS
  823.     IF U 0 < U 24 \>=
  824. OR
  825.     THEN {
  826. " yesterday" ""
  827. " tomorrow" } U
  828. SIGN 2 + GET
  829.     ELSE " today"
  830.     END + 'UT' STO
  831.   \>>
  832. \>>
  833. DEEP
  834. \<< 'SD' PURGE
  835.   IF 'SDAT' VTYPE 0
  836. <
  837.   THEN
  838. "Star data not exist"
  839. DOERR
  840.   ELSE 1 'n' STO {
  841. { "PICK"
  842.     \<< NEWS
  843.     \>> } { "CONT"
  844.     \<< 1.396971 t *
  845. .0003086 t t * * +
  846. .013056 t *
  847. .0000092 t t * * -
  848. 5.12362 .241614 t *
  849. .0001122 t t * * +
  850. + SDAT n GET OBJ\->
  851. DROP 'V' STO 'BV'
  852. STO "_pc " + 'R'
  853. STO \-> a b c A B
  854.       \<< B A c + SIN
  855. b * + '\Gb' STO A a B
  856. TAN A c + COS b * *
  857. - + '\Gl' STO n 10 +
  858. 'n' STO \Gl \Gb CONV
  859. '\Gd' STO '\Ga' STO
  860. AARS SDRV
  861.       \>>
  862.     \>> } { "BACK"
  863.     \<< 0 MENU
  864.     \>> } { "END"
  865.     \<< WIPE PURGE 2
  866. MENU
  867.     \>> } } TMENU
  868.   END
  869. \>>
  870. AARS
  871. \<< -16 CF ST 15 * \Ga
  872. - \-> h
  873.   \<< h COS \Gd COS \O/
  874. SIN * * \Gd SIN \O/ COS
  875. * - h SIN \Gd COS *
  876. \->V2 -16 SF V\-> F36
  877. \->HMS 'A' STO DROP h
  878. COS \Gd COS \O/ COS * *
  879. \O/ SIN \Gd SIN * +
  880. ASIN \->HMS 'a' STO
  881.     CASE n 10 ==
  882.       THEN -.8333
  883.       END n 9 ==
  884.       THEN \GD INV
  885. ASIN .7275 * .5667
  886. -
  887.       END n 9 \=/ n
  888. 10 \=/ AND
  889.       THEN -.5667
  890.       END
  891.     END
  892.     IF n 10 \=/
  893.     THEN \Ga 15 / GST
  894. - TZ \GlM 15 / - + DL
  895. + F24 'ET' STO
  896.     END SIN \Gd COS \O/
  897. COS * / \Gd TAN \O/ TAN
  898. * - DUP -16 CF
  899.     IF ABS 1 >
  900.     THEN DROP
  901.       IF a 0 >
  902.       THEN "Above"
  903.       ELSE "Below"
  904.       END "Horizon"
  905.     ELSE ACOS 15 /
  906. ET OVER - F24 \->HMS
  907. SWAP ET + F24 \->HMS
  908.     END 'SE' STO
  909. 'RI' STO
  910.   \>>
  911. \>>
  912. DAES
  913. \<< { '2451623.80984+
  914. 365242.37404*P+
  915. .05169*P^2-.00411*P
  916. ^3-.00057*P^4' '
  917. 2451716.56767+
  918. 365241.62603*P+
  919. .00325*P^2+.00888*P
  920. ^3-.0003*P^4' '
  921. 2451810.21715+
  922. 365242.01767*P-
  923. .11575*P^2+.00337*P
  924. ^3+.00078*P^4' '
  925. 2451900.05952+
  926. 365242.74049*P-
  927. .06223*P^2-.00823*P
  928. ^3+.00032*P^4' }
  929. \>>
  930. GEOC
  931. \<< -16 CF R \Gb SIN *
  932. R 2 ^ r 2 ^ + DUP 2
  933. R r \GlS \Gl - 180 +
  934. COS * * * SWAP OVER
  935. - \v/ '\GD' STO \Gb COS *
  936. - \v/ / ASIN 'b' STO
  937. \Gb COS \Gl COS R * * r
  938. \GlS COS * + R \Gb COS
  939. \Gl SIN * * r \GlS SIN
  940. * + \->V2 -16 SF V\->
  941. F36 'l' STO DROP
  942. \>>
  943. CALM
  944. \<< { '22640*SIN(G9)-
  945. 4586*SIN(G9-2*D)+
  946. 2370*SIN(2*D)+769*
  947. SIN(2*G9)-668*SIN(
  948. GS)-412*SIN(2*F9)-
  949. 212*SIN(2*G9-2*D)-
  950. 206*SIN(G9-2*D+GS)+
  951. 192*SIN(G9+2*D)+165
  952. *SIN(2*D-GS)+148*
  953. SIN(G9-GS)-125*SIN(
  954. D)-110*SIN(G9+GS)-
  955. 55*SIN(2*F9-2*D)-45
  956. *SIN(G9+2*F9)+40*
  957. SIN(G9-2*F9)-38*SIN
  958. (G9-4*D)+36*SIN(3*
  959. G9)-31*SIN(2*G9-4*D
  960. )+28*SIN(G9-2*D-GS)
  961. -24*SIN(2*D+GS)+19*
  962. SIN(G9-D)+18*SIN(D+
  963. GS)+15*SIN(G9+2*D-
  964. GS)+14*SIN(2*G9+2*D
  965. )+14*SIN(4*D)-13*
  966. SIN(3*G9-2*D)-11*
  967. SIN(G9+16*LS-18*L2)
  968. +10*SIN(2*G9-GS)+9*
  969. SIN(G9-2*F9-2*D)+9*
  970. COS(G9+16*LS-18*L2)
  971. -9*SIN(2*G9-2*D+GS)
  972. -8*SIN(G9+D)+8*SIN(
  973. 2*D-2*GS)-8*SIN(2*
  974. G9+GS)-7*SIN(2*GS)-
  975. 7*SIN(G9-2*D+2*GS)+
  976. 7*SIN(\GW)-6*SIN(G9-2
  977. *F9+2*D)-6*SIN(2*F9
  978. +2*D)-4*SIN(G9-4*D+
  979. GS)+4*T*COS(G9+16*
  980. LS-18*L2)-4*SIN(2*
  981. G9+2*F9)+4*T*SIN(G9
  982. +16*LS-18*L2)+3*SIN
  983. (G9-3*D)-3*SIN(G9+2
  984. *D+GS)-3*SIN(2*G9-4
  985. *D+GS)+3*SIN(G9-2*
  986. GS)+3*SIN(G9-2*D-2*
  987. GS)-2*SIN(2*G9-2*D-
  988. GS)-2*SIN(2*F9-2*D+
  989. GS)+2*SIN(G9+4*D)+2
  990. *SIN(4*G9)+2*SIN(4*
  991. D-GS)+2*SIN(2*G9-D)
  992. ' '18461*SIN(F9)+
  993. 1010*SIN(G9+F9)+
  994. 1000*SIN(G9-F9)-624
  995. *SIN(F9-2*D)-199*
  996. SIN(G9-F9-2*D)-167*
  997. SIN(G9+F9-2*D)+117*
  998. SIN(F9+2*D)+62*SIN(
  999. 2*G9+F9)+33*SIN(G9-
  1000. F9+2*D)+32*SIN(2*G9
  1001. -F9)-30*SIN(F9-2*D+
  1002. GS)-16*SIN(2*G9+F9-
  1003. 2*D)+15*SIN(G9+F9+2
  1004. *D)+12*SIN(F9-2*D-
  1005. GS)-9*SIN(G9-F9-2*D
  1006. +GS)-8*SIN(F9+\GW)+8*
  1007. SIN(F9+2*D-GS)-7*
  1008. SIN(G9+F9-2*D+GS)+7
  1009. *SIN(G9+F9-GS)-7*
  1010. SIN(G9+F9-4*D)-6*
  1011. SIN(F9+GS)-6*SIN(3*
  1012. F9)+6*SIN(G9-F9-GS)
  1013. -5*SIN(F9+D)-5*SIN(
  1014. G9+F9+GS)-5*SIN(G9-
  1015. F9+GS)+5*SIN(F9-GS)
  1016. +5*SIN(F9-D)+4*SIN(
  1017. 3*G9+F9)-4*SIN(F9-4
  1018. *D)-3*SIN(G9-F9-4*D
  1019. )+3*SIN(G9-3*F9)-2*
  1020. SIN(2*G9-F9-4*D)-2*
  1021. SIN(3*F9-2*D)+2*SIN
  1022. (2*G9-F9+2*D)+2*SIN
  1023. (G9-F9+2*D-GS)+2*
  1024. SIN(2*G9-F9-2*D)+2*
  1025. SIN(3*G9-F9)' '
  1026. 60.36298-3.27746*
  1027. COS(G9)-.57994*COS(
  1028. G9-2*D)-.46357*COS(
  1029. 2*D)-.08904*COS(2*
  1030. G9)+.03865*COS(2*G9
  1031. -2*D)-.03237*COS(2*
  1032. D-GS)-.02688*COS(G9
  1033. +2*D)-.02358*COS(G9
  1034. -2*D+GS)-.0203*COS(
  1035. G9-GS)+.01719*COS(D
  1036. )+.01671*COS(G9+GS)
  1037. +.01247*COS(G9-2*F9
  1038. )+.00704*COS(GS)+
  1039. .00529*COS(2*D+GS)-
  1040. .00524*COS(G9-4*D)+
  1041. .00398*COS(G9-2*D-
  1042. GS)-.00366*COS(3*G9
  1043. )-.00295*COS(2*G9-4
  1044. *D)-.00263*COS(D+GS
  1045. )+.00249*COS(3*G9-2
  1046. *D)-.00221*COS(G9+2
  1047. *D-GS)+.00185*COS(2
  1048. *F9-2*D)-.00161*COS
  1049. (2*D-2*GS)+.00147*
  1050. COS(G9+2*F9-2*D)-
  1051. .00142*COS(4*D)+
  1052. .00139*COS(2*G9-2*D
  1053. +GS)-.00118*COS(G9-
  1054. 4*D+GS)-.00116*COS(
  1055. 2*G9+2*D)-.0011*COS
  1056. (2*G9-GS)' }
  1057. \>>
  1058. CALS
  1059. { '9*COS(3.9+
  1060. 5507.55*p)+6*COS(
  1061. 1.73+5223.69*p)' '
  1062. 280*COS(3.199+
  1063. 84334.662*p)+102*
  1064. COS(5.422+5507.553*
  1065. p)+80*COS(3.88+
  1066. 5223.69*p)+44*COS(
  1067. 3.7+2352.87*p)+32*
  1068. COS(4+1577.34*p)' '
  1069. 6910*SIN(GS)+72*SIN
  1070. (2*GS)-17*T*SIN(GS)
  1071. -7*COS(GS-G4)+6*SIN
  1072. (L9-LS)+5*SIN(4*GS-
  1073. 8*G3+3*G4)-5*COS(2*
  1074. GS-2*G2)-4*SIN(GS-
  1075. G2)+4*COS(4*GS-8*G3
  1076. +3*G4)+3*SIN(2*GS-2
  1077. *G2)-3*SIN(G4)-3*
  1078. SIN(2*GS-2*G4)' '
  1079. 1.00014-.01675*COS(
  1080. GS)-.00014*COS(2*GS
  1081. )' }
  1082. CAL1
  1083. \<< { '84378*SIN(G1)+
  1084. 10733*SIN(2*G1)+
  1085. 1892*SIN(3*G1)-646*
  1086. SIN(2*F1)+381*SIN(4
  1087. *G1)-306*SIN(G1-2*
  1088. F1)-274*SIN(G1+2*F1
  1089. )-92*SIN(2*G1+2*F1)
  1090. +83*SIN(5*G1)-28*
  1091. SIN(3*G1+2*F1)+25*
  1092. SIN(2*G1-2*F1)+19*
  1093. SIN(6*G1)-9*SIN(4*
  1094. G1+2*F1)+8*T*SIN(G1
  1095. )+7*COS(2*G1-5*G2)' 
  1096. '24134*SIN(F1)+5180
  1097. *SIN(G1-F1)+4910*
  1098. SIN(G1+F1)+1124*SIN
  1099. (2*G1+F1)+271*SIN(3
  1100. *G1+F1)+132*SIN(2*
  1101. G1-F1)+67*SIN(4*G1+
  1102. F1)+18*SIN(3*G1-F1)
  1103. +17*SIN(5*G1+F1)-10
  1104. *SIN(3*F1)-9*SIN(G1
  1105. -3*F1)' '.39528-
  1106. .07834*COS(G1)-
  1107. .00795*COS(2*G1)-
  1108. .00121*COS(3*G1)-
  1109. .00022*COS(4*G1)' }
  1110. \>>
  1111. CAL2
  1112. \<< { '2814*SIN(G2)-
  1113. 181*SIN(2*F2)-20*T*
  1114. SIN(G2)+12*SIN(2*G2
  1115. )-10*COS(2*GS-2*G2)
  1116. +7*COS(3*GS-3*G2)' 
  1117. '12215*SIN(F2)+83*
  1118. SIN(G2+F2)+83*SIN(
  1119. G2-F2)' '.72335-
  1120. .00493*COS(G2)' }
  1121. \>>
  1122. CAL4
  1123. \<< { '38451*SIN(G3)+
  1124. 2238*SIN(2*G3)+181*
  1125. SIN(3*G3)-52*SIN(2*
  1126. F3)+37*T*SIN(G3)-22
  1127. *COS(G3-2*G4)-19*
  1128. SIN(G3-G4)+17*COS(
  1129. G3-G4)+17*SIN(4*G3)
  1130. -16*COS(2*G3-2*G4)+
  1131. 13*COS(GS-2*G3)-10*
  1132. SIN(G3-2*F3)-10*SIN
  1133. (G3+2*F3)+7*COS(GS-
  1134. G3)-7*COS(2*GS-3*G3
  1135. )-5*SIN(G2-3*G3)-5*
  1136. SIN(GS-G3)-5*SIN(GS
  1137. -2*G3)-4*COS(2*GS-4
  1138. *G3)+4*T*SIN(2*G3)+
  1139. 4*COS(G4)+3*COS(G2-
  1140. 3*G3)+3*SIN(2*G3-2*
  1141. G4)' '6603*SIN(F3)+
  1142. 622*SIN(G3-F3)+615*
  1143. SIN(G3+F3)+64*SIN(2
  1144. *G3+F3)' '1.53031-
  1145. .1417*COS(G3)-.0066
  1146. *COS(2*G3)-.00047*
  1147. COS(3*G3)' }
  1148. \>>
  1149. CAL5
  1150. \<< { '19934*SIN(G4)+
  1151. 5023*T+2511+1093*
  1152. COS(2*G4-5*G5)+601*
  1153. SIN(2*G4)-479*SIN(2
  1154. *G4-5*G5)-185*SIN(2
  1155. *G4-2*G5)+137*SIN(3
  1156. *G4-5*G5)-131*SIN(
  1157. G4-2*G5)+79*COS(G4-
  1158. G5)-76*COS(2*G4-2*
  1159. G5)-74*T*COS(G4)+68
  1160. *T*SIN(G4)+66*COS(2
  1161. *G4-3*G5)+63*COS(3*
  1162. G4-5*G5)+53*COS(G4-
  1163. 5*G5)+49*SIN(2*G4-3
  1164. *G5)-43*T*SIN(2*G4-
  1165. 5*G5)-37*COS(G4)+25
  1166. *SIN(2*L4)+25*SIN(3
  1167. *G4)-23*SIN(G4-5*G5
  1168. )-19*T*COS(2*G4-5*
  1169. G5)+17*COS(2*G4-4*
  1170. G5)+17*COS(3*G4-3*
  1171. G5)-14*SIN(G4-G5)-
  1172. 13*SIN(3*G4-4*G5)-9
  1173. *COS(2*L4)+9*COS(G5
  1174. )-9*SIN(G5)-9*SIN(3
  1175. *G4-2*G5)+9*SIN(4*
  1176. G4-5*G5)+9*SIN(2*G4
  1177. -6*G5+3*G6)-8*COS(4
  1178. *G4-10*G5)+7*COS(3*
  1179. G4-4*G5)-7*COS(G4-3
  1180. *G5)-7*SIN(4*G4-10*
  1181. G5)-7*SIN(G4-3*G5)+
  1182. 6*COS(4*G4-5*G5)-6*
  1183. SIN(3*G4-3*G5)+5*
  1184. COS(2*G5)-4*SIN(4*
  1185. G4-4*G5)-4*COS(3*G5
  1186. )+4*COS(2*G4-G5)-4*
  1187. COS(3*G4-2*G5)-4*T*
  1188. COS(2*G4)+3*T*SIN(2
  1189. *G4)+3*COS(5*G5)+3*
  1190. COS(5*G4-10*G5)+3*
  1191. SIN(2*G5)-2*SIN(2*
  1192. L4-G4)+2*SIN(2*L4+
  1193. G4)-2*T*SIN(3*G4-5*
  1194. G5)-2*T*SIN(G4-5*G5
  1195. )' '-4692*COS(G4)+
  1196. 259*SIN(G4)+227-227
  1197. *COS(2*G4)+30*T*SIN
  1198. (G4)+21*T*COS(G4)+
  1199. 16*SIN(3*G4-5*G5)-
  1200. 13*SIN(G4-5*G5)-12*
  1201. COS(3*G4)+12*SIN(2*
  1202. G4)+7*COS(3*G4-5*G5
  1203. )-5*COS(G4-5*G5)' '
  1204. 5.20883-.25122*COS(
  1205. G4)-.00604*COS(2*G4
  1206. )+.0026*COS(2*G4-2*
  1207. G5)-.0017*COS(3*G4-
  1208. 5*G5)-.00106*SIN(2*
  1209. G4-2*G5)-.00091*T*
  1210. SIN(G4)-.00084*T*
  1211. COS(G4)+.00069*SIN(
  1212. 2*G4-3*G5)-.00067*
  1213. SIN(G4-5*G5)+.00066
  1214. *SIN(3*G4-5*G5)+
  1215. .00063*SIN(G4-G5)-
  1216. .00051*COS(2*G4-3*
  1217. G5)-.00046*SIN(G4)-
  1218. .00029*COS(G4-5*G5)
  1219. +.00027*COS(G4-2*G5
  1220. )-.00022*COS(3*G4)-
  1221. .00021*SIN(2*G4-5*
  1222. G5)' }
  1223. \>>
  1224. CAL6
  1225. \<< { '23045*SIN(
  1226. G5)+5014*T-2689*COS
  1227. (2*G4-5*G5)+2507+
  1228. 1177*SIN(2*G4-5*G5)
  1229. -826*COS(2*G4-4*G5)
  1230. +802*SIN(2*G5)+425*
  1231. SIN(G4-2*G5)-229*T*
  1232. COS(G5)-153*COS(2*
  1233. G4-6*G5)-142*T*SIN(
  1234. G5)-114*COS(G5)+101
  1235. *T*SIN(2*G4-5*G5)-
  1236. 70*COS(2*L5)+67*SIN
  1237. (2*L5)+66*SIN(2*G4-
  1238. 6*G5)+60*T*COS(2*G4
  1239. -5*G5)+41*SIN(G4-3*
  1240. G5)+39*SIN(3*G5)+31
  1241. *SIN(G4-G5)+31*SIN(
  1242. 2*G4-2*G5)-29*COS(2
  1243. *G4-3*G5)-28*SIN(2*
  1244. G4-6*G5+3*G6)+28*
  1245. COS(G4-3*G5)+22*T*
  1246. SIN(2*G4-4*G5)-22*
  1247. SIN(G5-3*G6)+20*SIN
  1248. (2*G4-3*G5)+20*COS(
  1249. 4*G4-10*G5)+19*COS(
  1250. 2*G5-3*G6)+19*SIN(4
  1251. *G4-10*G5)-17*T*COS
  1252. (2*G5)-16*COS(G5-3*
  1253. G6)-12*SIN(2*G4-4*
  1254. G5)+12*COS(G4)-12*
  1255. SIN(2*G5-2*G6)-11*T
  1256. *SIN(2*G5)-11*COS(2
  1257. *G4-7*G5)+10*SIN(2*
  1258. G5-3*G6)+10*COS(2*
  1259. G4-2*G5)+9*SIN(4*G4
  1260. -9*G5)-8*SIN(G5-2*
  1261. G6)-8*COS(2*L5+G5)+
  1262. 8*COS(2*L5-G5)+8*
  1263. COS(G5-G6)-8*SIN(2*
  1264. L5-G5)+7*SIN(2*L5+
  1265. G5)-7*COS(G4-2*G5)-
  1266. 7*COS(2*G5)-6*T*SIN
  1267. (4*G4-10*G5)+6*T*
  1268. COS(4*G4-10*G5)+6*T
  1269. *SIN(2*G4-6*G5)-5*
  1270. SIN(3*G4-7*G5)-5*
  1271. COS(3*G4-3*G5)-5*
  1272. COS(2*G5-2*G6)+5*
  1273. SIN(3*G4-4*G5)+5*
  1274. SIN(2*G4-7*G5)+4*
  1275. SIN(3*G4-3*G5)+4*
  1276. SIN(3*G4-5*G5)+4*T*
  1277. COS(G4-2*G5)+3*T*
  1278. COS(2*G4-4*G5)+3*
  1279. COS(2*G4-6*G5+3*G6)
  1280. -3*T*SIN(2*L5)+3*T*
  1281. COS(2*G4-6*G5)-3*T*
  1282. COS(2*L5)+3*COS(3*
  1283. G4-7*G5)+3*COS(4*G4
  1284. -9*G5)+3*SIN(3*G4-6
  1285. *G5)+3*SIN(2*G4-G5)
  1286. +3*SIN(G4-4*G5)+2*
  1287. COS(3*G5-3*G6)+2*T*
  1288. SIN(G4-2*G5)+2*SIN(
  1289. 4*G5)-2*COS(3*G4-4*
  1290. G5)-2*COS(2*G4-G5)-
  1291. 2*SIN(2*G4-7*G5)+2*
  1292. COS(G4-4*G5)+2*COS(
  1293. 4*G4-11*G5)-2*SIN(
  1294. G5-G6)' '8297*SIN(
  1295. G5)-3346*COS(G5)+
  1296. 462*SIN(2*G5)-189*
  1297. COS(2*G5)+185+79*T*
  1298. COS(G5)-71*COS(2*G4
  1299. -4*G5)+46*SIN(2*G4-
  1300. 6*G5)-45*COS(2*G4-6
  1301. *G5)+29*SIN(3*G5)-
  1302. 20*COS(2*G4-3*G5)+
  1303. 18*T*SIN(G5)-14*COS
  1304. (2*G4-5*G5)-11*COS(
  1305. 3*G5)-10*T+9*SIN(G4
  1306. -3*G5)+8*SIN(G4-G5)
  1307. -6*SIN(2*G4-3*G5)+5
  1308. *SIN(2*G4-7*G5)-5*
  1309. COS(2*G4-7*G5)+4*
  1310. SIN(2*G4-5*G5)-4*T*
  1311. SIN(2*G5)-3*COS(G4-
  1312. G5)+3*COS(G4-3*G5)+
  1313. 3*T*SIN(2*G4-4*G5)+
  1314. 3*SIN(G4-2*G5)+2*
  1315. SIN(4*G5)-2*COS(2*
  1316. G4-2*G5)' '9.55774-
  1317. .53252*COS(G5)-
  1318. .01878*SIN(2*G4-4*
  1319. G5)-.01482*COS(2*G5
  1320. )+.00817*SIN(G4-G5)
  1321. -.00539*COS(G4-2*G5
  1322. )-.00524*T*SIN(G5)+
  1323. .00349*SIN(2*G4-5*
  1324. G5)+.00347*SIN(2*G4
  1325. -6*G5)+.00328*T*COS
  1326. (G5)-.00225*SIN(G5)
  1327. +.00149*COS(2*G4-6*
  1328. G5)-.00126*COS(2*G4
  1329. -2*G5)+.00104*COS(
  1330. G4-G5)+.00101*COS(2
  1331. *G4-5*G5)+.00098*
  1332. COS(G4-3*G5)-.00073
  1333. *COS(2*G4-3*G5)-
  1334. .00062*COS(3*G5)+
  1335. .00042*SIN(2*G5-3*
  1336. G6)+.00041*SIN(2*G4
  1337. -2*G5)-.0004*SIN(G4
  1338. -3*G5)+.0004*COS(2*
  1339. G4-4*G5)-.00028*T-
  1340. .00023*SIN(G4)+
  1341. .0002*SIN(2*G4-7*G5
  1342. )' }
  1343. \>>
  1344. CAL7
  1345. \<< { '19397*SIN(G6)+
  1346. 570*SIN(2*G6)-536*T
  1347. *COS(G6)+143*SIN(G5
  1348. -2*G6)+110*T*SIN(G6
  1349. )+102*SIN(G5-3*G6)+
  1350. 76*COS(G5-3*G6)-49*
  1351. SIN(G4-G6)+64*T-30*
  1352. T*COS(2*G6)+29*SIN(
  1353. 2*G4-6*G5+3*G6)+29*
  1354. COS(2*G6-2*G7)-28*
  1355. COS(G6-G7)+23*SIN(3
  1356. *G6)-21*COS(G4-G6)+
  1357. 20*SIN(G6-G7)+20*
  1358. COS(G5-2*G6)-19*COS
  1359. (G5-G6)+17*SIN(2*G6
  1360. -3*G7)+14*SIN(3*G6-
  1361. 3*G7)+13*SIN(G5-G6)
  1362. -24*T*COS(G6)-12*
  1363. COS(G6)+10*SIN(2*G6
  1364. -2*G7)-9*SIN(2*F6)-
  1365. 18*T*SIN(G6)+9*COS(
  1366. 2*G6-3*G7)+8*T*COS(
  1367. G5-2*G6)+7*T*COS(G5
  1368. -3*G6)-7*T*SIN(G5-3
  1369. *G6)+7*T*SIN(2*G6)+
  1370. 6*SIN(2*G4-6*G5+2*
  1371. G6)+6*COS(2*G4-6*G5
  1372. +2*G6)+5*SIN(G5-4*
  1373. G6)-4*SIN(3*G6-4*G7
  1374. )+4*COS(3*G6-3*G7)-
  1375. 3*COS(G7)-2*SIN(G7)
  1376. ' '2775*SIN(F6)+131
  1377. *SIN(G6-F6)+130*SIN
  1378. (G6+F6)' '19.21216-
  1379. .90154*COS(G6)-
  1380. .02488*T*SIN(G6)-
  1381. .02121*COS(2*G6)-
  1382. .00585*COS(G5-2*G6)
  1383. -.00508*T*COS(G6)-
  1384. .00451*COS(G4-G6)+
  1385. .00336*SIN(G5-G6)+
  1386. .00198*SIN(G4-G6)+
  1387. .00118*COS(G5-3*G6)
  1388. +.00107*SIN(G5-2*G6
  1389. )-.00103*T*SIN(2*G6
  1390. )-.00081*COS(3*G6-3
  1391. *G7)' }
  1392. \>>
  1393. CAL8
  1394. \<< { '3523*SIN(
  1395. G7)-50*SIN(2*F7)-43
  1396. *T*COS(G7)+29*SIN(
  1397. G4-G7)+19*SIN(2*G7)
  1398. -18*COS(G4-G7)+13*
  1399. COS(G5-G7)+13*SIN(
  1400. G5-G7)-9*SIN(2*G6-3
  1401. *G7)+9*COS(2*G6-2*
  1402. G7)-5*COS(2*G6-3*G7
  1403. )-4*T*SIN(G7)+4*COS
  1404. (G6-2*G7)+8*T*SIN(
  1405. G7)' '6404*SIN(F7)+
  1406. 55*SIN(G7+F7)+55*
  1407. SIN(G7-F7)-33*T*SIN
  1408. (F7)' '30.07175-
  1409. .25701*COS(G7)-
  1410. .00787*COS(2*L6-G6-
  1411. 2*L7)+.00409*COS(G4
  1412. -G7)-.00314*T*SIN(
  1413. G7)+.0025*SIN(G4-G7
  1414. )-.00194*SIN(G5-G7)
  1415. +.00185*COS(G5-G7)'
  1416. }
  1417. \>>
  1418. CAL9
  1419. \<< { '101577*SIN(G8)
  1420. +15517*SIN(2*G8)-
  1421. 3593*SIN(2*F8)+3414
  1422. *SIN(3*G8)-2201*SIN
  1423. (G8-2*F8)-1871*SIN(
  1424. G8+2*F8)+839*SIN(4*
  1425. G8)-757*SIN(2*G8+2*
  1426. F8)-285*SIN(3*G8+2*
  1427. F8)+454*T*SIN(G8)+
  1428. 218*SIN(2*G8-2*F8)+
  1429. 200*T*SIN(G8)' '
  1430. 57726*SIN(F8)+15257
  1431. *SIN(G8-F8)+14102*
  1432. SIN(G8+F8)+3870*SIN
  1433. (2*G8+F8)+1138*SIN(
  1434. 3*G8+F8)+472*SIN(2*
  1435. G8-F8)+353*SIN(4*G8
  1436. +F8)-144*SIN(G8-3*
  1437. F8)-119*SIN(3*F8)-
  1438. 111*SIN(G8+3*F8)' '
  1439. 40.74638-9.58235*
  1440. COS(G8)-1.16703*COS
  1441. (2*G8)-.22649*COS(3
  1442. *G8)-.04996*COS(4*
  1443. G8)' }
  1444. \>>
  1445. SDAT
  1446. { { 345.31156524
  1447. -59.3781705 26 -.16
  1448. .46 } { 69.78901344
  1449. -5.4674392 21 1.54
  1450. .85 } { 76.8297184
  1451. -31.1228063 280
  1452. -.03 .12 } {
  1453. 81.8577415
  1454. 22.8643432 13 .8
  1455. .08 } { 88.75414772
  1456. -16.0268472 95 1.85
  1457. .5 } { 104.9603567
  1458. -75.8237834 360 .15
  1459. -.72 } {
  1460. 104.0815749
  1461. -39.6052509 2.7 .01
  1462. -1.46 } {
  1463. 115.785452
  1464. -16.0196217 3.5 .42
  1465. .38 } { 113.2154745
  1466. 6.684069 11 1 1.14
  1467. } { 149.8287058
  1468. -.46469231 26 -.11
  1469. 1.35 } { 221.86971
  1470. -52.87803 110 .1
  1471. 1.41 } {
  1472. 221.64572335
  1473. -48.63873204 130
  1474. -.23 1.25 } {
  1475. 203.8410821
  1476. -2.0546734 79 -.23
  1477. .98 } { 233.7920649
  1478. -44.1374373 140
  1479. -.24 .61 } {
  1480. 204.2332698
  1481. 30.73620812 11 1.23
  1482. -.04 } { 239.479783
  1483. -42.594071 1.3 .68
  1484. 0 } { 249.7616937
  1485. -4.56997863 100
  1486. 1.83 .96 } {
  1487. 285.3154391
  1488. 61.7328474 8.1 0
  1489. .03 } { 301.775503
  1490. 29.3036472 5.1 .22
  1491. .77 } { 335.3285843
  1492. 59.9062706 560 .09
  1493. 1.25 } {
  1494. 333.8597432
  1495. -21.13538399 6.7
  1496. .09 1.16 } }
  1497. DAMO
  1498. { "Sun " "Mon "
  1499. "Tue " "Wed "
  1500. "Thu " "Fri "
  1501. "Sat " "Sun" "Jan "
  1502. "Feb " "Mar "
  1503. "Apr " "May "
  1504. "Jun " "Jul "
  1505. "Aug " "Sep "
  1506. "Oct " "Nov "
  1507. "Dec " }
  1508. NM
  1509. { "Mercury" "Venus"
  1510. "Mars" "Jupiter"
  1511. "Saturn" "Uranus"
  1512. "Neptune" "Pluto"
  1513. "Moon" "Sun"
  1514. "\Ga Eri" "\Ga Tau"
  1515. "\Gb Ori" "\Ga Aur"
  1516. "\Ga Ori" "\Ga Car"
  1517. "\Ga CMa" "\Ga CMi"
  1518. "\Gb Gem" "\Ga Leo"
  1519. "\Ga Cru" "\Gb Cru"
  1520. "\Ga Vir" "\Gb Cen"
  1521. "\Ga Boo" "\Ga Cen"
  1522. "\Ga Sco" "\Ga Lyr"
  1523. "\Ga Aql" "\Ga Gyg"
  1524. "\Ga PsA" }
  1525. NEWS
  1526. \<<
  1527. "Type a Star's number ?"
  1528. ":No:" INPUT OBJ\->
  1529. 'n' STO
  1530.   IF n 21 >
  1531.   THEN CLLCD
  1532. "Uncompanionable star"
  1533. 1 DISP 3 FREEZE
  1534.   ELSE CLLCD NM n
  1535. 10 + GET 1 DISP 3
  1536. FREEZE
  1537.   END
  1538. \>>
  1539. NOVA
  1540. \<< STD JG "  DST= "
  1541. DL + + "
  1542. TZ= " TZ
  1543. + "  LT= " LT \->HMS
  1544. + + + "
  1545. JD0 = " JD0
  1546. + + "
  1547. JD  = " JD +
  1548. + "
  1549. UT= " UT + +
  1550. "
  1551. " DT "  " DN + +
  1552. + + "
  1553. \Gl: " \GlM + +
  1554. "  \O/: " \O/ + + CLLCD
  1555. 1 DISP 3 FREEZE
  1556. \>>
  1557. WIPE
  1558. { t C ST GST T G4
  1559. G3 G2 GS LS L9 \GW p
  1560. r \Gl \Gb \GlS \GbS SD ET V
  1561. BV \Ga \Gd n A a RI SE
  1562. R PPAR DN }
  1563. ZZ
  1564. \<< + 1 \->GROB GOR
  1565. \>>
  1566. F24
  1567. \<< DUP 24 / IP 24 *
  1568. -
  1569.   IF DUP 0 <
  1570.   THEN 24 +
  1571.   END
  1572. \>>
  1573. F36
  1574. \<< DUP 360 / IP 360
  1575. * -
  1576.   IF DUP 0 <
  1577.   THEN 360 +
  1578.   END
  1579. \>>
  1580. \GD\165
  1581. .0016
  1582. \GD\Ge
  1583. .0001
  1584. \GlM
  1585. 6
  1586. \O/
  1587. 7
  1588. \Ge
  1589. 23
  1590. JD0
  1591. 4
  1592. JD
  1593. 5
  1594. JG
  1595. "Greetings "
  1596. UT
  1597. "Tyko v30 1992"
  1598. DT
  1599. "    by Kati Sinenmaa"
  1600. DL
  1601. 1
  1602. LT
  1603. 3
  1604. TZ
  1605. 2
  1606.  
  1607. END_RPL
  1608.