p = [0:.01:1];
L = (p.^2).*((1-p).^3);
plot(p,L);
xlabel('\pi')
ylabel('Likelihood L(\pi)')