OPTION ANGLE RADIANS
PRINT
PRINT "夏至と冬至： 北緯と日照時間  東経と日出日没時刻"
PRINT "        明石： 北緯３５度   東経１３５度"
PRINT "        札幌： 北緯４３度   東経１４１度"
PRINT "        那覇： 北緯２６度   東経１２７度"
PRINT
INPUT PROMPT "夏至ならば 1 を、冬至ならば 2 を入力  ":W
INPUT PROMPT "北緯？ ":xx
INPUT PROMPT "東経？ ":qq
PRINT
IF W=2 THEN LET W=-1
LET P=W*23.4*PI/180
IF xx>66.6 THEN
   PRINT "日照時間 ＝";12+W*12 ;"時間"
   GOTO 10
END IF
IF xx<-66.6 THEN
   PRINT "日照時間 ＝";12-W*12 ;"時間"
   GOTO 10
END IF
LET x=xx*PI/180
LET q=(qq-135)/15   !１５度あたり１時間の時刻差
LET Y=ASIN(TAN(P)*TAN(x))
LET Z=12+12*2*Y/PI
LET Z1=12-Z/2-q    ! 日の出
LET Z2=12+Z/2-q    ! 日の入り
PRINT "日照時間 ＝  ";INT(Z);"時";INT(MOD(Z,1)*60);"分間"
PRINT "日出時刻 ＝  ";INT(Z1);"時";INT(MOD(Z1,1)*60);"分"
PRINT "日没時刻 ＝  ";INT(Z2);"時";INT(MOD(Z2,1)*60);"分"
10 PRINT
END