pub fn bloch_term(y: f64) -> f64
Bloch correction L2(y), y = z α / β (see module docs).
L2(y)
y = z α / β