function val=f_binomial(data,theta)

y=data(:,1); n=data(:,2);

val=n.*y.*theta-n.*log(1+exp(theta));

