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