pub fn babylonian_sqrt(x: f32) -> f32