pub fn haversine_distance(lat1: f64, lon1: f64, lat2: f64, lon2: f64) -> f64
Calculate the distance between two points using the Haversine formula