2 'PC-125*,126*,136*
4 '座標変換プログラム　赤道座標→地平座標
6'
10 CLEAR : DIM Z(6)
20 USING "####.#"
30 PRINT "RA.DC>A.H"
40 INPUT "ﾈﾝ=";Y
50 INPUT "ﾂｷ=";Z(4)
60 INPUT "ﾆﾁ=";N
70 T=INT((Y-1901)*365.25+.1)
80 DATA 0,31,59,90,120,151
90 DATA 181,212,243,273,304,334
100 RESTORE 80
110 FOR I=1 TO Z(4)
120   READ R
130 NEXT I
140 Z(6)=Y-1900
150 IF INT (Y/4)*4<>Z(6) THEN 170
160 IF Z(4)=>3 THEN R=R+1
170 Z(3)=T+R+N
180 INPUT "TIME-H=";J
190 INPUT "TIME-M=";O
200 Y=(Z(3)+(J+O/60)/24-36160.5)/36525
210 '東経＆北緯
220 P=135 : Q=35　'明石市
230 INPUT "RA=";U
240 INPUT "DC=";K
250 Z(5)=(J+O/60)-9
260 Z(2)=2400.0513*Y+6.69736+Z(5)
270 Z(2)=Z(2)+2400
280 Z=Z(2)-INT(Z(2)/24)*24
290 H=(Z*15+P-U*15)
300 W=SIN H :G=COS H
310 V=SIN K :F=COS K
320 X=SIN Q :M=COS Q
330 A=-F*W
340 B=V*M-F*X*G
350 C=V*X+F*M*G
360 D=ATN(C/SQR(A^2+B^2))
370 E=ATN(A/B)
380 IF B<0 THEN LET E=E+180
390 IF E<0 THEN LET E=E+360
400 E=E+180
410 IF E>360 THEN LET E=E-360
420 PRIN "A=";E;"  H=";D
430 GOTO 180