estimate.rs

let alpha: f64 = match b {
    4 => 0.673,
    5 => 0.697,
    6 => 0.709,
    // b >= 7
    _ => 0.7213/(1.0 + 1.079/(m as f64))
};

let mut indicator: f64 = 0.0;
let base: f64 = 2.0;
for first_non_zero in &first_non_zero_by_experiment {
    indicator += base.powf(-(*first_non_zero as f64));
}
let m_multiplier = m as f64;
let mut estimate: f64 = (m_multiplier * m_multiplier * alpha) / indicator;