生成随机数（高斯分布）

高斯分布也称为正态分布（normal distribution）
常用的成熟的生成高斯分布随机数序列的方法由Marsaglia和Bray在1964年提出，C++版本如下：

#include   
#include   

double gaussrand()  
{  
    static double V1, V2, S;  
    static int phase = 0;  
    double X;  

    if ( phase == 0 ) {  
        do {  
            double U1 = (double)rand() / RAND_MAX;  
            double U2 = (double)rand() / RAND_MAX;  

            V1 = 2 * U1 - 1;  
            V2 = 2 * U2 - 1;  
            S = V1 * V1 + V2 * V2;  
        } while(S >= 1 || S == 0);  

        X = V1 * sqrt(-2 * log(S) / S);  
    } else  
        X = V2 * sqrt(-2 * log(S) / S);  

    phase = 1 - phase;  

    return X;  
}  

这样生成的高斯分布随机数序列的期望为0.0，标准差为1.0。若指定期望为E，标准差为V，则只需增加：
X = X * V + E;