ranpt.c

来自「numerical recipes in C source codes」· C语言 代码 · 共 12 行

extern long idum;

void ranpt(float pt[], float regn[], int n)
{
	float ran1(long *idum);
	int j;

	for (j=1;j<=n;j++)
		pt[j]=regn[j]+(regn[n+j]-regn[j])*ran1(&idum);
}
/* (C) Copr. 1986-92 Numerical Recipes Software "U,6VV'. */