pub fn cheb_sf(x: f64) -> f64
k < 0 => minimum mass above k k > 0 => maximum mass above k k = 0 => exact mass above k