Function c_eta 

pub fn c_eta(alpha: f64, beta: f64, smoothness_l: f64, eta: f64) -> f64

Geometric contraction coefficient c(eta) = (alpha-beta) - L eta (alpha+beta)^2 / 2.