Algorithm 1006

We present a computational procedure to evaluate the integral ∫ y x s p -1 e -μs ds for 0 ≤ x < y ≤ +∞,μ = ±1, p > 0, which generalizes the lower ( x =0) and upper ( y =+∞) incomplete gamma functions. To allow for large values of x , y , and p while avoiding under/overflow issues in the standard dou... View More