powered by
Get the total probability from a region of a normal pdf
pnorm.range.log(lower, upper, mu, sd)