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