unnorm_posterior = function(theta) {theta^y * (1 - theta)^(n-y)}
n = 10
y = rbinom(1,n,0.6)
curve(unnorm_posterior,from = 0,to = 1, main=paste("n = ",n,"; y = ",y))
n = 100
y = rbinom(1,n,0.6)
curve(unnorm_posterior,from = 0,to = 1, main=paste("n = ",n,"; y = ",y))
n = 1000
y = rbinom(1,n,0.6)
curve(unnorm_posterior,from = 0,to = 1, main=paste("n = ",n,"; y = ",y))