i_mmc = NewInput.MMC(lambda=(120), mu=(190), c=1, n=0, method=0)

o_mmc <- QueueingModel(i_mmc)

L = L(o_mmc)

p = RO(o_mmc)

Wq = Wq(o_mmc)

Lq = Lq(o_mmc)

W = W(o_mmc)

rbind(c('L','Lq','W','Wq','p','100/W'),c(L,Lq,W,Wq,p,100/W))
##      [,1]               [,2]               [,3]                
## [1,] "L"                "Lq"               "W"                 
## [2,] "1.71428571428571" "1.08270676691729" "0.0142857142857143"
##      [,4]                  [,5]                [,6]   
## [1,] "Wq"                  "p"                 "100/W"
## [2,] "0.00902255639097744" "0.631578947368421" "7000"

Response 1 - Avg time in sys : Histogram / Boxplot Response 2 - Avg time waiting in Queue : Boxplot Rotated