func Lnbeta(a float64, b float64) float64 {
return float64(C.gsl_sf_lnbeta(C.double(a), C.double(b)))
}
/* Logarithm of Beta Function
* Log[B(a,b)]
*
* a > 0, b > 0
*/
func LnBeta(a, b float64) (r float64) {
r = float64(C.gsl_sf_lnbeta(C.double(a), C.double(b)))
return
}