set.seed(12345)
round(datos<-rnorm(100, 25,5))
##   [1] 28 29 24 23 28 16 28 24 24 20 24 34 27 28 21 29 21 23 31 26 29 32 22
##  [24] 17 17 34 23 28 28 24 29 36 35 33 26 27 23 17 34 25 31 13 20 30 29 32
##  [47] 18 28 28 18 22 35 25 27 22 26 28 29 36 13 26 18 28 33 22 16 29 33 28
##  [70] 19 25 21 20 37 32 30 29 21 27 30 28 30 23 37 30 34 28 23 28 29 20 21
##  [93] 34 23 20 28 22 36 22 22
hist(datos,prob=TRUE, col="grey")
lines(density(datos))

shapiro.test(datos)
## 
##  Shapiro-Wilk normality test
## 
## data:  datos
## W = 0.9834, p-value = 0.242