set term postscript landscape enhanced  "Times-Roman" 14
set out 'MMK_compare.ps'
set data style steps
set size 1.0,1.0
#set ticscale 0.5
#set pointsize 0.4
#
set ytics nomirror
set xlabel 'Queue size L (server included)'
set xrange [1:200]
set title 'Loss probability evaluation M/M/1 vs. M/M/1/k'
#
set yrange [0.00000001:1]
set logscale y
set ylabel 'P_L' 0.0
a1=0.5
f05(x) = ((1-a1)*(a1**x))/(1-a1**(x+1)) 
g05(x) = a1**(x+1) 
a5=0.65
f06(x) = ((1-a5)*(a5**x))/(1-a5**(x+1)) 
g06(x) = a5**(x+1) 
a2=0.8
f08(x) = ((1-a2)*(a2**x))/(1-a2**(x+1)) 
g08(x) = a2**(x+1) 
a3=0.9
f09(x) = ((1-a3)*(a3**x))/(1-a3**(x+1)) 
g09(x) = a3**(x+1)
a4=0.95
f095(x) = ((1-a4)*(a4**x))/(1-a4**(x+1)) 
g095(x) = a4**(x+1) 
plot f05(x) t "M/M/1/L a=0.5" w l 1, \
     g05(x) t "M/M/1 a=0.5" w lp 1, \
     f06(x) t "M/M/1/L a=0.65" w l 5, \
     g06(x) t "M/M/1 a=0.65" w lp 5, \
     f08(x) t "M/M/1/L a=0.8" w l 2, \
     g08(x) t "M/M/1 a=0.8" w lp 2, \
     f09(x) t "M/M/1/L a=0.9" w l 3, \
     g09(x) t "M/M/1 a=0.9" w lp 3, \
     f095(x) t "M/M/1/L a=0.95" w l 4, \
     g095(x) t "M/M/1 a=0.95" w lp 4
########
set xrange [1:20]
set title 'Loss probability evaluation M/M/1 vs. M/M/1/k'
#
set out 'MMK_compare_zoom.ps'
set key bottom left
set yrange [0.0001:1]
plot f05(x) t "M/M/1/L a=0.5" w l 1, \
     g05(x) t "M/M/1 a=0.5" w lp 1, \
     f06(x) t "M/M/1/L a=0.65" w l 5, \
     g06(x) t "M/M/1 a=0.65" w lp 5, \
     f08(x) t "M/M/1/L a=0.8" w l 2, \
     g08(x) t "M/M/1 a=0.8" w lp 2, \
     f09(x) t "M/M/1/L a=0.9" w l 3, \
     g09(x) t "M/M/1 a=0.9" w lp 3, \
     f095(x) t "M/M/1/L a=0.95" w l 4, \
     g095(x) t "M/M/1 a=0.95" w lp 4
exit