minutes<-c(107, 116, 127, 126, 111, 90, 125, 124, 92, 121, 131, 137, 95, 101, 101, 111, 120, 114, 98, 106, 118, 82, 110, 115, 107, 121, 99, 100, 87)

m<-mean(minutes)

std<-sqrt(var(minutes))

hist(minutes,breaks=10,xlab="....",

main="Histogram of ....")

curve(dnorm(x, mean=m, sd=std),

col="darkblue", lwd=2, add=TRUE)