Она составлена в незапамятные времена (в начале 90-х) на доблестном ПК "Spectrum 48K"
Какое это было наслаждение!
Собственными руками написать программу по расчету восходов и заходов светил
Сейчас кажутся смешными и наивными 48 килобайт памяти и магнитофонные кассеты с записью программ
Но тогда это было даже необходимостью, поскольку АК уже было не достать, а вычислять на калькуляторе долго
Вот и составлял программы по восходам светил, по вычислению их координат и т.п.
Спасибо Вам за то, что позволили вспомнить о хорошо забытом старом!
Программа по восходам небольшая, поэтому вставляю ее листинг прямо в тело сообщения
Другие программы вместе с интерпретатором языка Бейсик для их использования интересующиеся могут найти во вложении
DECLARE FUNCTION ACS! (X!)
DECLARE FUNCTION ASN! (X!)
10 PRINT "Восход и заход светил для данной местности (по книге Duffet-Smith"
20 PRINT " Практическая астрономия с калькулятором)"
30 PRINT
31 PRINT "Введите прямое восхождение и склонение светила"
35 PRINT
40 INPUT "Часы ", A
41 INPUT "Минуты ", M
42 INPUT "Секунды ", T
43 INPUT "Градусы ", D
44 INPUT "Минуты", Z
45 INPUT "Секунды", Y
50 A = A + ((M + T / 60) / 60)
60 D = D + ((Z + Y / 60) / 60)
70 PRINT "Широта местности"
75 PRINT
76 PRINT "Долгота местности"
77 PRINT
80 PRINT "Прямое восхождение ="; A; " часов"
81 PRINT "Склонение ="; D; " градусов"
85 PRINT
90 P = 57.295779#
91 F = 54.166
100 D = D / P
101 L = 3.49
110 N = SIN(D) / COS(F)
111 R = ACS(N)
120 R = R * P
121 PRINT "Азимут восхода ="; R; " градусoв"
130 S = 360 - R
131 PRINT "Азимут захода ="; S; " градусов"
140 E = ACS(-TAN(F) * TAN(D))
141 E = E * P
145 H = E / 15
150 B = 24 + A - H
151 IF B > 24 THEN B = B - 24
160 G = A + H
161 IF G > 24 THEN G = G - 24
165 PRINT
170 PRINT "Местное звездное время восхода ="; B; "часов"
180 PRINT "Местное звездное время захода ="; G; "часов"
190 K = B - L
191 IF K > 24 THEN K = K - 24
192 IF K < 0 THEN K = K + 24
200 I = G - L
201 IF I > 24 THEN I = I - 24
202 IF I < 0 THEN I = I + 24
205 PRINT
210 PRINT "Введите день от января 0.0 текущего года и константу В"
211 PRINT
215 PRINT "Константа В для 1998 г.=17.369382, 1999 - 17.385297, 0 - 17.401211"
217 PRINT
220 INPUT "Количество дней =", W
221 INPUT "Константа В =", J
230 C = W * .0657098# - J
231 IF C < 0 THEN C = C + 24
235 PRINT
240 PRINT "То ="; C; " часов"
250 V = K - C
251 IF V < 0 THEN V = V + 24
260 X = I - C
261 IF X < 0 THEN X = X + 24
270 V = V * .99727
271 X = X * .99727
280 PRINT
300 PRINT "Время восхода по Гринвичу ="; V; " часов"
310 PRINT
315 PRINT
320 PRINT "Время захода по Гринвичу ="; X; "часов"
325 END
FUNCTION ACS (X)
ACS = 1.570796 - ASN(X)
END FUNCTION
FUNCTION ASN (X)
ASN = ATN(X / SQR(1 - X * X))
END FUNCTION
