Convolution of multiple probability density functions I have a series of tasks where when one task finishes the next task runs, until all of the tasks are done. I need to find the probability that everything will be finished at different points in time. How should I approach this? Is there a way to find this in polynomial time?
The pdfs for how long individual tasks will run have been found experimentally and are not guaranteed to follow any particular type of distribution.

I don't know if that is what you are looking for, but if $X_1,\ldots,X_n$ are i.i.d. random variables with mean $\mu$ and (finite) variance $\sigma^2$, then
$$
{\rm P}\bigg(\sum\limits_{i = 1}^n {X_i }  \le t \bigg) = {\rm P}\bigg(\frac{{\sum\nolimits_{i = 1}^n {X_i }  - n\mu }}{{\sigma \sqrt n }} \le \frac{{t - n\mu }}{{\sigma \sqrt n }}\bigg),
$$
and the random variable $\frac{{\sum\nolimits_{i = 1}^n {X_i }  - n\mu }}{{\sigma \sqrt n }}$ converges to the standard normal distribution as $n \to \infty$. Now, you can estimate the unknown parameters $\mu$ and $\sigma^2$ (assuming the variance is finite), so if $n$ is sufficiently large, you are actually done.
For further details see this (note the subsection Density functions). 
