def _cdf(self, x, mu):
        fac = sqrt(1.0/x)
        # Numerical accuracy for small `mu` is bad.  See #869.
        C1 = _norm_cdf(fac*mu/(x-mu))
        C1 += exp(1.0/mu) * _norm_cdf(-fac*(x+mu)/mu) * exp(1.0/mu)
        return C1
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/invgauss_gen._cdf/swapped operands of `/`