set.seed(12345) n<-10000; df_list<-c(5,10,15,30) par(mfrow=c(2,2)) for (i in 1:length(df_list)) { hist(rt(n, df=df_list[i], ncp=0), breaks=100, xlab="t", main=paste("df=", df_list[i])) }