Returns the value of the Beta function, defined in terms of the complete
gamma function, G, as: G(z)G(w)/G(z+w). The implementation follows Numerical
Recipes in C, section 6.1.