\begin{Schunk}
\begin{Sinput}
 x <- seq(0, 4, length.out = 200)
 y <- matrix(0, ncol = 5, nrow = 200)
 y[, 1] <- dt(x, df = 1)
 y[, 2] <- dt(x, df = 2)
 y[, 3] <- dt(x, df = 5)
 y[, 4] <- dt(x, df = 20)
 y[, 5] <- dt(x, df = 1000)
 matplot(x, y, type = "l", ylab = "probability density", col = "black")
 lines(x, dnorm(x), lty = 2, lwd = 3)
 text(0, 0.225, expression(t(nu == 1)), pos = 4)
 text(-0.2, 0.33, expression(t(nu == 2)), pos = 4)
 text(1, 0.25, pos = 4, expression(t(nu == 1000)))
 text(1.05, 0.23, pos = 4, "N(0,1)")
 legend("topright", legend = c(expression(nu == 1), expression(nu == 2), expression(nu == 
      5), expression(nu == 20), expression(nu == 1000), "N(0,1)"), lty = c(1:5, 
      2), lwd = c(1, 1, 1, 1, 1, 3))
\end{Sinput}
\end{Schunk}
