integrate x^(n-1) exp(-x-a/x)/Gamma(n) between 0 and infinity