a <- c(28,21,26,16,18,13,15,22,19,22)
sd(a)
sd(a)/length(a)
t.test(a, mu=23)