XProfan

Título de la ventana "   KÜRZESTER NAUTISCHER STARTKURS (Erdkugel R = 6371.221 km)"
'Sphärische Geometrie - angeregt de ´User0815´ (CL) CopyLeft 2018-03 P.Pájaro carpintero, sin Gewähr
Ventana de Estilo 24:font 2:Declarar lat0!,lon0!,lat1!,lon1!
var f!=pi()/180:var R! = 6371.221'km
:Proc ArcCos :Parámetros w!:declarar ac!:If w!=1:ac!=0:elseif w!=-1:ac!=pi()
más:ac!=Pi()/2-arctan(w!/Sqrt(1-Sqr(w!))):endif:volver ac!:ENDPROC

Proc Startkurs :parámetros lat0!,lon0!,lat1!,lon1!

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

ENDPROC

proc Dist :parámetros lat0!,lon0!,lat1!,lon1!

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

ENDPROC

Start:
Cls rgb(200+rnd(56),200+rnd(56),200+rnd(56))
imprimir "\n\n     Koordinaten des Ausgangspunktes"
imprimir "\n Geogr.Breite Phi° Äquator gegen Norden [+/-0...90, Lat.°]:",:input lat0!
imprimir "\n Geogr.Longitud (Lambda° Greenwich gg.West [+/-0..180, Lon.°]:",:input lon0!
imprimir "\n\n     Koordinaten des Zielpunktes"
imprimir "\n                                  Breitengrad [Latitude °]:",:input lat1!
imprimir "\n                                  Längengrad [Longitude °]:",:input lon1!
:if (lat0!=0) and (lon0!=0) and (lat1!=0) and (lon1!=0):font 0
:Imprimir "\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!
Imprimir "\n\n                     Kurs en Aufbruch: ",\
Startkurs(lat0!,lon0!,lat1!,lon1!)/f!,"  °"
imprimir   "\n                         Gesamtdistanz: ",\
Dist(lat0!,lon0!,lat1!,lon1!),"km"
waitinput:Goto "Start"

Título de la ventana "Großkreis-Distanzformel testen"
CLS:font 2
imprimir "Großkreis-Distanz:",Formato$("####.#############",\
GroßkreisDistanz(36.12,-86.67  ,  33.94,-118.4)),"km."
Waitinput' 2887.2599506071104 km
End

Proc ArcSin :Parámetros Angle!:Declarar Erg!

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

        Más :Erg!=ArcTan(Angle!/Sqrt(1-Sqr(Angle!)))
        EndIf:Volver Erg!

    ENDPROC

    Proc GroßkreisDistanz :Parámetros 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!
        volver hDist!

    ENDPROC