import mpmath


def var(mu=0, b=1):
    """
    Variance of the Laplace distribution.
    """
    if b <= 0:
        raise ValueError('b must be positive.')
    with mpmath.extradps(5):
        mu = mpmath.mpf(mu)
        b = mpmath.mpf(b)
        return 2*b**2