I am stuck with the following equation (1) that involves summation and integration both.
$$P_e = \sum_{i=0}^{M-r}b_{r,i}\int_0^\infty\!\!\left(\frac{1}{2}-\frac{\gamma(p,qz)}{2\Gamma(p)}\right)\frac{e^{-\mu_{r,i}\sqrt{z}}}{2\sqrt{z}}\,dz \qquad(1)$$
where, $b_{r,i}$, $\mu_{r,i}$ do not depend on $z$. $p,q$ are positive constants. $\gamma(\cdot,\cdot)$ is lower incomplete gamma function.
Any help in this regard is highly appreciated.