let gamma_p a x =
    assert (x >= 0. && a > 0.);
    if x < a +. 1.
    then gamma_0_x a x
    else 1. -. gamma_x_infty a x