/*  uniform.c   產生均勻分布的偽隨機數           */
/*  引入迭代公式Xj+1=(aXj+c)_(mod m)             */
/*      再取       Yj=Xj/m,  j=1,2,...,          */

double uniform(a,b,seed)
double a,b;
long int *seed;
{ double t;
    *seed=2045*(*seed)+1;                  /*取a=2045,c=1*//*m=2^20=1048576*/
    *seed=*seed-(*seed/1048576)*1048576;
   t=(*seed)/1048576.0;                    /*m=2^20=1048576*/
   t=a+(b-a)*t;                            /*通過變換Zj=a+(b-a)*Yj產生(a,b)區間上的隨機數*/
   return(t); 
}