なぜ不偏分散はn-1で割るのかをシミュレーションで！

for i=1:500

R=randn(10,1);    % 10件の正規乱数～N(0,1)を生成 (平均０、分散1 の正規乱数)

Bunsan(i)=sum((R-ones(10,1)*mean(R)).^2)/10;   % /9 にすると不偏分散になる

HeikinBunsan(i)=mean(Bunsan(1:i));

end

plot(HeikinBunsan)