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