library(MASS)
library(fitdistrplus)

#Simulation condition
k <- 1000 #Number of segment
n <- 10000 #Number of trial
runifmin <- 0
runifmax <- 6

#Random variable generation
s <- numeric(n)
for (i in 1:n) {
　x <- runif(k,runifmin,runifmax)
　s[i] <- sum(x)
}
truehist(s)
results <- summary(s)
print(results)
me <- mean(s)
st <- sd(s)
x<-seq(runifmin*k,runifmax*k,10)
lines(x,dnorm(x,me,st),lwd=2,lty=2,col="blue")