integral p * exp(- (sqrt((x-m)^2 + n))^2 / (2 * q^2)) dx from a to b