Let X_1,....,X_n be a random sample from a Bernoulli distribution, X_i\sim BIN(1, \theta) and assume a uniform prior \theta \sim UNIF(0,1). Derive a 100(1-\alpha)% Bayesian interval estimate of \theta.