Let $X \sim \Gamma(m, p)$ with a shape parameter $m$ and a scale parameter $p$ and $Y \sim \Gamma(m, q)$ with a shape parameter $m$ and a scale parameter $q$, and let $X$ and $Y$ be independent.
What will be the PDF and CDF of $X+Y$? How can I solve it in MatLab?