function [ Pb ] = erlang( A,k )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
num = (A.^k)./factorial(k);
sum = 0;
for i = 0:k
    sum = sum + (A.^i)./factorial(i);
end
Pb = num/sum;
end