def norm_ppf(q):
    """Gaussian percent point function (inverse of `cdf`) at q"""
    return special.ndtri(q)