clc; clearvars;
A = 1:50; % traffic intensity
lA = length(A);
k = 10:10:40;
lk = length(k);
Pb = zeros(lA,lk);
for i = 1:lk
    for j = 1:lA
        Pb(j,i) = erlang_r(A(j),k(i));       
    end
end

f = plot(Pb);
f(1).LineWidth = 1.5;
f(1).Marker = '*';
f(2).LineStyle = ':';
f(2).Marker = 's';
f(2).LineWidth = 1.5;
f(3).LineStyle = '-.';
f(3).Marker = 'd';
f(3).LineWidth = 1.5;
f(4).LineStyle = '--';
f(4).Marker = 'o';
f(4).LineWidth = 1.5;
grid on
title('Pb ')
xlabel('Traffic Intensity')
ylabel('Blocking Probability')