\begin{Schunk}
\begin{Sinput}
 x <- seq(from=00.01, to=10, by=0.1)
 y <- dgamma(x,shape=2,scale=1) 
 plot(x, y, type="l",xlab="x",ylab="",main="")
 polygon(x=c(0,x,10),y=c(0,y,0), col=gray(0.850))
\end{Sinput}
\end{Schunk}
