function Traffic_Density = ErLang_C_CalA(num_of_ch, Probablity_of_delay)
Traffic_Density = 1.0;
while Traffic_Density <= 800.0
sumValue = 1;
for loop = 1 : num_of_ch
item = 1;
for subcount = 0 : (loop - 1)
item = item * Traffic_Density./(loop - subcount);
end
sumValue = sumValue + item;
end
% Calculate A^Kcell / Kcell!
numValue = 1;
for count = 0 : (num_of_ch - 1)
numValue = numValue * (Traffic_Density./(num_of_ch - count));
end
P_kcell = numValue * (1 / (sumValue + numValue / (1 - Traffic_Density / num_of_ch)));
Probability = P_kcell / (1 - Traffic_Density / num_of_ch);
if (Probability >= Probablity_of_delay) && (Probability > 0)
Traffic_Density = Traffic_Density - 0.1;
break;
else
Traffic_Density = Traffic_Density + 0.1;
end
end
if (Traffic_Density >= 800)
Traffic_Density = 0;
end
end