Anyone know a fast and concise way of calculating the Beta $B(a,b)$ function for smallish (<10) real $a$ and $b$.

For integer $a$ and $b$ I have:

$B(a,b) = \prod\limits_{j=1}^b \frac{j}{a+j}$

which has tiny code and is also pretty fast. I've see some methods that rely on the Gamma or log Gamma function:

$\log B(a,b) = \log \Gamma(a) + \log \Gamma(b) - \log \Gamma(a+b)$

using approximations of $\log \Gamma$, but I was wondering if there was a more direct way.