erlang_b.m

把地理区域分成一个个小区

function erb=erlang_b(A,c)
num=A^c;
sum=0;
for k=0:c
    kfact=prod(1:k);
    term=(A^k)/kfact;
    sum=sum+term;
end
cfact=prod(1:c);
den=cfact*sum;
erb=num/den;