% Capacity vs. Traffic intensity given QoS
clc; clearvars; close all
a = 0.1:0.1:100; % traffic intensity
la = length(a);
Bc_obj = [0.2; 0.1; 0.01]; %blocking probability objective
lbp = length(Bc_obj);
n = zeros(la,lbp); %to store the number of channels the system needs given Bc

for bp = 1:lbp %recorrer los requerimiento de QoS planteados
    for j = 1:la %recorre el vector de intensidad de trafico
        ti = a(j);
        n_ch = 1; %proceso de biusqueda para determinar el n que cumpla la condicion
        Bc = erlang_r( ti,n_ch );
        if Bc < Bc_obj(bp)
            n(j,bp) = n_ch;
        else
            while Bc>=Bc_obj(bp)
                n_ch = n_ch + 1;
                Bc = erlang_r( ti,n_ch );
            end
        end
        n(j,bp) = n_ch;
    end
end

f = plot(a',n);
f(1).LineWidth = 1.5;
% f(1).Marker = '.';
f(2).LineWidth = 1.5;
f(3).LineWidth = 1.5;
grid on
title('Capacity vs. Traffic intensity')
xlabel({'Traffic intensity', 'a'}), ylabel({'Capacity', 'n'})
legend({['Bc_{obj} = ',num2str(Bc_obj(1))],['Bc_{obj} = ',num2str(Bc_obj(2))],...
    ['Bc_{obj} = ',num2str(Bc_obj(3))]}, 'Location','northeast')