Prob&stats

/// functie coin toss
n=40000;
u=rand(3,n);
y=u<0.5;
z=sum(y);
hist(z);

/// d
stairs(x,binocdf(x,2,0.5),"m");
stairs(x,binocdf(x,1,0.5),"m");