Function gammp 

pub fn gammp(a: f64, x: f64) -> f64

Regularized lower incomplete gamma P(a, x) = γ(a,x)/Γ(a), a > 0, x ≥ 0. Series for x < a+1, continued fraction (via Q) otherwise.