beta(a,b)

Returns the Beta function of reals a and b.
The Beta function is define as:

   beta(a,b)=gamma(a)*gamma(b)/gamma(a+b)

where gamma(a) returns the Gamma function of the real a.