\begin{Schunk}
\begin{Sinput}
  shape <- 2.2; scale <- 2.5
  xrange <- c(0.0001, 5*shape*scale)
  xseq <- seq( xrange[1], xrange[2], length.out = 200)
  dxseq <- dgamma(xseq, shape=shape, scale=scale)
  plot(dxseq ~ xseq, type="l", ylim=c(0,1.2* max(dxseq)), xlab="x",ylab=paste("dgamma(x, shape=",shape,", scale=",scale,")"), main=paste("PDF of Gamma(",shape,",",scale,")"))
\end{Sinput}
\end{Schunk}
