func GammaInc(a float64, x float64) float64 { return float64(C.gsl_sf_gamma_inc(C.double(a), C.double(x))) }
/* Non-normalized Incomplete Gamma Function * * Gamma(a,x) := Integral[ t^(a-1) e^(-t), {t,x,Infinity} ] * * x >= 0.0 * Gamma(a, 0) := Gamma(a) * */ func GammaInc(a, x float64) (r float64) { r = float64(C.gsl_sf_gamma_inc(C.double(a), C.double(x))) return }