Function normal_sf 

pub fn normal_sf(x: f64, mean: f64, sd: f64) -> Result<f64, EngineError>

Upper tail P(X > x), computed with erfc so large positive x keeps relative accuracy (never 1 - cdf).