n <- 1e2
arrivals <- cumsum(rexp(n, 1.8))
service <- rexp(n)
departures <- queue(
arrivals, service, servers = 2)
head(departures)
curve(ecdf(departures)(x) * n,
from = 0, to = max(departures),
xlab = "Time", ylab = "Number of customers")
curve(ecdf(arrivals)(x) * n,
from = 0, to = max(departures),
col = "red", add = TRUE)
Run the code above in your browser using DataLab