3rnd1.c

来自「1.书中所有的算法均用C语言描述」· C语言 代码 · 共 13 行

  double rnd1(r)
  double *r;
  { int m;
    double s,u,v,p;
    s=65536.0; u=2053.0; v=13849.0;
    m=(int)(*r/s); *r=*r-m*s;
    *r=u*(*r)+v; m=(int)(*r/s);
    *r=*r-m*s; p=*r/s;
    return(p);
  }