Kurswinkel Berechung

[robinWie]

Hallo, ich bin gerade dabei mir eine Unit zusammenzubasteln wo einige Funktionen zu Koordinaten drin sind. Momentan versuche ich eine Funktion zu schreiben, die einem den Kurswinkel (N=0° O=90° S=180° W=270°) von einer Koordinate zu einer anderen Ausgibt.
Im Internet hab hier ich auch eine Rechung dazu gefunden (weiter unten unter Beispielsrechnung) Aber beim Umsetzten in Delphi hab ich Probleme. Da kommt irgendwie nie das raus, was eigentlich rauskommt.
Es soll der Ankunftswinkel ausgerechnet werden.

Das hier ist mein Bisheriger Ansatz

Delphi-Quellcode:

function Direction(Lat1, Lat2, Lon1, Lon2: Extended): Single;
var
 Lt1: single;
 LT2: single;
 LG1: single;
 LG2: single;
 DLA: Extended;
 c: Extended;
 w: Extended;
begin
  Lt1 := DegToRad(Lat1);
  Lt2 := DegToRad(Lat2);
  Lg1 := DegToRad(Lon1);
  Lg2 := DegToRad(Lon2);
  DLA := DegToRad(Lon1 - Lon2);

  C := Sin(Lt1) * Sin(Lt2) + Cos(Lt1) * Cos(Lt2) * Cos(DLA);
  W := (Sin(Lt1) - Sin(Lt2) * cos(C)) / (cos(Lt2) * sin(C));
  Result := RadToDeg(W);
end;

Ich hab das mit den Radiant (mit den Delphi arbeitet) und den Grad soweit berücksichtigt, indem ich am Anfang konvertiere und am Ende. Aber wie gesagt kommt immer das falsche Raus

Danke im Voraus
Robin