func Gammastar(x float64) float64 { return float64(C.gsl_sf_gammastar(C.double(x))) }
/* Regulated Gamma Function, x > 0 * Gamma^*(x) = Gamma(x)/(Sqrt[2Pi] x^(x-1/2) exp(-x)) * = (1 + 1/(12x) + ...), x->Inf */ func GammaStar(x float64) (r float64) { r = float64(C.gsl_sf_gammastar(C.double(x))) return }