I have optimized the numerical evaluation for float and complex float values of the gamma_incomplete_regularized function and checked it out. The problem 15 in rtest_distrib.mac no longer fails. Dieter Kaiser