Constant FACTOR_BOUND 

const FACTOR_BOUND: f64 = 8.0;

Integration bound: |m| ≤ 8 captures the normal density to ~1e-15.