# ポアソン乱数（ポアソン定数 lambda）

function p_rand(lambda,     k)
{
   lambda = exp(lambda) * rand()
   k = 0
   while (lambda > 1) {
       lambda *= rand()
       k++
   }
   return k
}