\begin{Schunk}
\begin{Sinput}
  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}
