3rnds.c

来自「程序用于产生0到1之间均匀分布的随机数序列」· C语言 代码 · 共 15 行

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