% p_poisson( n, rate) - Poisson probability
%
function p = p_poisson( n, rate)

p = exp(-rate);

if( n > 0)
  p = p * prod(rate*ones(1,n)) / prod(1:n);
end
