pub fn ang_diff(x: f64, y: f64, e: &mut f64) -> f64
Compute the difference between two angles in degrees