\begin{Schunk}
\begin{Sinput}
 par(mfcol=c(4,2))
 for ( i in 1:8){
      shape <- c(0.5,1,2,3,5,10,15,30)
      scale <- 1
      xrange<-seq(from = 0, to = 50, by = 0.1)
      mainlabel<- paste("Gamma(shape=", shape[i],",scale=",scale,",)",sep="")
      xprob <- dgamma(xrange, shape = shape[i], scale = scale, log = F)
      plot(xrange, xprob, type = "l", main = mainlabel, xlab = "possible values of x", ylab = "probability of x")
  }
\end{Sinput}
\end{Schunk}
