service_time <- 2 num_customer <- 10 service_min <- 22 lamda <- 1/service_time prob <- 1 - pgamma(service_min, shape = num_customer, rate = lamda) prob
## [1] 0.3405106