\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, 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}
