sigma = (1/n) sqrt((1/n)S(x_0^1-x_1^1)^2)