01sub.c

来自「我们的大作业,包括runge-kutta法仿真,0-1分布,反变换的C语言原程序」· C语言 代码 · 共 16 行

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