\begin{Schunk}
\begin{Sinput}
 x <- seq(mu + -3.5 * sigma, mu + 3.5 * sigma, by = 0.01)
 xd <- dnorm(x, m = mu, s = sigma)
 plot(xd ~ x, type = "l", xlab = "x", ylab = "Probability Density of Normal")
\end{Sinput}
\end{Schunk}
