Survival function 1 - cdf(x), accurate in the upper tail.
1 - cdf(x)
double sf(double x) => 0.5 * erfc((x - mean) / (std * math.sqrt2));