\chapter{برنامه‌های کاربردی}
\begin{LTR}
\begin{verbatim}
 x=c(2.4, 2.42, 3.17, 3.75, 4.65, 4.95, 6.23, 6.68, 7.3(
 y=c(15.17,19.87,20.18,21.5,21.88,22.23,23.02,23.9,28.17,29.7(
 x.new=x*mean(x)/var(x(
 y.new=y*mean(y)/var(y(
 z=x.new
 l=function(alpha=5(
 -log(prod(dgamma(z,shape=alpha,rate=1(((
 library(stats4(
 summary(mle(l((
 \#output(z=x.new(
 Maximum likelihood estimation
 Call:
 mle(minuslogl = l(
 Coefficients:
       Estimate Std. Error
 alpha  6.40223  0.8107776
 -2 log L: 40.74708 
 \#output(z=y.new(
 Maximum likelihood estimation
 Call:
 mle(minuslogl = l(
 Coefficients:
       Estimate Std. Error
 alpha 29.85472   1.713405
 -2 log L: 61.30076 
\end{verbatim}
\end{LTR}