fn phi1<R: Scalar>(k: R, x: f64) -> R
φ(k, x) = (1 − e^{−k·x}) / k, by its series where |k·x| < 0.5.
φ(k, x) = (1 − e^{−k·x}) / k
|k·x| < 0.5