徐士良数学计算常用算法集（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				      { *r=u*(*r)+v; m=(int)(*r/s);
				        *r=*r-m*s; p[i]=*r/s;
				      }
				    return;
				  }