产生0-1之间均匀分布一个随机数.c

//产生0-1之间均匀分布一个随机数// /////////////////////////////////////////////////////////////// //传入参数:r--双精度实型

//////////////////////////////////////////////////////////////////////
  //             产生0-1之间均匀分布一个随机数               //
  /////////////////////////////////////////////////////////////
  //方法：参见《常用算法程序集》317页，徐士良。				 //
  //传入参数:r--双精度实型变量指针,指向的单元存放随机数种子值//                                        //
  //传出参数:p:产生0-1之间均匀分布一个随机数                 //
//////////////////////////////////////////////////////////////////////
double CAK630Doc::mrnd1(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);
  }