expdev.c

来自「数值算法」· C语言 代码 · 共 13 行

#include <math.h>

float expdev(idum)
int *idum;
{
	float dum,ran1();

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