fn exp_integral<R: Scalar>(k: R, v1: f64, v2: f64) -> R
∫_{v1}^{v2} e^{−k·v} dv = e^{−k·v1}·φ(k, v2 − v1).
∫_{v1}^{v2} e^{−k·v} dv = e^{−k·v1}·φ(k, v2 − v1)