3.2 产生0到1之间均匀分布的随机数序列 rnds.c

许士良常用算法程序集C语言,包括c++一些常用算法代码

void rnds(r,p,n)
double *r,p[];
int n;
{ 
	int i,m;
    double s,u,v;
    s=65536.0; u=2053.0; v=13849.0;
    for (i=0; i<=n-1; i++)
    { 
		*r=u*(*r)+v; m=(int)(*r/s);
        *r=*r-m*s; p[i]=*r/s;
    }
    return;
}