from numpy import *
from scipy import stats
from matplotlib.pyplot import *
from numpy.linalg import norm

t = arange(.01, 1, .01)

prior = .1*stats.norm.pdf(t-.5)+.9*stats.norm.pdf( 100*(t-.5))
figure(1)
plot(t,prior)
title("Prior")

N = 60

likelihood = (t**N)*((1-t)**(100-N))
figure(2)
plot(t, likelihood)
title("Likelihood")

posterior = likelihood*prior
figure(3)
plot(t, posterior)
title("Posterior")