一様乱数法による生成で、genrand_res53()が[0,1)の一様乱数に対応する。 int Poisson(double lambda){ int k=0; lambda=exp(lambda)*(genrand_res53()); while (lambda>1) { lambda*=(genrand_res53()); k++; } return k; } しかし、lambdaの値が大きくなる…