f = dfuni(0,1);
f2 = resize( convpf(f,f), 0,2);
f4 = resize( convpf(f2,f2), 0,4);
f8 = resize( convpf(f4,f4), 0,8);
f16 = resize( convpf(f8,f8), 0,16);
f32 = resize( convpf(f16,f16), 0,32);
f48 = resize( convpf(f32,f16), 0,48);
plotf(f48,'m')
F48 = getcf(f48);
pconv = pval(F48,25) - pval(F48,22)

rand('seed',0)
rand('seed')

ans =

   931316785


s1 = rand(8,1000);
rand('seed')

ans =

    65048263

s2 = rand(8,1000);
s3 = rand(8,1000);
s4 = rand(8,1000);
s5 = rand(8,1000);
s6 = rand(8,1000);
stot = s1+s2+s3+s4+s5+s6;
colsum = sum(stot);
psim = sum( colsum>=22 & colsum<=25 ) / length(colsum)

psim =

    0.5310

f48mod = [ f48(1,:) ; ((xvals(f48) - 24) / 2)' (pvals(f48)*2)' ];
plotf(f48mod,'m')
hold on
plotf( dfnorm(0,1),'g:' )
binom40 = mfbino(40,.2)

plotf(binom40)
hold off
m = mfberni(.2)

m2 = convpf(m,m)

m4 = convpf(m2,m2)

m8 = convpf(m4,m4)

m16 = convpf(m8,m8)

m32 = convpf(m16,m16)

m40 = convpf(m8,m32);
plotf(m40)
diff = binom40(:,2) - m40(:,2)

M40 = getcf(m40);
p10 = pval(M40,10)

rand('seed',0)
sm1 = rand(200,30);
sm1 = sm1<=.2;
row1 = sum(sm1');
tot1 = sum(row1>8)

tot1 =

    22

% This was my result for a sample of 1000 (5 of the above):
.132

ans =

    0.1320

diary off