expdev.txt

教你怎么用C++产各种各样的随机数

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