%
% Estudio del modelo de colas M/M/1/K
%
% 2006 Javier Valcarce García
%

p = 0:0.01:5;
k = [1 2 5 10 30];

KSIZE = length(k);
PSIZE = length(p);

B     = zeros(KSIZE, PSIZE);

clf
hold on

color  = ['b', 'g', 'r', 'm', 'y', 'k'];
marker = ['.', 'v', 'o', '+', '.', 'x'];

for i = 1:KSIZE
    
    B(i, :) = (p .^ k(i)) .* (1-p) ./ (1-p .^ (k(i) + 1));
    
    h = plot(p, B(i, :));
    %set(h, 'Marker', marker(i));
    set(h, 'Color', color(i));
    set(h, 'DisplayName', sprintf('K = %i', k(i)));
    
end

title('Probabilidad de bloqueo en el modelo M/M/1/K');
ylabel('Probabilidad de bloqueo (B)');
xlabel('Intensidad de tráfico ofrecido (\lambda / \mu) ');