ranpt.c

来自「可以用于科学计算的库」· 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);
}