clc; clearvars;
a = 1:0.1:100;
la = length(a);
n = zeros(la,1);
Bc_obj = 0.1;
for j = 1:la
    ti = a(j);
    n_ch = 1;
    Bc = erlang_r( ti,n_ch );
    if Bc < Bc_obj
        n(j) = n_ch;
    else
        while Bc>=Bc_obj
            n_ch = n_ch + 1;
            Bc = erlang_r( ti,n_ch );
        end
    end
    n(j) = n_ch;
end
f = plot(a',n);
f(1).LineWidth = 1.5; 
% f(1).Marker = '.';
grid on
title('Capacity vs. Traffic intensity')
xlabel({'Traffic intensity', 'a'}), ylabel({'Capacity', 'n'})
% legend({['n = ',num2str(n(1))],['n = ',num2str(n(2))],...
%     ['n = ',num2str(n(3))],['n = ',num2str(n(4))]}, 'Location','northeast')