XProfan

Titre de la fenêtre "   KÜRZESTER NAUTISCHER STARTKURS (Erdkugel R = 6371.221 km)"
'Sphärische Geometrie - angeregt de ´User0815´ (CL) CopyLeft 2018-03 P.Specht, sans Gewähr
Fenêtre Style 24:font 2:Déclarer lat0!,lon0!,lat1!,lon1!
var f!=pi()/180:var R! = 6371.221'km
:Proc ArcCos :Paramètres w!:declare ac!:Si w!=1:ac!=0:elseif w!=-1:ac!=pi()
d'autre: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))
imprimer "\n\n     Koordinaten des Ausgangspunktes"
imprimer "\n Geogr.Breite Phi° Äquator vers Norden [+/-0...90, Lat.°]:",:input lat0!
imprimer "\n Geogr.Longueur (Lambda° Greenwich gg.West [+/-0..180, Lon.°]:",:input lon0!
imprimer "\n\n     Koordinaten des Zielpunktes"
imprimer "\n                                  Breitengrad [Latitude °]:",:input lat1!
imprimer "\n                                  Längengrad [Longitude °]:",:input lon1!
:si (lat0!=0) and (lon0!=0) and (lat1!=0) and (lon1!=0):font 0
:Imprimer "\n\n SELBSTTEST: Leningrad -> SanFrancisco: s'il te plaît vérifier: 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!
Imprimer "\n\n                     Kurs chez Aufbruch: ",\
Startkurs(lat0!,lon0!,lat1!,lon1!)/f!,"  °"
imprimer   "\n                         Gesamtdistanz: ",\
Dist(lat0!,lon0!,lat1!,lon1!),"km"
waitinput:Goto "Start"

Titre de la fenêtre "Großkreis-Distanzformel testen"
CLS:font 2
imprimer "Großkreis-Distanz:",Format $("####.#############",\
GroßkreisDistanz(36.12,-86.67  ,  33.94,-118.4)),"km."
Waitinput' 2887.2599506071104 km
Fin

Proc ArcSin :Paramètres Angle!:Déclarer Erg!

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

        D'autre :Erg!=ArcTan(Angle!/Sqrt(1-Sqr(Angle!)))
        EndIf:Retour Erg!

    ENDPROC

    Proc GroßkreisDistanz :Paramètres 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