sd1 <- function (x)
{
	m <- mean(x);
	sqrt (sum ((x - m)^2) / length(x));
}