#### Ejercicio: Centro de llamadas con 5 operadores ####
library(queueing)

# Parámetros
lambda <- 29   # Tasa de llegada (llamadas por hora)
mu <- 6        # Tasa de servicio por operador (llamadas por hora)
s <- 5         # Número de servidores (operadores)

# Definir el modelo M/M/s
modelo_llamadas <- NewInput.MMC(lambda = lambda, mu = mu, c = s)

# Resolver el modelo
resultado <- QueueingModel(modelo_llamadas)

# Mostrar resultados principales
cat("Factor de utilización (rho):", resultado$RO, "\n")
## Factor de utilización (rho): 0.9666667
cat("Número promedio de clientes en la fila (Lq):", resultado$Lq, "\n")
## Número promedio de clientes en la fila (Lq): 26.61578
cat("Número promedio de clientes en el sistema (L):", resultado$L, "\n")
## Número promedio de clientes en el sistema (L): 31.44912
cat("Tiempo promedio en la fila (Wq):", resultado$Wq, "horas\n")
## Tiempo promedio en la fila (Wq): 0.9177856 horas
cat("Tiempo promedio en el sistema (W):", resultado$W, "horas\n")
## Tiempo promedio en el sistema (W): 1.084452 horas