Figure 4.9:

Histogram of 10,000 samples of \displaystyle y=\DOTSB \sum@ \slimits@ _{i=1}^{10} x_i .

Code for Figure 4.9

Text of the GNU GPL.

main.m

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
nsam = 10000;
nsum = 10;
x = rand(nsum,nsam);
y   = sum(x);

nbins = 20;

figure()
hist (y,nbins);
figure()
hist (x(1,:), nbins);

[ny, ysam] = hist(y,nbins);
[ysamb, nyb] = obar(ysam, ny);
[nx1, x1sam] = hist(x(1,:), nbins);
[x1samb, nx1b] = obar(x1sam, nx1);

data = [ysamb, nyb, x1samb, nx1b];

mean(x(1,:))
var(x(1,:))
mean(y)
var(y)

save cenlim.dat data