Computes the Beta Function, B(a,b) = \Gamma(a)\Gamma(b)/\Gamma(a+b) for a > 0, b > 0.
err contains an estimate of the absolute error in the value z.
This function is from the GNU Scientific Library, see http://www.gnu.org/software/gsl/ for documentation.
Package: gsl