rm(list=ls())
cat("\014") 

# ****************************************************************************
#                  FUNCIÓN PARA TABLA PECK & HAZELWOOD
# ****************************************************************************
PeckHazelwood<-function(N,X,M){
# Calcula el factor de eficiencia (F) y la probabilidad de esperar (D)
# Argumentos:
      # N: número total de clientes del sistema (N' en el texto)
      # X: factor de servicio X=Ts/(Ts+U)
      # M: número de canales de atención

  library(queueing)  
  
  U<-(1-X)/X
  lambda<-1/U
  mu<-1
  Ts<-1/mu
  
  Peck_Hazelwood <- NewInput.MMCKM(lambda, mu, c=M, k=N, m=N)
  model <- QueueingModel(Peck_Hazelwood)
  
  # F: Factor de eficiencia
  F<-(Ts+U)/(Ts+U+Wq(model))
  
  # Dnum: número de clientes que piden el servicio en el intervalo Dt
  # y tienen que esperar
  Dnum<-0                       
  for(n in M:(N-1)){(Dnum<-Dnum+model$Pn[n+1]*(N-n))}
  # Probabilidad de esperar= Dnum/J
  D<-Dnum/(N-L(model))

  
  N<-as.integer(N)
  M<-as.integer(M)
  output_p<-paste("N=",N,"  X=",X,"  M=",M,"  D=",round(D,3),"  F=",round(F,3))
  # descomentar la línea siguiente para impresión en consola de inputs/outputs
  print(output_p, quote=FALSE)  
  output_l<-list(N,M,X,round(D,3),round(F,3))
  # descomentar la línea siguiente para devolver lista de inputs/ouputs
  #return(output_l)             
}

# ****************************************************************************

# Ejemplos de uso de la función
#
# 1) Teoria de Colas - M.Miranda (pag135-136)
N<-10
X<-0.2
M<-2
PeckHazelwood(N,X,M)
#
#
# 2) Teoria de Colas - M.Miranda Ejemplo 4.3
N<-10
X<-0.110
M<-3
PeckHazelwood(N,X,M)
#
#
# 3)
N<-4
X<-0.2
M<-1
PeckHazelwood(N,X,M)