def _cdf(self, x, mu):
        trm1 = 1.0/mu - x
        trm2 = 1.0/mu + x
        isqx = 1.0/sqrt(x)
        return 1.0-_norm_cdf(isqx*trm1)-exp(2.0/mu)*_norm_cdf(-isqx*trm2)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/recipinvgauss_gen._cdf/original