谁能解以正态分布产生数据的算法

[fuladuo]

举例：利用该算法读取1-10之间的这十个数字，被读取的次数符合正态分布

[dcyu]

double _random(void)
{
int a;
double r;

a=rand()%32767;
r=(a+0.00)/32767.00;

return r;

}

double _sta(double mu,double sigma)
{
double _random(void);

int i;
double r,sum=0.0;

if(sigma<=0.0) { printf("Sigma<=0.0 in _sta!"); exit(1); }
for(i=1;i<=12;i++)
sum = sum + _random();
r=(sum-6.00)*sigma+mu;

return r;

}
函数_sta可以产生均值为mu，方差为sigma的正态分布。
主要是利用了大数定理，产生了12个（0，1）分布的随机数进行计算。

[dcyu]

还有一种算法是书上介绍的，但是我测试过，取的次数很多的话，会有溢出的可能。如下：
double _sta2(double mu,double sigma)
{
double r1,r2;

r1=_random();
r2=_random();

return sqrt(-2*log(r1))*cos(2*M_PI*r2)*sigma+mu ;

}