\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", xlim = c(-20,30))
      }
  }
\end{Sinput}
\end{Schunk}
