function val=g_beta(theta,par)

eta=exp(par(1))/(1+exp(par(1))); K=exp(par(2));

val=-betaln(K*eta,K*(1-eta))+K*eta*theta-K*log(1+exp(theta));

