/*   gauss1.c   產(chǎn)生高斯隨機數(shù)         */

/*由uniform.c產(chǎn)生（0，1）上的n個相互獨立的均勻分布隨機數(shù)，均值為1/2,方差為1/12*/
/*根據(jù)中心極限定理，當(dāng)n較大時近似有Z=sqrt(12/n)*((r1+r2+...+rn)-n/2) ~ N(0,1)*/
/*取n=12，則Z=(r1+r1+..+r12)-6 ~ N(0,1)*/
/*再通過變換X=mean+Z*sigma得到均值為mean方差為sigma^2的高斯分布隨機數(shù)X*/

#include "uniform.c"

double gauss1(mean,sigma,s)
double mean,sigma;
long int *s;
{
 int i;
 double x,y;
 double uniform();
 for(x=0,i=0;i<12;i++)
   x+=uniform(0.0,1.0,s);
 x=x-6.0;
 y=mean+x*sigma;
 return(y);
}