draw a large sample from the prior (the "Sampling" step)
p.prior=rbeta(N,a,b)