More work on optimising performance of binomial distribution call
Extract original impls of binomial and binomial_coef.  Restore parens around ** terms