简单几行C语言代码实现高斯分布

由于项目需要产生一个高斯分布，所以去网上查找，结果很多，但是都没有足够的注释或者原理讲解，所以大部分代码都看不懂，也没法确定代码的结果是否正确。因此想从原理上来解决这个问题。具体的理论推导请看 http://blog.sina.com.cn/s/blog_9ce5a1b501018c1g.html。 我这里只是简单的提供我测试这个原理的代码以及测试结果。测试是用MATLAB仿真的：

clear all;clc
N=500000;
mu   =0;  %高斯均值
theda=1; %高斯方差
%产生两个均匀分布的0~1的随机序列
t1=rand(1,N);
t2=rand(1,N);
%极坐标的两个随机变量分布序列
a = 2*pi*t1;            %a是极坐标的角度：变成了0~2*pi的均匀分布
r = sqrt(-2*log(t2));   %r是极坐标的距离：
%用极坐标(a,r)转换成笛卡尔坐标(x,y)，这里x,y都符合高斯分布
for i=1:N
    x(i) = mu+sqrt(theda)*r(i)*cos(a(i));
    y(i) = mu+sqrt(theda)*r(i)*sin(a(i));
end
x = 10*x;
mi = floor(min(x));
ma = ceil(max(x));
len= abs(mi)+abs(ma)+1;
array = zeros(1,len);
for i=1:N
    array(floor(x(i))-mi+1) = array(floor(x(i))-mi+1)+1;
end