library(matlib)
\[ px(j)=\binom{n}{j} p^{j}q^{n-j}for\ 0\leq j\leq n\ (binomial\ distribution)\ then \] \[ g(t)=\sum^{n}_{j=0} e^{tj}\binom{n}{j} p^{j}q^{n-j} \] \[ =\sum^{n}_{j=0} \binom{n}{j} (pe^{t})^{j}q^{n-j} \] \[ =(pe^{t}+q)^{n} \] Note that \[ \mu_{1} =g\prime (0)=n(pe^{t}+q)^{n-1}pe^{t}\mid_{t=0} =np,\] \[ \mu_{2} =g^{\prime \prime }(0)=n(n-1)p^{2}+np, \]
\[ so\ that\ \mu =\mu_{1} =np,and\ \sigma^{2} =\mu_{2} -\mu^{2}_{1} =np(1-p),as\ expected. \]
Let X have range [0,∞) and density function fX(x) = λe−λx (exponential density with parameter λ). In this case \[ \mu_{n} =\int^{\infty }_{0} x^{n}\lambda e^{-\lambda x}dx=\lambda (-1)^{n}\frac{d^{n}}{d\lambda^{n} } \int^{\infty }_{0} e^{-\lambda x}dx \] \[ g(t)=\sum^{\infty }_{k=0} \frac{\mu_{k} t^{k}}{k!} \] \[ =\sum^{\infty }_{k=0} [\frac{t}{\lambda } ]^{k}=\frac{\lambda }{\lambda -t} . \] Here the series converges only for |t| < λ. Alternatively, we have \[ g(t)=\int^{\infty }_{0} e^{tx}\lambda e^{-\lambda x}dx \] \[ =\frac{\lambda e^{(t-\lambda )x}}{t-\lambda } |^{\infty }_{0}=\frac{\lambda }{\lambda -t} \] Now we can verify directly that \[ \mu_{n} =g^{(n)}(0)=\frac{\lambda n!}{(\lambda -t)^{n+1}} |_{t=0}=\frac{n!}{\lambda^{n} } . \]