x <- c(2, 5, 34, 62, 45, 71)
mean(x)
## [1] 36.5
var(x)
## [1] 820.3
sqrt(var(x))
## [1] 28.64088
sum(x^2)
## [1] 12095
n <- length(x) (sum(x^2) - n * mean(x)^2) / (n-1)