function lb = log_binom_approx(n,k) 

    lb = (gammaln(n+1) - gammaln(k+1) - gammaln(n-k+1));

end