# File lib/standard/facets/math/beta.rb, line 9
  def self.beta(x, y)
    #exp(lgamma(x).first + lgamma(y).first - lgamma(x+y).first)
    tgamma(x) * tgamma(y) / tgamma(x + y)
  end