17 lines
650 B
Text
17 lines
650 B
Text
10 DEF FN r(x)=x*PI/180
|
|
20 DEF FN d(x)=x*180/PI
|
|
30 INPUT "Enter latitude (degrees): ";latitude
|
|
40 INPUT "Enter longitude (degrees): ";longitude
|
|
50 INPUT "Enter legal meridian (degrees): ";meridian
|
|
60 PRINT "Latitude: ";latitude
|
|
70 PRINT "Longitude:";longitude
|
|
80 PRINT "Legal meridian: ";meridian
|
|
90 PRINT '" Sun Dial"
|
|
100 PRINT "Time hour angle hour line ang."
|
|
110 PRINT "________________________________"
|
|
120 FOR h=6 TO 18
|
|
130 LET hra=15*h-longitude+meridian-180
|
|
140 LET hla=FN d(ATN (SIN (FN r(latitude))*TAN (FN r(hra))))
|
|
150 IF ABS (hra)>90 THEN LET hla=hla+180*SGN (hra*latitude)
|
|
160 PRINT h;" ";hra;" ";hla
|
|
170 NEXT h
|