func LnBeta(x, y float64) float64 {
	z, _ := math.Lgamma(x + y)
	x, _ = math.Lgamma(x)
	y, _ = math.Lgamma(y)
	return x + y - z
}
// LnBeta computes the natural logarithm of the beta function.
LnBeta computes the natural logarithm of the beta function.
// LnBeta computes the natural logarithm of the beta function .
func LnBeta(x, y float64) float64 {
	z, _ := math.Lgamma(x + y)
	x, _ = math.Lgamma(x)
	y, _ = math.Lgamma(y)
	return x + y - z
}
func LnBeta ( x , y float64 ) float64 { z , _ := math . Lgamma ( x + y ) x , _ = math . Lgamma ( x ) y , _ = math . Lgamma ( y ) return x + y - z }
LnBeta computes the natural logarithm of the beta function.
LnBeta computes the natural logarithm of the beta function .