pub fn phi(x: f64) -> f64
Phi (CDF for Normal Dist)
$$\Phi(x) = \frac{1}{2}\left[1 + \text{erf}\left(\frac{x}{\sqrt{2}}\right) \right]$$