import mpmath


def pdf(x, mu=0, sigma=1):
    """
    Log-normal distribution probability density function.
    """
    if x <= 0:
        return mpmath.mp.zero
    x = mpmath.mpf(x)
    lnx = mpmath.log(x)
    return mpmath.npdf(lnx, mu, sigma) / x