Python script that performs haversine calculations

haversine_calculation/__init__.py

@@ -0,0 +1 @@
+# This file makes haversine_calculation a Python package

haversine_calculation/haversine.py

@@ -0,0 +1,31 @@
+import math
+
+
+def haversine(lat1, lon1, lat2, lon2):
+    """
+    Calculate the great circle distance in kilometers between two points
+    on the earth (specified in decimal degrees)
+    """
+    # Convert decimal degrees to radians
+    lat1, lon1, lat2, lon2 = map(math.radians, [lat1, lon1, lat2, lon2])
+
+    # Haversine formula
+    dlat = lat2 - lat1
+    dlon = lon2 - lon1
+    a = (
+        math.sin(dlat / 2) ** 2
+        + math.cos(lat1) * math.cos(lat2) * math.sin(dlon / 2) ** 2
+    )
+    c = 2 * math.asin(math.sqrt(a))
+    r = 6371  # Radius of earth in kilometers. Use 3956 for miles.
+    return c * r
+
+
+# Example usage
+if __name__ == "__main__":
+    # Coordinates of New York City and London
+    lat1, lon1 = 40.7128, -74.0060  # New York City
+    lat2, lon2 = 51.5074, -0.1278  # London
+
+    distance = haversine(lat1, lon1, lat2, lon2)
+    print(f"Distance between New York City and London: {distance:.2f} km")