expdev.c

Numerical Recipes Software 提供的算法子程序集

#include <math.h>

float expdev(idum)
long *idum;
{
	float ran1();
	float dum;

	do
		dum=ran1(idum);
	while (dum == 0.0);
	return -log(dum);
}
/* (C) Copr. 1986-92 Numerical Recipes Software . */