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
}