pub fn robinson_g(eps: f64) -> f64
Robinson’s fit g(ε) = 3.4008 ε^(1/6) + 0.40244 ε^(3/4) + ε.
g(ε) = 3.4008 ε^(1/6) + 0.40244 ε^(3/4) + ε