ranpt.c

适合大型数值计算代码 现在网络上已经找不到了 购买需要20$

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);
}