gammainc の結果です。
lower の対数です (lower がオーバーフローする箇所でも有限です)。
lower
upper の対数です。
upper
下側不完全ガンマ関数 gamma(a, x) = integral from 0 to x of t^(a-1) e^-t dt です。
上側不完全ガンマ関数 Gamma(a, x) = integral from x to Inf of t^(a-1) e^-t dt です。
gammainc の結果です。