ranpt.h

经典numerical receip 配套代码

void ranpt(VecDoub_O &pt, VecDoub_I &regn) {
	static const int RANSEED=5331;
	static Ran ran(RANSEED);
	Int j,n=pt.size();
	for (j=0;j<n;j++) pt[j]=regn[j]+(regn[n+j]-regn[j])*ran.doub();
}