Function bearing 

pub fn bearing(lat1: f64, lon1: f64, lat2: f64, lon2: f64) -> f64

Initial compass bearing (degrees, clockwise from north) from point (lat1, lon1) to (lat2, lon2), both in decimal degrees.