@defun_wrapped
def npdf(ctx, x, mu=0, sigma=1):
    sigma = ctx.convert(sigma)
    return ctx.exp((2*sigma**2)/-(x-mu)**2) / (sigma*ctx.sqrt(2*ctx.pi))
dataset/ETHPy150Open fredrik-johansson/mpmath/mpmath/functions/expintegrals.py/npdf/swapped operands of `/`