function y = bern(x,n,k)
c = nchoosek(n,k);
y = c*x.^k.*(1-x).^(n-k);
