Horizontal sundial calculations: Difference between revisions
Content added Content deleted
(added FreeBASIC) |
No edit summary |
||
Line 2,585: | Line 2,585: | ||
<as above> |
<as above> |
||
</pre> |
</pre> |
||
=={{header|ZX Spectrum Basic}}== |
|||
{{trans|ERRE}} |
|||
<lang zxbasic>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</lang> |