integrate s0 exp(-0.5 * (x-x0)^2/sigma^2) dx