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')