ポアソン分布の乱数を得るには？

λがnum*pであるようなポアソン分布の乱数を得るには、

double lambda=num*p;
int k=0;
double l=exp(-lambda);
double m=1;
do
{
    k++;
    m*=rand()/double(RAND_MAX);
}while(m>=l);
int po=k-1;

poにポアソン分布の乱数が入る。

最終更新日 2009-06-08

本文中のFC4はFedora ProjectのFedora Core 4を、FC5はFedora Core 5を、FC6はFedora Core 6をopenSUSEはNovellのSUSE Linux OSSを表します。Fedora7以降は、単にFedora7、Fedora8、Fedora9、Fedora10、Fedora11、Fedora12、Fedora13、Fedora14、Fedora15と表示しています。Ubuntuは、必要に応じて20.04、21.04のようにバージョン番号をつけて区別しています。