Matlab生成正态分布随机数问题

题目是
利用均匀分布 U（–1，1）,使用中心极限定理产生正态分布的随机数1000个,其中每次用于产生正态分布随机数的均匀分布样本容量 n >50

n个均匀分布 U（–1，1）数的平均数期望和方差为0和(1-(-1))^2/(12*n)=1/(3*n)的正态分布
所以

a=zeros(1000,1);
n = 60;    %n取60>50
for i = 1:1:size(a,1)
x = 2*rand(n,1)-1;
a(i) = sum(x)/n/sqrt(1/(3*n));   %产生第i个随机数
end
histogram(a)   %画一下分布