pub fn bakhshali_sqrt(x: f32) -> f32