#include #include int ix=1,iy=1; main() { int i,n=100; double urnd(void); for(i=1;i<=10000;i++) urnd(); for(i=1;i<=n;i++) printf("%lf\n",urnd()); } /* ================================================ */ float urnd(void) { int kx,ky; double rn; /* Input: ix, iy: Seeds Output: rn: Uniform Random Draw U(0,1) */ kx=ix/53668; ix=40014*(ix-kx*53668)-kx*12211; ky=iy/52774; iy=40692*(iy-ky*52774)-ky*3791; rn=(double)(ix-iy)/2147483563.; rn-=(int)rn; if( rn<0.) rn++; return rn; }