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