expdev.cpp

来自「C++数值算法原书代码」· C++ 代码 · 共 14 行

#include <cmath>
#include "nr.h"
using namespace std;

DP NR::expdev(int &idum)
{
	DP dum;

	do
		dum=ran1(idum);
	while (dum == 0.0);
	return -log(dum);
}