expdev.cpp

来自「工程算法 这是一个很有用的工程数值算法集锦」· C++ 代码 · 共 11 行

double expdev(int & idum)
{
	double temp;
    double dum;
	dum=0.0;
    while (dum == 0)
        dum = ran1(idum);
    temp = -log(dum);
	return temp;
}