pub fn two_one_diff(a1: f64, a0: f64, b: f64) -> [f64; 3]