expdev.txt

来自「教你怎么用C++产各种各样的随机数」· 文本 代码 · 共 10 行

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