\begin{Schunk}
\begin{Sinput}
 par(mfrow=c(4,2))
 for (i in 1:4){
      for (j in 1:2) {
         mu <- 3*j
         sigma <- 2*i
         xrange <- seq(from=mu-3*sigma,to=mu+3*sigma,by=0.2)
         mainlabel <- paste("N(",mu,",",sigma*sigma,",)",sep="")
         xprob <- dnorm(xrange, mean=mu, sd=sigma, log = FALSE)
         plot(xrange, xprob, type="l", main=mainlabel, xlab="possible values of x", ylab="probability of x")
      }
  }
\end{Sinput}
\end{Schunk}
