plot(function(p) dbeta(p,4,3), xlab=expression(theta), ylab="Density")
lines(c(0,1), c(1,1), lty=2)
legend("topleft", c("Prior", "Posterior"), lty=c(2,1))