% posterior for the data: 7000x heads and 1000x tails
close all;
x = 0:0.001:1;
figure;
area(x,betapdf(x,5+7000,5+1000));
axis([-0.25 1.25 0 5]); xlabel('\theta');ylabel('p(\theta)'); title('posterior');