In my work I came across the following integral which stems from basically computing the statistics of an output signal of a nonlinear system given white noise at the input, so the integral is (if I did not make a mistake) represent a probability.
$$\int_{0}^{\infty} \frac{\lambda_i^{N_i} x^{N_i-1} e^{-\lambda_i x}}{(N_i-1)!} \prod_{j=1,j\neq i}^{N} (e^{-\lambda_j x} \sum_{m=0}^{N_j-1}\frac{(\lambda_j x)^m}{m!} ) \prod_{l=N+1}^M (1-e^{-\lambda_l x} \sum_{n=0}^{N_l - 1} \frac{(\lambda_n x)^n}{n!}) dx $$
$N_i$ are natural numbers in the range of 1 to 8, $M$ is a natural in the range of 20, $N$ equals 8. $\lambda_i$ are positive real values. There is actually another sum across $i$, however, due to linearity of the integral it is of no interest here. That is where the i index comes from.
I had a look at the residuum theorem, but nothing seems to exist for this complicated case which could help solving this integral. One can formally compute the products involving the polynomials and then integrate each sum, however, this gets very difficult to tract afterwards if I want to evaluate this as the number of summands, without any relation between the resulting coefficients, is rather large.
Anyone got an idea how to solve this one?