rgsl::function

type function<T> = fn(x: f64, p: &mut T) -> f64;