clear;

n = 1000;	// n回コイントスを行う
m = 1200;	// n回のコイントスをN回繰り返す

p = 0.5;

X = bool2s(rand(m,n) < p);
hatS = 1 / n * sum(X,'c');
Z = n / sqrt(n * p * (1 - p)) * (hatS - p);

// *** グラフのプロット ***
// 規格化された正規分布関数
function y = gauss(x)
  y = 1 / sqrt(2 * %pi) * exp(- 1 / 2 * x .^ 2);
endfunction
x = linspace(-4,4);
plot(x,gauss(x));

// ヒストグラム
histplot(100,Z);