pub fn poisson_process(tmax: f64, lambda: f64) -> Array1<f64>
Simulate a Poisson process with constant intensity $\lambda$.