pub fn variable_poisson<F>( tmax: f64, lambda: &F, max_lambda: f64, ) -> TimeProcessResultwhere F: Fn(f64) -> f64 + Send + Sync,
Simulate a Poisson process with variable intensity.