amackey

nmin <- 120;
time <- nmin * 60;
lambda <- 10/60; # 10 passengers every minute, on average
arrivals <- rpois(time, lambda);
population <- sum(arrivals);

nelev <- 4; # number of elevators
nfloor <- 7; # number of floors
elevators <- matrix(1, ncol=nelev, nrow=time)