pub fn binomial_stable(n: u64, k: u64) -> CoreResult<f64>
Stable computation of binomial coefficients