pub fn inhomogeneous_poisson( rate_func: fn(f64) -> f64, duration_s: f64, dt: f64, seed: u64, ) -> Vec<f64>
Inhomogeneous Poisson via thinning (Lewis & Shedler 1979).