function p=erlangB(A,N)
  rec=1;
  for x=1:N
    rec=1+x*rec/A;
  endfor;
  p=1/rec;
endfunction