XProfan

WindowTitle "   KÜRZESTER NAUTISCHER STARTKURS (Erdkugel R = 6371.221 km)"
'Sphärische Geometrie - angeregt von ´User0815´ (CL) CopyLeft 2018-03 P.Specht, ohne Gewähr
WindowStyle 24:font 2:Declare lat0!,lon0!,lat1!,lon1!
var f!=pi()/180:var R! = 6371.221'km
:Proc ArcCos :Parameters w!:declare ac!:If w!=1:ac!=0:elseif w!=-1:ac!=pi()
else:ac!=Pi()/2-arctan(w!/Sqrt(1-Sqr(w!))):endif:return ac!:EndProc

Proc Startkurs :parameters lat0!,lon0!,lat1!,lon1!

    var cosg!=sin(lat0!)*sin(lat1!)+cos(lat0!)*cos(lat1!)*cos(lon1!-lon0!)
    return arccos(cosg!)

endproc

proc Dist :parameters lat0!,lon0!,lat1!,lon1!

    return R!*StartKurs(lat0!,lon0!,lat1!,lon1!)

endproc

Start:
Cls rgb(200+rnd(56),200+rnd(56),200+rnd(56))
print "\n\n     Koordinaten des Ausgangspunktes"
print "\n Geogr.Breite Phi° Äquator gegen Norden [+/-0...90, Lat.°]:",:input lat0!
print "\n Geogr.Länge (Lambda° Greenwich gg.West [+/-0..180, Lon.°]:",:input lon0!
print "\n\n     Koordinaten des Zielpunktes"
print "\n                                  Breitengrad [Latitude °]:",:input lat1!
print "\n                                  Längengrad [Longitude °]:",:input lon1!
:if (lat0!=0) and (lon0!=0) and (lat1!=0) and (lon1!=0):font 0
:Print "\n\n SELBSTTEST: Leningrad -> SanFrancisco: Bitte prüfen: 79.74°, 8867 km ???"
:lat0!=59.9:lon0!=-30.3:lat1!=37.8:lon1!=122.4:font 2:endif
lat0!=lat0!*f!:lon0!=lon0!*f!:lat1!=lat1!*f!:lon1!=lon1!*f!
Print "\n\n                     Kurs bei Aufbruch: ",\
Startkurs(lat0!,lon0!,lat1!,lon1!)/f!,"  °"
print   "\n                         Gesamtdistanz: ",\
Dist(lat0!,lon0!,lat1!,lon1!),"km"
waitinput:Goto "Start"

WindowTitle "Großkreis-Distanzformel testen"
CLS:font 2
print "Großkreis-Distanz:",Format$("####.#############",\
GroßkreisDistanz(36.12,-86.67  ,  33.94,-118.4)),"km."
Waitinput' 2887.2599506071104 km
End

Proc ArcSin :Parameters Angle!:Declare Erg!

    If Abs(Angle!)=1:Erg!=Angle!*Pi()/2

        Else :Erg!=ArcTan(Angle!/Sqrt(1-Sqr(Angle!)))
        EndIf:Return Erg!

    EndProc

    Proc GroßkreisDistanz :Parameters n1!,e1! , n2!,e2!'Dezimalgrad

        var D2R!=ArcTan(1)/45:var diam!=2*6372.8:var Lg1m2!=(e1!-e2!)*D2R!
        var Lt1!=n1!*D2R!:var Lt2!=n2!*D2R!:var dz!=sin(Lt1!)-sin(Lt2!)
        var dx!=cos(Lg1m2!)*cos(Lt1!)-cos(Lt2!):var dy!=sin(Lg1m2!)*cos(Lt1!)
        var hDist!=ArcSin(Sqrt(sqr(dx!)+sqr(dy!)+sqr(dz!))/2)*diam!
        return hDist!

    Endproc