pub fn slab_escape(a: f64, b: f64) -> (f64, f64, f64)
g(a, b), ∂g/∂a and ∂g/∂b for optical depths a, b of 0 or more.
g(a, b)
∂g/∂a
∂g/∂b
a
b