By trying to find a marginal distribution I came accross integration of the product series. For the sake of generality, lets assume the integral is of following form: $$\int \prod_{k=1}^{n}\left ( x+a_{k} \right )^{b_{k}}dx.$$ $a_{k}$ is a real coefficient and $b_{k}$ is positive integers. Is there any method, that could be used to integrate this analyticaly? It is of no problem to calculate it by some numerical method. But the first thing I would like to try is to find its analytical expression.

Since the derivative of such a product has a similar form $$\left(\sum_{k=1}^n\frac{b_k}{x+a_k} \right)\prod_{k=1}^{n}\left ( x+a_{k} \right )^{b_{k}},$$ One might hope that something similar is true of some antiderivative. (I am implicitly assuming that the $b_k$ are positive integers, although the formula is valid more generally. However there are cases that will involve some $\ln x+a_k$ . I will continue to neglect these.) If we imagine that we have some one antiderivative $$F(x)=\int \prod_{k=1}^{n}\left ( x+a_{k} \right )^{b_{k}}dx$$ and let $F_k=F(x)F(a_k),$ then $F_k(x)$ is also an antiderivative and is divisible by $(x+a_k)^{b_k+1}$ That is not particularly deep. We can use integration by parts to obtain an answer without any explicit multiplying although marching through to the end will not be very pretty. It is more satisfying to know that we can do it than to actually do so.
It would be only slightly more complicated (but also only slightly more satisfying) to have products of factors $(m_kx+a_k)^{b_k}$. A large table of integrals might have formulas such as this one. 


The magic words are "Lauricella Hypergeometric Function". 

