Horizontal sundial calculations: Difference between revisions

← Older edit

Horizontal sundial calculations (view source)

Revision as of 20:53, 9 December 2023

3,903 bytes added , 5 months ago

m

→‎{{header|Wren}}: Minor tidy

PureFox

9,476

edits

Revision as of 02:55, 4 February 2023 (view source) Peak (talk \| contribs) (→‎{{header\|jq}}) ← Older edit		Latest revision as of 20:53, 9 December 2023 (view source) PureFox (talk \| contribs) m (→‎{{header\|Wren}}: Minor tidy)
(7 intermediate revisions by 4 users not shown)
Line 190: HR= +6; HRA= +90.500; HLA= +84.225 </pre> =={{header\|Arturo}}== {{trans\|Nim}} <syntaxhighlight lang="arturo">degToRad: function [deg]-> deg * pi // 180 radToDeg: function [rad]-> rad * 180 // pi lat: to :floating input "Enter latitude => " lng: to :floating input "Enter longitude => " med: to :floating input "Enter legal meridian => " print "" slat: sin degToRad lat print " sine of latitude: " ++ to :string .format:".3f" slat print " diff longitude: " ++ to :string .format:".3f" lng-med print "" print "Hour, sun hour angle, dial hour line angle from 6am to 6pm" loop (neg 6)..6 'h [ hra: med + (to :floating 15h) - lng hla: radToDeg atan slat tan degToRad hra print "HR=" ++ (to :string .format:"3d" h) ++ "; " ++ "HRA=" ++ (to :string .format:"7.3f" hra) ++ "; " ++ "HLA=" ++ (to :string .format:"7.3f" hla) ]</syntaxhighlight> {{out}} <pre>Enter latitude => -4.95 Enter longitude => -150.5 Enter legal meridian => -150 sine of latitude: -0.086 diff longitude: -0.500 Hour, sun hour angle, dial hour line angle from 6am to 6pm HR= -6; HRA=-89.500; HLA= 84.225 HR= -5; HRA=-74.500; HLA= 17.283 HR= -4; HRA=-59.500; HLA= 8.334 HR= -3; HRA=-44.500; HLA= 4.847 HR= -2; HRA=-29.500; HLA= 2.795 HR= -1; HRA=-14.500; HLA= 1.278 HR= 0; HRA= 0.500; HLA= -0.043 HR= 1; HRA= 15.500; HLA= -1.371 HR= 2; HRA= 30.500; HLA= -2.910 HR= 3; HRA= 45.500; HLA= -5.018 HR= 4; HRA= 60.500; HLA= -8.671 HR= 5; HRA= 75.500; HLA=-18.451 HR= 6; HRA= 90.500; HLA= 84.225</pre> =={{header\|AutoHotkey}}== Line 294 ⟶ 345: =={{header\|BASIC}}== ==={{header\|~~Minimal~~ANSI BASIC}}=== {{trans\|ALGOL-68\|ANSI BASIC has a possibility to use degrees in calculations.}} ~~{{works with\|QBasic}}~~ {{works with\|Decimal BASIC}} ~~{{trans\|ALGOL-68}}~~ <syntaxhighlight lang="~~gwbasic~~basic"> 10100 REM Horizontal sundial calculations 110 INPUT PROMPT "Enter latitude => ": Lat ~~20 DEF FNM(X) = INT(X1000+0.5)/1000~~ 30120 ~~PRINT~~INPUT PROMPT "Enter ~~latitude~~longitude => ";: Lng 130 INPUT PROMPT "Enter legal meridian => ": Ref ~~40 INPUT L~~ ~~50 PRINT "Enter longitude";~~ ~~60 INPUT L1~~ ~~70 PRINT "Enter legal meridian";~~ ~~80 INPUT R~~ ~~90 PRINT~~ ~~100 LET P = 4ATN(1)~~ ~~110 LET S1 = SIN(LP/180)~~ ~~120 PRINT " sine of latitude:"; S1~~ ~~130 PRINT " diff longitude: "; FNM(L1-R)~~ 140 PRINT 150 OPTION ANGLE DEGREES ~~150 PRINT "Hour, sun hour angle, dial hour line angle from 6am to 6pm"~~ 160 ~~FOR~~LET HSlat = ~~-6 TO 6~~SIN(Lat) 170 PRINT " sine of latitude: "; Slat ~~170 LET A1 = 15H~~ 180 PRINT " diff longitude: "; Lng - Ref ~~180 REM Correct for longitude difference:~~ 190 PRINT ~~190 LET A1 = A1-(L1-R)~~ 200 PRINT "Hour, sun hour angle, dial hour line angle from 6am to 6pm" ~~200 LET A2 = ATN(S1TAN(A1P/180))180/P~~ 210 ~~PRINT~~FOR ~~"HR~~Hour ="; H;-6 TO 6 220 LET HourAngle = 15 Hour ~~220 PRINT TAB(9); "HRA ="; FNM(A1);~~ 230 LET HourAngle = HourAngle - (Lng - Ref) ! correct for longitude difference ~~230 PRINT TAB(24); "HLA ="; FNM(A2)~~ 240 LET HourLineAngle = ATN(Slat * TAN(HourAngle)) ~~240 NEXT H~~ 250 PRINT USING "HR=###; HRA=####.###; HLA=####.###": Hour, HourAngle, HourLineAngle ~~250 END~~ 260 NEXT Hour 270 END </syntaxhighlight> {{out}} <pre> Enter latitude => -4.95 Enter longitude => -150.5 Enter legal meridian => -150 sine of latitude: -8.62863657979234E-2 diff longitude: -.5 Hour, sun hour angle, dial hour line angle from 6am to 6pm HR= -6; HRA= -89.500; HLA= 84.225 HR= -5; HRA= -74.500; HLA= 17.283 HR= -4; HRA= -59.500; HLA= 8.334 HR= -3; HRA= -44.500; HLA= 4.847 HR= -2; HRA= -29.500; HLA= 2.795 HR= -1; HRA= -14.500; HLA= 1.278 HR= 0; HRA= .500; HLA= -.043 HR= 1; HRA= 15.500; HLA= -1.371 HR= 2; HRA= 30.500; HLA= -2.910 HR= 3; HRA= 45.500; HLA= -5.018 HR= 4; HRA= 60.500; HLA= -8.671 HR= 5; HRA= 75.500; HLA= -18.451 HR= 6; HRA= 90.500; HLA= 84.225 </pre> ==={{header\|BASIC256}}=== Line 360 ⟶ 428: next hora end subroutine</syntaxhighlight> ==={{header\|BBC BASIC}}=== {{works with\|BBC BASIC for Windows}} <syntaxhighlight lang="bbcbasic"> INSTALL @lib$+"FNUSING" INPUT "Enter latitude (degrees) : " latitude INPUT "Enter longitude (degrees) : " longitude INPUT "Enter legal meridian (degrees): " meridian PRINT '" Time", "Sun hour angle", "Dial hour line angle" FOR hour = 6 TO 18 hra = 15hour - longitude + meridian - 180 hla = DEG(ATN(SIN(RAD(latitude)) TAN(RAD(hra)))) IF ABS(hra) > 90 hla += 180 * SGN(hra * latitude) PRINT FNusing("##.##", hour), FNusing(" ####.### ", hra), FNusing(" ####.###", hla) NEXT hour</syntaxhighlight> {{out}} (note the correct '''negative''' value for time 18:00) <pre>Enter latitude (degrees) : -4.95 Enter longitude (degrees) : -150.5 Enter legal meridian (degrees): -150.0 Time Sun hour angle Dial hour line angle 6.00 -89.500 84.225 7.00 -74.500 17.283 8.00 -59.500 8.334 9.00 -44.500 4.847 10.00 -29.500 2.795 11.00 -14.500 1.278 12.00 0.500 -0.043 13.00 15.500 -1.371 14.00 30.500 -2.910 15.00 45.500 -5.018 16.00 60.500 -8.671 17.00 75.500 -18.451 18.00 90.500 -95.775</pre> ==={{header\|Chipmunk Basic}}=== {{works with\|Chipmunk Basic\|3.6.4}} <syntaxhighlight lang="qbasic">10 CLS 20 DEF FNR(x) = x * PI/180 30 DEF FND(x) = x * 180/PI 40 INPUT " Enter latitude (degrees): ";latitude 50 INPUT " Enter longitude (degrees): ";longitude 60 INPUT "Enter legal meridian (degrees): ";meridian 70 PRINT 80 PRINT " Time Sun hour angle Dial hour line angle" 90 FOR h = 6 TO 18 100 hra = 15h-longitude+meridian-180 110 hla = FN D(ARCTAN(SIN(FN R(latitude))TAN(FN R(hra)))) 120 IF ABS(hra) > 90 THEN hla = hla+180SGN(hralatitude) 130 PRINT USING "##.##";h;TAB (10) 140 PRINT USING "####.###";hra;TAB (27) 150 PRINT USING "####.###";hla 160 NEXT h 170 END</syntaxhighlight> {{out}} <pre>Same as FreeBASIC entry.</pre> ==={{header\|FreeBASIC}}=== {{trans\|BBC BASIC}} <syntaxhighlight lang="freebasic">' version 04-11-2016 ' compile with: fbc -s console #Macro deg2rad (x) (x) * Atn(1) / 45 #EndMacro #Macro rad2deg (x) (x) * 45 / Atn(1) #EndMacro ' ------=< MAIN >=------ Dim As Double latitude, longitude, meridian, hra, hla Dim As ULong h Input " Enter latitude (degrees): ", latitude Input " Enter longitude (degrees): ", longitude Input "Enter legal meridian (degrees): ", meridian Print Print " Time Sun hour angle Dial hour line angle" For h = 6 To 18 hra = h * 15 - longitude + meridian - 180 hla = rad2deg(Atn(Sin(deg2rad(latitude)) * Tan(deg2rad(hra)))) If Abs(hra) > 90 Then hla += 180 * Sgn(hra * latitude) Print Using "##.## ####.### ####.###"; h; hra; hla Next ' empty keyboard buffer While InKey <> "" : Wend Print : Print "hit any key to end program" Sleep End</syntaxhighlight> {{out}} <pre> Enter latitude (degrees): -4.95 Enter longitude (degrees): -150.5 Enter legal meridian (degrees): -150 Time Sun hour angle Dial hour line angle 6.00 -89.500 84.225 7.00 -74.500 17.283 8.00 -59.500 8.334 9.00 -44.500 4.847 10.00 -29.500 2.795 11.00 -14.500 1.278 12.00 0.500 -0.043 13.00 15.500 -1.371 14.00 30.500 -2.910 15.00 45.500 -5.018 16.00 60.500 -8.671 17.00 75.500 -18.451 18.00 90.500 -95.775</pre> ==={{header\|FutureBasic}}=== <syntaxhighlight lang="futurebasic">window 1 def fn rad2deg( theta as double ) as double = theta * 180 / pi def fn deg2rad( theta as double ) as double = theta * pi / 180 local fn SolarHourAngle( latitude as double, longitude as double, meridian as double ) long hour double hra, hla, t CFStringRef ap print "Latitude = "; latitude; chr$(13); "Longitude = "; longitude; chr$(13); "Meridian = "; meridian print : print "sine of latitude: "; sin(latitude * pi / 180 ); chr$(13); " diff longitude: "; longitude - meridian print : print "Time", "Sun hour angle", "Dial hour line angle" for hour = 6 to 18 hra = ( 15 * hour ) - longitude + meridian - 180 hla = fn rad2deg( atn( sin( fn deg2rad( latitude ) ) * tan( fn deg2rad( hra ) ))) if abs( hra ) > 90 then hla = hla + 180 * sgn( hra * latitude ) if hour > 12 then t = hour - 12 : ap = @" a.m." else t = hour : ap = @" p.m." print using "##"; t; ap, using "####.##"; hra, using "####.###"; hla next hour end fn fn SolarHourAngle( -4.95, -150.5, -150.0 ) HandleEvents</syntaxhighlight> Output: <pre> Latitude = -4.95 Longitude = -150.5 Meridian = -150 sine of latitude: -0.0862863658 diff longitude: -0.5 Time Sun hour angle Dial hour line angle 6 p.m. -89.50 84.225 7 p.m. -74.50 17.283 8 p.m. -59.50 8.334 9 p.m. -44.50 4.847 10 p.m. -29.50 2.795 11 p.m. -14.50 1.278 12 p.m. 0.50 -0.043 1 a.m. 15.50 -1.371 2 a.m. 30.50 -2.910 3 a.m. 45.50 -5.018 4 a.m. 60.50 -8.671 5 a.m. 75.50 -18.451 6 a.m. 90.50 -95.775 </pre> ==={{header\|GW-BASIC}}=== {{trans\|ALGOL 68}} {{works with\|PC-BASIC\|any}} <syntaxhighlight lang="qbasic"> 10 ' Horizontal sundial calculations 20 PRINT "Enter latitude => "; 30 INPUT LAT 40 PRINT "Enter longitude => "; 50 INPUT LNG 60 PRINT "Enter legal meridian => "; 70 INPUT REF 80 PRINT 90 LET PI = 4 * ATN(1) 100 LET SLAT = SIN(LAT * PI / 180) 110 PRINT " sine of latitude: "; USING "#.##^^^^"; SLAT 120 PRINT " diff longitude: "; USING "####.###"; LNG - REF 130 PRINT 140 PRINT "Hour, sun hour angle, dial hour line angle from 6am to 6pm" 150 FOR H% = -6 TO 6 160 LET HRA = 15 * H% 170 LET HRA = HRA - (LNG - REF): ' correct for longitude difference 180 LET HLA = ATN(SLAT * TAN(HRA * PI / 180)) * 180 / PI 190 PRINT "HR="; USING "+##"; H%; 200 PRINT "; HRA="; USING "+###.###"; HRA; 210 PRINT "; HLA="; USING "+###.###"; HLA 220 NEXT H% 230 END </syntaxhighlight> {{out}} <pre> Enter latitude => ? -4.95 Enter longitude => ? -150.5 Enter legal meridian => ? -150 sine of latitude: -.86E-01 diff longitude: -0.500 Hour, sun hour angle, dial hour line angle from 6am to 6pm HR= -6; HRA= -89.500; HLA= +84.225 HR= -5; HRA= -74.500; HLA= +17.283 HR= -4; HRA= -59.500; HLA= +8.334 HR= -3; HRA= -44.500; HLA= +4.847 HR= -2; HRA= -29.500; HLA= +2.795 HR= -1; HRA= -14.500; HLA= +1.278 HR= +0; HRA= +0.500; HLA= -0.043 HR= +1; HRA= +15.500; HLA= -1.371 HR= +2; HRA= +30.500; HLA= -2.910 HR= +3; HRA= +45.500; HLA= -5.018 HR= +4; HRA= +60.500; HLA= -8.671 HR= +5; HRA= +75.500; HLA= -18.451 HR= +6; HRA= +90.500; HLA= +84.225 </pre> ==={{header\|IS-BASIC}}=== <syntaxhighlight lang="is-basic">100 PROGRAM "SunDial.bas" 110 OPTION ANGLE RADIANS 120 TEXT 80 130 LET DR=PI/180:LET RD=180/PI 140 INPUT PROMPT "Enter latitude: ":LAT 150 INPUT PROMPT "Enter longitude: ":LNG 160 INPUT PROMPT "Enter legal meridian: ":REF 170 LET S=SIN(LATDR) 180 PRINT :PRINT "Sine of latitude :";S 190 PRINT "Diff longitude: ";LNG-REF 200 PRINT " Hour,",,"sun hour angle, dial hour line angle from 6am to 6 pm" 210 FOR H=6 TO 18 220 LET HRA=15H-LNG+REF 230 LET HLA=ATN(STAN(HRADR))RD 240 PRINT "HR = ";H,"HRA =";HRA,"HLA =";HLA 250 NEXT</syntaxhighlight> ==={{header\|Liberty BASIC}}=== Based on Algol & BBC BASIC versions. Note Liberty BASIC works in radians. <syntaxhighlight lang="lb">global pi pi = 3.14159265 input "Enter latitude (degrees) : "; latitude ' -4.95 input "Enter longitude (degrees) : "; longitude ' -150.5 input "Enter legal meridian (degrees): "; meridian ' -150.0 print print "Time Sun hour angle Dial hour line angle" for hour = 6 TO 18 hra = 15 hour - longitude + meridian - 180 hla = rad2deg(atn(sin(deg2rad(latitude)) * tan(deg2rad(hra)))) if abs(hra) > 90 then hla = hla + 180 * sgn(hra * latitude) print using( "##.##", hour), using("####.### ", hra), using("####.###", hla) next hour end function rad2deg(theta) rad2deg = theta * 180 / pi end function function deg2rad(theta) deg2rad = theta * pi / 180 end function function sgn(x) if x > 0 then sgn = 1 else sgn = -1 end function</syntaxhighlight> {{out}} <pre> Enter latitude (degrees) : -4.95 Enter longitude (degrees) : -150.5 Enter legal meridian (degrees): -150.0 Time Sun hour angle Dial hour line angle 6.00 -89.500 84.225 7.00 -74.500 17.283 8.00 -59.500 8.334 9.00 -44.500 4.847 10.00 -29.500 2.795 11.00 -14.500 1.278 12.00 0.500 -0.043 13.00 15.500 -1.371 14.00 30.500 -2.910 15.00 45.500 -5.018 16.00 60.500 -8.671 17.00 75.500 -18.451 18.00 90.500 -95.775 </pre> ==={{header\|Microsoft Small Basic}}=== {{trans\|ALGOL-68}} <syntaxhighlight lang="microsoftsmallbasic"> TextWindow.Write("Enter latitude => ") lat = TextWindow.ReadNumber() TextWindow.Write("Enter longitude => ") lng = TextWindow.ReadNumber() TextWindow.Write("Enter legal meridian => ") ref = TextWindow.ReadNumber() sLat = Math.Sin(Math.GetRadians(lat)) TextWindow.WriteLine("") TextWindow.Write(" sine of latitude: ") TextWindow.WriteLine(Math.Round(sLat * 10000) / 10000) TextWindow.Write(" diff longitude: ") TextWindow.WriteLine(lng - ref) TextWindow.WriteLine("") TextWindow.WriteLine("Hour, sun hour angle, dial hour line angle from 6am to 6pm") For hour = -6 To 6 hourAngle = 15 * hour hourAngle = hourAngle - (lng - ref) ' correct for longitude difference hourLineAngle = math.GetDegrees(Math.ArcTan(sLat * Math.Tan(Math.GetRadians(hourAngle)))) TextWindow.Write("HR=") TextWindow.CursorLeft = 3 + (3 - Text.GetLength(hour)) TextWindow.Write(hour) TextWindow.Write("; HRA=") TextWindow.CursorLeft = 12 + (6 - Text.GetLength(hourAngle)) TextWindow.Write(hourAngle) TextWindow.Write("; HLA=") TextWindow.CursorLeft = 24 + (4 - Text.GetLength(Math.Floor(hourLineAngle))) TextWindow.Write(Math.Round(hourLineAngle * 1000) / 1000) TextWindow.WriteLine("") EndFor </syntaxhighlight> {{out}} Enter latitude => -4.95 Enter longitude => -150.5 Enter legal meridian => -150 sine of latitude: -0.0863 diff longitude: -0.5 Hour, sun hour angle, dial hour line angle from 6am to 6pm HR= -6; HRA= -89.5; HLA= 84.225 HR= -5; HRA= -74.5; HLA= 17.283 HR= -4; HRA= -59.5; HLA= 8.334 HR= -3; HRA= -44.5; HLA= 4.847 HR= -2; HRA= -29.5; HLA= 2.795 HR= -1; HRA= -14.5; HLA= 1.278 HR= 0; HRA= 0.5; HLA= -0.043 HR= 1; HRA= 15.5; HLA= -1.371 HR= 2; HRA= 30.5; HLA= -2.91 HR= 3; HRA= 45.5; HLA= -5.018 HR= 4; HRA= 60.5; HLA= -8.671 HR= 5; HRA= 75.5; HLA= -18.451 HR= 6; HRA= 90.5; HLA= 84.225 ==={{header\|Minimal BASIC}}=== {{works with\|QBasic}} {{trans\|ALGOL-68}} <syntaxhighlight lang="gwbasic"> 10 REM Horizontal sundial calculations 20 DEF FNM(X) = INT(X1000+0.5)/1000 30 PRINT "Enter latitude"; 40 INPUT L 50 PRINT "Enter longitude"; 60 INPUT L1 70 PRINT "Enter legal meridian"; 80 INPUT R 90 PRINT 100 LET P = 4ATN(1) 110 LET S1 = SIN(LP/180) 120 PRINT " sine of latitude:"; S1 130 PRINT " diff longitude: "; FNM(L1-R) 140 PRINT 150 PRINT "Hour, sun hour angle, dial hour line angle from 6am to 6pm" 160 FOR H = -6 TO 6 170 LET A1 = 15H 180 REM Correct for longitude difference: 190 LET A1 = A1-(L1-R) 200 LET A2 = ATN(S1TAN(A1P/180))180/P 210 PRINT "HR ="; H; 220 PRINT TAB(9); "HRA ="; FNM(A1); 230 PRINT TAB(24); "HLA ="; FNM(A2) 240 NEXT H 250 END </syntaxhighlight> ==={{header\|PureBasic}}=== {{trans\|ALGOL 68}} <syntaxhighlight lang="purebasic">If OpenConsole() Define.f lat, slat, lng, ref Define.i h Print("Enter latitude => "): lat=ValF(Input()) Print("Enter longitude => "): lng=ValF(Input()) Print("Enter legal meridian => "): ref=ValF(Input()) PrintN("") slat=Sin(lat2#PI/360) PrintN(" sine of latitude: "+StrF(slat,3)) PrintN(" diff longitude: "+StrF((lng-ref),3)+#CRLF$) PrintN("Hour, sun hour angle, dial hour line angle from 6am to 6pm") For h=-6 To 6 Define.f hra, hla hra=15h hra=hra-(lng-ref) hla=ATan(slatTan(hra2#PI/360))360/(2#PI) PrintN("HR="+RSet(Str(h),3)+"; HRA="+RSet(StrF(hra,3),7)+"; HLA="+RSet(StrF(hla,3),7)) Next EndIf</syntaxhighlight> {{out}} <pre>Enter latitude => -4.95 Enter longitude => -150.5 Enter legal meridian => -150 sine of latitude: -0.086 diff longitude: -0.500 Hour, sun hour angle, dial hour line angle from 6am to 6pm HR= -6; HRA=-89.500; HLA= 84.225 HR= -5; HRA=-74.500; HLA= 17.283 HR= -4; HRA=-59.500; HLA= 8.334 HR= -3; HRA=-44.500; HLA= 4.847 HR= -2; HRA=-29.500; HLA= 2.795 HR= -1; HRA=-14.500; HLA= 1.278 HR= 0; HRA= 0.500; HLA= -0.043 HR= 1; HRA= 15.500; HLA= -1.371 HR= 2; HRA= 30.500; HLA= -2.910 HR= 3; HRA= 45.500; HLA= -5.018 HR= 4; HRA= 60.500; HLA= -8.671 HR= 5; HRA= 75.500; HLA=-18.451 HR= 6; HRA= 90.500; HLA= 84.225 </pre> ==={{header\|Run BASIC}}=== <syntaxhighlight lang="runbasic">global pi pi = 22 / 7 print "Enter latitude (degrees) : "; :input latitude ' -4.95 print "Enter longitude (degrees) : "; :input longitude ' -150.5 print "Enter legal meridian (degrees): "; :input meridian ' -150.0 print print "Time Sun hour angle Dial hour line angle" for hour = 6 TO 18 hra = (15 hour) - longitude + meridian -180 hla =rad2deg( atn( sin( deg2rad( latitude)) tan( deg2rad( hra)))) if abs( hra) >90 then hla =hla +180 sgn( hra latitude) print using( "##", hour);" ";using("####.##", hra);" ";using("####.###", hla) next hour function rad2deg( theta) rad2deg =theta 180 /pi end function function deg2rad( theta) deg2rad =theta pi /180 end function function sgn( x) if x >0 then sgn =1 else sgn =-1 end function end</syntaxhighlight> {{out}} <pre> Enter latitude: -4.95 Enter longitude: -150.5 Enter legal meridian: -150 Time Sun hour angle Dial hour line angle 6 -89.50 84.606 7 -74.50 17.316 8 -59.50 8.342 9 -44.50 4.850 10 -29.50 2.796 11 -14.50 1.279 12 0.50 -0.043 13 15.50 -1.371 14 30.50 -2.911 15 45.50 -5.021 16 60.50 -8.680 17 75.50 -18.488 18 90.50 -96.224</pre> ==={{header\|True BASIC}}=== Line 383 ⟶ 927: NEXT hora END</syntaxhighlight> ==={{header\|XBasic}}=== {{trans\|ALGOL-68}} {{works with\|Windows XBasic}} <syntaxhighlight lang="xbasic"> PROGRAM "sundial" VERSION "0.0001" IMPORT "xma" DECLARE FUNCTION Entry() FUNCTION Entry() lat! = SINGLE(INLINE$("Enter latitude => ")) lng! = SINGLE(INLINE$("Enter longitude => ")) ref! = SINGLE(INLINE$("Enter legal meridian => ")) PRINT slat! = SIN(lat! $$PI / 180.0) PRINT " sine of latitude: "; FORMAT$("#.##^^^^", slat!) PRINT " diff longitude: "; FORMAT$("#.###", lng! - ref!) PRINT PRINT "Hour, sun hour angle, dial hour line angle from 6am to 6pm" FOR hour@ = -6 TO 6 hourAngle! = 15 * hour@ hourAngle! = hourAngle! - (lng! - ref!) ' correct for longitude difference hourLineAngle! = ATAN(slat! * TAN(hourAngle! * $$PI / 180.0)) * 180.0 / $$PI PRINT "HR="; FORMAT$("###", hour@); PRINT "; HRA="; FORMAT$("####.###", hourAngle!); PRINT "; HLA="; FORMAT$("####.###", hourLineAngle!) NEXT hour@ END FUNCTION END PROGRAM </syntaxhighlight> {{out}} <pre> Enter latitude => -4.95 Enter longitude => -150.5 Enter legal meridian => -150 sine of latitude: -8.63E-02 diff longitude: -0.500 Hour, sun hour angle, dial hour line angle from 6am to 6pm HR= -6; HRA= -89.500; HLA= 84.225 HR= -5; HRA= -74.500; HLA= 17.283 HR= -4; HRA= -59.500; HLA= 8.334 HR= -3; HRA= -44.500; HLA= 4.847 HR= -2; HRA= -29.500; HLA= 2.795 HR= -1; HRA= -14.500; HLA= 1.278 HR= 0; HRA= 0.500; HLA= -0.043 HR= 1; HRA= 15.500; HLA= -1.371 HR= 2; HRA= 30.500; HLA= -2.910 HR= 3; HRA= 45.500; HLA= -5.018 HR= 4; HRA= 60.500; HLA= -8.671 HR= 5; HRA= 75.500; HLA= -18.451 HR= 6; HRA= 90.500; HLA= 84.225 </pre> ==={{header\|Yabasic}}=== Line 414 ⟶ 1,015: end sub</syntaxhighlight> ==={{header\|ZX Spectrum Basic}}=== =={{~~header~~trans\|~~BBC BASIC~~ERRE}}== <syntaxhighlight lang="zxbasic">10 DEF FN r(x)=xPI/180 ~~{{works with\|BBC BASIC for Windows}}~~ 20 DEF FN d(x)=x180/PI ~~<syntaxhighlight lang="bbcbasic"> INSTALL @lib$+"FNUSING"~~ 30 INPUT "Enter latitude (degrees): ";latitude 40 INPUT "Enter ~~latitude~~longitude (degrees) : " ~~latitude~~;longitude 50 INPUT "Enter ~~longitude~~legal meridian (degrees) : " ~~longitude~~;meridian 60 PRINT "Latitude: ";latitude ~~INPUT "Enter legal meridian (degrees): " meridian~~ 70 PRINT "Longitude:";longitude 80 PRINT "Legal meridian: ";meridian ~~PRINT '" Time", "Sun hour angle", "Dial hour line angle"~~ 90 PRINT '" Sun Dial" 100 PRINT "Time hour angle ~~FOR~~ hour ~~= 6 TO~~line 18ang." 110 PRINT "________________________________" ~~hra = 15hour - longitude + meridian - 180~~ 120 FOR h=6 TO 18 ~~hla = DEG(ATN(SIN(RAD(latitude)) TAN(RAD(hra))))~~ 130 LET hra=15h-longitude+meridian-180 ~~IF ABS(hra) > 90 hla += 180 SGN(hra * latitude)~~ 140 LET hla=FN d(ATN (SIN (FN r(latitude))TAN (FN r(hra)))) ~~PRINT FNusing("##.##", hour), FNusing(" ####.### ", hra), FNusing(" ####.###", hla)~~ 150 IF ABS (hra)>90 THEN LET hla=hla+180SGN (hralatitude) ~~NEXT hour</syntaxhighlight>~~ 160 PRINT h;" ";hra;" ";hla ~~{{out}}~~ 170 NEXT h</syntaxhighlight> ~~(note the correct '''negative''' value for time 18:00)~~ ~~<pre>Enter latitude (degrees) : -4.95~~ ~~Enter longitude (degrees) : -150.5~~ ~~Enter legal meridian (degrees): -150.0~~ ~~Time Sun hour angle Dial hour line angle~~ ~~6.00 -89.500 84.225~~ ~~7.00 -74.500 17.283~~ ~~8.00 -59.500 8.334~~ ~~9.00 -44.500 4.847~~ ~~10.00 -29.500 2.795~~ ~~11.00 -14.500 1.278~~ ~~12.00 0.500 -0.043~~ ~~13.00 15.500 -1.371~~ ~~14.00 30.500 -2.910~~ ~~15.00 45.500 -5.018~~ ~~16.00 60.500 -8.671~~ ~~17.00 75.500 -18.451~~ ~~18.00 90.500 -95.775</pre>~~ =={{header\|C}}== Line 802 ⟶ 1,385: =={{header\|EasyLang}}== <syntaxhighlight lang="text">~~func getn s$ . v .~~ func ~~write~~getn s$ . write s$ ~~v = number input~~ v = number input ~~print v~~ print v return v . ~~call~~lat = getn "Enter latitude: " ~~lat~~ ~~call~~lng = getn "Enter longitude: " ~~lng~~ ~~call~~merid = getn "Enter legal meridian: " ~~merid~~ slat = sin lat diff = lng - merid Line 818 ⟶ 1,403: print "Hour\tSun hour angle\tDial hour line angle" for h = -6 to 6 hra = 15 h - diff hla = atan2 (slat * sin hra) cos hra print h + 12 & "\t" & hra & "\t\t" & hla . ~~.</syntaxhighlight>~~ </syntaxhighlight> =={{header\|ERRE}}== Line 1,129 ⟶ 1,715: end program SunDial</syntaxhighlight> ~~=={{header\|FreeBASIC}}==~~ ~~{{trans\|BBC BASIC}}~~ ~~<syntaxhighlight lang="freebasic">' version 04-11-2016~~ ~~' compile with: fbc -s console~~ ~~#Macro deg2rad (x)~~ ~~(x) * Atn(1) / 45~~ ~~#EndMacro~~ ~~#Macro rad2deg (x)~~ ~~(x) * 45 / Atn(1)~~ ~~#EndMacro~~ ~~' ------=< MAIN >=------~~ ~~Dim As Double latitude, longitude, meridian, hra, hla~~ ~~Dim As ULong h~~ ~~Input " Enter latitude (degrees): ", latitude~~ ~~Input " Enter longitude (degrees): ", longitude~~ ~~Input "Enter legal meridian (degrees): ", meridian~~ ~~Print~~ ~~Print " Time Sun hour angle Dial hour line angle"~~ ~~For h = 6 To 18~~ ~~hra = h * 15 - longitude + meridian - 180~~ ~~hla = rad2deg(Atn(Sin(deg2rad(latitude)) * Tan(deg2rad(hra))))~~ ~~If Abs(hra) > 90 Then hla += 180 * Sgn(hra * latitude)~~ ~~Print Using "##.## ####.### ####.###"; h; hra; hla~~ ~~Next~~ ~~' empty keyboard buffer~~ ~~While InKey <> "" : Wend~~ ~~Print : Print "hit any key to end program"~~ ~~Sleep~~ ~~End</syntaxhighlight>~~ ~~{{out}}~~ ~~<pre> Enter latitude (degrees): -4.95~~ ~~Enter longitude (degrees): -150.5~~ ~~Enter legal meridian (degrees): -150~~ ~~Time Sun hour angle Dial hour line angle~~ ~~6.00 -89.500 84.225~~ ~~7.00 -74.500 17.283~~ ~~8.00 -59.500 8.334~~ ~~9.00 -44.500 4.847~~ ~~10.00 -29.500 2.795~~ ~~11.00 -14.500 1.278~~ ~~12.00 0.500 -0.043~~ ~~13.00 15.500 -1.371~~ ~~14.00 30.500 -2.910~~ ~~15.00 45.500 -5.018~~ ~~16.00 60.500 -8.671~~ ~~17.00 75.500 -18.451~~ ~~18.00 90.500 -95.775</pre>~~ ~~=={{header\|FutureBasic}}==~~ ~~<syntaxhighlight lang="futurebasic">window 1~~ ~~def fn rad2deg( theta as double ) as double = theta * 180 / pi~~ ~~def fn deg2rad( theta as double ) as double = theta * pi / 180~~ ~~local fn SolarHourAngle( latitude as double, longitude as double, meridian as double )~~ ~~long hour~~ ~~double hra, hla, t~~ ~~CFStringRef ap~~ ~~print "Latitude = "; latitude; chr$(13); "Longitude = "; longitude; chr$(13); "Meridian = "; meridian~~ ~~print : print "sine of latitude: "; sin(latitude * pi / 180 ); chr$(13); " diff longitude: "; longitude - meridian~~ ~~print : print "Time", "Sun hour angle", "Dial hour line angle"~~ ~~for hour = 6 to 18~~ ~~hra = ( 15 * hour ) - longitude + meridian - 180~~ ~~hla = fn rad2deg( atn( sin( fn deg2rad( latitude ) ) * tan( fn deg2rad( hra ) )))~~ ~~if abs( hra ) > 90 then hla = hla + 180 * sgn( hra * latitude )~~ ~~if hour > 12 then t = hour - 12 : ap = @" a.m." else t = hour : ap = @" p.m."~~ ~~print using "##"; t; ap, using "####.##"; hra, using "####.###"; hla~~ ~~next hour~~ ~~end fn~~ ~~fn SolarHourAngle( -4.95, -150.5, -150.0 )~~ ~~HandleEvents</syntaxhighlight>~~ ~~Output:~~ ~~<pre>~~ ~~Latitude = -4.95~~ ~~Longitude = -150.5~~ ~~Meridian = -150~~ ~~sine of latitude: -0.0862863658~~ ~~diff longitude: -0.5~~ ~~Time Sun hour angle Dial hour line angle~~ ~~6 p.m. -89.50 84.225~~ ~~7 p.m. -74.50 17.283~~ ~~8 p.m. -59.50 8.334~~ ~~9 p.m. -44.50 4.847~~ ~~10 p.m. -29.50 2.795~~ ~~11 p.m. -14.50 1.278~~ ~~12 p.m. 0.50 -0.043~~ ~~1 a.m. 15.50 -1.371~~ ~~2 a.m. 30.50 -2.910~~ ~~3 a.m. 45.50 -5.018~~ ~~4 a.m. 60.50 -8.671~~ ~~5 a.m. 75.50 -18.451~~ ~~6 a.m. 90.50 -95.775~~ ~~</pre>~~ =={{header\|Go}}== Line 1,296 ⟶ 1,773: 5 75.500 -18.451 6 90.500 -95.775 ~~</pre>~~ ~~=={{header\|GW-BASIC}}==~~ ~~{{trans\|ALGOL 68}}~~ ~~{{works with\|PC-BASIC\|any}}~~ ~~<syntaxhighlight lang="qbasic">~~ ~~10 ' Horizontal sundial calculations~~ ~~20 PRINT "Enter latitude => ";~~ ~~30 INPUT LAT~~ ~~40 PRINT "Enter longitude => ";~~ ~~50 INPUT LNG~~ ~~60 PRINT "Enter legal meridian => ";~~ ~~70 INPUT REF~~ ~~80 PRINT~~ ~~90 LET PI = 4 * ATN(1)~~ ~~100 LET SLAT = SIN(LAT * PI / 180)~~ ~~110 PRINT " sine of latitude: "; USING "#.##^^^^"; SLAT~~ ~~120 PRINT " diff longitude: "; USING "####.###"; LNG - REF~~ ~~130 PRINT~~ ~~140 PRINT "Hour, sun hour angle, dial hour line angle from 6am to 6pm"~~ ~~150 FOR H% = -6 TO 6~~ ~~160 LET HRA = 15 * H%~~ ~~170 LET HRA = HRA - (LNG - REF): ' correct for longitude difference~~ ~~180 LET HLA = ATN(SLAT * TAN(HRA * PI / 180)) * 180 / PI~~ ~~190 PRINT "HR="; USING "+##"; H%;~~ ~~200 PRINT "; HRA="; USING "+###.###"; HRA;~~ ~~210 PRINT "; HLA="; USING "+###.###"; HLA~~ ~~220 NEXT H%~~ ~~230 END~~ ~~</syntaxhighlight>~~ ~~{{out}}~~ ~~<pre>~~ ~~Enter latitude => ? -4.95~~ ~~Enter longitude => ? -150.5~~ ~~Enter legal meridian => ? -150~~ ~~sine of latitude: -.86E-01~~ ~~diff longitude: -0.500~~ ~~Hour, sun hour angle, dial hour line angle from 6am to 6pm~~ ~~HR= -6; HRA= -89.500; HLA= +84.225~~ ~~HR= -5; HRA= -74.500; HLA= +17.283~~ ~~HR= -4; HRA= -59.500; HLA= +8.334~~ ~~HR= -3; HRA= -44.500; HLA= +4.847~~ ~~HR= -2; HRA= -29.500; HLA= +2.795~~ ~~HR= -1; HRA= -14.500; HLA= +1.278~~ ~~HR= +0; HRA= +0.500; HLA= -0.043~~ ~~HR= +1; HRA= +15.500; HLA= -1.371~~ ~~HR= +2; HRA= +30.500; HLA= -2.910~~ ~~HR= +3; HRA= +45.500; HLA= -5.018~~ ~~HR= +4; HRA= +60.500; HLA= -8.671~~ ~~HR= +5; HRA= +75.500; HLA= -18.451~~ ~~HR= +6; HRA= +90.500; HLA= +84.225~~ </pre> Line 1,450 ⟶ 1,874: HR= 5pm HRA=75.5, HLA=-18.45099922256532 HR= 6pm HRA=90.5, HLA=-95.77516739863968</pre> ~~=={{header\|IS-BASIC}}==~~ ~~<syntaxhighlight lang="is-basic">100 PROGRAM "SunDial.bas"~~ ~~110 OPTION ANGLE RADIANS~~ ~~120 TEXT 80~~ ~~130 LET DR=PI/180:LET RD=180/PI~~ ~~140 INPUT PROMPT "Enter latitude: ":LAT~~ ~~150 INPUT PROMPT "Enter longitude: ":LNG~~ ~~160 INPUT PROMPT "Enter legal meridian: ":REF~~ ~~170 LET S=SIN(LATDR)~~ ~~180 PRINT :PRINT "Sine of latitude :";S~~ ~~190 PRINT "Diff longitude: ";LNG-REF~~ ~~200 PRINT " Hour,",,"sun hour angle, dial hour line angle from 6am to 6 pm"~~ ~~210 FOR H=6 TO 18~~ ~~220 LET HRA=15H-LNG+REF~~ ~~230 LET HLA=ATN(STAN(HRADR))RD~~ ~~240 PRINT "HR = ";H,"HRA =";HRA,"HLA =";HLA~~ ~~250 NEXT</syntaxhighlight>~~ =={{header\|J}}== Line 1,807 ⟶ 2,213: 5 PM +75.500 -18.451 6 PM +90.500 -95.775 ~~</pre>~~ ~~=={{header\|Liberty BASIC}}==~~ ~~Based on Algol & BBC BASIC versions. Note Liberty BASIC works in radians.~~ ~~<syntaxhighlight lang="lb">global pi~~ ~~pi = 3.14159265~~ ~~input "Enter latitude (degrees) : "; latitude ' -4.95~~ ~~input "Enter longitude (degrees) : "; longitude ' -150.5~~ ~~input "Enter legal meridian (degrees): "; meridian ' -150.0~~ ~~print~~ ~~print "Time Sun hour angle Dial hour line angle"~~ ~~for hour = 6 TO 18~~ ~~hra = 15 hour - longitude + meridian - 180~~ ~~hla = rad2deg(atn(sin(deg2rad(latitude)) * tan(deg2rad(hra))))~~ ~~if abs(hra) > 90 then hla = hla + 180 * sgn(hra * latitude)~~ ~~print using( "##.##", hour), using("####.### ", hra), using("####.###", hla)~~ ~~next hour~~ ~~end~~ ~~function rad2deg(theta)~~ ~~rad2deg = theta * 180 / pi~~ ~~end function~~ ~~function deg2rad(theta)~~ ~~deg2rad = theta * pi / 180~~ ~~end function~~ ~~function sgn(x)~~ ~~if x > 0 then sgn = 1 else sgn = -1~~ ~~end function</syntaxhighlight>~~ ~~{{out}}~~ ~~<pre>~~ ~~Enter latitude (degrees) : -4.95~~ ~~Enter longitude (degrees) : -150.5~~ ~~Enter legal meridian (degrees): -150.0~~ ~~Time Sun hour angle Dial hour line angle~~ ~~6.00 -89.500 84.225~~ ~~7.00 -74.500 17.283~~ ~~8.00 -59.500 8.334~~ ~~9.00 -44.500 4.847~~ ~~10.00 -29.500 2.795~~ ~~11.00 -14.500 1.278~~ ~~12.00 0.500 -0.043~~ ~~13.00 15.500 -1.371~~ ~~14.00 30.500 -2.910~~ ~~15.00 45.500 -5.018~~ ~~16.00 60.500 -8.671~~ ~~17.00 75.500 -18.451~~ ~~18.00 90.500 -95.775~~ </pre> Line 1,962 ⟶ 2,318: 5 75.5 -18.451 6 90.5 84.2248</pre> ~~=={{header\|Microsoft Small Basic}}==~~ ~~{{trans\|ALGOL-68}}~~ ~~<syntaxhighlight lang="microsoftsmallbasic">~~ ~~TextWindow.Write("Enter latitude => ")~~ ~~lat = TextWindow.ReadNumber()~~ ~~TextWindow.Write("Enter longitude => ")~~ ~~lng = TextWindow.ReadNumber()~~ ~~TextWindow.Write("Enter legal meridian => ")~~ ~~ref = TextWindow.ReadNumber()~~ ~~sLat = Math.Sin(Math.GetRadians(lat))~~ ~~TextWindow.WriteLine("")~~ ~~TextWindow.Write(" sine of latitude: ")~~ ~~TextWindow.WriteLine(Math.Round(sLat * 10000) / 10000)~~ ~~TextWindow.Write(" diff longitude: ")~~ ~~TextWindow.WriteLine(lng - ref)~~ ~~TextWindow.WriteLine("")~~ ~~TextWindow.WriteLine("Hour, sun hour angle, dial hour line angle from 6am to 6pm")~~ ~~For hour = -6 To 6~~ ~~hourAngle = 15 * hour~~ ~~hourAngle = hourAngle - (lng - ref) ' correct for longitude difference~~ ~~hourLineAngle = math.GetDegrees(Math.ArcTan(sLat * Math.Tan(Math.GetRadians(hourAngle))))~~ ~~TextWindow.Write("HR=")~~ ~~TextWindow.CursorLeft = 3 + (3 - Text.GetLength(hour))~~ ~~TextWindow.Write(hour)~~ ~~TextWindow.Write("; HRA=")~~ ~~TextWindow.CursorLeft = 12 + (6 - Text.GetLength(hourAngle))~~ ~~TextWindow.Write(hourAngle)~~ ~~TextWindow.Write("; HLA=")~~ ~~TextWindow.CursorLeft = 24 + (4 - Text.GetLength(Math.Floor(hourLineAngle)))~~ ~~TextWindow.Write(Math.Round(hourLineAngle * 1000) / 1000)~~ ~~TextWindow.WriteLine("")~~ ~~EndFor~~ ~~</syntaxhighlight>~~ ~~{{out}}~~ ~~Enter latitude => -4.95~~ ~~Enter longitude => -150.5~~ ~~Enter legal meridian => -150~~ ~~sine of latitude: -0.0863~~ ~~diff longitude: -0.5~~ ~~Hour, sun hour angle, dial hour line angle from 6am to 6pm~~ ~~HR= -6; HRA= -89.5; HLA= 84.225~~ ~~HR= -5; HRA= -74.5; HLA= 17.283~~ ~~HR= -4; HRA= -59.5; HLA= 8.334~~ ~~HR= -3; HRA= -44.5; HLA= 4.847~~ ~~HR= -2; HRA= -29.5; HLA= 2.795~~ ~~HR= -1; HRA= -14.5; HLA= 1.278~~ ~~HR= 0; HRA= 0.5; HLA= -0.043~~ ~~HR= 1; HRA= 15.5; HLA= -1.371~~ ~~HR= 2; HRA= 30.5; HLA= -2.91~~ ~~HR= 3; HRA= 45.5; HLA= -5.018~~ ~~HR= 4; HRA= 60.5; HLA= -8.671~~ ~~HR= 5; HRA= 75.5; HLA= -18.451~~ ~~HR= 6; HRA= 90.5; HLA= 84.225~~ =={{header\|МК-61/52}}== Line 2,615 ⟶ 2,914: 6:00 PM 90.5 84.225 ~~</pre>~~ ~~=={{header\|PureBasic}}==~~ ~~{{trans\|ALGOL 68}}~~ ~~<syntaxhighlight lang="purebasic">If OpenConsole()~~ ~~Define.f lat, slat, lng, ref~~ ~~Define.i h~~ ~~Print("Enter latitude => "): lat=ValF(Input())~~ ~~Print("Enter longitude => "): lng=ValF(Input())~~ ~~Print("Enter legal meridian => "): ref=ValF(Input())~~ ~~PrintN("")~~ ~~slat=Sin(lat2#PI/360)~~ ~~PrintN(" sine of latitude: "+StrF(slat,3))~~ ~~PrintN(" diff longitude: "+StrF((lng-ref),3)+#CRLF$)~~ ~~PrintN("Hour, sun hour angle, dial hour line angle from 6am to 6pm")~~ ~~For h=-6 To 6~~ ~~Define.f hra, hla~~ ~~hra=15h~~ ~~hra=hra-(lng-ref)~~ ~~hla=ATan(slatTan(hra2#PI/360))360/(2#PI)~~ ~~PrintN("HR="+RSet(Str(h),3)+"; HRA="+RSet(StrF(hra,3),7)+"; HLA="+RSet(StrF(hla,3),7))~~ ~~Next~~ ~~EndIf</syntaxhighlight>~~ ~~{{out}}~~ ~~<pre>Enter latitude => -4.95~~ ~~Enter longitude => -150.5~~ ~~Enter legal meridian => -150~~ ~~sine of latitude: -0.086~~ ~~diff longitude: -0.500~~ ~~Hour, sun hour angle, dial hour line angle from 6am to 6pm~~ ~~HR= -6; HRA=-89.500; HLA= 84.225~~ ~~HR= -5; HRA=-74.500; HLA= 17.283~~ ~~HR= -4; HRA=-59.500; HLA= 8.334~~ ~~HR= -3; HRA=-44.500; HLA= 4.847~~ ~~HR= -2; HRA=-29.500; HLA= 2.795~~ ~~HR= -1; HRA=-14.500; HLA= 1.278~~ ~~HR= 0; HRA= 0.500; HLA= -0.043~~ ~~HR= 1; HRA= 15.500; HLA= -1.371~~ ~~HR= 2; HRA= 30.500; HLA= -2.910~~ ~~HR= 3; HRA= 45.500; HLA= -5.018~~ ~~HR= 4; HRA= 60.500; HLA= -8.671~~ ~~HR= 5; HRA= 75.500; HLA=-18.451~~ ~~HR= 6; HRA= 90.500; HLA= 84.225~~ </pre> Line 3,120 ⟶ 3,371: HR = 6; HRA = 90.500; HLA = 84.225 </pre> ~~=={{header\|Run BASIC}}==~~ ~~<syntaxhighlight lang="runbasic">global pi~~ ~~pi = 22 / 7~~ ~~print "Enter latitude (degrees) : "; :input latitude ' -4.95~~ ~~print "Enter longitude (degrees) : "; :input longitude ' -150.5~~ ~~print "Enter legal meridian (degrees): "; :input meridian ' -150.0~~ ~~print~~ ~~print "Time Sun hour angle Dial hour line angle"~~ ~~for hour = 6 TO 18~~ ~~hra = (15 * hour) - longitude + meridian -180~~ ~~hla =rad2deg( atn( sin( deg2rad( latitude)) tan( deg2rad( hra))))~~ ~~if abs( hra) >90 then hla =hla +180 sgn( hra latitude)~~ ~~print using( "##", hour);" ";using("####.##", hra);" ";using("####.###", hla)~~ ~~next hour~~ ~~function rad2deg( theta)~~ ~~rad2deg =theta 180 /pi~~ ~~end function~~ ~~function deg2rad( theta)~~ ~~deg2rad =theta pi /180~~ ~~end function~~ ~~function sgn( x)~~ ~~if x >0 then sgn =1 else sgn =-1~~ ~~end function~~ ~~end</syntaxhighlight>~~ ~~{{out}}~~ ~~<pre>~~ ~~Enter latitude: -4.95~~ ~~Enter longitude: -150.5~~ ~~Enter legal meridian: -150~~ ~~Time Sun hour angle Dial hour line angle~~ ~~6 -89.50 84.606~~ ~~7 -74.50 17.316~~ ~~8 -59.50 8.342~~ ~~9 -44.50 4.850~~ ~~10 -29.50 2.796~~ ~~11 -14.50 1.279~~ ~~12 0.50 -0.043~~ ~~13 15.50 -1.371~~ ~~14 30.50 -2.911~~ ~~15 45.50 -5.021~~ ~~16 60.50 -8.680~~ ~~17 75.50 -18.488~~ ~~18 90.50 -96.224</pre>~~ =={{header\|Rust}}== Line 3,493 ⟶ 3,693: {{trans\|Go}} {{libheader\|Wren-fmt}} <syntaxhighlight lang="~~ecmascript~~wren">import "io" for Stdin, Stdout import "./fmt" for Fmt var getNum = Fn.new { \|prompt\| Line 3,689 ⟶ 3,889: o_ft: db "HR= %3d; ",9," HRA=%7.3f; ",9," HLA= %7.3f", 10, 0</syntaxhighlight> ~~=={{header\|XBasic}}==~~ ~~{{trans\|ALGOL-68}}~~ ~~{{works with\|Windows XBasic}}~~ ~~<syntaxhighlight lang="xbasic">~~ ~~PROGRAM "sundial"~~ ~~VERSION "0.0001"~~ ~~IMPORT "xma"~~ ~~DECLARE FUNCTION Entry()~~ ~~FUNCTION Entry()~~ ~~lat! = SINGLE(INLINE$("Enter latitude => "))~~ ~~lng! = SINGLE(INLINE$("Enter longitude => "))~~ ~~ref! = SINGLE(INLINE$("Enter legal meridian => "))~~ ~~PRINT~~ ~~slat! = SIN(lat! $$PI / 180.0)~~ ~~PRINT " sine of latitude: "; FORMAT$("#.##^^^^", slat!)~~ ~~PRINT " diff longitude: "; FORMAT$("#.###", lng! - ref!)~~ ~~PRINT~~ ~~PRINT "Hour, sun hour angle, dial hour line angle from 6am to 6pm"~~ ~~FOR hour@ = -6 TO 6~~ ~~hourAngle! = 15 * hour@~~ ~~hourAngle! = hourAngle! - (lng! - ref!) ' correct for longitude difference~~ ~~hourLineAngle! = ATAN(slat! * TAN(hourAngle! * $$PI / 180.0)) * 180.0 / $$PI~~ ~~PRINT "HR="; FORMAT$("###", hour@);~~ ~~PRINT "; HRA="; FORMAT$("####.###", hourAngle!);~~ ~~PRINT "; HLA="; FORMAT$("####.###", hourLineAngle!)~~ ~~NEXT hour@~~ ~~END FUNCTION~~ ~~END PROGRAM~~ ~~</syntaxhighlight>~~ ~~{{out}}~~ ~~<pre>~~ ~~Enter latitude => -4.95~~ ~~Enter longitude => -150.5~~ ~~Enter legal meridian => -150~~ ~~sine of latitude: -8.63E-02~~ ~~diff longitude: -0.500~~ ~~Hour, sun hour angle, dial hour line angle from 6am to 6pm~~ ~~HR= -6; HRA= -89.500; HLA= 84.225~~ ~~HR= -5; HRA= -74.500; HLA= 17.283~~ ~~HR= -4; HRA= -59.500; HLA= 8.334~~ ~~HR= -3; HRA= -44.500; HLA= 4.847~~ ~~HR= -2; HRA= -29.500; HLA= 2.795~~ ~~HR= -1; HRA= -14.500; HLA= 1.278~~ ~~HR= 0; HRA= 0.500; HLA= -0.043~~ ~~HR= 1; HRA= 15.500; HLA= -1.371~~ ~~HR= 2; HRA= 30.500; HLA= -2.910~~ ~~HR= 3; HRA= 45.500; HLA= -5.018~~ ~~HR= 4; HRA= 60.500; HLA= -8.671~~ ~~HR= 5; HRA= 75.500; HLA= -18.451~~ ~~HR= 6; HRA= 90.500; HLA= 84.225~~ ~~</pre>~~ =={{header\|XPL0}}== Line 3,851 ⟶ 3,994: <as above> </pre> ~~=={{header\|ZX Spectrum Basic}}==~~ ~~{{trans\|ERRE}}~~ ~~<syntaxhighlight lang="zxbasic">10 DEF FN r(x)=xPI/180~~ ~~20 DEF FN d(x)=x180/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=15h-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+180SGN (hralatitude)~~ ~~160 PRINT h;" ";hra;" ";hla~~ ~~170 NEXT h</syntaxhighlight>~~