Suppose $X, Y$ are independent random variables with $X\sim NB(r, p)$ and $Y\sim NB(s, q)$. From a previous post, I understand that when the success probabilities are equal, $p = q$, then, $$ X+Y \sim NB(r + s, p). $$
What is the distribution of $X+Y$ in cases when $p \neq q$? If there is no analytical solution, I'd still be interested in an approximation.