Suppose $X$ has a gamma distribution. Prove that its moments are given by
$$
E\left[X^n\right]=\left[\prod_{k=1}^{n-1}\left(1+k C_X^2\right)\right] E[X]^n \quad n=2,3, \cdots .
$$
Since the squared coefficient of variation is given by
$$
C_X^2=\frac{1}{\beta},
$$
this means that we can write the above formula as
$$
E\left[X^n\right]=\frac{\beta(\beta+1)(\beta+2) \cdots(\beta+n-1)}{\alpha^n}, \quad n=1,2,3, \cdots .
$$
This result also implies that for the Erlang- $k$ random variable, we have
$$
E\left[X^n\right]=\left(1+\frac{1}{k}\right)\left(1+\frac{2}{k}\right) \cdots\left(1+\frac{n-1}{k}\right) E[X]^n, \quad n=1,2,3, \cdots .
$$