public static void Main()
var φ1 = lat1 * Math.PI/180;
var φ2 = lat2 * Math.PI/180;
var Δφ = (lat2-lat1) * Math.PI/180;
var Δλ = (lon2-lon1) * Math.PI/180;
var a = Math.Sin(Δφ/2) * Math.Sin(Δφ/2) +
Math.Cos(φ1) * Math.Cos(φ2) *
Math.Sin(Δλ/2) * Math.Sin(Δλ/2);
var c = 2 * Math.Atan2(Math.Sqrt(a), Math.Sqrt(1-a));
Console.WriteLine(d / 1000);