20/07/11 16:37:49.32 GlB5Fub6.net
# simulation for N coin flips with n(+) sequentical heads
sim <- function(N=1000,n=10,p=0.5,greater=TRUE){ # n+ sequential heads after N flips (greater=T)
x=rbinom(N,1,p)
(y=rle(x))
i=which(y$values==1)
ifelse(greater,max(y$lengths[i])>=n,max(y$lengths[i])==n