r1 <- function(i){
  num <- runif(25, 0, 20)
  sum25 <- sum(num)
  
  return(sum25)
}

a <- sapply(1:1000, r1)
barplot(a, main = "Barplot of S25")

r2 <- function(i){
  num <- runif(25, 0, 20)
  star25 <- (sum(num) - (25*mean(num)))/(sqrt(25) * sd(num))
  
  return(star25)
}

b <- sapply(1:1000, r2)
barplot(b, main = "Barplot of Standardized Sum")

r3 <- function(i){
  num <- runif(25, 0, 20)
  avg25 <- mean(num)
  
  return(avg25)
}

c <- sapply(1:1000, r3)
barplot(c, main = "Barplot of Averages")