pub fn rel_equal(one:f64, other: f64, diff: f64)->bool{
    if other != 0.0{
        (1.0-one/other).abs()<diff
    } else { one.abs()< diff }
}